Материал: Переходные процессы и спектры. Матвеев Б.В

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

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

Источник: https://studfile.net/preview/16563501/