Операционное исчисление
.pdf
4.2. Решение задачи Коши для обыкновенных дифференциальных уравнений с постоянными коэффициентами
Операционный метод значительно упрощает вычисления по сравнению со стандартными методами, в частности, при этом осо-
бенно эффективен, когда f t – разрывная функция.
Пример 24. Решить задачу Коши:
24.1. x x 1, x 0 1,
24.3. x 2x t t3 , 3
x 0 x 0 0,
24.2. x x 2cost,
x(0) 0, |
|
1, |
x (0) |
24.4. x 3x 2x 2e3t , x(0) 1, x (0) 3,
24.5. x x 6x 3(cos3t sin3t), |
24.6. x x 6x 0, |
|
x 0 0, x 0 3, |
x 0 15, x 0 2, x 0 56. |
|
24.1. Пусть функция x t имеет изображение X p . |
Тогда по |
|
теореме о дифференцировании |
оригинала получим |
x t |
pX p x 0 pX p 1. Применим преобразование |
Лапласа |
|
к обеим частям уравнения. Выпишем получившееся операторное
уравнение pX p 1 X p |
|
1 |
. Откуда получим X p |
1 |
. Та- |
|||||||||||
|
|
p |
||||||||||||||
|
|
|
|
|
|
|
|
p |
|
|
|
|
|
|
||
ким образом x t 1. |
|
|
|
|
|
|
|
|
|
|
|
|||||
24.2. Перейдем от оригиналов к изображениям |
|
|
|
|
|
|||||||||||
|
|
|
x t X p , |
x t pX p x 0 pX p , |
|
|
|
|||||||||
|
|
t |
|
2 |
|
|
0 p |
2 |
X p 1, |
|
|
p |
|
|
||
x |
p |
|
|
|
cost p2 1 . |
|||||||||||
|
|
X p px 0 x |
|
|||||||||||||
61
Запишем уравнение для изображений
p2X p 1 X p p22 p 1.
Решим уравнение для изображений
X p |
2p |
|
1 |
|
. |
|
p2 1 2 |
p2 1 |
|||||
|
|
|
||||
По теореме о дифференцировании изображения найдем оригинал первого слагаемого
p22p1 2 p21 1 t sint .
Следовательно, решение имеет вид
x t t sint sint t 1 sint.
24.3. Пусть x t X p . Перейдем в уравнении к изображениям
|
2 |
|
|
|
|
|
|
|
1 |
|
|
1 3! |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
p |
|
|
|
0 2X p p2 3 p4 . |
|||||||||||
|
X p px 0 x |
||||||||||||||
Так как x 0 x 0 0, то |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
p2 2 X |
|
p2 2 |
, |
X p |
|
1 |
. |
|
|
|
|||
|
|
|
4 |
4 |
|
|
|
||||||||
|
|
|
|
|
p |
|
|
|
|
p |
|
|
|
||
62
Найдя оригинал по данному изображению, получим решение задачи Коши
|
|
X p |
|
1 |
3! |
|
|
t3 |
|
|
|
t3 |
x t |
. |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
p4 |
3!p4 |
3! |
|
6 |
|
|
||||||||||||||||||||||
24.4. Перейдем от оригиналов к изображениям |
|
|
|
|
|
|||||||||||||||||||||||||
|
x t X p , |
x t pX p 1, |
|
|||||||||||||||||||||||||||
|
|
t |
|
|
|
2 |
X p p 3, |
|
|
3t |
|
|
|
1 |
|
|
||||||||||||||
x |
p |
|
|
e |
|
|
|
|
|
p 3. |
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
Запишем уравнение для изображений |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
p2X p p 3 3pX p 3 2X p |
2 |
. |
||||||||||||||||||||||||||||
|
|
p 3 |
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Решим уравнение для изображений |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
p2 3p 2 X p |
|
|
2 |
|
|
|
|
p, |
|
|
|
|
|
|||||||||||||||
|
|
|
p |
3 |
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
X p |
|
|
|
|
|
p2 |
3p 2 |
|
|
|
|
|
|
|
1 |
|
|
. |
|
|||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
p 3 p2 3p 2 |
p 3 |
|
||||||||||||||||||||||||||
Найдем оригинал для функции X p |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
X p |
1 |
|
e3t |
|
x t . |
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
p 3 |
|
|
|
|
|
|
|||||||||||||||||||||
24.5. Пусть x t |
X p . Перейдем в уравнении к изображе- |
|||||||||||||||||||||||||||||
ниям
p2X p px 0 p 0 pX p x 0 6X p
3 p2p 9 p23 9 .
63
Так как x 0 0, |
x 0 3, получим |
|
|
|
p2X p 3 pX p 6X p |
3 p 3 |
, |
||
p2 9 |
||||
|
|
|
||
X p p23 9 sin3t x t .
24.6. Перейдем от оригиналов к изображениям
x t |
X p , |
x t pX p 15, |
||||||||
|
x t p2X p 15p 2, |
|||||||||
x t p3X p 15p2 2p 56. |
||||||||||
Решим уравнение для изображений |
|
|
|
|||||||
p3 p2 6p X p 15p2 13p 36, |
X p |
15p2 13p 36 |
. |
|||||||
|
||||||||||
|
|
|
|
|
|
|
|
|
p p 3 p 2 |
|
Функция X p |
правильная рациональная несократимая дробь, |
|||||||||
для которой точки |
p1 0, p2 3, p2 |
2 |
являются простыми полю- |
|||||||
сами. Так как |
|
|
|
|
|
|
|
|
|
|
P p 15p2 13p 36, |
Q p p3 p2 6p, |
|||||||||
|
|
Q p 3p2 2p 6, |
||||||||
то, для p1 0 |
|
|
|
|
|
|
|
|
|
|
|
|
P p |
|
p 0 |
36 |
6 , |
||||
|
|
|
||||||||
|
|
Q p |
|
|
6 |
|||||
|
|
|
|
|
|
1 |
|
|
|
|
64
для |
p2 3 |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
P |
|
p |
|
|
|
|
|
60 |
4 , |
|||
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
Q' p |
|
p 3 |
15 |
|||||||||
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
для |
p3 2 |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
P |
|
p |
|
|
|
|
|
50 |
5, |
|||
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|||||||||
|
|
Q' p |
|
p 2 |
10 |
||||||||||
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
ипо второй теореме разложения получим
xt 6 5e 2t 4e3t .
Замечание. Во многих практических задачах правая часть дифференциального уравнения задается графически. В этом случае алгоритм решения не изменяется, а для нахождения изображения оригинала, заданного графиком, используются теорема запаздывания и методы из разд. 2.2.
Пример 25. Решить задачу Коши для дифференциального уравнения с правой частью, заданной графически:
f(t)
25.1. |
x x f t , |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
x(0) x (0) 0, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
|
|
||
|
|
|
0 |
|
2 |
|
t |
||||
|
|
2 |
f(t) |
|
|
|
|
|
|
|
|
25.2. |
x 4x f t , |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x(0) x (0) 0, |
|
|
|
|
|
|
|
|
|
t |
|
|
0 |
|
1 |
2 |
|
|
|
|||
65
25.1. Перейдем от оригиналов к изображениям
x t X p , |
x t pX p , |
x t p2X p , |
f t t 2 t 4 1p e 2p e 4p .
Решим уравнение для изображений
p2 1 X p 1p e 2p e 4p ,
X p p p12 4 e 2p e 4p .
Так как
p p12 1 1p p2p 1 1 cost ,
то
x t 1 cos t 2 t 2 1 cos t 4 t 4 .
Решение задачи Коши можно представить в аналитическом виде:
|
0, |
|
t 2, |
|
1 cos t 2 , |
|
2 t 4, |
x t |
|
||
|
|
, |
t 4. |
cos t 4 cos t 2 |
|||
25.2. Перейдем от оригиналов к изображениям |
|||
x t X p , |
x t pX p , |
x t p2X p , |
|
66
ft 2t t 2t t 1 4 2t t 1 4 2t t 2
2t t 4 t 1 t 1 2 t 2 t 2
|
|
2 |
|
|
4 |
|
|
|
|
2 |
e 2p |
2 |
|
|
1 2e p e 2p . |
||||||||||||
|
|
|
|
|
e p |
|
|
|
|
|
|||||||||||||||||
p2 |
p2 |
p2 |
|
p2 |
|||||||||||||||||||||||
Решим уравнение для изображений |
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
p2 4 X p |
2 |
1 2e p e 2p , |
|
|
||||||||||||||||||||
|
|
|
p2 |
|
|
||||||||||||||||||||||
|
|
|
X p |
2 |
|
1 2e p e 2p . |
|
|
|||||||||||||||||||
|
|
|
p2 p2 4 |
|
|
||||||||||||||||||||||
Так как |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
1 |
1 |
|
|
1 |
|
|
|
|
|
1 |
1 |
|
|
||||||
|
|
|
|
|
2 |
|
|
|
|
|
|
|
2 t |
2sin2t |
, |
||||||||||||
|
p2 p2 |
4 |
|
p2 |
|
p2 |
4 |
|
|||||||||||||||||||
то |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x t |
1 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
1 |
|
|
1 |
|||||
2 |
t |
2 |
|
sint |
t t |
|
2 |
sin t 1 t |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
1 |
|
|
|
|
1 |
sin t |
2 |
|
t 2 . |
|
|
|||||||||||
|
|
|
|
2 |
|
t 2 |
2 |
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Замечание. Если в задаче Коши за начальный момент времени взято t t0 0, то вводят новую переменную t t0 . Тогда 0
при t t0 .
Пример 26. Решить задачу Коши |
|
|
|
|
|
|
|
|||
|
x x t, |
|
|
x x 2sint, |
|
|
|
|||
26.1. |
x 1 0, |
26.2. |
|
|
0, |
|
|
1. |
||
|
x 1 1, |
|
x |
|
x |
2 |
|
|||
|
|
|
|
|
2 |
|
|
|
|
|
67
26.1. Положим |
t 1, x t x 1 z . Тогда уравнение |
||
и начальные условия примут вид |
|
||
z z 1, |
z 0 1, |
z 0 0. |
|
Перейдем от оригиналов к изображениям |
|||
z Z p , |
z pZ p 1, |
z p2Z p p, |
|
1 p12 1p .
Запишем уравнение для изображений
p2Z p p pZ p 1 p12 1p .
Решая операторное уравнение и переходя к оригиналам, получим
Z p p13 1p 1 22 z .
Возвращаясь к исходной переменной t , получим решение задачи Коши
|
x t |
1 |
t 1 2 |
. |
|
|||
|
2 |
|
|
|
||||
|
|
|
|
|
|
|
|
|
26.2. Положим t |
|
, |
|
|
|
|
z . Тогда уравне- |
|
2 |
x t x |
|
||||||
|
|
|
|
|
|
2 |
|
|
ние и начальные условия примут вид |
|
|
|
|
||||
z z 2sin |
|
|
|
z 0 0, |
z 0 1. |
|||
|
, |
|||||||
|
|
|
2 |
|
|
|
|
|
68
Перейдем от оригиналов к изображениям |
|
|
|
|
||||||
z Z p , |
z pZ p , |
z p2Z p 1, |
||||||||
2sin |
|
|
|
2cos |
|
2p |
. |
|||
|
|
|
|
|
||||||
p2 |
1 |
|||||||||
|
|
|
2 |
|
|
|
||||
Запишем уравнение для изображений
p2Z p 1 Z p p22p1.
Решим уравнение для изображений
|
|
Z p |
1 |
|
|
|
2p |
. |
|
|
|
|
|
p2 1 |
p2 1 2 |
|
|
||||||
|
|
|
|
|
|
|
|
||||
Переходя к оригиналам получим |
|
|
|
|
|||||||
1 |
sin , |
|
p |
sin cos |
1 |
sin , |
|||||
|
p2 1 |
p2 |
1 2 |
|
2 |
||||||
z 1 sin .
Возвращаясь к исходной переменной t , получим решение исходной задачи Коши
|
t |
|
|
|
|
|
|
|||
x t 1 |
2 |
sin t |
2 |
|
t 1 |
2 |
cost . |
|||
|
|
|
|
|
|
|
|
|||
69
4.3.Общее решение дифференциального уравнения
спостоянными коэффициентами
Операционное исчисление позволяет найти не только частное, но и общее решение дифференциального уравнения
an x n an 1x n 1 ... a1x a0x f (t).
Для этого достаточно положить x(k )(0) Ck ,k 0,...,n 1. Решив
задачу Коши с произвольными начальными условиями, мы получим общее решение уравнения. Подставляя в полученное общее реше-
ние конкретные значения для Ck можно находить решения различных задач Коши для заданного уравнения.
Пример 27. Найти общие решения дифференциальных уравнений
27.1. |
x 4x 0, |
|
|
27.2. x 2x x et . |
|
|||||||
27.1. |
Выберем |
произвольные |
начальные |
условия: x 0 C1, |
||||||||
x 0 C2 . Пусть |
x t X p . Тогда |
x t pX p C1, |
||||||||||
x t |
p2X p pC1 C2. Для задачи Коши x 4x 0, |
x 0 |
||||||||||
C1, |
x 0 C2 операторное уравнение иметь вид |
|
||||||||||
|
|
p2X p |
C p C |
2 |
4X p 0. |
|
||||||
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
C1p C2 |
C1 cos2t |
C2 |
|
|
|||||
|
|
X p |
|
|
|
sin2t x t . |
|
|||||
|
|
p2 4 |
|
2 |
|
|||||||
Решение задачи зависит от двух произвольных постоянных C1 и C2 , является общим решением данного однородного уравнения. Общее решение дифференциального уравнения x 4x 0 можно также записать в виде x t c1 cos2t c2 sin2t , где c1,c2 - произвольные постоянные.
70
