книги / Обыкновенные дифференциальные уравнения
..pdfкомплексных корней характеристического уравнения (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 e− x + 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 |
( n−1) |
|
2 |
( n−2) |
|
|
1 ( n−1) |
|
|
|
|||
(uυ ) |
( n−1) |
= |
u |
|
υ + |
Cn−1u |
υ |
+′ |
Cn−1u |
υ |
+ ′′ |
+... υ u |
( n−1) |
, |
|
|
||
|
|
|
( n−1) |
|
1 |
( n−2) |
|
2 |
( n−3) |
|
|
|
|
……..
(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]+ ... + |
υ (n−1) |
|
L |
|
[u]+ |
υ (n) |
L |
[u] |
. |
|||||||||||||
|
|
(n −1)! |
|
n! |
||||||||||||||||||||||
|
|
1! |
1 |
|
2! |
2 |
|
n−1 |
|
|
|
n |
||||||||||||||
Здесь введены обозначения: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
L(y) = y(n) + a1y(n−1) + a2 y(n−2) +...+ an−1y′+ an y, |
, |
|
|
|
||||||||||||||||||||||
1 |
= |
ny(n−1) |
+ |
|
|
− |
1 |
+ |
|
+ |
|
n−2 |
y′ |
+ |
n−1 |
|
|
|
|
|||||||
L (y) |
|
|
(n |
|
1)a y(n−2) |
|
|
... |
|
2a |
|
|
a |
y |
|
|
|
|||||||||
L (y) = n(n −1)y(n−2) |
+ (n −1)(n − 2)a y(n−3) |
+... + 2 1 a |
n−2 |
y , |
|
|||||||||||||||||||||
2 |
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
…
Ln−1 ( 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 |
xm−1L [ekx ]+ |
m(m −1) |
xm−2L [ekx |
]+...+Cm−1xL |
[ekx ]+ L |
|
[ekx ] |
||||||||
|
|
|
||||||||||||||
1 |
|
1 |
|
1 2 |
2 |
|
m |
m−1 |
|
m |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
||||||
или, поскольку, в силу формулы (4.46) |
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
L (ekx ) = ekxF (r) (k) , |
|
|
|
|
|
|
|
|||||
получим: |
|
r |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
kx |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
m |
e |
= |
|
|
|
|
|
|
|
||
|
|
|
|
L x |
|
|
|
|
|
|
(4.47) |
|||||
= ekx {xm F (k ) + Cm1 xm−1F ′(k ) + Cm2 xm−2 F ′′(k ) + ... + Cmm−1xF m−1 (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 −kp−1)x = 0. |
(4.52) |
Но тождество (4.52) невозможно, так как e(k p −k p−1)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) + C3e−x. |
|
|||||
Пример 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 + qm−1xm−1 + ... + 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 |
xm−1F′α( |
+) |
|
...+ |
F (mα) ( |
+)} |
|
||||||
m |
|
|
|
m |
|
|
|
|
|
|
|
|
|
(4.58) |
+qm−1{x |
m−1 |
F (α )+ |
1 |
|
m−2 |
F′α( |
+) |
|
...+ |
F |
(m−1) |
+)} ... |
||
|
Cm x |
|
|
α ( |
|
|||||||||
|
... + q1{xF (α )+ |
F′α( |
)}+ |
q0 Fα ( |
). |
|
|
Приравнивая выражение (4.58) многочлену Pm(x) и отождествляя коэффициенты при одинаковых степенях x, получим m + 1 уравнений с m + 1 неизвестными qm, qm–1, ..., q0:
95
qm F (α ) = pm |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
qm−1F (α ) + qmC1m F ′(α ) = pm−1 |
|
|
|
|
|
|
|
|||||||||
q |
m−2 |
F (α ) + q |
C1 |
|
|
F ′(α |
) + q |
m |
C 2 F ′′(α |
) = |
p |
m−2 |
||||
|
|
m−1 |
m−1 |
|
|
m |
|
|
|
|||||||
... |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
q |
m−r |
F (α ) + q |
m−r |
C1 |
|
|
F ′(α |
) + ... + q |
m |
C r |
F (r) (α ) |
|||||
|
|
+1 |
m−r +1 |
|
|
|
m |
|
|
|
||||||
... |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
q |
0 |
F (α ) + q F ′(α ) |
+ q |
2 |
F ′′(α |
) + q F ′′′(α ) + ... + q |
m |
F (m) |
||||||||
|
|
1 |
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(4.59) |
|
|
|
= pm−r |
|
|
|
|
|
|
|
|
|
|
|
(α ) = 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+1≠ 0.
II. Пусть теперь α является корнем характеристического уравнения
кратности r ≥ 1.
Тогда F(α ) = F′(α ) = ... = F(r−1) (α ) = 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+r−1 |
+... + q xr ). |
(4.60) |
m |
m |
m−1 |
0 |
|
Подставляя это выражение в уравнение (4.55) и требуя, чтобы (4.60) было решением уравнения, мы приходим к условию
|
|
|
e−α xL(xrQ (x)eα x )= P (x). |
(4.61) |
||||
|
|
|
|
|
|
m |
m |
|
Снова вычисляем левую часть, |
пользуясь формулой (4,47) и помня, что |
|||||||
F(α ) = F′(α ) = ... = F(r−1) (α ) = 0 , |
F (r) (α ) ≠ 0 . Имеем: |
|
||||||
e−α xL(xrQ (x)eα x )= q |
{xmCr |
|
F |
(r) (α ) + Cr+1 |
xm−1F(r+1) (α ) + ... + F(m+r) |
(α )}+ |
||
|
m |
m |
m+r |
|
m+r |
|
|
|
+ q |
{Cr |
xm−1F(r) (α ) + Cr+1 |
xm−2F(r+1) (α ) + ... + F(m+r−1) (α )}+ ... |
|
||||
m−1 m+r−1 |
|
m+r−1 |
|
|
||||
... + 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 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
qm−1Cmr +r −1F (r ) (α ) + qmCmr ++1r F (r +1) (α ) = pm−1 |
|
|
|
|
|
||||||||||||||
... |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
C r |
|
F (r ) |
(α ) + q |
C r +1 |
|
F (r +1) |
(α ) + ... + q |
|
C r +l |
F (r +l ) (α ) = p |
|
|
||||
q |
|
|
|
|
|
|
(4.62) |
||||||||||||
|
m−l |
m+r −1 |
|
|
|
|
m−l +1 m+r −l +1 |
|
|
|
m |
m+r |
|
m−l |
|
||||
... |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
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 e− x |
|
|
|
|
|
|
|
|
|
|
|
|
|
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, поэтому |
|
|
|
|
|
|||
|
|
|
= e−x (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 = e−x (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) + a1xn−1y(n−1) +...+ an−1xy′ + an y = 0 , |
(4.65) |
где a1, a2, …, an – действительные постоянные. |
|
|
Пусть x |
(0;+ ∞ ). В этом случае уравнение (4.65) удовлетворяет всем |
условиям теоремы Коши – Пикара. Покажем, что замена независимой x = et преобразует уравнение (4.65) в уравнение с постоянными коэффициентами.
100