Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Математика для инженеров(теория)I том

.pdf
Скачиваний:
128
Добавлен:
18.02.2016
Размер:
4.09 Mб
Скачать

После того, как число 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¢ =

×et =

+ z,

y¢¢ =

et = ç

 

 

+

÷et , исход-

 

 

 

 

 

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¢¢

 

et

 

 

 

 

 

 

=

 

. Отсюда,

в силу соотношения

1+ t

d( y′) = y′′dx ,

т.е.

dt = y′′dx , имеем

(1+ t)dt = et 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 , если

существуют

постоянные числа

α12 ,...,α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 - yy

2

) = 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

 

1

 

 

 

 

2

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В

 

силу

того,

 

 

что

 

y y

- yy

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

- yy

2

= 0 . Допустим, что

y ¹ 0

на I .

 

 

1

2

 

1

 

1

 

Тогда,

на

основании

последнего равенства,

можно

напи-

537

 

y y

- yy

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)называется характеристическим.

Таким образом, решение однородного дифференциального уравнения свелось к решению алгебраического уравнения. Этот метод решения называется методом Эйлера.

Пусть λ12 ,...,λn корни уравнения (3), причем среди них могут быть и кратные. Возможны следующие случаи.

1. λ12 ,...,λ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 - = 0 . Находим его корни: λ

= 0,λ

= -1,λ

 

= 3 . Так

1

2

3

 

как они действительные и различные, то общее решение

имеет вид

y= C1 + C2ex + 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= C1ex + C2 xex + 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 - 4 + 3 - 2 + λ - 2 = 0 .

542