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

VM_pidr

.pdf
Скачиваний:
21
Добавлен:
19.03.2015
Размер:
2.73 Mб
Скачать

y

1

= y

10

, y

= y

 

 

,

y

2

= y

20

 

, y

= y.

 

 

 

 

 

 

 

 

1

 

10

 

 

 

 

 

 

 

2

 

 

 

20

 

 

 

 

 

Тоді С1 y10 + С2 y20

= y0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С

y

+ С

2

y

 

= y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

10

 

 

 

 

20

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

Розв’яжемо цю систему відносно С1 іС2 , одержимо:

 

 

 

 

 

С1

=

y

0

y

 

yy

20

=

y

0

y

yy

20 ,

 

 

 

 

 

 

 

 

20

0

 

 

20

 

0

 

 

 

 

 

 

 

y

10

y

 

y

y

20

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

20

10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С2

=

y

10

y

y

y

 

0

=

y

10

y

yy

0 .

 

 

 

 

 

 

 

0

10

 

 

0

 

10

 

 

 

 

 

 

 

 

y

10

y

 

y

y

 

20

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

20

10

 

 

 

 

 

 

 

 

 

 

 

Оскільки

0

0 , бо система розв’язків

y1 і y2 фундаменталь-

на, то для С1 і С2 справді знайдемо потрібні значення.

Таким чином, з розв’язку y = С1 y1 + С2 y2 можна знайти будь-

який частинний розв’язок, тобто розв’язок, що відповідає будь-яким початковим умовам, а це означає, що розв’язок y = С1 y1 + С2 y2 є

загальним розв’язком рівняння (7.24).

3.2Лінійні неоднорідні диференціальні рівняння другого порядку

Означення. Лінійним неоднорідним диференціальним рівнянням другого порядку називається рівняння :

y''+ p( x )y'+ g( x )y = f ( x ) ,

(7.27)

де функція f(x) називається правою частиною рівняння.

Рівняння y''+ p( x )y'+ g( x )y = 0 , яке одержується із рівнян-

ня(7.27), коли f(x)=0, називається однорідним рівнянням, що відповідає рівнянню (7.27).

ТЕОРЕМА ( про структуру розв’язку неоднорідного диференціального рівняння). Загальний розв’язок неоднорідного ди-

ференціального рівняння другого порядку дорівнює сумі загального розв’язку відповідного лінійного однорідного диференціального рівняння і будь-якого частинного розв’язку лінійного неоднорідного диференціального рівняння другого порядку.

Надалі будемо користуватись позначенням:

y=yз.о +yч.н,

(7.28)

де у - загальний розв’язок лінійного неоднорідного диференціального рівняння другого порядку; yз.о. - загальний розв’язок відповід-

371

ного однорідного диференціального рівняння; yч.н. - будь-який час-

тинний розв’язок лінійного неоднорідного рівняння. Доведення. Знайдемо похідні:

 

 

 

 

y′ = y

 

+ y

, y′′ = y′′

 

 

+ y′′

.

 

 

 

 

 

 

 

 

з.о.

 

ч.н.

 

 

з.о.

 

ч.н.

 

y, y′′ в рівнян-

 

Підставимо

y = yз.о. + yч.н.

і знайдені похідні

 

ня(7.27), матимемо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y′′

+ y′′

+ p( x )( y

+ y

 

) + g( x )( y

з.о.

+ y

ч.н.

) =

f ( x ) або

з.о.

ч.н.

 

 

з.о.

 

ч.н.

 

 

 

 

 

 

 

 

 

 

y′′

+ p( x )y

+ g( x )y

з.о.

+ y′′

+ р( x )y

 

+ g( x )y

ч.н.

= f ( x ) .

з.о.

 

з.о.

 

 

 

 

 

ч.н.

 

 

ч.н.

 

 

 

 

 

 

 

Оскільки

yз.о. є загальним розв’язком однорідного рівняння

y''+ p( x )y'+ g( x )y = 0 , то в останній рівності

 

 

 

 

 

y''з.о. + p( x )y'з.о. + g( x )yз.о.

= 0 .

 

 

 

 

 

 

 

 

 

 

 

 

А з того, що

yч.н. - частинний розв’язок рівняння (7.27), мати-

мемо тотожність: y''ч.н. + p( x )y'ч.н. + g( x )yч.н.

= f ( x ) .

 

 

Отже, ми довели, що функція

y = yз.о. + yч.н. - розв’язок дифе-

ренціального рівняння (7.27)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Легко показати,

що кожний розв’язок рівняння (7.27), який

задовольняє будь-які початкові умови (з певної області їх визначення), можна знайти з функції y = yз.о. + yч.н. . Отже, ця функція є за-

гальним розв’язком. Довести це твердження можна аналогічно до доведення такого ж твердження для відповідного однорідного рівняння, яке ми привели вище, тому тут на цьому доведенні не зупиняємось.

3.3. Метод варіації довільних сталих

Розглянемо метод варіації довільних сталих (метод Лагранжа) для знаходження частинного розв’язку неоднорідного диференціального рівняння.

Суть цього методу полягає в наступному: щоб знайти частинний розв’язок неоднорідного лінійного рівняння (7.27), досить у ви-

раз для загального розв’язку

y=C1y1(x)+C2y2(x)

 

відповідного одно-

рідного рівняння замість сталих C1 і C2 підставити функції незалеж-

ної змінної

x,

похідні від яких С1

( x )

і

С2

 

задовольняють таку

 

 

 

 

 

 

 

( x )

 

 

систему алгебраїчних рівнянь:

+ С2

 

 

2

 

=

 

 

 

С1

( x )y

1

( x )

 

 

( x )

0 ,

 

 

 

 

( x )y

 

 

 

 

С( x )y( x ) +

C( x )y( x )

= f ( x ).

 

 

1

1

 

 

2

 

2

 

 

 

372

Доведемо це твердження.

Запишемо розв’язок диференціального рівняння у вигляді:

y=C1y1(x)+C2y2(x)

Знайдемо похідну:

y′ = Cy

1

( x ) + Cy

2

( x ) + C

1

( x )y

( x ) + C

2

( x )y

( x ) .

1

2

 

1

 

2

 

Ми хочемо визначити дві функції C1(x) і C2(x).

Одне співвідношення між ними ми можемо вибрати довільним. Поставимо вимогу, щоб C1(x) і C2(x) задовольняли рівність

C( x )y

1

( x ) + C( x )y

2

( x ) = 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

y

 

 

1

 

 

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

тоді

=

C

 

 

 

1

 

 

+

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( x )y( x )

 

 

 

 

( x )y( x ) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Знайдемо другу похідну:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y′′ =

C

( x )y( x ) +

C( x )y( x ) + C

1

( x )y′′( x ) + C

2

( x )y′′( x ).

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

 

 

2

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

2

 

 

 

Підставимо значення y,y′,y″в диференціальне рівняння (7.27),

будемо мати

 

+ р1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C( x )y′′( x )

y

( x )

+ р2

y

1

( x ))

+

C

2

( x )( y′′( x )

+ р1

y( x )

+

1

 

 

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

+ p

2

y

2

 

( x )) + ( C( x )y( x ) + C

( x )y

( x ) = f ( x ).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Оскільки y1(x) та y2(x) є розв’язками однорідного рівняння, то

вирази :

 

 

 

 

 

 

 

 

 

y′′( x ) + p y

( x ) + p

 

 

 

 

( x ) = 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

y

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y′′( x ) + p y

( x ) + p

2

y

2

( x ) = 0 .

 

 

 

 

 

 

 

 

Звідси

 

 

1

 

1

 

2

 

 

+

 

 

2

 

 

1

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

f ( x ) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C( x )y( x )

 

 

C( x )y( x )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отже, для того щоб функція

 

y = C1 ( x )y1 ( x ) + C2 ( x )y2 ( x ) ,

яка

 

 

задовольняє

 

умову

 

 

 

 

 

 

1

 

 

 

 

1

( x )

+

 

2

 

2

( x )

=

0

,

була

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

( x )y

 

 

 

C( x )y

 

 

 

розв’язком рівняння (7.27), необхідно і достатньо, щоб

 

 

 

 

 

виконувалась рівність

 

 

1

 

 

 

 

 

 

 

1

 

 

+

 

 

2

 

 

 

 

 

 

2

=

f ( x ) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

: C

( x )y( x )

 

 

C( x )y( x )

 

 

 

 

 

 

 

Таким чином дістаємо систему

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

x ) y

1

( x ) + C ( x ) y

2

( x ) = 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1(

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

( x ) y

 

( x ) + C

( x ) y

( x ) = f ( x ),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

 

 

 

 

1

та

 

 

 

2

 

2

 

 

 

 

 

 

 

2

 

 

система

має

єдиний

з якої визначаємо

 

 

 

 

 

 

 

 

 

 

C

 

 

 

 

 

 

 

 

Ця

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C( x )

 

 

 

 

 

 

( x ) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

розв’язок, бо за нашим припущенням y1(x) та y2(x) -лінійно незалежні розв’язки однорідного рівняння.

Нехай

C

( x )

= ϕ 1

( x )

і C

( x )

= ϕ

2

( x ) , тоді, інтегруючи,

 

1

 

 

2

 

 

 

 

одержимо

C1 ( x ) = ϕ1 ( x )dx +

 

1 , C2 ( x ) = ϕ2 ( x )dx +

 

2 , де

С

С

С1 , С2 - довільні сталі. Підставимо знайдені С1 ( x ) і С2 ( x ) в

373

співвідношення y = С1 ( x )y1 ( x ) + С2 ( x )y2 ( x ) і отримаємо загальний розв’язок диференціального рівняння (7.27).

Приклад. Знайти загальний розв’язок рівняння y′′ − y= x x

Розв’язування. Знайдемо загальний розв’язок цього рівняння методом варіації довільних сталих. Знайдемо загальний розв’язок

відповідного однорідного рівняння

 

y′′ −

y

= 0 ,

 

y′′

=

1

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

y

Інтегруючи обидві частини рівняння, одержимо

ln

 

y

 

= ln

 

x

 

+ ln

 

C

 

; y′ = Cx , звідки y = C1 x2 + C2 .

 

 

 

 

 

 

Тепер припустимо, що C1

 

іC2 є функціями від x і складаємо

систему для знаходження

 

1

 

і

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C( x )

 

 

C( x ) :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C( x )x2 + C( x ) = 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2C( x )x = x.

 

 

 

 

 

 

 

 

 

 

 

Звідки C

( x )

=

1

,

C ( x )

= −

 

1

 

 

x 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

2

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В результаті інтегрування, маємо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C1 ( x ) =

 

 

 

 

+ C 1 , C2

( x ) = −

 

 

 

 

+ C 2 , де C 1 ,C 2 - довільні сталі.

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Підставляємо C1 ( x ) і C2 ( x ) в загальний розв’язок відповід-

ного однорідного рівняння

y = C1 ( x )x2 + C2 ( x ) ,

 

 

 

 

 

 

одержимо загальний розв’язок даного рівняння:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y =

 

1 x 2

+

 

2

+

x 3

x 3

 

або y =

 

1 x 2

+

 

2 +

x 3

.

 

 

 

С

С

С

С

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

3.4. Лінійні однорідні диференціальні рівняння другого

 

 

 

 

 

 

 

 

 

порядку з сталими коефіцієнтами

 

 

 

 

 

 

Означення. Рівняння вигляду

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y''+ py'+ gy = 0 ,

 

 

 

 

 

(7.29)

де p,g- сталі числа, називається лінійним однорідним диференціальним рівнянням другого порядку з сталими коефіцієнтами.

Знайдемо загальний розв’язок цього рівняння. За теоремою про структуру загального розв’язку лінійного однорідного диференціального рівняння другого порядку маємо: y=C1y1(x)+C2y2(x),

де y1(x),y2(x) - лінійно незалежні частинні розв’язки.

374

Частинний розв’язок рівняння (7.29) будемо шукати у вигляді y = ekx , k - стала , яку треба знайти.

Знайдемо y′ = kekx , y′′ = k 2ekx і підставимо в рівняння (7.29),

маємо: ekx ( k 2 + pk + g ) = 0. Оскільки ekx 0

завжди, то

k2 + pk + g = 0 .

(7.30)

Рівняння (7.30) називається характеристичним рівнянням для диференціального рівняння (7.29).

Позначимо через k1 ,k2 - корені характеристичного рівняння. При розв’язанні характеристичного рівняння, яке є квадратним рів-

нянням відносно k , можливі три випадки:

k1 k2 ,

1)

Корені характеристичного рівняння дійсні, різні:

( D > 0 ) .

 

2)

Корені характеристичного рівняння дійсні, рівні:

k1 = k2 ,

( D = 0 ) .

 

3)

Корені характеристичного рівняння комплексні числа

k1,2 = α ± βi , ( D < 0 ) .

Розглянемо кожен випадок окремо.

1)Корені характеристичного рівняння дійсні, різні: k1 k2 . Частинними розв’язками рівняння (7.29) в цьому випадку бу-

дуть функції y1 = ek1 x та

y2 = ek2 x .

 

Ці розв’язки лінійно незалежні, бо

 

 

y2

=

ek2 x

= e( k2 k1 )x const .

 

 

y1

 

 

 

 

ek1 x

 

Отже, загальний розв’язок рівняння (7.29) запишеться

 

y = C1ek1 x + C2ek2 x .

(7.31)

Приклад 1. Розв’язати диференціальне рівняння

 

 

y′′ + y′ − 2 y = 0 .

 

Розв’язування. Складаємо характеристичне рівняння:

k 2 + k 2 = 0. Тут корені k1 = 1;k2 = −2 -дійсні, різні.

Загальний

розв’язок диференціального рівняння запишеться: y = C1ex + C2e2 x . 2) Корені характеристичного рівняння дійсні, рівні: k1 = k2 . Один частинний розв’язок може бути y1 = ek1 x .

375

Розв’язок ek2 x лінійно залежний з першим ek1 x , тому його розглядати в якості другого частинного розв’язку не можна. Будемо

шукати другий частинний розв’язок у вигляді y2 = u( x )ek1 x , де u( x )- невідома функція, яку можна знайти.

 

 

Знайдемо

y

та

y′′ :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

= uek1 x + k

uek1 x = ek1 x ( u′ + k

u ),

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

1

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

y′′

=

u′′ek1 x

+

2k

1

uek1 x

+

k 2uek1 x

=

ek

1 x ( u′′

+

2k

1

u

+

k

2u )

і підста

-

2

 

 

 

 

 

 

1

 

 

 

 

 

 

1

 

вимо в рівняння (7.29):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ek1 x ( u′′ + ( 2k1 + р)u′ + ( k12 + pk1 + q )u ) = 0.

 

 

 

Оскільки

 

 

k1 - кратний корінь характеристичного рівняння

(7.30), то k12 + pk1

+ q = 0. Крім того k1 = k2 = −

p

 

або 2k1 = − p ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

2k1 + p = 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отже, для того, щоб знайти u( x ), треба розв’язати рівняння

 

 

 

ek1 x u′′ = 0

 

 

або

u′′ = 0.Інтегруючи

 

двічі,

одержимо

u = Аx + B. Якщо покласти A = 1,B = 0 , то u = x.

 

 

 

 

 

 

 

Таким чином, в якості другого частинного розв’язку будемо

брати функцію

 

y2

= xek1 x .Розв’язки

y1 та

y2

 

будуть лінійно неза-

лежними, бо

y2

= x const .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отже, загальний розв’язок диференціального рівняння (7.29)

буде :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = C1ek1 x + C2 xek1 x = ek1 x ( C1 + C2 x ) .

 

(7.32)

 

 

Приклад 2. Розв’язати диференціальне рівняння

 

 

 

 

 

 

 

 

 

 

 

 

 

y′′ − 4 y′ + 4 y = 0.

 

 

 

 

 

 

 

 

Розв’язування. Характеристичне рівняння k 2 4k + 4 = 0 , корені якого k1 = k2 = 2 - дійсні, рівні.

Загальний розв’язок диференціального рівняння буде:

y= C1e2 x + C2 xe2 x = e2 x ( C1 + C2 x ) .

3)Корені характеристичного рівняння комплексні числа. Якщо дискримінант характеристичного рівняння D < 0 , то

дійсних коренів це рівняння не має.Введемо поняття комплексного

376

числа. Будь-яке комплексне число можна представити у вигляді: z = α + βi де i = 1 або i2 = −1.

При цьому α - називається дійсною частиною комплексного числа; β - уявною частиною комплексного числа .

Числа z = α + βi та z = α − βi називаються комплексно

спряженими числами.

Повернемось тепер до розв’язання характеристичного рівняння (7.30), дискримінант якого менший нуля. Тоді його корені будуть комплексно спряженими числами, тобто k1,2 = α ± βi .

У випадку, коли k1 = α + βi , розв’язок диференціального

рівняння буде y1 = e( α + βi )x = eαx eiβx .

Якщоалгебраїчнерівнянняздійснимикоефіцієнтамимаєкорінь k1 = α + βi , то воно має і спряжений з ним корінь k2 = α − βi ,

а тому y2 = e( α −βi )x = eαx eiβx .

За формулами Ейлера : e± xi = cos x ± i sin x ,запишемо функції y1 та y2 у вигляді:

y1 = eαx (cosβx + i sinβx ) = eαx cosβx + ieαx sinβx , y2 = eαx (cosβx i sinβx ) = eαx cosβx ieαx sinβx .

Далі використаємо таку теорему.

Теорема. Якщо диференціальне рівняння (7.29) має розв’язком функцію u( x ) + iv( x ) , то кожна з функцій u( x ) і

v( x ) є розв’язком рівняння (7.29).

Доведення. Дійсно, підставимо функцію u( x ) + iv( x ) в

рівняння (7.29), маємо

( u( x ) + iv( x ))′′ + p( u( x ) + iv( x ))′ + q( u( x ) + iv( x )) 0 або ( u′′ + pu′ + qu ) + i( v′′ + pv′ + qv ) 0 .

Але комплексна функція дорівнює тотожно нулю тоді і тільки

тоді, коли рівні нулю її дійсна і уявна частини, тобто: u′′ + pu′ + qu = 0 ,

v′′ + pv′ + qv = 0.

А це означає, що u( x ) і v( x ) є розв’язками рівняння (7.29).

377

На основі доведеної вище теореми, частинними розв’язками

будуть функції

 

1 = eαx cosβx ,

 

2 = eαx sinβx ,

причому вони бу-

y

y

дуть лінійно незалежними, бо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

αx cosβx

 

 

 

 

 

 

 

 

 

 

y

1

=

= ctgβx const .

 

 

 

 

 

 

 

 

 

2

e

αx sinβx

 

 

 

 

 

 

 

 

y

 

 

 

Отже, загальний розв’язок рівняння(7.29) буде

y = С1

 

1 + С2

 

2 = С1eαx cosβx + С2eαx sinβx

або

y

y

 

 

 

 

y = eαx ( С1 cosβx + С2 sinβx ) .

(7.33)

Приклад 3. Знайти частинний розв’язок диференціального рів-

няння y′′ + 2 y′ + 5 y = 0 ,

що задовольняє таким початковим умовам

y x = 0 = 0 , yx = 0 = 1 .

Розв’язування. Шукаємо спочатку загальний розв’язок. Для цього складаємо характеристичне рівняння і розв’язуємо його:

k 2 + 2k + 5 = 0 , D = 4 20 = −16 < 0 ,

 

k1,2

= 2 ± − 16 = 2 ± 16 ( 1 ) = 2 ± 4i = −1 ± 2i -

 

2

2

2

корені комплексно спряжені числа. Загальний розв’язок запишеться: y = ex ( C1сos2 x + C2 sin 2 x ) .

Знайдемо частинний розв’язок, використовуючи початкові умови :

0 = e0 ( C1 cos( 2 0 ) + C2 sin( 2 0 )), C1 = 0 ;

 

 

 

y′ = ex 2C2сos2 x exC2 sin 2 x; 1 = 2C2 , C2

=

1

,

 

 

2

 

звідки y = 1 ex sin 2 x.

2

3.5 Лінійні неоднорідні диференціальні рівняння другого порядку з сталими коефіцієнтами

Означення. Рівняння

y''+ py'+ gy = f ( x ),

(7.34)

де p, g - сталі числа, f ( x ) 0 , називається

лінійним

неоднорідним диференціальним рівнянням другого порядку з сталими коефіцієнтами.

За теоремою про структуру розв’язку лінійного неоднорідного

378

рівняння другого порядку загальний розв’язок рівняння (7.34) є сумою загального розв’язку відповідного однорідного рівняння і будьякого частинного розв’язку неоднорідного рівняння. Це твердження записано формулою (7.28): y = yз.о. + yч.н. .

Загальний розв’язок однорідного рівняння ми детально розглянули вище. Тепер перейдемо до знаходження частинного розв’язку неоднорідних рівнянь із спеціальною правою частиною, розв’язок яких можна знайти не вдаючись до інтегрування.

 

 

 

 

1)

Нехай

 

права

частина рівняння

(7.34)

має

вигляд

f ( x ) = Pn ( x )eαx

 

, де Pn ( x ) - многочлен n-го степеня. Тут можливі

два випадки:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а)

α

- не є коренем

характеристичного

рівняння, тоді

y

ч

.н.

= Q ( x )eαx ,

де Q

n

( x ) - многочлен n-го степеня із неозначе-

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ними коефіцієнтами.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б)

α

- є

 

коренем

характеристичного

рівняння кратності

r ( r = 1 або r = 2 ) , тоді

y

ч.н.

= xrQ ( x )eαx .

 

 

 

 

 

 

 

 

 

 

Зауваження. Якщо

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

f ( x ) = Pn ( x ) ,

то вважаємо , що α = 0 і

перевіряємо , чи 0 є коренем характеристичного рівняння.

 

 

 

 

 

Приклад 1. Знайти загальний розв’язок рівняння

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 y′′ − y′ − y = 4 xe2 x .

 

 

 

 

 

(7.35)

 

 

 

 

Розв’язування. Загальний

розв’язок

шукаємо

у

вигляді

y = yз.о. + yч.н. .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Спочатку знайдемо

загальний

розв’язок

yз.о.

відповідного

однорідного рівняння:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 y′′ − y′ − y = 0; 2k 2 k 1 = 0; k1 = 1,k2

= −

1

,

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yз.о. = C1e x + C2e 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Частинний розв’язок

yч.н.

неоднорідного рівняння шукаємо у

вигляді

правої

частини

рівняння,

а

саме

y

= ( Ax + B )e2 x ,

оскільки

α = 2

 

 

 

 

 

 

 

 

 

 

 

 

 

ч.н.

 

 

 

не є коренем характеристичного рівняння. Тут

потрібно знайти неозначені коефіцієнти

A та В . Для цього знайде-

мо ( yч.н. )

та ( yч.н. )′′ , маємо:

 

 

 

 

 

 

 

 

 

( y

 

)′ = 2e2 x ( Ax + B ) + e2 x A;

 

 

 

 

 

 

 

 

 

 

 

ч.н.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

379

( yч.н. )′′ = 4e2 x ( Ax + B ) + 2e2 x A + 2e2 x A = 4e2 x ( Ax + B ) + 4 Ae2 x

Підставляємо yч.н. , ( yч.н. )та ( yч.н. )′′ у рівняння (7.35), одержимо

8Axe2 x + 8Be2 x + 8Ae2 x 2Axe2 x 2Be2 x Ae2 x Axe2 x Be2 x = 4xe2 x .

Розділимо ліву і праву частини рівняння на e2 x , маємо :

8 Ax + 8B + 8 A 2 Ax 2B A Ax B = 4 x ,

5 Ax + 5B + 7 A = 4 x .

Відомо, що многочлени рівні між собою тоді і тільки тоді, коли рівні їх коефіцієнти при однакових степенях x лівої і правої частини, тобто

 

 

 

5 A = 4

 

, звідки

 

A =

4

 

; B = −

28

.

 

 

 

7 A + 5 B = 0

 

 

 

 

 

 

 

 

5

 

25

 

Отже,

 

 

 

 

 

 

4

 

 

 

 

28

 

2 x

і

загальний розв’язок буде

 

 

yч.н.

 

=

 

 

 

x

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

25

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

x

 

 

4

 

 

 

 

28

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

2 x

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = C

1e

 

+ C

2 e

 

 

 

+

 

 

 

x

 

 

 

e

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

25

 

 

 

 

 

 

 

 

 

 

 

 

Приклад 2. Знайти розв’язок диференціального рівняння

 

 

 

 

 

 

 

 

 

y′′ − 2 y′ + y = 1 + x ,

(7.36)

який задовольняє початковим умовам y

 

x =0 = 2; y

 

x =0 = −3 .

 

 

 

Розв’язування.

Тут характеристичне рівняння k 2 2k + 1 = 0

має дійсні, рівні корені

 

k1 = k2 = 1 ,

 

тому загальний розв’язок

відповідного однорідного рівняння буде

 

yз.о. = ( C1 + C2 x )e x .

Права частина рівняння (7.36) має вигляд Pn ( x ) = 1 + x , причому α = 0 не є коренем характеристичного рівняння, тому частин-

ний

розв’язок

рівняння

будемо

шукати

у

вигляді

yч.н. = Ax + B.Продиференціювавши

yч.н. ,

підставимо у рівняння

(7.36), маємо :

 

 

 

 

 

 

2 A + Аx + B = 1 + x , звідки

A = 1;

2 A + B = 1; B = 3.

 

 

Частинним

розв’язком

заданого

рівняння

є

функція

yч.н.

= x + 3, а його загальним розв’язком функція

 

 

y = ( C1 + C2 x )e x + ( x + 3 ).

Знайдемо тепер частинний розв’язок рівняння (7.36), який задовольняє заданим початковим умовам. Для цього знайдемо

y′ = ( C1 + C2 x )e x + С2e x + 1 , тоді

380

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]