Дифференциальные уравнения. Учебное пособие-1
.pdf
Подставив |
y ч.н. |
|
и |
ее |
|
|
производные |
в |
исходное |
уравнение, |
получим |
|||||||||||
ex Ax 2 4Ax Bx 2A 2B |
3ex Ax 2 2Ax Bx B |
2ex Ax 2 Bx xex . |
||||||||||||||||||||
Сокращаем на ex и приводим подобные члены: 2 Ax 2 A B x . |
||||||||||||||||||||||
Приравниваем коэффициенты при одинаковых степенях x в |
левой и |
|||||||||||||||||||||
правой частях тождества, находим: A 0,5, |
B 1. |
|
|
|
|
|
||||||||||||||||
|
x2 |
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Имеем y ч.н. |
|
|
x e |
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Общее решение уравнения имеет вид |
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
y y0.0. y |
|
x |
|
|
2 x |
x2 |
|
|
x |
. |
|
|
|||||||
|
|
|
ч.н. C e |
|
C e |
|
|
|
x e |
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
1 |
|
|
2 |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Указание к задаче 14.1. |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
Найти общее решение уравнения |
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
y 2 y |
5y 17cos2x . |
|
|
|
|
|
||||||
Решение. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Составим |
характеристическое |
|
уравнение |
и |
найдем |
его корни |
||||||||||||||||
K2 2K 5 0, |
K |
1 |
4 1 2i, |
1, |
2, тогда |
|
||||||||||||||||
|
|
|
|
1,2 |
|
|
|
|
e x C cos2x C sin2x . |
|
|
|
|
|
||||||||
|
|
|
|
|
y |
0.0. |
|
|
|
|
|
|||||||||||
|
i 2i |
|
1 |
|
|
|
|
2 |
|
|
|
|
|
|
|
|||||||
Так как |
|
не |
|
совпадает |
с |
корнями |
характеристического |
|||||||||||||||
уравнения, |
то |
y ч.н. Acos2x Bsin2x . |
Подставляя |
y ч.н. |
в |
исходное |
||||||||||||||||
уравнение, получаем тождество |
|
|
||
4Acos2x 4Bsin2x 2 2Asin2x 2Bcos2x 5 Acos2x Bsin2x 17cos2x 0. |
||||
Приравнивая к |
нулю суммарные коэффициенты при cos2x и sin2x , |
|||
получаем: cos2x: 4 A 4B 5A 17; |
sin 2x: 4B 4 A 5B 0 . |
|
||
Решая систему, |
находим: |
A 1, B 4 , тогда y ч.н. cos2x 4sin2x . |
Для |
|
общего решения получаем выражение |
|
|||
y C e x cos2x C e x |
sin 2x cos2x 4 sin 2x . |
|
||
|
1 |
2 |
|
|
Указание к задаче 14.2. |
|
|
||
Найти общее решение уравнения |
|
|||
|
y 4 y 4 y e2 x cos x . |
|
||
Решение. |
|
|
уравнение K2 4K 4 0 |
|
Составим |
характеристическое |
или |
||
K 2 2 0 . Корень K 2 |
кратности 2, поэтому y0.0. C1e2 x C2 xe2 x . Так |
|||
как i 2 i |
не совпадает с корнями характеристического уравнения, |
|||
|
|
31 |
|
|
положим y ч.н. Ae2 x cos x Be2 x sin x . Подставляя y ч.н. в исходное
уравнение и сокращая на e2 x , будем иметь:
3Acos x 4Asin x 3Bsin x 4Bcos x 4 2Acos x Asin x 2Bsin x Bcos x
4 Acos x Bsin x cos x 0.
Приравнивая к нулю суммарные коэффициенты при cos x и sin x , получим A 1, B 0, y ч.н. e2 x cos x . Общее решение имеет вид y C1e2 x C2 xe2 x e2 x cos x .
Указание к задаче 15.
Найти общее решение уравнения
y 2 y y ex ex sin x .
Решение. |
|
|
|
|
Характеристическое уравнение K2 2K 1 0 имеет корень |
K 1, |
|||
кратности 2. Значит, y |
0.0 |
C ex |
C xex . |
|
|
1 |
2 |
|
|
Для отыскания частного решения исходного уравнения найдем |
||||
частные решения двух уравнений: |
|
|||
|
|
y 2 y y ex , |
(1.18) |
|
|
|
y 2 y y ex sin x . |
(1.19) |
|
Правая часть уравнения (1.18) есть экспонента ex , значит, в формуле (1.16)
имеем |
1, 0, |
s 0, m 0 . |
Число |
i 1 |
является |
корнем |
|
характеристического |
уравнения |
кратности |
2, значит, |
частное решение |
|||
y(1) ч.н. |
уравнения (1.18) ищем в виде |
y(1) ч.н. Ax2ex . Подставив |
y(1) ч.н. |
||||
в уравнение (1.18), |
получим A 1 , т. |
е. |
y(1) ч.н. 1 x2ex . Правая часть |
||||
|
|
|
2 |
|
2 |
|
|
уравнения (1.19) равна: f (x) ex sin x , т. е. 1, 1, |
s 0, m 0 ; число |
||||||
i 1 i не является корнем характеристического уравнения. |
Значит, |
||||||
частное решение y(2) ч.н. уравнения (1.19) ищем в виде |
|
|
|||||
|
|
y(2) ч.н. ex Acos x Bsin x . |
|
|
|||
Подставляя в (1.19), имеем:
ex 2 Asin x 2B cos x 2ex B A sin x B A cos xex Acos x Bsin x ex sin x.
Сокращаем на ex , приводим подобные слагаемые:
sin x 2 A 2B 2 A B cos x 2B 2B 2 A A sin x 0 .
32
Приравнивая к нулю коэффициенты при sin x, cos x, получим A 0 , B 0 , y(2) ч.н. ex sin x . Общее решение исходного уравнения:
yC1ex C2 xex ex sin x 12 x 2ex .
2.5.Линейные уравнения второго порядка с переменными коэффициентами. Метод вариации произвольных постоянных
Признак. Уравнение имеет вид |
|
y P1 (x) y P2 (x) y Q(x) . |
(1.20) |
Метод решения. Согласно методу вариации постоянных общее решение уравнения (1.20) следует искать в виде y C1 (x) y1 C2 (x) y2 , где
y1 , y2 |
– фундаментальная система |
решений |
соответствующего |
|||
однородного уравнения |
|
|
|
|
||
|
|
y P1 (x) y P2 (x) y 0 . |
|
(1.21) |
||
Функции C1 (x) и C2 (x) определяются из системы: |
|
|
||||
|
|
y1C1 y2C2 0, |
y1C1 y2C2 Q(x) . |
(1.22) |
||
Сначала находим C1 , |
C2 , затем интегрированием C1 |
C2 . |
|
|||
Заметим, |
что при всяком фиксированном |
x система (1.22) |
представляет |
|||
собой систему из двух линейных уравнений относительно неизвестных C1 C1 (x) и C2 C2 (x) . Определитель матрицы системы (1.22) равен
определителю Вронского системы функций y1 ( x), y2 ( x) . Так как функции y1 ( x), y2 ( x) образуют фундаментальную систему решений однородного
уравнения (1.21), то определитель Вронского отличен от 0 при любом x . Значит, линейная система (1.22) имеет единственное решение.
Указание к задаче 16.
Найти решение задачи Коши:
y y sin x 5
2 cos x 1
2 ,
|
|
|
|||
y |
|
1, |
|
|
0. |
2 |
y |
2 |
|||
|
|
|
|
||
Решение.
Соответствующее однородное уравнение y y 0 .
Его характеристическое уравнение |
K2 1 0 имеет корни K |
i , и |
|
1,2 |
|
общее решение однородного уравнения имеет вид y0.0. C1 cos x C2 sin x ,
33
где |
|
|
|
|
|
|
y1 cos x, |
y2 |
sin x . |
|
|
|
Общее: |
|
решение |
|
исходного |
|
|
|
неоднородного |
|||||||||||||||||||||||||||||||||||||||||||||||||||||
уравнения ищем в виде y C1 (x) cos x C2 (x) sin x . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Составляем систему уравнений (1.22) для C1 , |
C2 : |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
cos xC1 (x) sin xC2 |
(x) O |
|
|
|
|
|
|
|
|
|
|
|
|
(1.23) |
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5 2 |
|
|
|
1 2 |
. |
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
sin xC1(x) cos xC2 (x) sin x |
|
|
|
cos x |
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||
Решаем систему (1.23) методом Крамера: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
cos x |
|
|
sin x |
|
cos2 x sin2 |
x 1 0; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||
|
|
sin x |
|
|
|
cos x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
1 |
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
sin x |
|
sin x 1,5 cos x 0,5 ; |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||
|
|
|
sin x |
|
2,5 |
cos x |
|
0,5 |
|
|
cos x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
2 |
|
|
|
|
cos x |
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
sin x 2,5 |
cos x 0,5 ; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||
|
|
|
sin x |
|
|
sin x |
|
2,5 |
cos x |
|
0,5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
получаем: C1 (x) |
1 |
sin x 1,5 cos x 0,5 |
|
|
|
|
|
1 |
|
|
; |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||
|
sin3 x cos x |
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
C2 (x) |
|
2 |
|
sin x 2,5 cos x 0,5 |
|
cos x |
|
|
. Интегрируя, находим |
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
sin5 x |
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
C (x) |
|
|
|
|
|
|
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
dx |
|
|
|
|
|
dtgx |
|
|
|
t tgx |
|
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
sin3 x cos x |
|
|
|
sin3 x cos x cos2 x |
|
|
|
|
tg3 x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
cos4 x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
dt |
|
|
|
|
|
t 1,5dt |
t 0,5 |
D |
|
2 |
D |
|
|
2 |
|
|
D |
2 |
|
|
|
ctgx D ; |
|
||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
t 3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,5 |
|
|
1 |
|
|
|
|
t |
1 |
|
|
|
|
tgx |
|
|
|
1 |
|
|
|
|
|
|
1 |
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
C2 (x) |
|
|
|
|
cos x |
dx |
|
cos x cos4 x |
|
|
dx |
|
|
1 |
|
|
|
dx |
|
|
|
|
|
|
t tgx |
|
|
|||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
sin5 |
|
|
|
|
|
|
sin5 x |
|
|
|
cos2 x |
|
tg5 x cos2 |
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
dt |
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
2 |
|
|
|
D2 |
|
|
|
|
|
2 |
|
|
|
|
D2 |
|
2 |
|
|
|
ctg 3 x D2 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
|
|
|
t 5 |
|
3 |
t 3 |
|
3 |
tg3 x |
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
34
Общее решение исходного уравнения есть
y C1 y1 C2 y2 |
2 |
ctgx D1 |
|
|
2 |
ctg |
3 |
x D2 |
|
cos x |
3 |
|
sin x |
||||||
|
|
|
|
|
|
|
|
|
D1 cos x D2 sin x 2
ctgx cos x 23 ctgx
ctgx sin x
D1 cos x D2 sin x 43
ctgx cos x.
Витоге получили
|
y D cos x D |
|
sin x 4 |
ctgx cos x , |
|
|
(1.24) |
||||||||||
|
1 |
|
|
|
2 |
|
|
|
|
3 |
|
|
|
|
|
|
|
где D1 , D2 |
– произвольные числа. |
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Теперь требуется найти такие значения D1 , D2 , |
при которых функция |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(1.24) удовлетворяла бы начальным условиям: y |
1, |
|
0 . |
|
|||||||||||||
y |
|
||||||||||||||||
Подставляя в (1.24) x , |
|
|
|
|
|
|
|
|
|
|
2 |
|
2 |
|
|
|
|
y 1, получаем D 1. |
|
|
|
|
|
||||||||||||
|
2 |
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Находим производную функции (1.24): |
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y C1 y1 C2 y2 C1 y1 C2 y2 C1 y1 C2 y2 C1 y1 C2 y2 |
|
|
|||||||||||||||
2 |
|
|
|
|
2 |
|
ctg |
3 |
x D |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
||||||||
ctgx D cos x |
3 |
|
|
sin x |
|
|
|
|
|||||||||
|
1 |
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
||
2 |
|
|
|
|
2 |
ctg3x D |
|
|
|
|
|
|
|||||
ctgx D sin x |
3 |
cos x. |
|
|
|
|
|||||||||||
|
1 |
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|||
Подставляя |
в последнее |
выражение |
|
x |
|
и приравнивая |
его |
к |
0, |
||||||||
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
D1 0, D2 |
1 |
|
|
получаем: |
0 D1 0D2 , |
D1 0 . |
|
Подставляя |
значения |
в |
|||||||||||
выражение (1.24), получаем решение задачи: y sin x 4 |
ctgx cos x . |
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
2.6. Некоторые приближенные методы решения дифференциальных уравнений
А) Приближенный метод решения с помощью рядов
Метод может быть использован для получения приближенных решений уравнения F(x, y, y ,..., y (n) ) 0 , где коэффициенты уравнения и
35
искомая функция таковы, что в окрестности рассматриваемой точки (или на заданном интервале) они могут быть представлены степенными рядами.
Метод решения.
а) Все функции, входящие в выражение для коэффициентов, необходимо разложить в ряд Тейлора по степеням (x x0 ) . Искомое
решение |
y(x) |
также |
представляется |
степенным |
рядом |
y a0 a1 (x x0 ) a2 (x x0 )2 |
|
|
|
||
... ak (x x0 )k |
с неопределенными |
||||
k 0
коэффициентами a0 , a1 , a2 ,... . Подставляя ряды в уравнение, собирая члены с одинаковыми степенями (x x0 ) k , k 0,1,2,... и приравнивая к нулю коэффициенты перед (x x0 ) k , k 0,1,2,..., получим систему уравнений, связывающих между собой a0 , a1 , a2 ,... . Для задачи Коши с условиями в точке x x0 коэффициенты a0 , a1 , a2 ,..., an 1 , оставшиеся
произвольными после удовлетворения уравнения, определяются из начальных условий
y(x0 ) y0(0) , |
y (x0 ) y0(1) , |
y (x0 ) y0(2) ,..., |
y (n 1) (x0 ) y0(n 1) , где |
y0(0) , |
||||||||
y (1) , y |
(2) ,... – заданные числа, при этом ak |
y0(k ) |
, k 0,1,..., n 1. |
|
||||||||
|
|
|||||||||||
0 |
0 |
|
|
|
|
|
|
k! |
|
|
||
б) |
Если уравнение имеет вид, |
разрешенный относительно старшей |
||||||||||
производной |
y |
(n) |
|
|
(n 1) |
) , |
|
то |
коэффициенты |
ak , |
||
|
f (x, y, y ,..., y |
|
|
|||||||||
k 0,1,..., n 1 |
для задачи Коши в точке x x0 |
определяются из начальных |
||||||||||
условий по указанным выше формулам, а остальные коэффициенты ak , k n , могут быть найдены дифференцированием этого уравнения (k n) раз:
ak |
|
y (k ) (x0 ) |
|
1 |
|
|
d k n f |
|
. |
|
|
||||||||
k! |
k! |
dx k n |
|
||||||
|
|
|
|
|
x x0 |
||||
Замечание. Ряды Тейлора для коэффициентов уравнения и для полученного решения y(x) могут оказаться сходящимися не при всех x из
заданного интервала, а только лишь в некоторой окрестности точки x0 . Тогда решение y(x) можно использовать для исследования локального поведения интегральных кривых только в окрестности этой точки.
Пример 2.5. Найти четыре первых отличных от нуля члена разложения в степенной ряд решения дифференциального уравнения,
удовлетворяющего начальному условию: y xy 2 e2 x , y(0) 1.
36
Решение. а) Представим функцию e2x |
рядом Тейлора |
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
e |
2x |
1 |
2x |
2x |
2 |
|
4 |
x |
3 |
... . |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|||||||||||||||
Искомое решение y(x) также |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
зададим |
|
степенным |
рядом с |
|||||||||||||||||||||||||
неопределенными коэффициентами ( x0 0 ): |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
y a |
0 |
|
a x a |
2 |
x 2 a |
3 |
x3 |
a |
4 |
x 4 |
... . |
|
|
|
|
|
||||||||||
Подставляя e2x |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
и y(x) |
|
в дифференциальное уравнение, получим |
|||||||||||||||||||||||||||
a |
2a |
2 |
x 3a |
3 |
x 2 |
... x(a |
2 2a |
0 |
a x a 2 x 2 |
2a |
0 |
a |
2 |
x 2 |
...) |
||||||||||||||
1 |
|
|
|
|
|
|
|
|
|
0 |
|
|
|
1 |
|
|
|
1 |
|
|
|
|
|||||||
|
|
|
|
|
|
|
(1 2x |
2x |
2 |
|
|
4 |
x |
3 |
...) . |
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Приравнивая коэффициенты при одинаковых степенях x в левой и правой частях равенства, имеем
a 1, 2a |
2 |
a 2 2 , 3a |
3 |
2a |
0 |
a 2 . |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
1 |
|
|
|
0 |
|
|
|
|
|
1 |
|
|
|
|
|
|
1 , a |
|
4 . |
|
||||||||
Так как |
y(0) 1, то a |
0 |
1. Тогда a |
1, a |
2 |
3 |
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
2 |
|
3 |
|
||||
Первые четыре члена разложения для y(x) |
|
|
|
|
|
|
||||||||||||||||||||||
имеют вид |
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
y 1 x |
1 |
x |
2 |
|
4 |
x |
3 |
.... |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
2 |
|
3 |
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
б) Так как уравнение имеет вид, разрешенный относительно |
|||||||||||||||||||||||||||
производной, коэффициенты ak |
можно определить дифференцированием |
|||||||||||||||||||||||||||
этого |
уравнения. |
|
Имеем |
|
|
|
y xy 2 e2 x , |
|
y y 2 2xyy 2e2 x , |
|||||||||||||||||||
y 4 yy 2xy 2 |
2xyy 4e2 x . |
|
|
|
|
|
|
|
|
y |
(0) 8 . |
Учитывая, |
что |
|||||||||||||||
Тогда |
y(0) 1, |
y (0) 1, |
|
y (0) 1, |
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
ak |
|
y (k ) (0) |
, получим заданное выше разложение. |
|
|
|
|
|||||||||||||||||||||
|
|
k! |
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Б) Метод малого параметра |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
Метод можно использовать для построения приближенных решений |
|||||||||||||||||||||||||||
уравнения F (x, y,..., y (n) ) 0 , |
при |
условии, что уравнение или условия |
||||||||||||||||||||||||||
для |
искомой функции y(x) |
содержат |
(или |
позволяют |
ввести) малый |
|||||||||||||||||||||||
параметр . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
y(x) |
Метод решения. Коэффициенты уравнения, а также искомое решение |
|||||||||||||||||||||||||||
|
представляется в виде разложения в ряд (как правило, степенной) по |
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
малому |
параметру |
: |
y(x) f k |
( ) yk (x) , |
|
1. |
Во |
многих задачах |
||||||||||||||||||||
f0 ( ) 1, а |
y0 (x) |
|
|
|
|
|
k 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
– |
известное |
|
возмущаемое |
решение. |
Уравнения |
для |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
37 |
|
|
|
|
|
|
|
|
|
|
|
|
yk (x) получим, приравнивая к нулю суммарные коэффициенты для членов одинакового порядка.
Пример 2.6. Построить приближенное решение нелинейного д. у. y y 16 y3 , 1, y O(1) .
Решение. Будем искать y(x) в виде ряда по степеням : y y0 (x) y1 (x) .... Подставляя это разложение в исходное уравнение,
получим
( y0 y1 ...) ( y0 y1 ...) 6 ( y03 3y02 y1 ...) .
Приравнивая к нулю коэффициенты при одинаковых степенях в левой и правой частях уравнения, получим систему уравнений для определения y0 , y1, ... .
y y0 |
0 , y1 y1 16 y03 , ... . |
(1.25) |
Решение этой линейной системы значительно проще, чем решение исходного нелинейного уравнения. Находя сначала общее решение уравнения для y0 , а затем для y1, получим
|
|
|
|
|
|
y0 a0 sin(x 0 ) , |
|
|
|
|
|
|
|
(1.26) |
||||||||||
|
|
|
|
|
|
|
|
a3 |
|
|
|
|
|
a3 |
|
|
|
|
|
|
|
|
||
y |
a sin(x |
1 |
) |
|
0 |
x cos(x |
0 |
) |
0 |
|
sin 3(x |
0 |
) , |
|
|
|||||||||
|
|
|
|
|
|
|||||||||||||||||||
|
1 |
1 |
|
|
|
16 |
|
|
|
|
|
|
192 |
|
|
|
|
|
|
|||||
где ak , k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
– произвольные постоянные. Приближенное решение задачи |
|||||||||||||||||||||||
имеет вид |
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
y a0 sin(x 0 ) |
a1 sin(x 1 ) |
a0 |
|
x cos(x 0 ) |
|
a0 |
sin 3(x 0 ) |
.... |
||||||||||||||||
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
16 |
|
|
|
|
|
192 |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Пример 2.7. Построить первые два члена разложения функции |
y(x) , |
|||||||||||||||||||||||
удовлетворяющей д. у. |
y sin y 0 , для малых y |
( y o(1) ). |
|
|
|
|
||||||||||||||||||
Решение. Учитывая, что y мало |
( y O( ) ), будем |
искать |
||||||||||||||||||||||
приближенное решение в виде y y |
0 |
(x) |
3 y (x) .... Разложим sin y по |
|||||||||||||||||||||
степеням |
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
1 y3 ) 3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
: sin y y |
|
( y |
... . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
0 |
1 |
6 |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
38
Подставляя разложение для y в уравнение, получим
( y0 y1 3 ...) [ y0 ( y1 16 y03 ) 3 ...] 0 .
Приравнивая к нулю коэффициенты перед одинаковыми степенями , получим для y0 , y1,... систему уравнений (1.25), общее решение которой
имеет вид (1.26).
Пример 2.8. y cos(x y) , 1, |
y 0(1) , |
y(x) ? . |
|
|
||||||||||
Решение. Будем искать решение в виде |
y y0 y1 .... Тогда, |
|||||||||||||
учитывая, что cos(x y) cos x y sin x ..., получим |
|
|
|
|
||||||||||
|
|
y0 y1 ... cos x sin x( y0 y1 2 ...) ... . |
|
|
||||||||||
Отсюда y0 cos x , y1 y0 sin x . |
Общее решение последней системы |
|||||||||||||
уравнений |
|
|
|
|
|
|
x |
|
1 sin 2x C |
|
|
|||
y |
0 |
(x) sin x C , |
y (x) C cos x |
|
2 |
, |
||||||||
|
||||||||||||||
|
1 |
1 |
1 |
|
2 |
|
4 |
|
|
|
||||
где C1 , C2 – произвольные постоянные. Тогда |
|
|
|
|
|
|||||||||
y(x) (sin x C ) [C cos x |
x |
1 sin 2x C |
2 |
] O( 2 ) . |
||||||||||
|
||||||||||||||
|
|
1 |
1 |
2 |
4 |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
||||
3. Системы дифференциальных уравнений
Укажем два способа решения систем дифференциальных уравнений. А) Метод Эйлера применим к системе линейных однородных
уравнений с постоянными коэффициентами ( aij const )
|
|
a |
a |
|
... |
|
|
a |
|
|
|
|
|
||||
|
|
11 |
|
12 |
|
|
|
|
|
1n |
|
|
|
|
|||
dY |
AY , |
a21 |
a22 |
... |
|
|
a2n |
|
|
|
|||||||
dx |
A |
|
... |
... |
|
|
|
... |
, Y |
|
|||||||
|
... |
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
an2 |
... |
|
|
|
|
|
|
|
|
|
||||
|
|
an1 |
|
|
ann |
|
|
|
|
||||||||
или, в развернутом виде |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
dy1 |
a |
|
y |
a |
|
y |
|
... a |
|
y , |
|
||||
|
|
dx |
11 |
1 |
12 |
|
2 |
|
|
1n |
|
n |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dy2 |
a21 y1 |
a22 y2 ... a2n yn , |
|||||||||||||
|
|
|
|||||||||||||||
|
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
................................................. |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dyn |
a |
|
y |
a |
n2 |
y |
2 |
|
... a |
nn |
y |
. |
|||
|
|
|
|
n1 1 |
|
|
|
|
|
n |
|
||||||
|
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y1 y2
... , yn
(1.27)
(1.28)
39
Общее решение системы (1.27) имеет вид
|
y (i) |
|
|
||
|
|
1 |
|
|
|
n |
|
(i) |
, i 1,2,..., n , |
(1.29) |
|
Y CiY (i) , Y (i) y2 |
|
||||
i 1 |
|
... |
|
|
|
|
|
(i) |
|
|
|
|
yn |
|
|
|
|
где Ci – произвольные постоянные, а векторы Y (i) образуют линейно-
независимую систему решений уравнений (1.27). Для отыскания Y (i) следует найти корни r (i) (i 1, n ) характеристического уравнения
a11 r |
a12 |
a21 |
a22 r |
... |
... |
an1 an2
Для каждого из корней необходимо алгебраических уравнений относительно i
... a1n
... a2n
... ...
... ann
решить
( i 1, n ):
0 . |
(1.30) |
систему линейных
(a |
r) |
1 |
a |
|
1 |
... a |
|
n |
0, |
|
||||||||
|
11 |
|
|
|
|
|
12 |
|
|
|
1n |
|
|
|
||||
a |
21 1 |
(a22 |
r) 2 |
... a2n n |
0, |
(1.31) |
||||||||||||
|
.................................................., |
|||||||||||||||||
|
|
|||||||||||||||||
a |
|
1 |
a |
n2 |
|
2 |
... (a |
nn |
r) |
n |
0. |
|
||||||
|
n1 |
|
|
|
|
|
|
|
|
|
|
|
||||||
Если все корни характеристического уравнения (1.30) различные, то
каждому действительному корню r (i) |
в выражении (1.29) соответствует |
||||||||||||||
решение |
|
|
|
|
|
|
|
|
|
(i) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
Y (i) er |
(i) |
x |
|
(i) |
, |
|
|
|
|
||||
|
|
|
|
2 |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
.... |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(i) |
|
|
|
|
|
а каждому комплексному корню r (i) |
|
|
|
|
n |
|
|
|
|
|
|||||
|
– два решения вида |
|
|||||||||||||
|
|
|
|
(i) |
|
|
|
|
|
|
|
|
|
(i) |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
1 |
|||
(i) |
Re e |
r(i) x |
(i) |
, Y |
(i) |
Im e |
r(i) x |
|
|
(i) |
|||||
Y |
|
|
2 |
|
|
|
2 . |
||||||||
|
|
|
|
.... |
|
|
|
|
|
|
|
|
.... |
||
|
|
|
|
(i) |
|
|
|
|
|
|
|
|
|
(i) |
|
|
|
|
n |
|
|
|
|
|
|
|
|
n |
|||
40
