340 ':1. О-. Диффереrщирльные уравнения
или, в ма1т!.qн:>й фор:. е,
X(t) = A(l)X\t),
где
|
|
В обJ1асти непрерьшности иоэффицпентсв aij (t), |
i, j |
1, |
|
|
. , |
п, си |
стема (12) удовлетворяет условиям тсоремь.t существования и единствен |
ности решения задачи :Коши. |
|
|
|
|
|
|
|
|
|
|
|
|
называется |
|
|
Фундаментальной систе.мой решет.tй системы (12) |
совоиупность(k)произвольных |
|
п линейно независимых решений Xk (t) = |
- |
( |
Х |
(k) ( |
t |
) |
, Х2 |
|
( |
t |
|
(k) ( |
) |
) |
т, k - 1, 2, |
. |
. . |
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
) , . . . , Xn t |
|
|
|
|
|
|
|
|
|
|
|
|
|
- |
|
Если Xk(t), |
k = 1, 2, . . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, |
п, |
- |
фундаментальнаяn система решений |
системы (12), то общее решение имеат вид X(t) |
|
|
C |
|
|
), |
.где |
|
С1 , |
|
|
2:: k=Xk(t. |
|
|
|
С2 |
, . . ., Сп - произвольные постоянные. |
|
|
|
|
= k=1 |
|
|
|
|
|
|
|
|
|
|
Интегрирование систеыы (12) обычно пр•>Dодится методом исилюче |
ния (см. |
пример 3) . |
|
|
|
|
|
|
|
|
|
енциальных урав |
|
|
|
ий |
|
|
|
Р |
е |
шит |
ь |
с |
истемы линей |
ных диффер |
нен |
: |
|
|
|
|
|
|
|
|
|
. |
|
|
|
, n. |
|
|
|
|
|
|
|
|
|
10.427. |
dy |
|
|
у |
|
|
dz |
|
|
|
2у |
|
|
z |
|
|
|
|
|
|
|
|
|
|
|
|
|
- = - - |
+ xz, |
|
|
- = -- + |
- . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dx |
|
|
х |
|
|
|
|
dx |
|
|
|
х3 |
|
|
х |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dy |
= -у + zx, |
|
|
dz |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
10.428. |
х dx |
|
х2 dx |
= -2у + zx. |
|
|
|
|
|
|
|
|
|
|
|
|
10.429. х = -,t |
iJ = _ :t: . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
10.430. |
. |
|
|
|
2 |
у. = у + |
t + 2 |
х. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
х = - -х, |
--t |
|
|
|
|
|
|
н о э ф ф и ц и е н т а |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В частном случае с и с т е м с |
п о с т о я н н ы м и |
м и, |
ногда матрица A(t) в правой части (13) |
|
не зависит от |
t, |
|
для отыс |
иания фундаментальной системы решений Xk (t), |
|
k |
= 1, 2, |
. |
., |
п, |
могут |
быть использованы методы линейной алгебры. |
|
|
|
|
|
|
|
|
|
|
|
|
|
Из хараитеристичесиого уравнения |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
det. (А - |
>..Е) = О |
|
|
|
|
|
|
|
|
|
|
|
(14) |
на |
ходятся различные иорни Л |
1 , Л2, |
|
|
Лs |
и для всяиого иорня |
|
(с уче |
том его нратности) определяется соответствующее, |
ему частное. |
решение>.. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
• |
• • |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
342 |
Гл. 10. Дифференциальные уравнения |
Для нахождения частного решения нонстанты С1, С2, Сз определяем из следующей системы:
х (о) ( ) с, Н)
отнуда С1 |
= |
1, |
С2 |
2, С3 |
|
|
(::3=: ) |
|
решения получаем |
|
= |
X(t) |
|
|
|
|
= |
x |
(t) |
|
|
|
|
Ш + с, О)
= |
( - g : |
g:) ' |
|
7С1 + С2 |
+ 5Сз |
Онончательно для исномого частного
et + (2) e4t + (15) еы. t>
|
б) Л - номпленсный норень нратности 1. |
Тогда норнем харанте |
ристичесного уравнения (14) является танже сопряженное с ). число Л. |
Вместо номпленсных частных решений X("J (t) и x(XJ (t) следует взять |
действительные частные решения xf"J(t) |
= |
Re X("J (t) |
и |
X "J(t) |
= |
= |
lm X(Л)(t). |
= |
|
|
|
|
П р и м е р 11. Найти общее решение системы |
|
|
|
X1 (t) |
= |
Х1 + Х2, |
|
|
|
|
|
±2(t) |
-2х1 + Зх2. |
|
|
|
<J Харантеристичесное уравнение |
|
|
|
|
|
имеет номпленсно сопряженные норни >.1, |
2 |
2 ± i. |
Для нахождения соб |
ственного вентора, соответствующего норню Л |
|
2 |
+ i, получаем систему |
|
|
( |
- 1 |
i)y ,\) |
+ |
= |
|
у2(,\) - |
о |
' |
|
|
|
--2у,\) |
+ |
|
|
|
= |
Полагая |
(,\) |
|
(1 - i)y ,\) |
= |
|
|
|
(,\) |
|
1 + |
i, |
= |
о . |
|
у1 |
= 1, находим у2 |
|
|
. |
т. е. |
|
|
§ 3.
Системы дифференциальных уравнений
Отсюда пара действительных частных решений имеет следующий вид: |
|
1 |
(t) |
- |
R |
е (( |
|
|
) е |
|
|
) |
|
( |
|
|
e2t |
cos t |
|
|
|
) |
- |
( |
|
cos t |
) |
е |
2t |
|
х(Л) |
|
|
|
1 |
|
(2+i)t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
- |
I |
|
|
l + i |
|
|
|
|
- |
|
|
2 |
t(cos t - sin t) |
|
- |
|
cos t |
- sin t |
|
|
|
' |
|
|
|
|
|
|
|
|
|
|
|
|
e |
|
|
|
) |
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
2t |
. t |
|
|
|
|
|
|
|
. |
|
t |
|
|
|
х2(Л) |
(t) - |
m (( |
1 |
+ i |
) е(2+i)t) |
- |
( |
2 |
е SШ |
|
|
|
) |
- ( |
sш |
|
е |
2t· |
|
|
|
|
|
|
|
|
|
|
|
|
e |
t(cos t + sin t) |
|
|
|
cos t + sin t |
|
|
|
|
Оиончательно (см. формулу (15)) |
получаем общее решение |
|
|
|
|
|
|
х |
(t) |
= |
с |
( |
cos t |
|
|
2t |
+ |
С2 |
( |
|
sin t |
) |
е |
2t |
= |
|
|
|
|
|
|
|
|
|
|
1 |
|
|
cos t - sin t) е |
|
|
|
cos t + sin t |
|
|
|
|
|
|
|
|
|
· |
!> |
|
|
|
|
|
|
|
|
|
|
|
|
= ((C1 |
+ C2) cost + |
(C2 - C1) sint) е2t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
С1 |
cos t + С2 |
sin t |
|
|
|
|
|
|
|
в) ,\ - иорень иратности r 2. Соответствующее этому иорню ре шение системы (13) ищетсн в виде веюора
|
а |
(1) |
|
(a.n. . |
|
|
1 |
|
|
(!) |
|
2 |
|
(!} |
|
|
иоэффициенты иоторого a j) , |
i = 1 |
, |
. " п ; j = 1 |
, |
. . " r, определнютсн |
из системы линейных уравнений, получающейсн приравниванием ноэф |
фициентов при одинаиовых степенпх t |
в результате подстановни вектора |
(16) |
исходную систему (13). |
|
|
|
|
системы |
|
|
Пвр и м е р 12. Найти общее решение. |
|
|
|
1 |
|
4 |
± 1 |
(t) |
= |
2х1 - х2, |
|
|
|
|
:i:2 |
|
1 |
|
+ 6х2. |
|
|
|
|
|
|
(t) = 4х1 |
|
|
<] Харантеристичесиое уравнение |
|
|
|
|
|
|
|
|
2 |
- |
|
-1 |
|
= (Л - 4)2 = О |
|
|
|
|
|
Л |
6 - ,\ |
|
|
|
имеет норень ,\ = 4 кратности r = 2. Поэтому ищем решение системы в виде
344 |
Гл. 10. Дифференциальные уравнения |
Подставляем это выражение в исходную систему и соиращаем на e4 t :
Приравнивая коэффициенты при одинаиовых степеннх t, получаем:
/31 + 2а1 + а2 = О,
/32 - 4а1 - 2а2 = О,
2/31 + /32 = о,
-2/32 - 4/31 = о.
Полагая а1 = С1 и /31 = С2 , имеем /32 = - 2С2 и а2 = -2С1 - С2 . Таиим образом, общее решение системы имеет вид
Решить |
следующие |
|
системы линейных |
диф ференциальных |
уравнений |
с |
постоннными коэффициентами. |
Там, где даны на |
чальные условия, кроме общего решения, найти соответствующее |
частное решение: |
|
|
|
|
|
|
|
|
|
|
10.431. |
х |
= у, |
у = -2х + Зу. |
|
|
|
|
|
|
10.432. |
х |
= х + Зу, |
у = -х + 5у; |
х (О) = З, |
у (О) = 1 . |
|
10.433. |
х |
= Зх - |
2у, у = 4х + 7у; х (О) = 1 , |
у(О) = О. |
|
|
|
у = 5х - 6у. |
|
|
|
|
|
10.434. |
i; |
= 2х - 5у; |
|
|
|
|
|
|
10.435. |
х |
= х - 4у, |
у = х - Зу. |
|
|
|
|
|
10.436. |
х |
= |
-х |
+ |
2у, |
у = -2х - 5у; х (О) |
= |
О, у(О) |
= 1 . |
|
10.437. |
х |
= у, |
у = z, |
i = х; |
х (О) |
= у (О) |
= z(O) = 1 . |
|
10.438. |
х = у + z , |
у = z + х, i |
= х + у; |
х (О) = у (О) |
= 2 , |
z(O) = - 1 . |
i; |
= х - 2 у - z , у = - х + у + z , z = х - z . |
|
|
10.439. |
|
|
10.440. |
х |
= 5 х + 2у - Зz, у = 4х + 5у - 4z, |
z = 6 х + 4у - 4z. |
5. Линейные неоднородные системы. Нормальная линейная неодно- |
родная система дифференциальных уравнений имеет вид |
|
|
|
±1 |
= al l (t)x1 + a12 (t)x2 |
+ . . . + aln (t)xn + fi (t) , |
|
(17) |
|
±2 |
= |
а21 (t)x1 + a22 (t)x2 |
+ . . . + a2n (t)x,,; + /2 (t) , |
|
§ 3. Системы дифференциальных уравнений |
345 |
где по ирайней мере одна из фуниций fk(t) не равна тождественно нулю. В матричной форме система (17) имеет вид
где F(t) = |
|
|
|
X(t) = A(t)X(t) + F(t), |
(18) |
(f1 (t), |
f |
(t), |
. . ., fn(t)) т. Интегрирование системы (17) |
мож |
|
|
2 |
|
|
|
но проводить методом исилючения (см. пример 3), однаио иногда предпо |
чтительнее найти предварительно решение Xo(t) соответствующей (18) |
однородной системы |
= A(t)X(t) |
(19) |
X(t) |
и иаиое-либо частное решение X(t) системы (18) . Тогда общее решение системы (18) имеет вид
X(t) = X0(t) + X(t). |
|
k = 1, 2, |
(20) |
Если известна фундаментальная система Xk(t), |
. . ., п, ре |
шений однородной системы (19), то общее решение X(t) можно найти |
методом вариации произвольных постоянных. |
Именно, полагая |
n |
|
|
|
(21) |
k=1 |
|
|
|
X(t) = L:·Ck |
(t)Xk(t), |
|
|
|
определяем фуниции Ck (t) подстановиой (21) в систему (18) . |
Учитывая |
при этом равенства |
k = 1, |
2, . . . |
, п, |
|
Xk(t) - A(t)Xk (t) = о, |
|
приходим и системе уравнений относительно Ck (t): |
|
|
n |
|
|
|
(22) |
L Ck(t)Xk(t) = F(t). |
|
|
k=1 |
|
|
|
|
Из этой системы находим Ck (t) = 'Pk(t) и, интегрируя, получаем фуни
ции Ck (t) с точностью до произвольных постоянных. Подставляя их в
(21), получаем исиомое общее решение неоднородной системы (18) . Пр и м е р 13. Зная фундаментальную систему решений 
однородной системы
:Х1 :Х2
346 |
Гл. 10. Дифференциальные уравнения |
найти общее решение неоднородной системы
Х1 = 6Х1 + Х2 + t,1.
±2 = .5х1 + 2х2 +
<J Воспользуемся методом вариации произвольных постоянных. Для
фующий |
C\ (t) |
и |
(\(t) |
составим систему вида (22) |
|
|
(1)1 |
е |
7t |
+ . |
(t) |
5 |
t |
= |
(1t) · |
|
|
|
С. |
1 (t) |
|
C2 |
(-1) е |
|
|
|
|
|
|
|
|
Найдя |
с. (t) |
ы |
6 |
|
-7t |
|
C1 (t) = - i(425 |
|
+ 1 |
|
|
t + 492 ) |
е-7t |
и проинтегрировав, |
получим |
|
|
|
= --е |
|
,
+С1 , 
Таним образом, общее решение системы запишется в виде
i, j |
Если ноэффициенты aij (t) системы (17) |
постоянны, т. е. a;j ( |
t) |
= |
aij, |
= |
1, |
|
|
а |
фующии fi(t) имеют вид произведений |
|
|
|
|
|
|
(P(t) cos {Зt + Q(t) sin {Зt)e°'t, |
X(t) |
|
(23) |
|
Следует. . " n, |
|
|
|
|
|
|
|
|
|
|
|
|
где P(t) и Q(t) - многочлены, то частное решение |
X(t) можно найти |
методом неопределенных ноэффициентов, записав |
|
|
|
в виде, аналогич |
ном |
(23), с учетом совпадения или несовпадения чисел о : ±i{З с норнями |
харантеристичес1юго уравнения. |
|
наибольшая степень многочленов |
|
|
и |
|
иметь в виду, что если k |
P(t) |
Q(t) |
в |
(23) и |
Л |
|
i{З - |
- |
|
|
|
|
|
харантеристиче- |
|
|
о : + |
норень нратности |
|
сного уравнения, |
то |
частное |
решение X(t) |
ищется в |
виде |
|
|
|
|
|
|
|
|
|
|
= |
|
|
r |
|
|
|
|
')'1otk+l |
|
20t |
I |
. |
? . . |
|
/notk+I |
|
+')'11tk +
2.1 .-:.
+/nl tk +
'. ' . ' . : :
·· · + 'Yn,
§3. Системы дифференциальных уравнений
П р и м е р 14. Найти частное решение системы
|
|
|
Х1 = -Х2 + t2 , |
|
|
|
<J |
|
|
Х2 = Х1 + et . |
= 1 |
|
|
Та:к :каи хараитеристичесиое уравнение |
с |
= О имеет иорни |
Л1 |
, 2 |
= |
±i, ищем частное решение системы1-/в |
уммы многочлена |
|
виде |
|
второй степени и фуниции вида Det : |
|
|
|
Подставшш эти фуниции в заданную систему, получим равенства
2A1 t + В1 + D1 et = -A2t2 - B2t - С2 - D2et + t2 , 2A2t + В2 + D2et = Ai t2 + B1t + С1 + D2et + et .
ПриравниваR иоэффициенты при одинаиовых степенRх t и при et , полу чим систему
D1 |
= -D2 , |
l - A2 = О, |
D2 |
= D1 + 1, |
Ai = 0. |
В |
2 |
= С1 = О, |
А |
2 |
= 1, В1 = 2, |
С |
= -2, D |
2 |
|
|
|
|
|
2 |
|
исиомое частное решение имеет вид |
|
|
|
П р и м е р 15. Найти общее решение системы
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X (t) = AX(t) + F(t) , |
где А = |
2 |
- 1 |
|
и F = |
t + 1 |
|
е |
Зt |
. |
1 |
4 |
|
2t |
|
|
|
<J Хараитеристичесиое( ) |
уравнение( ) |
|
|
|
|
1 |
2 - |
л |
- 1 |
= Л |
2 |
- 6Л + 8 + 1 = (Л - 3)2 = О |
|
1 |
4 _ Л 1 |
|
имеет иорень Л = 3· : кратност( и 2). Общее решение однородной системы
ищем в виде X0 (t) = at + (3б e3t , подставив иоторое в однородную it +
348 |
Гл. 10. Дифференциальные уравнения |
систему и со1 ращая на e3t , имеем
Получим систему
З(аt + .В) + .В = 2(at + .В) - (1t + 6) ,
З(1t + 6) + б = at + .В + 4(1t + 6) ,
из ноторой следуют два независимых соотношения а = -1 и .В + а = -6. Полагая а = С1 и .В = С2 , имеем / = -С1 и 6 = - С1 - С2 , т. е.
Таи нан F(t) содержит множитель e3t, причем Л = 3 - норень харанте ристичесного уравнения нратности 2, то ищем частное решение в виде
|
|
2 |
A1 t + В1 |
|
з . |
( |
а не в виде t |
|
2 |
|
2) |
|
|
A |
t + В |
е |
t |
Подставив |
X((t) |
|
в заданную)систему и сонратив на e3t , получаем ма |
тричное равенство |
|
|
|
|
|
ноторое можно записать в .виде равенств
A1 t3 |
+ B1 t2 |
+ D1 t + 3A1 t2 |
+ 2B1t + D1 = -A2 t3 - B2 t2 - D2t + t + 1 , |
-A2 t3 |
- B2t2 |
- D2t + ЗА2t2 |
+ 2B2t + D2 = Ai ti + B1 t2 + Di t + 2t. |
Сравнивая ноэффициенты при одинановых степенях t, получаем систему
уравнений
Ai + А2 = О, В1 + ЗА1 + В2 = О, Di + 2В1 + D2 = 1 ,
D1 = 1,
|
Ai + А2 |
= |
О, |
В1 + В2 |
- ЗА2 |
= О, |
Di |
+ D2 |
- 2В2 |
= -2, |
|
|
D2 |
= |
О. |