0.04 |
|
|
|
|
|
0.02 |
|
|
|
|
|
i(t) |
|
|
|
|
|
0 |
|
|
|
|
|
0.02 |
|
4 |
|
4 |
|
0 |
|||||
1 10 |
t |
2 10 |
|
||
|
|
|
|
Рис. 2.4. Переходной процесс при R 2
2.2. Решение дифференциальных уравнений операторным методом
Если дифференциальное уравнение, которое необходимо решить, имеет постоянные коэффициенты, то для этих целей можно использовать преобразование Лапласа. Используя его свойства, переводят дифференциальное уравнение в соответствующее алгебраическое, аналитически решив которое и применив к нему операцию обратного преобразования Лапласа, можно получить решение дифференциального уравнения в общем виде.
Возьмем в качестве примера для решения дифференциальное уравнение (1.11), описывающее процессы в последовательном колебательном контуре
d2i t |
|
R |
|
di t |
|
1 |
i t 0. |
|
dt2 |
|
|
|
|
||||
|
L dt LC |
|||||||
Приступая к решению уравнения, получим, прежде всего, его изображение с помощью прямого преобразования Лапласа. Для этого используем специальный оператор
«laplace» (Лапласа) панели Symbolic (символьные): ● laplace, ● →.
60
В левом маркере этого оператора должно быть определено преобразуемое выражение без знака равенства, в правом
– переменная, по которой осуществляется преобразование:
d2 |
i t |
R |
|
d |
i t |
1 |
i t laplace, t S (S laplace i t ,t,s |
||||||||
dt2 |
L dt |
|
|||||||||||||
|
|
|
|
|
|
L C |
|
|
|||||||
i 0 ) |
|
t 0 |
R |
(S laplace i t ,t,s i 0 ) |
1 |
laplace i t ,t,s |
|||||||||
|
|||||||||||||||
|
|
||||||||||||||
|
d |
|
|
|
|
L |
|
|
L C |
||||||
|
i t . |
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
||||||||
|
|
|
dt |
|
|
|
|
|
|
|
|||||
Далее необходимо упростить это громоздкое выражение. Для этого переменную laplace i t ,t,s заменим на Z . Начальное условие для первой производной
t 0
d i t dt
запишем в виде i 0.
В результате получим выражение
|
R |
1 |
|
|
|
S Z i0 |
|
Z. |
|
S S Z i0 i 0 |
|
|
||
|
L |
L C |
||
После проведенных замен уравнение следует разрешить относительно переменной Z, т.е. использовать специальный оператор символьного решения уравнений solve пане-
ли Symbolic (Символьные) и принять, что i0 0, а i 0 E .
L
При таких начальных условиях после применения оператора solve получаем, что
S S Z |
E |
|
R |
S Z |
1 |
|
E |
|
|
|
|
|
Z solve,Z |
S2LC rSC 1 |
C. |
||
L |
L |
L C |
||||||
61
Теперь, чтобы получить решение дифференциального уравнения, к полученному решению-изображению нужно применить операцию обратного преобразования Лапласа. Сделать это можно при помощи оператора invlaplace панели
Symbolic (Символьные), т. е.
E |
invlaplase,S |
|
S2 L C r S C 1 C |
simplify |
. |
При расчете обратного преобразования Лапласа желательно одновременно провести и упрощение получаемого громоздкого выражения, что можно сделать, задействовав одновременно с оператором invlaplace и специальный оператор
Simplify (Упростить).
Найденная функция i t после применения обратного преобразования Лапласа будет иметь следующий вид
1
4 L R2 C 2
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
R |
|
|
|
|
|
|
|
L2C |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
i t 2E Cexp |
|
|
|
|
t |
L |
|
|
|
|
|
|
|
|||||||
2 |
|
|
|
|
4 L R2 |
C |
||||||||||||||
|
|
|
|
|
|
|
L |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
||||
|
|
1 |
4 L R2 |
C |
|
|
|
|
|
|
|
|
||||||||
|
2 |
|
|
|
|
|
||||||||||||||
sin |
|
|
|
|
|
|
|
|
|
|
|
|
|
t . |
|
|
|
|
||
|
2 |
L2 C |
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Такое решение дифференциального уравнения является пусть и записанным не в самой удобной форме, но полностью соответствует поставленной задаче и дает такие же результаты, как и при применении процедуры Given odesolve, как это было сделано в пп.1.5.1 для решения дифференциального уравнения второго порядка.
После соответствующего ввода численных значений будем иметь:
|
|
C : 5 10 9 |
|
|
|
E : 100 |
L : 0.05 |
R : |
L C 0,5 |
||
62
|
|
|
|
|
|
|
|
|
|
|
4 L R2C |
||||
|
|
|
|
|
1 R |
|
|
|
|
|
|
L2 |
C |
|
|
|
|
|
|
|
L |
|
|
|
|||||||
i t : 2 E exp |
|
t |
|
4 L R2 C |
|
||||||||||
|
|
|
|
2 L |
|
|
|
|
|
|
|
||||
|
|
1 |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
sin |
|
4 L R2 C 2 |
|
|
|
|
|
|
|||||||
|
2 |
|
L2 |
C |
|
|
|
t |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t : 0, 0.000001.. 0.0002 |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
0.04 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0.02 |
|
|
|
|
|
|
|
|
|
|
|
|
i(t) |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0.02 |
|
|
|
|
4 |
|
|
|
|
4 |
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
||||
|
|
|
|
|
1 10 |
t |
|
|
2 10 |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 2.5. График переходного процесса |
|
||||||||||||
Этот |
|
результат |
(рис. |
2.5) |
|
|
соответствует |
условию |
|||||||
R 2 , где |
L. |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
C |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
63 |
|
|
|
|
|
|
2.3. Применение процедур численного дифференцирования и интегрирования к вычислениям методом интеграла Дюамеля
Вычисление интеграла Дюамеля (1.53) предусматривает применение процедур численного дифференцирования и интегрирования. При необходимости можно предварительно осуществить процедуру символьного дифференцирования, которое может быть использовано для подавляющего числа функций, в отличие от того же символьного интегрирования, имеющего весьма ограниченное применение.
Приступая к вычислению производной какой-либо функции, следует найти соответствующий оператор. Оператор производной расположен на панели Calculus (Вычисления) и, помимо того вводится сочетанием клавиш [Shift] + [/], т.е.
d . d
Если в результате дифференцирования должна быть получена функция производной, следует обратиться к возможностям символьного вычисления. В качестве оператора вывода следует использовать оператор символьного вывода
(«→» – Evaluate Symbolically).
При символьном дифференцировании можно оперировать функциями нескольких переменных и функциями с параметрами. Оператор дифференцирования может сочетаться с любым вычислительным оператором.
Пример 2.7. Символьное дифференцирование
d |
ln ex sin x |
|
1 cos x |
ln exp x sin x |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, |
||
dx |
|
|
|
x |
|
|
|
x |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
x2 |
|||||
|
|
d a |
|
1 |
|
|
4 |
|
|
|
|
|||||
|
|
|
|
|
|
|
a |
|
. |
|
|
|
||||
dx |
x3 |
30 |
|
|
|
|
|
|||||||||
x 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
64