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

книги / Обыкновенные дифференциальные уравнения

..pdf
Скачиваний:
1
Добавлен:
12.11.2023
Размер:
2.7 Mб
Скачать

комплексных корней характеристического уравнения (4.38) соответствует два действительных частных решения уравнения (4.34) вида (4.43).

Пример 6.

y′′′ + y = 0.

 

Характеристическое уравнение

k 3 + 1 = 0

его корни

k

= −1

,

 

 

 

,

1

 

k2,3

= 1

± i

3

. Следовательно, общее решение имеет вид

 

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

3 + C

 

3

 

 

 

 

 

 

 

 

y = C ex + e 2

 

cos x

sin x

 

 

 

 

 

 

 

C

2

 

 

 

 

 

 

 

 

1

 

 

 

2

3

 

2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

II. Пусть среди корней уравнения (4.38) существуют кратные.

В этом случае количество различных в ряду (4.39) чисел будет меньше n, и, соответственно, число линейно независимых частных решений вида (4.40) будет меньше n; этих решений недостаточно для получения общего решения. Для получения недостающих решений изучим выражение линейного оператора L от произведения двух функций uυ .

По формуле Лейбница имеем:

(uυ )

( n )

= u

 

υ +

Cnu

υ +′

Cn u

υ +

′′

+...

Cnu′ +

υ u

( n )

,

 

 

 

( n )

 

1

( n1)

 

2

( n2)

 

 

1 ( n1)

 

 

 

(uυ )

( n1)

=

u

 

υ +

Cn1u

υ

+′

Cn1u

υ

+ ′′

+... υ u

( n1)

,

 

 

 

 

 

( n1)

 

1

( n2)

 

2

( n3)

 

 

 

 

……..

(uυ )′′= uυ′′+ 2uυ ′+′ υ u ′′,

(uυ )′= uυ+′

υu ,

(uυ )= uυ ,

 

 

где C nk =

n!

 

– число сочетаний из n элементов по k элементов.

 

 

 

 

 

k!(n k )!

Умножая первую строку на 1, вторую на a1, …, последнюю на a n и складывая получим:

L(uυ ) = υ L[u]+

υ ′

L [u]+

υ ′′

L [u]+ ... +

υ (n1)

 

L

 

[u]+

υ (n)

L

[u]

.

 

 

(n 1)!

 

n!

 

 

1!

1

 

2!

2

 

n1

 

 

 

n

Здесь введены обозначения:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L(y) = y(n) + a1y(n1) + a2 y(n2) +...+ an1y′+ an y,

,

 

 

 

1

=

ny(n1)

+

 

 

1

+

 

+

 

n2

y

+

n1

 

 

 

 

L (y)

 

 

(n

 

1)a y(n2)

 

 

...

 

2a

 

 

a

y

 

 

 

L (y) = n(n 1)y(n2)

+ (n 1)(n 2)a y(n3)

+... + 2 1 a

n2

y ,

 

2

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

Ln1 ( y) = n(n 1)...2 y′ + (n 1)(n 2)...1 a1 y, Ln ( y) = n(n 1)...2 1 y.

(4.44)

(4.45)

91

Операторы Lr = 1, 2, …, n составлены из L по правилу, аналогичному правилу дифференцирования многочлена, только роль показателей играют указатели порядка производной.

Формула (4.44) применима к любому линейному оператору. Если, в частности, коэффициенты a1, a2, …, an есть постоянные, то каждому оператору Lr соответствует характеристический многочлен Fr(k), причем

легко видеть, что Fr(k) есть r

 

производная

по k многочлена

 

F(k),

соответствующего оператору L:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F (k) = F (r) (k) .

 

 

 

 

(4.46)

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

Вычислим теперь

выражение

(4.44), если

 

u = e kx , υ

= x m , где

 

m

целое неотрицательное число.

 

 

 

 

 

 

 

 

 

 

 

Мы получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L[xmekx ]= xmL[ekx ]+

m

xm1L [ekx ]+

m(m 1)

xm2L [ekx

]+...+Cm1xL

[ekx ]+ L

 

[ekx ]

 

 

 

1

 

1

 

1 2

2

 

m

m1

 

m

 

 

 

 

 

 

 

 

 

 

 

 

или, поскольку, в силу формулы (4.46)

 

 

 

 

 

 

 

 

 

 

 

 

L (ekx ) = ekxF (r) (k) ,

 

 

 

 

 

 

 

получим:

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kx

 

 

 

 

 

 

 

 

 

 

 

 

 

m

e

=

 

 

 

 

 

 

 

 

 

 

 

L x

 

 

 

 

 

 

(4.47)

= ekx {xm F (k ) + Cm1 xm1F (k ) + Cm2 xm2 F ′′(k ) + ... + Cmm1xF m1 (k ) + F m (k )} .

Пусть теперь k1 есть корень характеристического уравнения (4.38) кратности m1; тогда, как известно,

F (k1) = 0; F (k1) = 0; ... F (m1 1) (k1) = 0; F (m1 ) (k1) 0 .

Если в выражении (4.47) взять показатель m при x меньшим, чем m1, то все члены в скобке правой части обратятся в нуль; следовательно, мы получаем m1 частных решений дифференциального уравнения (4.34), соответствующих корню k1:

e k1x , xe k1x , x 2 e k1x ,...., x m1 1e k1x.

(4.48)

Аналогично, если имеются другие корни характеристического уравнения, k 2

кратности m2, …, kp кратности mp, mi

1, причем m1 + m2 +…+ mp = n, и все

k r уже различны, то им будут соответствовать частные решения:

 

ek2 x , xek2 x , x2ek2 x , ..., xm2 1ek2 x ,

 

 

 

 

 

 

 

 

 

 

 

 

 

...

 

 

, xe p

 

, x2e

 

 

, ..., x

 

 

(4.49)

e

p

x

x

p

x

m

p e p .

 

k

 

k

k

 

 

1 k x

 

Совокупность решений (4.48) и (4.49) даст в общем случае кратных корней n частных решений. Остается доказать, что они линейно независимы, т.е. образуют фундаментальную систему.

92

Допустим, что между этими решениями существует тождественно линейное соотношение

p

 

 

 

 

 

 

mr 1 ) kr x

p

kr x

 

,

 

(A

(r)

+ A

(r)

x + ... + A

(r)

x

= P (x)e

= 0

(4.50)

 

 

 

e

 

 

0

1

mr 1

 

 

r

 

 

 

 

r=1

 

 

 

 

 

 

 

r=1

 

 

 

 

где коэффициенты A(jr) – постоянные. Без ограничения общности можно предположить, что в многочлене Pp(x) по крайней мере один коэффициент отличен от нуля. Разделим обе части этого соотношения на ek1x :

P (x) +

p P (x)e(kr k1)x = 0

.

1

r

 

r=1

 

Дифференцируя последнее (предполагаемое) тождество m1 раз по x, мы вместо первого многочлена получим нуль, а все многочлены, стоящие множителями при показательных функциях, заменяются новыми многочленами тех же степеней, и мы получим новое тождество:

p

 

Qr (x)e(kr k1)x = 0.

(4.51)

r=2

 

Очевидно, что Q p ( x) не равно нулю тождественно. Сумма

(4.51)

содержит уже p – 1 слагаемых; продолжая тот же процесс, мы придем, наконец, к тождеству

Rp (x)e(kp kp1)x = 0.

(4.52)

Но тождество (4.52) невозможно, так как e(k p k p1)x ≠ 0 , а многочлен R p ( x), будучи той же степени, что и Pp ( x ), имеет по крайней мере один

коэффициент отличным от нуля и не тождественно равен нулю.

Доказав линейную независимость частных решений (4.48) и (4.49), мы можем написать общее решение уравнения (4.34) в случае кратных корней в виде

p

y = Gr (x)ekr x = 0 , (4.53)

r=1

где Gr (x) есть многочлен степени mr 1 с произвольными коэффициентами. Число произвольных постоянных в выражении (4.53)

m1 + m 2 + ... + m p = n,

т.е. порядку уравнения, как и должно быть.

В случае комплексных кратных корней характеристического уравнения выражение (4.53) неудобно, так как оно окажется комплексной функцией действительного переменного x. Заметим, что комплексные корни войдут попарно сопряженными с одинаковой кратностью. Если корень k1 = α + β i

93

(β ≠ 0) имеет кратность m1 , то сопряженный корень k2 = α − β i имеет ту же кратность. Соответствующая корню k1 совокупность решений (4.48) будет:

e (α

i ) x , xe (α +β i ) x ,..., x m1 1e (α +β i ) x .

Отделим в этих

выражениях

действительные части от мнимых

и получим, таким образом, 2m1 решений:

 

 

 

 

 

e

α

x

 

α

x

cosβ

x, ..., x

m1

1α

x

cosβ

x,

 

 

cosβ x, xe

 

 

e

 

e

α x

α

x

sinβ

x, ..., x

m1

1α

x

sinβ

x.

 

 

sin β x, xe

 

 

e

 

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

Пример 7.

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

Характеристическое уравнение:

 

k 3 k 2 k + 1 = 0

его корни:

 

 

,

k1 = k2 =1, k3 = −1. Общее решение:

y = ex (C1 + C2 x) + C3ex.

 

Пример 8.

 

 

 

 

 

 

 

y(4)

+

8y′′

+

16y

0

 

 

 

 

 

= .

 

 

Характеристическое уравнение: k 4 + 8k 2 + 16 = 0

или (k2 + 4)2 = 0, его

корни: k1 = k2 = 2i, k3 = k4 = – 2i.

Общее

решение:

y = (C1 + C2)cos2x +

+(C3 + C4x)sin2x.

4.5.Линейные неоднородные уравнения с постоянными коэффициентами и правой частью специального вида

Рассмотрим уравнение вида

 

 

 

Ly = f (x).

(4.54)

Общее решение этого уравнения имеет, как известно, вид

y =

 

+ y*.

y

Метод нахождения общего решения y однородного уравнения (4.34) указан выше. Частное решение y* уравнения (4.54) может быть найдено методом вариации произвольных постоянных. Однако, если правая часть f(x) имеет специальный вид, то y* определяется без интегрирования, а лишь с помощью

алгебраических операций.

Начнем с теоремы, относящейся к любому линейному уравнению с правой частью.

Теорема. Пусть y1* и y2* соответственно частные решения уравнений:

Ly = f1(x), Ly = f2 (x) .

94

 

Тогда

функция

y* = y

*

+ y

*

 

есть

частное

решение

уравнения

 

 

 

 

 

1

 

2

 

 

 

 

 

 

Ly = f1(x) + f2 (x) . Действительно L( y1* + y2* ) = Ly1* + Ly2*.

 

 

L( y1* + y2* ) =

 

Но,

по

условию,

Ly1* = f1 ( x),

Ly2* = f2 ( x),

откуда

= f1(x) + f2 (x), что и требовалось доказать.

 

 

 

 

 

Рассмотрим

теперь

линейное

уравнение

с

постоянными

коэффициентами и правой частью специального вида:

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

r ≥ 0

 

 

 

Ly = P (x)eα

r x

 

P (x)

– многочлены степени

.

 

 

 

r =1

r

, где

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В силу только что доказанной теоремы достаточно уметь находить

частное решение уравнения вида

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ly = Pm (x)eα

x ,

 

 

(4.55)

где

P (x) = p xm + ...+ p

 

 

 

 

 

 

 

 

 

 

m

m

 

0 есть многочлен степени m 0.

 

 

 

Рассмотрим два случая:

I. α не является корнем характеристического уравнения (4.38): F(α ) 0 .

Докажем, что в этом случае существует частное решение того же вида, что и правая часть, а именно:

y* = Qm (x)eα x ,

(4.56)

где Qm (x) = qmxm + qm1xm1 + ... + q0 .

 

Рассматривая коэффициенты qm, qm–1, ..., + q0 как

неизвестные,

покажем, что их можно определить так, чтобы выполнялось следующее тождество по x:

L (Qm (x)e

α x

)

 

 

α

x

или

 

 

= Pm (x)e

 

 

e

α x

L (Qm

 

 

α

x

) = Pm

(x) .

(4.57)

 

(x)e

 

Вычислим левую часть, применяя формулу (4.47):

 

 

e

−α x

 

 

 

 

α x

) =

 

 

 

 

 

 

 

L(Qm (x)e

 

 

 

 

= q {xm F (α )+

C1

xm1F′α(

+)

 

...+

F (mα) (

+)}

 

m

 

 

 

m

 

 

 

 

 

 

 

 

 

(4.58)

+qm1{x

m1

F (α )+

1

 

m2

F′α(

+)

 

...+

F

(m1)

+)} ...

 

Cm x

 

 

α (

 

 

... + q1{xF (α )+

F′α(

)}+

q0 Fα (

).

 

 

Приравнивая выражение (4.58) многочлену Pm(x) и отождествляя коэффициенты при одинаковых степенях x, получим m + 1 уравнений с m + 1 неизвестными qm, qm–1, ..., q0:

95

qm F (α ) = pm

 

 

 

 

 

 

 

 

 

 

 

 

 

qm1F (α ) + qmC1m F (α ) = pm1

 

 

 

 

 

 

 

q

m2

F (α ) + q

C1

 

 

F (α

) + q

m

C 2 F ′′(α

) =

p

m2

 

 

m1

m1

 

 

m

 

 

 

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

mr

F (α ) + q

mr

C1

 

 

F (α

) + ... + q

m

C r

F (r) (α )

 

 

+1

mr +1

 

 

 

m

 

 

 

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

0

F (α ) + q F (α )

+ q

2

F ′′(α

) + q F ′′′(α ) + ... + q

m

F (m)

 

 

1

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.59)

 

 

= pmr

 

 

 

 

 

 

 

 

 

 

(α ) = p0

 

 

 

 

Поскольку по условию α не является корнем характеристического уравнения, F(α )0. Система дает возможность последовательно вычислить

qm, qm–1, ..., q0:

 

 

 

qm =

 

pm

 

,

 

 

 

 

 

 

F (α

)

 

 

 

 

 

 

 

 

 

 

 

qm 1

=

1

(pm 1 qm C1m F (α ))=

pm 1

C1m

qm

F (α ) и т.д.

 

 

[F (α )]2

 

 

F (α )

 

F (α )

 

Таким

образом, мы находим

 

искомое

частное решение (4.56)

(разрешимость системы (4.59) относительно qm, qm–1, ..., q0 можно сразу усмотреть из того, что ее определитель (F(α ))m+10.

II. Пусть теперь α является корнем характеристического уравнения

кратности r ≥ 1.

Тогда F(α ) = F(α ) = ... = F(r1) (α ) = 0 ,

F (r) (α ) 0 . Формула

(4.47)

показывает, что в этом случае L(eα x xm ) есть произведение eαx

на многочлен

степени m r. Чтобы получить в результате

подстановки в

левую

часть

уравнения eαx, умноженное на многочлен степени m, естественно искать частное решение в этом случае в виде

y

* = xrQ

(x)eα x = eα x (q

xm+r + q xm+r1

+... + q xr ).

(4.60)

m

m

m1

0

 

Подставляя это выражение в уравнение (4.55) и требуя, чтобы (4.60) было решением уравнения, мы приходим к условию

 

 

 

e−α xL(xrQ (x)eα x )= P (x).

(4.61)

 

 

 

 

 

 

m

m

 

Снова вычисляем левую часть,

пользуясь формулой (4,47) и помня, что

F(α ) = F(α ) = ... = F(r1) (α ) = 0 ,

F (r) (α ) 0 . Имеем:

 

e−α xL(xrQ (x)eα x )= q

{xmCr

 

F

(r) (α ) + Cr+1

xm1F(r+1) (α ) + ... + F(m+r)

(α )}+

 

m

m

m+r

 

m+r

 

 

+ q

{Cr

xm1F(r) (α ) + Cr+1

xm2F(r+1) (α ) + ... + F(m+r1) (α )}+ ...

 

m1 m+r1

 

m+r1

 

 

... + q

{Cr xF(r) (α ) + F(r+1) (α )}+ q F(r) (α ).

 

 

1

r+1

 

 

 

 

0

 

 

96

Подставляя это выражение в равенство (4.61) и приравнивая после этого коэффициенты при одинаковых степенях x в обеих частях равенства (4.61), снова получаем систему m + 1 уравнений для определения qm, qm–1, ...,

q0:

 

r

 

 

 

(r )

(α )

= pm

 

 

 

 

 

 

 

 

 

 

Cm+r qm F

 

 

 

 

 

 

 

 

 

 

 

 

qm1Cmr +r 1F (r ) (α ) + qmCmr ++1r F (r +1) (α ) = pm1

 

 

 

 

 

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C r

 

F (r )

(α ) + q

C r +1

 

F (r +1)

(α ) + ... + q

 

C r +l

F (r +l ) (α ) = p

 

 

q

 

 

 

 

 

 

(4.62)

 

ml

m+r 1

 

 

 

 

ml +1 m+r l +1

 

 

 

m

m+r

 

ml

 

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q0 F

(r )

(α ) + q1F

(r +1)

(α

) + ... + qm F

(r +m)

(α )

= p0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определитель системы (4.62):

Cr + Cr + − ... Cr [F (r) (α )]m+1 ≠ 0 ,

m r m r 1 r

поэтому все неизвестные qi (i = 0, 1, 2, …, m) определяются однозначно, и мы получаем решение вида (4.60).

Итак, мы приходим к следующему результату: частное решение линейного уравнения с постоянными коэффициентами и с правой частью

специального вида Pm (x)eα x может быть найдено в виде xrQm(x)eα x , где r ≥ 0, есть кратность корня α характеристического уравнения, Qm есть многочлен той же степени, что Pm .

На практике для нахождения частного решения обычно записывают его в форме (4.56) или (4.60) с неопределенными коэффициентами у многочлена Qm (x) . Подставляя это выражение частного решения в заданное уравнение, сокращая на eαx и приравнивая коэффициенты при одинаковых степенях x, получают систему линейных уравнений для этих коэффициентов. По доказанному, эта система всегда имеет определенные решения.

Пример 9. Решить уравнение

y′′′ + y′′ = x2 +1+ 3xex .

Согласно теореме 3 (см. стр. 76) будем искать частные решения двух

уравнений: y′′′ + y′′ = x2 +1 и y′′′ + y′′ = 3xex .

Характеристическое

уравнение,

очевидно, имеет вид k 3 + k 2 = 0 , его корни:

k1 = k2 = 0 , k3 = −1 .

Рассмотрим

сначала первое уравнение; в правой части нет показательного множителя, следовательно, α = 0; но нуль есть двукратный корень характеристического уравнения; поэтому мы, согласно изложенному выше, должны искать частное решение в виде y1* = x2Q2 ( x) = a2 x4 + a1 x3 + a0 x2 . Находим:

y1*= 4a2 x3 + 3a1x2 + 2a0 x,

97

 

 

 

 

 

 

y*

= 12a x2

+ 6a x + 2a ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

2

 

 

 

 

 

1

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y′′′*

=

24a

2

x

+

6a x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

1 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя в уравнение, получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

24a x + 6a +12a x2

+ 6a x + 2a = x2

+1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

1

 

 

2

 

 

 

 

1

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Приравнивая коэффициенты при равных степенях x, находим систему

уравнений: 12a2 = 1;

 

24a2 + 6a1 = 0;

 

6a1 + 2a0= 1.

Из

нее

 

определяем

коэффициенты:

a =

1

 

a =

1

a =

3

 

 

следовательно,

Y =

1

 

4

1

 

3

+

3

2

 

 

 

;

 

;

 

;

 

 

x

 

 

 

x

 

 

x

.

 

 

 

 

 

 

 

 

 

 

 

 

2

 

12

 

1

3

0

 

2

 

 

 

 

 

 

1

12

 

 

 

3

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Переходим ко второму уравнению; здесь α = 1 не является корнем характе-

ристического уравнения. Ищем частное решение в виде

y*

= ex (b x + b )

.

 

2

 

 

 

 

 

 

1

0

Находим:

 

 

 

 

 

 

 

 

 

 

y*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= ex (b x + b + b )

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

1

 

0

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y*

= ex (b x + b + 2b ) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

1

 

0

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y*(3) = ex (b x + b + 3b ) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

1

 

0

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

подставляя

 

в

уравнение

и

сокращая

 

на

 

ex,

получаем

 

 

 

b1x + b0 +3b1 +

+ b1x + b0 +2b1 = 3x.

 

 

Приравниваем

коэффициенты: 2b1 = 3;

 

2b0 + 5b1 = 0,

 

 

3

 

 

 

15

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

3

 

 

 

15

 

 

 

 

 

откуда b =

 

;

b = −

 

 

.

Искомое частное решение: Y2 = e

 

 

 

 

x

 

 

 

 

 

. Общее

2

4

2

4

 

 

1

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

решение данного уравнения:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = C ex

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

1

 

 

1

 

 

 

 

3

 

 

 

15

 

 

 

 

 

 

 

 

 

 

 

+ C

2

+ C

3

x + +

 

 

 

x

2 + −

 

x3

+

 

 

 

 

x

4 + e x

 

x

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

3

 

 

12

 

 

 

2

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C1, C2, C3 – произвольные постоянные.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим уравнение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ly = eα x (P(1) cosβ x + P(2) sinβ x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.63)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m1

 

 

 

m2

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку

 

 

 

eiβ x + e

iβ x

 

 

 

 

 

 

eiβ x e

iβ x

 

 

 

уравнение

 

 

 

 

 

 

 

может

cos β x =

 

 

 

 

 

 

 

 

 

 

, sin β x =

 

 

 

 

 

 

 

,

(4.63)

 

 

 

 

 

2

 

 

 

 

 

 

2i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

быть записано в виде

 

Ly =

1

e (α

 

 

i) x {P (1) iP (2) }+

1

e (α

−β i) x

 

{P (1)

+ iP (2) }

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

m1

 

 

m 2

2

 

 

 

 

 

 

 

 

 

 

 

m1

 

 

m 2

.

Обозначим

m = max{m1;m2}. По предыдущему,

мы должны искать частное

решение в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e(α

i)xQ(1)

(x) + e(α

−β i)xQ(2) (x) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.64)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

если

α ± iβ

не являются корнями характеристического уравнения;

если же

α iβ

есть корни характеристического уравнения кратности r, то выражение

(4.64) должно быть умножено на xr. Легко видеть, что из системы уравнений

98

(4.59) и (4.62) для определения коэффициентов многочлена Qm(1)

соответствующая система для Qm(2) получается переходом к комплексным сопряженным значениям коэффициентов уравнений; следовательно,

коэффициенты многочлена Qm(2) окажутся комплексными сопряженными

с соответствующими коэффициентами многочлена

Qm(1) . Поэтому, отделяя

действительную часть от мнимой, мы получим: если

Qm(1) (x) = Qm* (x) + iQm**(x) ,

то Qm(2) (x) = Qm* (x) iQm**(x) . Подставляя эти многочлены в выражение (4.64) и переходя снова от показательных функций к тригонометрическим, находим искомое частное решение:

y* = eα x (2Q

(1) (x)cosβ x 2Q(2) (x)sin β x).

 

 

 

m

m

 

 

Это выражение уже не содержит комплексных величин.

 

 

В случае, если

α ± iβ

есть корни

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

уравнения

кратности r, предыдущее выражение надо умножить на xr.

 

 

Итак, имеем окончательно: частное

решение линейного

уравнения

с постоянными коэффициентами и с правой частью специального вида

 

f (x) = eα x {P(1) (x) cosβ x + P(2) (x)sin β x}

 

 

 

 

m

m

 

 

 

 

1

2

 

 

можно найти в форме

xr eα x (Qm(1) (x) cosβ x + Qm(2) (x)sin β x) , где Qm(1) и Qm(2)

многочлены степени

m = max{m1;m2}, а r

0 есть кратность корня α ±

iβ

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

 

 

 

На практике опять пишут многочлены Qm(1) и Qm(2) с неопределенными

коэффициентами, подставляют в уравнение и приравнивают коэффициенты в

обеих частях при выражениях вида xl cosβ

x и xl sinβ x (l = r, r +1,..., r + m).

Пример 10. Найти общее решение уравнения

 

 

 

 

y′′ + 2y′ + 5y = 2cos x.

 

Рассмотрим однородное уравнение y′′ + 2 y′ + 5y = 0.

 

Составим

характеристическое уравнение

k 2

+ 2k + 5 = 0

оно имеет корни

 

,

k1,2 = −1 ± 2i, поэтому

 

 

 

 

 

 

 

 

= ex (C cos2x + C sin2x).

 

 

y

 

 

 

 

 

1

 

2

 

 

Найдем y*. В нашем случае

f (x) = 2cos x. Сравним с

f (x) в уравнении

(4.63):

α = 0; β = 1. Числа

α ± β

i = ±i

не являются корнями

характеристического уравнения. Поэтому

 

 

 

 

y* = Acos x + Bsin x .

99

Подставляя y* в исходное уравнение и приравнивая коэффициенты при

cos x и при sin x, получим A =

2

,

B =

1

, следовательно,

 

 

5

 

5

 

y* = 2 cosx + 1 Bsin x. 5 5

Таким образом, общее решение уравнения таково:

y = y* + y = ex (C1 cos 2x + C2 sin 2x) + 2 cos x + 1 sin x. 5 5

Пример 11. Решить уравнение

y′′ + 4 y = cos 2x.

Составим характеристическое уравнение

k 2 + 4 = 0.

Его корни k1,2 = ±2i. Поэтому y = C1 cos 2x + C2 sin 2x. Правая часть имеет вид

f (x) = cos2x , т.е. α

= 0,β = 2 . Числа α

± β i = ±2i являются простыми корнями

характеристического уравнения. Поэтому

 

 

 

 

 

y* = x(Acos2x + Bsin2x).

 

 

 

 

Подставим y* в исходное уравнение и приравняем коэффициенты при

cos 2x и sin 2x .

Получим

A = 0,

B =

1

. Следовательно

y* =

1

x sin 2x.

 

 

Окончательно

 

 

4

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

y = C

cos2x + C

sin 2x +

1

xsin 2x

.

 

 

 

 

 

 

 

 

 

1

2

4

 

 

 

 

 

 

 

 

 

 

 

 

4.6.Уравнения, приводящиеся к уравнениям

спостоянными коэффициентами

Рассмотрим теперь некоторые типы линейных уравнений с переменными коэффициентами, которые могут быть приведены к уравнениям с постоянными коэффициентами с помощью замены независимого переменного. Отметим, что уравнение при такой замене остается линейным (подробно этот вопрос изложен в книге [1], глава V).

I. Уравнение Эйлера.

 

xn y(n) + a1xn1y(n1) +...+ an1xy′ + an y = 0 ,

(4.65)

где a1, a2, …, an – действительные постоянные.

 

Пусть x

(0;+ ∞ ). В этом случае уравнение (4.65) удовлетворяет всем

условиям теоремы Коши – Пикара. Покажем, что замена независимой x = et преобразует уравнение (4.65) в уравнение с постоянными коэффициентами.

100

Соседние файлы в папке книги