Математика для инженеров(теория)I том
.pdfПосле того, как число k найдено, необходимо сделать замену переменных x = et , y = zekt , где z = z(t) – новая неизвестная функция,
а t – новая независимая переменная. Получаем уравнение, в которое не входит независимая переменная x . Порядок такого уравнения понижается одним из ранее рассмотренных способов (см.п. 3).
Пример 8. Решить уравнение x4 y¢¢ + (xy¢ - y)3 = 0.
Решение. Покажем сначала, что это уравнение
обобщенное однородное. Найдем число k. Будем считать x, y, y′, y′′ величинами соответственно 1,k,(k −1) , (k − 2) -ой
степени и, после сравнения степеней всех слагаемых в
левой части данного уравнения, получим |
4 + k − 2 = 3k , |
|||||||||||||
откуда k =1. Теперь сделаем подстановку |
|
|
x = et |
, y = zet . |
||||||||||
|
dy |
|
dz |
|
|
dy¢ |
æ |
d |
2 |
z |
|
dz |
ö |
|
Так как y¢ = |
×e−t = |
+ z, |
y¢¢ = |
e−t = ç |
|
|
+ |
÷e−t , исход- |
||||||
|
|
|
|
|
2 |
|
||||||||
|
dt |
|
dt |
|
dt |
ç |
|
|
|
dt |
÷ |
|||
|
|
|
è dt |
|
|
|
ø |
|||||||
æ |
d |
2 |
z |
|
dz |
ö |
|
ное уравнение примет вид e4t ç |
|
|
+ |
÷e |
|||
|
|
2 |
|
||||
ç |
dt |
|
dt |
÷ |
|||
è |
|
|
|
ø |
|||
или
−t |
æ |
t æ dz |
ö |
|
t ö3 |
||
|
+ çe |
ç |
|
+ z ÷ |
- ze |
÷ |
= 0 |
|
|
||||||
|
è |
è dt |
ø |
|
ø |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d 2 z |
|
|
dz |
|
æ dz ö3 |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+ |
|
+ ç |
|
÷ = 0 . |
|
|
|
|
(3) |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
dt2 |
|
|
dt |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
è dt ø |
|
|
|
|
|
|||||||||
Допустим теперь |
dz |
|
= u и примем z за независимую переменную. |
|||||||||||||||||||||||||||||
dt |
||||||||||||||||||||||||||||||||
|
d2 z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Тогда |
= |
|
du |
|
= |
|
du dz |
= |
|
du |
u . После подстановки в (3) будем |
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
dt2 |
|
dt |
|
|
dz dt |
|
dz |
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
иметь |
du |
u + u + u3 = 0 или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
dz |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
du |
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+1+ u2 = 0 . |
|
|
|
|
(4) |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dz |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
получим u = tg(C1 - z) . |
|
|||||||||
Интегрируя |
уравнение |
(4), |
За- |
|||||||||||||||||||||||||||||
менив u на |
dz |
будем иметь |
|
dz |
= tg(C - z) или |
|
|
dz |
= dt . |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
1 |
- tg(z - C1) |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
Проинтегровав |
|
еще раз, получим: |
ln | sin(z - C1) | +t = ln | C2 | . |
Воз- |
||||||||||||||||||||||||||||
вращаясь к переменным x и y по формулам t = ln x, z = |
y |
, находим |
||||||||||||||||||||||||||||||
x |
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
общий интеграл исходного уравнения в виде
533
|
|
ln |
sin |
æ y |
- C |
ö |
+ ln x = ln | C |
|
| |
|||
|
|
ç |
|
÷ |
2 |
|||||||
|
|
|
|
|
|
x |
1 |
|
|
|||
|
|
|
|
|
|
è |
|
ø |
|
|
|
|
|
æ |
|
|
B |
ö |
(A = C1, B = ±C2 ) . □ |
|
|||||
или |
y = xç A + arcsin |
|
|
÷ |
|
|||||||
x |
|
|||||||||||
|
è |
|
|
ø |
|
|
|
|
|
|
||
В некоторых случаях найти решение в виде явной или неявной
функции непросто, однако удается получить решение в параметриче-
ской форме.
Пример 9. Решить уравнение |
¢¢ |
¢ |
y′ |
=1 . |
|
|
||||||
y (1 |
+ y )e |
|
|
|
||||||||
Решение. |
Вводя |
параметр t |
по |
формуле |
y′ = t , |
из |
||||||
уравнения получаем y¢¢ |
|
e−t |
|
|
|
|
|
|
||||
= |
|
. Отсюда, |
в силу соотношения |
|||||||||
1+ t |
||||||||||||
d( y′) = y′′dx , |
т.е. |
dt = y′′dx , имеем |
(1+ t)dt = e−t dx . |
Следова- |
||||||||
тельно, |
x = tet + C . |
|
|
|
|
|
|
|
|
|
||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
Из уравнения y′ = t |
находим |
|
|
|
|
|
|
|||||
|
|
y = òtdx + C2 = òt(1+ t)et dt + C2 = (t2 - t +1)et |
+ C2 . |
|
||||||||
Итак, параметрические уравнения общего решения |
||||||||||||
и |
м |
е |
ю |
|
|
т |
|
|
в |
|
и |
д |
x = tet + C1 , y = (t2 - t +1)et + C2 . □
§ 2. Линейные дифференциальные уравнения высших порядков
10. Основные определения. Линейным дифференциальным уравнением n-го порядка называют уравнение
|
y(n) + pn−1(x)y(n−1) + ... + p1(x)y¢ + p0 (x)y = f (x) , |
(1) |
|||
где |
функции |
pi (x),i = |
|
|
|
0,n -1 и f (x) непрерывны на не- |
|||||
котором |
интервале |
I = (a;b) . |
|
||
Линейным дифференциальным оператором n-го порядка называют
выражение
L |
= |
d n |
+ p |
(x) |
d n−1 |
+ ... + p (x) |
d |
+ p (x) . |
(2) |
|
|
|
|||||||
n |
|
dxn |
n−1 |
|
dxn−1 |
1 |
dx |
0 |
|
|
|
|
|
|
|
|
534
Тогда линейное дифференциальное уравнение, с ис- пользованием обозначения (2), примет вид
|
|
|
|
Ln[ y] = f (x) . |
|
|
|
(1’) |
||
|
Нетрудно убедиться, что оператор Ln удовлетворяет |
|||||||||
у |
с |
л |
о |
в |
|
и |
я |
|
м |
: |
|
а) однородности: |
Ln[cy] = cLn [ y] ; |
|
|
|
|
||||
|
б) аддитивности: Ln[ y1 + y2 ] = Ln[y1] + Ln[ y2 ]. |
|
|
|
||||||
|
Уравнение (1) называют линейным неоднородным, если |
f (x) ¹ 0 . |
||||||||
Если |
f (x) º 0 , т.е. уравнение имеет вид |
|
|
|
||||||
|
|
|
|
Ln[ y] = 0 , |
|
|
|
|
(3) |
|
его называют линейным однородным. |
|
|
|
|
|
|||||
|
Задача Коши формулируется так: найти решение |
|||||||||
уравнения (1), удовлетворяющее условиям |
|
|
|
|||||||
|
|
|
¢ |
¢ |
,..., y |
(n−1) |
(n−1) |
, |
|
(4) |
|
|
y(x0 ) = y0 , y (x0 ) = y0 |
|
(x0 ) = y0 |
|
|||||
называемыми начальными условиями. |
|
|
|
|
|
|||||
|
Общим решением уравнения (1) называется функция |
|||||||||
y = ϕ(x,C1,...,Cn ) , |
зависящая от n произвольных постоянных |
|||||||||
и удовлетворяющая условиям: 1) при любых постоянных C1,...,Cn эта
функция является решением (1); 2) каковы бы ни были на- чальные условия (4), можно подобрать такие значения
C10 ,...,Cn0 постоянных C1,...,Cn , что решение y = ϕ (x,C10 ,...,Cn0 )
уравнения (1) будет удовлетворять начальным условиям (4).
20. Свойства линейных однородных уравнений. Пусть
n
y1, y2 ,..., yn – решения уравнения (3). Тогда y = åCk yk также решение
k=1
уравнения (3) при произвольных постоянных Ck ,k =1,...,n .
Действительно, |
по условию Ln[yk ] = 0 . Имеем: |
|
|
||||||||||||||
L [ y] = L |
é |
n |
|
ù |
= |
n |
|
] = |
n |
|
|
] = |
n |
C |
|
×0 = 0 . |
|
ê |
C y |
|
ú |
L [C y |
|
C L [ y |
|
å |
|
||||||||
n |
n |
ê |
å k |
k |
ú |
|
å n k |
k |
|
å k n |
k |
|
|
k |
|
||
|
|
ëk=1 |
|
û |
|
k=1 |
|
|
k=1 |
|
|
|
k=1 |
|
|
|
|
Систему функций |
y1(x), y2 (x),..., yn (x) |
называют линейно |
|||||||||||||||
зависимой |
на |
|
I , если |
существуют |
постоянные числа |
||||||||||||
α1,α2 ,...,αn , хотя бы одно из которых отлично от нуля, что выполняется:
535
α1y1(x) +α2 y2 (x) + ... +αn yn (x) = 0 x (a;b) . |
(5) |
Если равенство (5) имеет место только |
при |
α1 = α2 = ... = αn = 0 , то система функций y1(x), y2 (x),..., yn (x)
называется линейно независимой на I.
Пример 1. Показать, что система функций 1, x, x2 ли- нейно независима на любом интервале I.
Решение. Действительно, равенство (5) в данном случае имеет вид
α +α |
2 |
x +α |
3 |
x2 = 0 x (a;b) , |
1 |
|
|
||
что эквивалентно α1 = α2 = α3 = 0 , т.к. многочлен обра- |
||||
щается в тождественный нуль на любом интервале I тогда |
||||
и только тогда, когда все его коэффициенты равны нулю
(п.10.8.20 ). □
Рассмотрим систему n |
функций y1(x), y2 (x),..., yn (x) не- |
||||||
прерывных вместе со своими производными |
до (n −1) -го |
||||||
порядка включительно на I . |
|
|
|
|
|
||
Определителем |
Вронского |
или |
вронскианом |
||||
W [ y1, y2 ,..., yn ] |
этой системы функций называется определи- |
||||||
тель |
|
|
y1 |
y2 ... |
yn |
|
|
|
|
|
|
||||
|
|
|
y′ |
y′ ... |
y′ |
|
|
|
W [y1, y2 ,..., yn ] = |
1 |
2 |
n |
. |
||
|
K |
K ... |
K |
||||
|
|
|
y(n−1) |
y(n−1) ... |
y(n−1) |
|
|
|
|
|
1 |
2 |
n |
|
|
Утверждение 1. Если функции y1, y2 ,..., yn |
линейно за- |
||||||
висимые на |
интервале |
I , |
то определитель |
Вронского |
|||
W [y1, y2 ,..., yn ] |
равен нулю на этом интервале. |
|
|
|
|||
Упражнение 1. Докажите утверждение 1.
Установим некоторые основные свойства линейных однородных уравнений, ограничиваясь при доказательствах уравне- ниями второго порядка.
Утверждение 2. Если определитель Вронского, со- ставленный для решений y1 и y2 дифференциального
уравнения y′′ + py′ + qy = 0 , не равен нулю при x = x0 I , то он не обращается в нуль ни при одном значении x из этого интервала.
Доказательство. По условию y1 и y2 – решения ука- занного уравнения, значит
536
ì |
y¢¢+ py¢ |
+ qy = 0, |
|
í |
1 |
1 |
1 |
¢¢ |
¢ |
+ qy2 = 0. |
|
îy2 |
+ py2 |
||
|
Умножая первое уравнение на y2 , а второе – на y1 , |
|||||||||||||||||||||||||||||||||||||
п |
|
|
о |
|
|
л |
|
|
|
|
|
|
|
у |
|
|
|
|
|
|
|
ч |
|
|
|
|
|
|
|
|
|
и |
|
|
м |
|||
|
|
|
|
|
|
|
ìy¢¢y |
2 |
+ py¢y |
2 |
+ qy y |
2 |
= 0, |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
í |
1 |
|
|
|
1 |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
¢¢ |
|
|
|
|
¢ |
|
|
+ qy2 y1 = 0. |
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
î y2 y1 |
+ py2 y1 |
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
Вычитая из второго уравнения системы первое, будем |
|||||||||||||||||||||||||||||||||||||
и |
|
|
|
м |
|
|
|
|
|
|
|
|
|
е |
|
|
|
|
|
|
|
|
|
|
|
|
т |
|
|
|
|
|
ь |
|||||
|
|
|
|
|
|
y y′′ - y′′y |
2 |
+ p(y′ |
y - y′y |
2 |
) = 0 . |
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
1 |
2 |
|
|
|
1 |
|
|
|
|
2 |
1 |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
В |
|
силу |
того, |
|
|
что |
|
y y′ |
- y′y |
2 |
= W[y , y |
2 |
] = W (x) , |
|
а |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
2 |
1 |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
||||||
y y′′ - y |
2 |
y′′ = W ′(x) , |
запишем |
|
W ′ + pW = 0 |
|
|
или |
W ′ |
= - p . |
Интег- |
|||||||||||||||||||||||||||
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
1 |
2 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
W |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
рируя |
последнее |
|
уравнение, |
получим |
|
|
ln |
|
W |
|
= - ò pdx + ln |
|
C |
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x0 |
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
− ò pdx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
или W = Ce x0 |
. При x = x |
|
имеем: |
W (x ) = C = W . |
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
− ò pdx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
Формулу вида W =W e x0 |
|
|
называют формулой Лиу- |
||||||||||||||||||||||||||||||||||
вилля. |
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Из этой формулы и следует утверждение 2. □ |
|
|
|
|
|
|||||||||||||||||||||||||||||||||
|
Для уравнения (3) имеет место формула Остроградского-Лиувилля: |
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
− ò |
pn−1(t)dt |
|
|
|
|
|
|||||
|
|
|
W [y , y |
|
,..., y |
n |
] = W [y , y |
|
,..., y |
n |
] |
|
|
|
×e x0 |
|
|
. |
|
|
|
|
|
|||||||||||||||
|
|
|
|
1 |
2 |
|
|
|
|
|
|
|
1 |
|
2 |
|
|
|
|
x |
=x0 |
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
Утверждение 3. |
|
Если y1 |
и |
|
y2 – решения уравнения |
||||||||||||||||||||||||||||||||
y′′ + py′ + qy = 0 , |
линейно |
|
независимые |
|
|
на интервале |
I , |
|
то |
|||||||||||||||||||||||||||||
определитель Вронского W , составленный из этих реше- ний, не обращается в нуль ни в одной точке этого интер- вала.
Доказательство. Допустим, что W[ y1, y2 ] = 0 в некото- рой точке x0 Î I . Тогда, по утверждению 1, определитель Вронского W[ y1, y2 ] будет равен нулю во всех точках ин-
тервала |
I , |
т.е. y y′ |
- y′y |
2 |
= 0 . Допустим, что |
y ¹ 0 |
на I . |
||
|
|
1 |
2 |
|
1 |
|
1 |
|
|
Тогда, |
на |
основании |
последнего равенства, |
можно |
напи- |
||||
537
|
y y′ |
- y′y |
2 |
|
сать |
1 2 |
1 |
= 0 или |
|
|
y2 |
|
||
|
|
|
|
|
|
|
1 |
|
|
y2 = C = const , т.е. решения y1
æ |
|
y2 |
ö¢ |
= 0 . Отсюда следует |
|
ç |
÷ |
||||
y |
|||||
è |
|
ø |
|
||
1 |
|
||||
y1 |
и |
y2 |
линейно зависимы, что |
||
противоречит предположению об их линейной независи- мости. □
Упражнение 2. Проверить утверждение 3 в случае y1 = 0.
30. Структура общего решения линейного однородного уравнения. Систему функций y1, y2 ,..., yn , являющихся ли-
нейно независимыми решениями уравнения (3), называют
фундаментальной системой решений этого уравнения. Для такой системы вронскиан W (x) =W [y1(x), y2 (x),..., yn (x)] ¹ 0 .
Теорема 1. Если y1, y2 ,..., yn – фундаментальная система реше-
ний однородного дифференциального уравнения (3), то его
общее решение имеет вид
y = C1 y1 + C2 y2 + ... + Cn yn , |
(6) |
где Ci , i =1,n , – произвольные постоянные.
Доказательство. Выражение (6) является решением уравнения (3) на основании свойства линейного диффе-
р е н ц и а л ь н о г о |
|
о п е р а т о р а . |
|||||||||||||||
Докажем, что если |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
y = y0 , |
y¢ = y0¢,..., y(n−1) = y0(n−1) |
при |
x = x0 , |
(7) |
||||||||||||
то C1,C2 ,...,Cn можно подобрать таким образом, |
что |
||||||||||||||||
(6) будет |
удовлетворять |
условиям |
(7). Подставляя в |
(6) |
|||||||||||||
x = x0 и обозначая yi (x0 ) = yi0 ,i = |
1, n |
, |
получим систему |
ли- |
|||||||||||||
нейных алгебраических уравнений |
|
|
|
|
|
|
|
|
|||||||||
|
ìC y |
+ C y |
20 |
+ ... + C y |
n0 |
= y , |
|
|
|
|
|||||||
|
ï 1 |
10 |
|
2 |
|
|
n |
|
|
0 |
|
|
|
|
|
||
|
ïC y¢ |
+ C y¢ |
+ ... + C y¢ |
= y¢ |
, |
|
|
|
|
||||||||
|
ï 1 |
10 |
|
2 |
20 |
|
|
n |
n0 |
|
|
0 |
|
|
|
|
|
|
|
¢¢ |
|
|
¢¢ |
|
|
|
¢¢ |
|
|
¢¢ |
, |
|
|
|
(8) |
|
íC1 y10 |
+ C2 y20 |
+ ... + Cn yn0 |
= y0 |
|
|
|
||||||||||
|
ï............................................. , |
|
|
|
|
|
|
||||||||||
|
ï |
(n−1) |
|
|
(n−1) |
|
|
|
|
(n−1) |
|
(n−1) |
|
|
|||
|
ï |
|
|
|
|
|
|
|
. |
|
|||||||
|
îC1 y10 |
|
+ C2 y20 |
|
+ ... + Cn yn0 |
|
= y0 |
|
|||||||||
Так как главным определителем системы (8) является определитель |
|||||||||||||||||
Вронского |
при |
x = x0 |
и |
W (x0 ) ¹ 0 , |
то эта |
система имеет |
|||||||||||
538
единственное решение относительно неизвестных C1,C2 ,...,Cn , при котором функция (6) удовлетворяет условиям (7). Итак, дока- зано, что если y1, y2 ,..., yn – система линейно независимых
решений уравнения (3), то функция (6) является решением этого уравнения, и любое решение этого уравнения можно получить из формулы (6) соответствующим выбором по- стоянных C1,C2 ,...,Cn . Таким образом равенство (6) опреде-
ляет общее решение линейного однородного уравнения
(3). □
Сформулированные выше свойства позволяют решить следующую задачу: по данной системе n линейно незави- симых функций y1(x), y2 (x),..., yn (x), непрерывных вместе со
своими производными до n-го порядка включительно на (a;b) , построить линейное однородное дифференциальное
уравнение n-го порядка, решениями которого являются данные функции. Искомым дифференциальным уравнени-
ем является следующее равенство: |
|
|
|
||
|
y |
y1 |
... |
yn |
|
|
|
||||
|
y′ |
y′ |
... |
y′ |
|
W[ y, y1,..., yn ] = |
|
1 |
|
n |
= 0 . |
........ ........... |
... |
........... |
|||
|
y(n) |
y(n) |
... |
y(n) |
|
|
|
1 |
|
n |
|
§ 3. Линейные однородные дифференциальные уравнения
n-го порядка с постоянными коэффициентами
Линейным однородным дифференциальным уравнением с посто-
янными коэффициентами называется уравнение
y(n) + a y(n−1) |
+ a y(n−2) |
+ ... + a |
n−1 |
y′ + a |
n |
y = 0 , |
(1) |
1 |
2 |
|
|
|
|
где ai , i =1, n , – действительные числа.
Будем искать частные решения этого уравнения в ви-
де |
|
y = eλx , |
(2) |
где λ – некоторое число.
539
Подставляя функцию (2) и ее производные в (1), по-
лучим
|
|
λneλx + λn−1a eλx + λn−2a |
eλx + ... + λa eλx + a eλx = 0 |
|
|||||
|
|
|
1 |
2 |
|
|
n−1 |
n |
|
|
или eλx (λn + a λn−1 + ... + a |
λ + a ) = 0 . |
|
|
|||||
|
|
|
1 |
n−1 |
n |
|
|
|
|
|
Поскольку eλx ¹ 0 , то (2) |
будет решением уравнения |
|||||||
( |
1 |
) |
, |
|
е |
|
с |
л |
и |
|
|
|
λn + a λn−1 + ... + a λ + a |
n |
= 0 , |
|
(3) |
||
|
|
|
1 |
|
n−1 |
|
|
|
|
т.е., если λ будет корнем уравнения (3). Уравнение
(3)называется характеристическим.
Таким образом, решение однородного дифференциального уравнения свелось к решению алгебраического уравнения. Этот метод решения называется методом Эйлера.
Пусть λ1,λ2 ,...,λn – корни уравнения (3), причем среди них могут быть и кратные. Возможны следующие случаи.
1. λ1,λ2 ,...,λn – вещественные и различные. Тогда фун- даментальная система решений уравнения (1) имеет вид
eλ1x ,eλ2x ,..., eλnx , |
(4) |
||
а общим решением этого уравнения будет |
|||
y = C eλ1x + C |
eλ2x + ... + C |
eλnx , |
|
1 |
2 |
n |
|
где |
C1,C2 ,...,Cn – произвольные постоянные. |
П р и м е р 1 . Н а й т и о б щ е е р е ш е н и е у р а в н е н и я |
|
y′′′ - 2y′′ - 3y′ = 0 |
|
|
|
. |
Решение. Составляем характеристическое |
уравнение |
|||
λ3 - 2λ2 - 3λ = 0 . Находим его корни: λ |
= 0,λ |
= -1,λ |
|
= 3 . Так |
1 |
2 |
3 |
|
|
как они действительные и различные, то общее решение
имеет вид
y= C1 + C2e−x + C3e3x . □
2.Корни характеристического уравнения вещественные, но среди
них есть кратные. Пусть, например, λ1 = λ2 = ... = λk , т.е. λ1 является k-кратным корнем уравнения (3), а все остальные
(n − k) корней различные. |
Фундаментальная система реше- |
|||||||
ний уравнения (1) в этом случае имеет вид |
|
|
||||||
eλ1x , xeλ1x ,..., xk−1eλ1x ,eλk +1x ,...,eλnx , |
|
(5) |
||||||
а общее решение |
|
|
|
eλk +1x + ... + C |
eλnx . |
|||
y = C eλ1x + C |
2 |
xeλ1x + ... + C |
k |
xk−1eλ1x + C |
||||
1 |
|
|
|
k+1 |
n |
|
||
540
Упражнение 1. Показать, что корню λ1 кратности k будет соответствовать линейно независимая система пер-
вых k решений из (5). |
|
|
|
|
|
|
|
Пример 2. |
Найти общее решение уравнения y′′′ + 2y′′ + y′ = 0 . |
||||||
Решение. |
Характеристическое |
уравнение имеет вид |
|||||
λ3 + 2λ2 + λ = 0 . |
Отсюда λ |
= λ |
= −1, λ |
|
= 0 . |
Корни |
действи- |
|
1 |
2 |
3 |
|
λ = −1 – |
|
|
тельные, причем один из них, |
а именно |
двукрат- |
|||||
ный, поэтому общее решение имеет вид: |
|
|
|||||
y= C1e−x + C2 xe−x + C3 . □
3.Среди корней характеристического уравнения (3)
есть комплексные. Пусть для определенности |
λ1 = α + iβ , |
||
λ2 = α − iβ , λ3 =ν + iδ , λ4 =ν − iδ , а остальные корни веще- |
|||
ственные и различные. Поскольку коэффициенты |
ai , i = |
|
, |
1,n |
|||
уравнения (3) вещественные, то комплексные корни этого уравнения попарно сопряженные. Согласно (2), будем
иметь
y |
= eλ1x |
= e(α +iβ )x = eα x (cos β x + isin β x) , |
|
|
|||||||
1 |
|
|
|
|
|
|
|
|
|
|
|
y2 = eα x (cos β x − isin β x) , |
|
|
|
|
|
|
|
||||
y |
= eν x (cosδ x + isinδ x) , |
|
|
y |
4 |
= eν x (cosδ x − isinδ x) , |
|||||
3 |
|
|
|
|
|
|
|
|
|
|
|
y = eλ5x |
, ..., y |
n |
= eλnx. |
|
|
|
|
|
|
|
|
5 |
|
|
|
|
|
|
|
|
|
|
|
Фундаментальная система решений в этом случае за- |
|||||||||||
пишется: |
|
|
|
|
|
eν x sinδ x , eλ5x , ..., eλnx , (6) |
|||||
eα x cos β x , |
eα x sin β x , |
eν x cosδ x , |
|||||||||
а общее решение |
|
|
|
|
|
|
|
||||
|
|
|
y = C eα x cos β x + C |
eα x sin β x + C eν x cosδ x + |
|||||||
|
|
|
|
1 |
2 |
|
|
|
|
3 |
|
|
|
|
|
+ C eν x sinδ x + C eλ5x + ... + C eλnx. |
|
||||||
|
|
|
|
4 |
|
5 |
|
|
|
n |
|
П р и м е р 3 . Н а й т и о б щ е е р е ш е н и е у р а в н е н и я |
|||||||||||
y′′′ + 4y′′ +13y′ = 0 |
|
|
|
|
|
|
|
. |
|||
Решение. |
Характеристическое |
уравнение |
|||||||||
λ3 + 4λ2 +13λ = 0 |
имеет |
корни |
λ |
= 0,λ |
= −2 − 3i,λ = −2 + 3i . |
||||||
Общее решение: |
|
|
1 |
|
2 |
3 |
|||||
|
|
|
|
|
|
|
|||||
|
|
|
|
y = C1 + C2e−2x cos3x + C3e−2x sin 3x . |
□ |
||||||
541
4. Пусть λ1 = α + iβ является k -кратным корнем уравнения (3)
æ k £ |
n |
ö |
, |
λ = α - iβ также будет k -кратным корнем и пусть |
|
÷ |
|||
ç |
2 |
|
2 |
|
è |
ø |
|
|
остальные корни вещественны и различны. Тогда фунда- ментальная система решений уравнения (1) имеет вид
eα x cos β x, eα x sin β x, xeα x cos β x, xeα x sin β x,...,
(7)
xk−1eα x cos β x, xk−1eα x sin β x, eλ2k +1x ,...,eλnx.
Значит общее решение однородного дифференциаль- ного уравнения (1) запишется
y = C1eα x cos β x + C2eα x sin β x + C3xeα x cos β x + C4 xeα x sin β x + ... + +C2k−1xk−1eα x cos β x + C2k xk−1eα x sin β x + C2k+1eλ2k +1x + ... + Cneλnx.
Пример 4. Найти общее решение уравнения yV - 2yIV + 2y¢¢¢ - 4y¢¢ + y¢ - 2y = 0 .
Решение. Характеристическое уравнение имеет вид
λ5 - 2λ4 + 2λ3 - 4λ2 + λ - 2 = 0 .
542
