Бабич Диференціальні рівняння
.pdfякою б не була точка (x , y , y0 , . . . , y(n−1)) D, у деякому околi точки x iсну¹ ¹диний розв'язок рiвняння0 0 0 (2.10.2), який задовольня¹ умови (2.1.3)0.
Як i для рiвнянь першого порядку, задачу вiдшукання розв'язку рiвняння (2.1.1) або (2.1.2), який задовольня¹ початковi умови (2.1.3), називають задачею Кошi.
öiþЗагальним розв'язком рiвняння n-го порядку (2.1.2) називають функ-
|
|
|
|
y = ϕ(x, C1, C2, . . . , Cn), |
|
(2.1.4) |
||||||||||
залежну вiд |
x |
òà |
n |
сталих |
C1 |
, |
C2 |
, . . . , |
Cn |
|
|
|
|
|
|
|
íÿ (2.1.2) çà |
|
|
|
|
, якщо вона ¹ розв'язком рiвнян- |
|||||||||||
|
будь-яких значень сталих |
|
||||||||||||||
початкових умов (2.1.3) iснують ¹динiC1 |
, |
C2 |
, . . . , |
Cn |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
значення сталихi якщо за будь-яких1 |
|||||||
такi, що функцiя |
|
|
|
|
|
|
|
|
|
|
C0 |
C0 |
C0 |
|||
|
|
|
|
|
|
|
|
|
|
1 , |
2 , . . . , |
n |
умови.
задовольня¹ заданi початковi
Загальний розв'язок може бути записаний i в неявнiй формi:
Φ(x, y, C1, . . . , Cn) = 0.
Функцiя y = ϕ(x, C10, C20, . . . , Cn0), отримана iз загального розв'язку y = ϕ(x, C1, C2, . . . , Cn) рiвняння n-го порядку за певних значень сталих
цьогоC |
0 |
|
C |
0 |
, . . . , |
C |
|
C |
0 |
, назива¹ться частинним розв'язком |
рiвнянняC , C . |
2 |
n = |
n |
|||||||
1 = |
1 |
2 = |
|
|
|
|
|
Для побудови розв'язку задачi Кошi для рiвняння n-го порядку загальний розв'язок рiвняння пiдпорядковують початковим умовам:
|
|
y(x0) = ϕ(x, C1, C2, . . . , Cn) x=x0 |
= y0, |
|
|
|
|
|||||||||||||
|
|
|
∂ϕ(x, C1, C2 |
, . . . , Cn) |
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
y0 |
(x0) = |
|
|
|
|
|
|
|
|
|
|
= y0 , |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
∂x |
|
|
|
|
x=x0 |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂ ϕ(x, C1, C2, . . . , Cn) |
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(2.1.5) |
|
|
|
y00(x0) = |
|
2 |
|
|
|
|
|
|
= y00 |
, |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
x=x |
|
0 |
|
|
|
|
|||||
|
|
|
|
|
|
∂x |
|
|
|
|
0 |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . , |
|
|
|||||||||||||||||
|
|
|
|
n 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(n 1) |
|
∂ − |
ϕ(x, C1, C2, . . . , Cn) |
|
|
|
|
|
(n 1) |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y |
− |
(x |
) = |
|
|
|
|
|
|
|
|
|
|
|
= y |
− . |
|
|
||
|
|
n 1 |
|
|
|
|
|
|
|
|
||||||||||
|
|
|
0 |
|
|
|
|
|
|
|
x=x |
|
|
0 |
|
|
||||
|
|
|
|
|
|
∂x − |
|
|
|
|
|
0 |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
рiвнянь вiдносно |
|
íåâi- |
|||||
Рiвностi (2.1.5) являють собою систему |
n |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
домих C , C , . . . , C
ма¹ ¹диний1 2ðîçâ'ÿçîên. За виконання умов твердження 2.1.1 ця система
C10, C20, . . . , Cn0. Розв'язком задачi Кошi ¹ функцiя
y = ϕ(x, C10, C20, . . . , Cn0).
1Тут ми ма¹мо на увазi будь-якi точки (x0, y0, y00 , . . . , y0(n−1)) D, тобто точки, у яких виконуються вимоги твердження 2.1.1.
21
Для видiлення частинного розв'язку рiвняння порядку вище першо-
го додатковi умови можуть бути заданi не в однiй точцi |
x = x0, à íà |
кiнцях деякого промiжку |
[a, b], тобто при x = a i x = b. Таку задачу називають крайовою задачею. Зазначимо, що на вiдмiну вiд задачi Кошi за виконання умов твердження 2.1.1 крайова задача може мати або ¹диний розв'язок, або безлiч розв'язкiв, або взагалi не мати розв'язку.
Перейдемо до розгляду окремих випадкiв диференцiальних рiвнянь, iнтегровних у квадратурах.
2.2.Рiвняння y(n) = f(x)
Iнтегруючи за змiнною x рiвняння
y(n) = f(x), |
(2.2.1) |
отрима¹мо:
Z Z Z
y(n)(x) dx = f(x) dx + C1, àáî y(n−1) = f(x) dx + C1,
äå C
Iнтегруючи1 довiльнаотриманестала. рiвняння, матимемо:
ZZ
y(n−2) = f(x) dx dx + C1x + C2.
Остаточно, пiсля n-кратного iнтегрування знайдемо розв'язок рiвняння (2.2.1):
ZZ
y(x) = · · · f(x) dx . . . dx + C1xn−1 + C2xn−2 + . . . + Cn−1x + Cn,
|{z}
n ðàçiâ
äå C1, C2, . . . , Cn äîâiëüíi ñòàëi.
П р и к л а д 1. Для диференцiального рiвняння y00 = 0 побудувати загальний розв'язок та знайти частиннi розв'язки, якi задовольняють такi додатковi
умови: |
à) y(0) = 1, y0(0) = 2; |
á) y(0) |
= 1, y(1) = 4; |
|
â) y0(0) = 1, y0(1) = 1; |
ã) y0(0) |
= 1, y0(1) = 2. |
Р о з в ' я з а н н я. Буду¹мо загальний розв'язок рiвняння. Iнтегруючи задане рiвняння, отриму¹мо
y0 = C1.
22
Iнтегру¹мо отримане рiвняння ще раз i отриму¹мо загальний розв'язок:
y = C1x + C2.
Розв'язу¹мо задачу а:
y(0) |
= |
C1x + C2 |
x=0 = C2 = 1, |
||||||
|
|
|
|
|
|
|
x=0 |
|
|
0 |
(0) |
= |
1 |
|
|
x0 |
= C1 |
= 2. |
|
y |
(C |
x + C2) |
|
|
|||||
|
|
|
|
|
|
|
|
|
Отже, розв'язком задачi а ¹ функцiя y = 2x + 1. Розв'язу¹мо задачу б :
y(0) |
= |
C1x + C2 |
x=0 = C2 = 1, |
|
||||
|
|
|
|
|
x=1 |
|
|
|
y(1) |
= |
1 |
|
|
= C1 |
+ C2 |
= 4. |
|
(C |
x + C2) |
|
||||||
|
|
|
|
|
|
|
|
|
Знаходимо, що C1 = 3, C2 = 1, i задача б ма¹ ¹диний розв'язок y = 3x + 1. Буду¹мо розв'язок задачi в:
y0(0) = (C1x + C2)0x = C1 = 1.
x=0
y0(1) = (C1x + C2)0x = C1 = 1.
x=1
Як бачимо, крайовi умови y0(0) = 1, y0(1) = 1 задовольняють такi розв'язки
заданого диференцiального рiвняння: y = x+C2, äå C2 довiльна стала, тобто крайова задача в ма¹ безлiч розв'язкiв.
Очевидно, що крайова задача г не ма¹ розв'язку, оскiльки система
y0(0) = (C1x + C2)0x = C1 = 1.
x=0
y0(1) = (C1x + C2)0x = C1 = 2.
x=1
¹ несумiсною.
23
2.3.Рiвняння, що допускають зниження порядку
Вкажемо два окремi випадки, коли порядок рiвняння можна знизити. 1. Рiвняння виду F (x,y (k), . . . ,y (n)) = 0. Уведемо нову змiнну
z = y(k). Тодi порядок рiвняння знизиться на k одиниць:
F (x, z, z0, . . . , z(n−k)) = 0.
Якщо знайдемо загальний розв'язок останнього рiвняння
z = ϕ(x, C1, . . . , Cn−k),
òî y визнача¹ться iз рiвняння
y(k) = ϕ(x, C1, . . . , Cn−k).
Рiвняння такого типу розглянутi в попередньому параграфi.
2. Рiвняння виду F (y,y 0, . . . ,y (n)) = 0. Якщо рiвняння не мiстить незалежно¨ змiнно¨ x, то вiзьмемо за незалежну змiнну y i покладемо
y0 = p(y),
äå p(y) невiдома функцiя. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Òîäi |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y00 |
= |
dy0 |
= |
dp(y) |
= |
|
dp |
p, |
|
|
|
|
|
|
|||||||
dx |
|
dx |
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
dy |
|
|
|
|
|
|
|||||||
|
|
d dp |
|
|
d dp |
d2p |
|
dp |
|
2 |
|||||||||||
y000 |
= |
|
|
|
|
p = |
|
|
|
|
p p = |
|
p2 |
+ |
|
|
p, |
||||
dx |
dy |
|
dy |
dy |
dy2 |
dy |
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Звiдси видно, що порядок рiвняння знизиться на одиницю.
Якщо буде знайдено розв'язок рiвняння, отриманого в результатi замiни,
p = ϕ(y, C1, . . . , Cn−1),
òî íåâiäîìà y ¹ розв'язком рiвняння першого порядку iз вiдокремлюваними змiнними
y0 = ϕ(y, C1, . . . , Cn−1).
24
П р и к л а д 1. Знайти частинний розв'язок рiвняння y00 + xy0 = x за додаткових умов y(0) = 2, y0(0) = 1.
Р о з в ' я з а н н я. Рiвняння не мiстить змiнно¨ y. Вводимо нову невiдому
çìiííó z = y0. Враховуючи, що y00 = z0, перепишемо рiвняння в нових позна- ченнях:
Отримане рiвняння ¹ рiвнянням першого порядку з вiдокремлюваними змiн- |
|||||||||||||
íèìè: |
|
|
|
|
|
dz |
|
|
|
|
|
||
|
|
|
|
|
|
|
= x(1 − z). |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dx |
|
|
|
|
||
Буду¹мо його загальний розв'язок: |
|
|
|
|
|
||||||||
|
dz |
|
dz |
|
|
|
|
x2 |
2 |
|
|||
|
|
= x dx; |
Z |
|
= Z |
x dx; |
ln(1 − z) = − |
|
+ ln C1; |
z = 1 − C1e−x |
/2. |
||
|
1 − z |
1 − z |
2 |
Враху¹мо додаткову умову y0(0) = 1, àáî z(0) = 1:
z(0) = 1 − C1 = 1.
Звiдси знаходимо, що C1 = 0, à z(x) = 1.
Îñêiëüêè z = y0, то для визначення невiдомо¨ y ма¹мо таке рiвняння:
y0 = 1.
Загальним розв'язком отриманого рiвняння ¹ функцiя y = x + C2. Значення довiльно¨ стало¨ C2 знаходимо з умови y(0) = 2:
y(0) = 0 + C2 = 2; C2 = 2.
Отже, розв'язком задачi ¹ функцiя
П р и к л а д 2. Знайти загальний розв'язок рiвняння y000 − y00 = 1.
Р о з в ' я з а н н я. Рiвняння не мiстить y i y0. Çàìiíèìî y00 = z. У результатi рiвняння набуде вигляду:
Буду¹мо загальний розв'язок отриманого рiвняння:
z0 − z = 1; |
dx = z + 1; |
z + 1 = dx; |
Z |
z + 1 = Z |
dx; ln(z + 1) = x + ln C1; |
|
dz |
dz |
|
dz |
|
z = C1ex − 1.
Враховуючи, що z = y00, запишемо рiвняння для визначення y:
y00 = C1ex − 1.
Iнтегруючи двiчi останн¹ рiвняння, знаходимо y:
y0 = C1ex − x + C2; y = C1ex − x2 + C2x + C3.
2
25
Отже, задане рiвняння ма¹ загальний розв'язок
y = C1ex + C2x + C3 − x2 .
2
П р и к л а д 3. Знайти розв'язок диференцiального рiвняння yy00 − y02 = 0 за додаткових умов y(0) = 1, y0(0) = 2.
Р о з в ' я з а н н я. Рiвняння не мiстить незалежно¨ змiнно¨ x. Çàìiíèìî
y0 = p(y), y00 = p ddyp.
Перепишемо задане рiвняння в нових позначеннях:
yp dy |
− p2 = 0, àáî p |
y dy − p |
= 0. |
||
|
dp |
|
|
dp |
|
Останн¹ рiвняння розпада¹ться на два рiвняння
p = 0 òà y |
dp |
− p = 0. |
dy |
Розв'язком першого з них y = const нехту¹мо, оскiльки вiн не задовольня¹
початкову умову y0(0) = 2.
Буду¹мо розв'язок другого рiвняння:
y |
dp |
− p = 0; |
dp |
= |
p |
; |
dp |
= |
dy |
; ln p = ln y + ln C1; p = C1y. |
dy |
dy |
y |
p |
y |
Îñêiëüêè p = y0, то для визначення y ма¹мо таке рiвняння:
y0 = C1y.
Знаходимо його розв'язок:
dy |
= C1y; |
dy |
= C1 dx; ln y = C1x + ln C2; y = C2eC1x. |
|
dx |
y |
|||
|
|
Значення сталих C1 i C2 знайдемо з початкових умов:
|
|
y(0) |
= |
C2eC1x x=0 = C2 = 1, |
|||
|
|
y0(0) |
= |
|
C1x |
x=0 = C2C1 = 2. |
|
|
|
C2C1e |
|||||
Знаходимо, що |
C1 |
, |
= 1, à |
розв'язком |
задано¨ задачi Кошi ¹ функцiя |
||
|
= 2 C2 |
|
|
|
y = e2x.
26
 ï ð à â è
Розв'язати рiвняння:
17. y00 = xex. Початковi умови y(0) = 0, y0(0) = 1.
18. yIV = sin2 x. Початковi умови y(0) = 1, y0(0) = y00(0) = y000(0) = 0.
19. xy00 + y0 = 0. Початковi умови y(1) = 2, y0(1) = 3. 20. y02 + yy00 = 1. Крайовi умови y(0) = 1, y(1) = 2.
21. (x − 1)y000 − y00 = 0. Початковi умови y(2) = 0, y0(2) = 1, y00(2) = 0.
22. y00 = y0ey. Початковi умови y(0) = 0, y0(0) = 1. 23. y3y00 + 1 = 0. Початковi умови y(1) = 1, y0(1) = 0.
24. (1+sin x)y000 = cos x y00. Початковi умови y(0) = 0, y0(0) = 1, y00(0) = 0.
25. |
yIV + y000 |
= x. |
|||
26. |
(x + 1)y000 |
+ y00 = (x + 1). Початковi умови y(0) = y0(0) = y00(0) = 0. |
|||
27. |
y000 |
tg x = y00. |
|||
28. |
y00 |
+ 8 sin y cos3 y = 0. Початковi умови y(0) = 0, y0(0) = 2. |
|||
29. |
y000 |
tg x = y00 + 1. |
|||
|
xy000 + y00 |
= p |
|
. |
|
30. |
x |
27
Òåìà 3
Лiнiйнi рiвняння вищих порядкiв
3.1.Основнi поняття
Означення 13. Лiнiйним диференцiальним рiвнянням n-го порядку називають рiвняння виду
y(n) + a1(x)y(n−1) + . . . + an−1(x)y0 + an(x)y = f(x), |
(3.1.1) |
де функцi¨ a1(x), . . . , an−1(x), an(x) i f(x) визначенi й неперервнi на деякому промiжку a < x < b.
ßêùî f(x) ≡ 0, рiвняння (3.1.1) називають о д н о р i д н и м, якщо f(x) ≡/ 0 í å î ä í î ð i ä í è ì.
Позначимо лiву частину рiвняння (3.1.1) через L[y]:
L[y] ≡ y(n) + a1(x)y(n−1) + . . . + an−1(x)y0 + an(x)y |
(3.1.2) |
i назвемо L[y] диференцiальним оператором.
Безпосередньою перевiркою можна переконатися, що оператор L[y] ма¹ двi такi важливi властивостi:
òîäi1) нехай y1 i y2 будь-якi функцi¨, що мають n неперервних похiдних,
L[y1 + y2] = L[y1] + L[y2]; |
(3.1.3) |
2) нехай y n разiв диференцiйовна функцiя, а α стала, тодi
L[αy] = αL[y]. |
(3.1.4) |
28
Розглянемо деякi властивостi лiнiйних однорiдних рiвнянь
y(n) + a1(x)y(n−1) + . . . + an−1(x)y0 + an(x)y = 0, |
(3.1.5) |
або у позначеннях (3.1.2)
L[y] = 0.
Очевидно, що функцiя y(x) = 0 ¹ розв'язком однорiдного рiвняння (3.1.5).
Означення 14. Ðîçâ'ÿçîê y = 0 лiнiйного однорiдного рiвняння називають т р и в i а л ь н и м р о з в ' я з к о м.
Нас цiкавитимуть розв'язки однорiдного рiвняння (3.1.5), в i д м i н - н i вiд тривiального, i далi мова йтиме саме про такi розв'язки.
Iз властивостей (3.1.3) i (3.1.4) оператора L випливають два такi твердження:
Твердження 3.1. 1. ßêùî y(x) який-небудь розв'язок однорiдно-
го рiвняння (3.1.5), то функцiя C · y(x), äå C довiльна стала, теж ¹ його розв'язком.
Твердження 3.1. 2. ßêùî y1 i y2 розв'язки рiвняння (3.1.5), то
íÿ (3.1.5). |
, äå |
|
, |
|
довiльнi сталi, також ¹ розв'язком рiвнян- |
|
C1y1 + C2y2 |
|
C1 |
|
C2 |
|
|
Í à ñ ë i ä î ê. ßêùî |
y1, y2, . . . , yn частиннi розв'язки рiвняння (3.1.5), |
|||||
то вираз |
|
|
|
|
||
|
|
|
|
|
C1y1 + C2y2 + . . . + Cnyn |
(3.1.6) |
також ¹ його розв'язком (тут C1, C2, . . . ,
Вираз (3.1.6) ¹ розв'язком рiвняння n-го порядку i мiстить n довiльних сталих. Виника¹ питання: чи не ¹ вираз (3.1.6) загальним розв'язком рiвняння (3.1.5)? Вiдповiдь на це питання буде ствердною, якщо частиннi
ðîçâ'ÿçêè y1, y2, . . . , yn ¹ ë i í i é í î í å ç à ë å æ í è ì è.
Означення 15. Система функцiй {y1, y2, . . . , yn} назива¹ться лiнiйно незалежною на промiжку (a, b), ÿêùî ðiâíiñòü
C1y1 + C2y2 + . . . + Cnyn = 0 |
(3.1.7) |
викону¹ться тотожно на (a, b) òîäi é òiëüêè òîäi, êîëè
C1 = C2 = . . . = Cn = 0.
29
Означення 16. Систему n лiнiйно незалежних розв'язкiв рiвняння |
|
|
(3.1.5) називають ф у н д а м е н т а л ь н о ю системою розв'язкiв цього |
|
|
рiвняння. |
|
|
Твердження 3.1. 3. Для будь-якого лiнiйного однорiдного рiвняння |
|
|
iсну¹ фундаментальна система розв'язкiв. |
|
- |
Твердження 3.1. 4. ßêùî (3.1.5), то йогоутворюютьзагальнийфундаментальрозв'язок |
||
ну систему розв'язкiв рiвняння y1, y2, . . . , yn |
|
|
визнача¹ться формулою |
y(x) |
|
|
|
|
y = C1y1 + C2y2 + . . . + Cnyn. |
(3.1.8) |
Останн¹ твердження опису¹ структуру загального розв'язку однорiдного рiвняння (3.1.5). Вiдповiдь на питання щодо структури розв'язку лiнiйного н е о д н о р i д н о г о рiвняння (3.1.1) да¹ таке твердження:
Твердження 3.1. 5. Загальний розв'язок y(x) неоднорiдного рiвня-
ння (3.1.1) явля¹ собою суму загального розв'язку y(x) вiдповiдного одно-
рiдного рiвняння (3.1.5) i якого-небудь частинного розв'язку y (x) неоднорiдного рiвняння (3.1.1), тобтî
y = y + y . |
(3.1.9) |
Перейдемо до побудови розв'язкiв лiнiйних диференцiальних рiвнянь. При цьому ми обмежимось окремим випадком рiвняннями зi сталими
коефiцi¹нтами, тобто випадком, коли a1, a2, . . . , an
3.2.Однорiднi рiвняння зi сталими коефiцi¹нтами
Розглянемо лiнiйне однорiдне рiвняння
y(n) + a1y(n−1) + . . . + an−1y0 + any = 0, |
(3.2.1) |
äå a1, . . . , an−1, an дiйснi числа. Його розв'язок шукатимемо у виглядi
y = eλx, äå λ поки що невiдоме число. Пiсля пiдстановки цi¹¨ функцi¨ у рiвняння (3.2.1) отрима¹мо таку рiвнiсть:
eλx λn + a1λn−1 + . . . + an−1λ + an = 0.
Îñêiëüêè eλx 6= 0, то з останньо¨ рiвностi ма¹мо такий висновок: функцiя y = eλx буде розв'язком рiвняння (3.2.1), якщо λ ¹ коренем рiвняння
λn + a1λn−1 + . . . + an−1λ + an = 0. |
(3.2.2) |
30