Лекции по обыкновенным дифференциальным уравнениям Ржавинская Белякова Жаркова
.pdfæç y1(t)
Пусть Y (t) = ç ...
çè yn (t)
ö
÷÷ - решение системы (4.6), функции fi (t, y1(t),..., yn (t)), i =1,..., n имеют
÷
ø
непрерывные частные производные по всем аргументам до порядка n – 1 включительно. Подставляем Y(t) в (4.6), получим n тождеств:
ì dy1 |
= f (t, y (t),..., y |
n |
(t)) |
||||||
ï |
dt |
1 |
1 |
|
|
||||
ï |
|
|
... |
|
|
|
|
||
í |
dyn |
|
|
|
|
|
|
||
ï |
= f |
n |
(t, y (t),..., y |
n |
(t)) |
||||
|
|||||||||
ï |
dt |
|
1 |
|
|
||||
î |
|
|
|
|
|
|
|
в частности,
|
|
|
|
|
|
|
|
dy1 = f |
(t, y (t),..., y |
n |
(t)) . |
(4.12) |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
dt |
|
1 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Продифференцируем (4.12) по t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
d |
2 y |
|
¶f |
|
|
¶f dy |
|
|
|
|
|
|
|
|
¶f |
|
|
dy |
n |
|
|||||||||||
|
|
|
|
|
|
1 |
= |
1 |
+ |
1 |
|
1 |
+ ...+ |
|
|
1 |
|
|
|
|
= |
|||||||||||||
|
|
|
|
dt2 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
¶t |
|
|
¶y1 dt |
|
|
|
|
|
|
|
|
¶yn dt |
|||||||||||||||||
|
|
|
|
|
∂f1 |
|
n |
∂f1 |
|
dyk |
|
|
|
∂f1 |
|
|
|
|
|
n |
∂f1 |
|
|
|||||||||||
|
|
= |
+ |
å |
|
|
= |
+ å |
fk |
|||||||||||||||||||||||||
|
|
∂yk dt |
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
∂t |
|
k=1 |
|
|
|
|
∂t |
|
|
|
|
|
k =1 |
∂yk |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d 2 y |
|
¶f |
|
n |
|
¶f |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
1 |
= |
|
1 + å |
|
1 |
|
|
|
fk . |
|
|
|
|
(4.13) |
|||||||||||
|
|
|
|
|
|
|
|
dt |
2 |
|
¶yk |
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
¶t |
k =1 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
Обозначим правую часть в (4.13) через F2(t, y1, …, yn): |
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
¶f |
|
|
|
|
n |
|
|
¶f |
|
|
|
|
|
|
||
|
|
|
|
|
F2 (t, y1,K, yn ) = |
|
|
|
1 |
+ |
å |
|
1 |
fk . |
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
¶t |
|
|
k =1 |
¶y |
k |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
Дифференцируем (4.13) по t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
d3 y |
|
¶F |
|
|
n |
|
¶F dy |
|
|
¶F |
|
|
|
n |
¶F |
||||||||||||||||||
|
|
1 |
|
= |
|
|
2 |
+ å |
2 |
|
|
|
k |
|
= |
|
|
|
|
2 |
+ |
å |
|
2 fk . |
||||||||||
|
dt |
3 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
¶t |
|
k =1 |
¶yk dt |
|
|
|
|
¶t |
|
k =1 |
¶yk |
83
PDF created with pdfFactory Pro trial version www.pdffactory.com
Обозначим
F3(t, y1,K, yn ) = |
¶F |
n |
¶F |
fk |
|
2 |
+ å |
¶y |
2 |
||
|
¶t |
k =1 |
k |
|
|
|
|
|
|
и т.д. n – 1 раз.
Получили n – 1 тождество
ì |
|
|
dy1 = f (t, y (t),..., y |
n |
(t)) |
|||
ï |
|
|
dt |
1 |
1 |
|
|
|
ï |
|
|
|
|
|
|
|
|
|
d 2 y |
|
|
|
|
|
||
ï |
|
|
1 |
= F (t, y (t),..., y |
n |
(t)) |
||
|
|
|
||||||
í |
|
dt2 |
2 |
1 |
|
|
||
|
|
... |
|
|
|
|||
ï |
|
|
|
|
|
|
|
|
ïd n−1y |
|
|
|
|
|
|||
ï |
|
|
1 |
= Fn−1 |
(t, y1(t),..., yn (t)) |
|||
|
|
n−1 |
||||||
ï |
dt |
|||||||
î |
|
|
|
|
|
|
и еще одно уравнение
d n y1 = Fn (t, y1(t),..., yn (t)) . dtn
(4.14)
(4.15)
Систему (4.14) можно рассматривать как систему, неявно задающую функции y2(t), y3(t), …, yn(t). Предположив, что в рассматриваемой области якобиан
D( f1, F2, F3,..., Fn−1) ¹ 0 , систему (4.14) можем разрешить относительно y2, y3, …, yn, выра-
D( y2, y3,..., yn )
|
dy |
|
d n−1y |
|
|
зив их через переменные t, y , |
1 |
,..., |
1 |
. Подставив эти выражения в уравнение (4.15), |
|
dtn−1 |
|||||
1 |
dt |
|
|
получим дифференциальное уравнение порядка n относительно y1(t):
d n y |
|
dy |
|
d n−1y |
|
1 |
= F(t, y , |
1 |
,..., |
1 |
) . |
dtn |
|
||||
1 |
dt |
|
dtn−1 |
Справедливо следующее утверждение, которое мы приводим без доказательства.
Теорема 4.2. Пусть дана система дифференциальных уравнений
dyi = f |
(t, y ,..., y |
n |
), i =1,...,n . |
(4.16) |
|
dt |
i |
1 |
|
|
|
|
|
|
|
|
Функции fi (t, y1,..., yn ) имеют непрерывные частные производные по всем аргументам в некоторой области
G :t0 - a £ t £ t0 + a; yi0 - bi £ y £ yi0 + bi , i =1,..., n ;
J = D( f1, F2 , F3,..., Fn−1) ¹ 0 в области G.
D(y2 , y3,..., yn )
Тогда следствием системы (4.16) является уравнение
84
PDF created with pdfFactory Pro trial version www.pdffactory.com
d n y |
|
dy |
|
d n−1y |
|
|
1 |
= F(t, y , |
1 |
,..., |
1 |
), i =1,..., n . |
(4.17) |
dtn |
|
|||||
1 |
dt |
|
dtn−1 |
|
Обратное утверждение. По любому решению y1(t) уравнения (4.17) можно подобрать
æç y1(t) n – 1 функцию y2(t), …, yn(t), такую, что вектор-функция Y (t) = ç ...
çè yn (t)
системы (4.16).
(Доказательство см., например, [1, гл. 3, § 2].)
ö
÷÷ является решением
÷
ø
ì dx |
= x - y . |
|
В примере 4.2 íï dt |
||
ï |
dy |
= y - 4x |
î dt |
|
В обозначениях теоремы 4.2: f1 = x – y, f2 = y – 4x и J = ¶f2 = -4 ¹ 0 |
Þ |
систему (4.7) |
||
|
|
|
¶x |
Теор. 4.2 |
|
можно |
свести к одному уравнению второго порядка относительно y(t) |
(аналогично |
|||
~ |
= |
¶f |
= -1 ¹ 0 и систему (4.7) можно свести к одному уравнению относительно x(t)). |
||
J |
1 |
||||
|
|
¶y |
|
|
|
4.3. Системы линейных дифференциальных уравнений
Определение 4.2. Система дифференциальных уравнений называется линейной, если уравнения системы линейны относительно всех компонент неизвестной вектор- функции и их производных, т.е. имеют вид
dyi |
|
n |
|
|
+ f (t), i =1,...,n . |
(4.18) |
= |
å |
a (t)y |
|
|||
dt |
|
i j |
j |
i |
|
|
|
|
|
|
j=1
Обозначим
æ y (t) ö |
|
||
ç |
1 |
÷ |
, |
Y (t) = ç |
L |
÷ |
|
ç |
|
÷ |
|
è yn (t)ø |
|
dY dt
çdy1(t)
çdt
=ç L
çdyn (t)
çè dtæ
ö |
æ |
f (t) ö |
|
|
÷ |
|
|||
÷÷ , |
ç |
1 |
÷ |
; A =||ai j (t)||, i, j =1,K,n . |
F = ç |
L |
÷ |
||
÷ |
ç |
|
÷ |
|
÷ |
è |
fn (t)ø |
|
|
ø |
|
|
|
|
Тогда систему (4.18) можно записать кратно в матричном виде
dY |
= AY + F . |
(4.19) |
dt |
|
|
85
PDF created with pdfFactory Pro trial version www.pdffactory.com
|
æ |
0 ö |
|
|
Определение 4.3. Если в системе (4.19) |
ç |
|
÷ |
, система (4.19) называется сис- |
F = çL÷ |
||||
|
ç |
0 |
÷ |
|
|
è |
ø |
|
темой линейных однородных дифференциальных уравнений.
Рассмотрим сначала систему линейных однородных дифференциальных уравнений.
|
|
|
|
|
dY = AY . |
(4.20) |
|
|
|
|
|
dt |
|
Пусть C(1) |
- совокупность всех вектор-функций, компоненты которых непрерывны |
|||||
[a, b] |
|
|
|
|
|
|
вместе со своими производными первого порядка на отрезке [a, b], C[a,b] - |
совокупность |
|||||
всех вектор-функций, компоненты которых непрерывны на отрезке [a, b]. |
|
|||||
Обозначим L[Y ] = dY - AY , L :C(1) |
] |
® C |
[a,b] |
- дифференциальный оператор. |
|
|
|
[a,b |
|
|
|
dt
Утверждение. L - линейный оператор.
В самом деле, пусть
|
|
|
æ |
(1) |
(1) |
|
ç |
(t) = |
ç |
||
Y1(t)ÎC[a,b], Y2 |
(t)ÎC[a,b] , Y1 |
ç |
|
|
|
|
ç |
|
|
|
è |
y1(1) (t)ö÷
L ÷ , Y2
÷
yn(1) (t)÷ø
æç y1(2) (t
(t) = ç L
ç
çè yn(2) (t
)ö÷ ÷ .
÷
)÷ø
Имеем
éêæç y1(1) (t)
L[Y1 +Y2 ] = Lêç êç
êç y(1) (t)
ëè n
æ d(y(1) + y(2) )ö |
|
|||
ç |
1 |
1 |
÷ |
æ |
|
dt |
|||
ç |
|
÷ |
ç |
|
= ç |
|
... |
÷ |
- Aç |
ç d(yn(1) + yn(2) )÷ |
çç |
|||
ç |
|
|
÷ |
è |
|
|
|||
ç |
|
dt |
÷ |
|
è |
|
ø |
|
+ y1(2) (t)ö÷ùú
... ÷ú = ÷ú
+ y(2) (t)÷ú n øû
(1) |
|
(2) |
ö |
(y1 |
+ y1 |
)÷ |
|
|
|
|
÷ |
(1) |
... |
(2) |
÷ . |
(yn |
+ yn |
÷ |
|
)ø |
Привлекая свойства производных, а также свойство дистрибутивности умножения матриц, приходим к равенству:
|
|
æ dy(1) |
ö |
æ dy(2) |
ö |
|
|
|
||
|
|
ç |
1 |
÷ |
ç |
1 |
÷ |
|
|
|
L[Y1 |
+Y2 |
ç |
dt |
÷ |
ç |
dt |
÷ |
- AY1 |
- AY2 |
|
] = ç ... |
÷ |
+ ç ... |
÷ |
= |
||||||
|
|
ç |
(1) |
÷ |
ç |
(2) |
÷ |
|
|
|
|
|
ç dyn |
÷ |
ç dyn |
÷ |
|
|
|
||
|
|
ç |
dt |
÷ |
ç |
dt |
÷ |
|
|
|
|
|
è |
ø |
è |
ø |
|
|
|
= L[Y1] + L[Y2 ] .
Пусть λ - произвольное действительное число, Y (t)ÎC[(a1), b] .
86
PDF created with pdfFactory Pro trial version www.pdffactory.com
Рассмотрим L[λY ] .
Используя свойства производных, имеем:
|
|
|
æ d(ly ) ö |
|
|
|
|
||||
æ ly (t) ö |
ç |
|
1 |
÷ |
æ ly |
ö |
|
||||
dt |
|
||||||||||
ç |
1 |
÷ |
ç |
|
÷ |
ç |
1 |
÷ |
= |
||
L[lY ] = Lç ... |
÷ |
= ç |
|
... |
÷ |
- Aç ... |
÷ |
||||
ç |
|
÷ |
ç d(lyn ) ÷ |
ç |
|
÷ |
|
||||
èlyn (t)ø |
ç |
|
|
|
÷ |
èlyn ø |
|
||||
|
dt |
|
|
||||||||
|
|
|
è |
|
ø |
|
|
|
|
æ dy (t) ö |
|
|
|
||||
ç |
|
1 |
÷ |
æ y |
ö |
|
|
dt |
|
||||||
ç |
|
÷ |
ç 1 |
÷ |
= l L[Y ]. |
||
= lç ... |
÷ |
- l Aç ... |
÷ |
||||
ç dyn (t) ÷ |
ç |
÷ |
|
||||
ç |
|
|
|
÷ |
è yn ø |
|
|
|
dt |
|
|
||||
è |
|
ø |
|
|
|
Откуда заключаем, что L - линейный оператор.
Используя введенные обозначения, систему (4.20) запишем в следующем виде:
L[Y ] = 0. |
(4.21) |
Теорема 4.3. Пусть Y1(t), Y2(t), …, Yk (t) - решения системы линейных однородных дифференциальных уравнений L[Y ] = 0. Тогда при любых λ1, λ2, …, λk вектор-функция Y(t) = λ1Y1(t) + λ2Y2(t) + … + λkYk(t) является решением системы L[Y ] = 0.
Доказательство. Рассмотрим L[Y ] = L[λ1Y1(t) + λ2Y2(t) +…+ λkYk(t)].
L[Y] |
= |
λ1L[Y1(t)] + λ2L[Y2(t)] +…+ λk L[Yk(t)]. |
|
L − линейный |
|
Далее, так как Y1(t), Y2(t), …, Yk(t) - решения системы (4.21), L[Y1] = 0, L[Y2] = 0,…, L[Yk] = 0, и, таким образом, L[Y ] = 0 Þ Y(t) = λ1Y1(t) + + λ2Y2(t) +…+ λkYk(t) - решение системы
(4.21).
Итак, любая линейная комбинация решений системы линейных однородных дифференциальных уравнений является решением этой системы.
Теорема 4.4. Если система линейных однородных дифференциальных уравнений (4.21) с действительными коэффициентами aij(t) имеет решением комплекснозначную векторфункцию Y(t) = U(t) + i V(t), то U(t) и V(t) - решения системы (4.21).
Доказательство слово в слово повторяет доказательство аналогичного утверждения для одного дифференциального уравнения порядка n и мы его опустим.
Определение 4.4. Система вектор-функций Y1(t), Y2(t), …, Yk(t) называется линейно зависимой на отрезке [a, b], если найдутся такие числа α1,…, αk, не все равные нулю одновременно, что "t Î[a, b]
a1Y1(t) + ... + akYk (t) = 0 .
87
PDF created with pdfFactory Pro trial version www.pdffactory.com
Система вектор-функций Y1(t), Y2(t), …, Yk(t) называется линейно независимой на отрезке [a, b], если из равенства
a1Y1(t) + ...+ akYk (t) = 0
следует, что α1 = … = αk = 0.
Приведем без доказательства следующее утверждение.
Теорема 4.5 (о структуре общего решения системы линейных однородных диффе-
ренциальных уравнений). Пусть дана система линейных однородных дифференциальных уравнений
dY |
= A(t)Y , |
(4.22) |
|
dt |
|||
|
|
где A = ||ai j (t)||, i, j =1,K,n , функции ai j (t) непрерывны на отрезке [a, b], и пусть Y1(t), Y2(t),
…, Yn(t) - линейно независимые на отрезке [a, b] решения системы (4.22).
Тогда общее решение системы линейных однородных дифференциальных уравнений (4.22) имеет вид
Yобщ = С1Y1(t) +K+ СnYn (t) ,
где С1, …, Сn - произвольные действительные числа.
Пример 4.3.
ì dx |
= y |
|
|
ï |
dt |
. |
(4.23) |
í |
|||
ï |
dy |
= -x |
|
î dt |
|
Частные решения: |
x1(t) = cos t ; y1(t) = -sin t ; x2 (t) = sin t ; y2 (t) = cost (убеждаемся непо- |
||||
средственно подстановкой). В векторной форме |
|||||
|
æ |
cost ö |
|
æ sin t ö |
|
Y (t) = ç |
÷ |
; Y (t) = ç |
÷ . |
||
1 |
ç |
÷ |
2 |
ç |
÷ |
|
è |
- sin t ø |
|
ècost ø |
Вектор-функции Y1(t) и Y2(t) - линейно независимы на любом отрезке [a, b]. В самом деле, пусть λ1 Y1(t) + λ2 Y2(t) = θ, или
ì |
l cost + l |
2 |
sin t = 0 |
(4.24) |
|
í |
1 |
cost = 0. |
|||
î |
- l sin t + l |
2 |
|
||
1 |
|
|
|
88
PDF created with pdfFactory Pro trial version www.pdffactory.com
Определитель системы линейных однородных алгебраических уравнений (4.24)
D = |
cos t |
sin t |
= cos2 t + sin2 t = 1 ¹ 0 система (4.24) имеет единственное решение λ1 = λ2 = |
|
- sin t |
cost |
|
0 и вектор-функции Y1(t) и Y2(t) линейно независимы на любом отрезке [a, b]. Тогда по теореме 4.5 общее решение системы (4.23)
Y |
= С |
æ |
cost ö |
+ С |
æ sin t ö |
|
ç |
÷ |
ç |
÷ . |
|||
общ |
1 |
ç |
÷ |
|
2 ç |
÷ |
|
|
è |
- sint ø |
|
ècost ø |
Определение 4.5. Любые n линейно независимых частных решений системы линей-
ных однородных дифференциальных уравнений называются фундаментальной системой решений.
В примере 4.3 вектор-функции Y1 |
æ |
cost ö |
и Y2 |
æ sint ö |
- фундаментальная система |
|
ç |
÷ |
ç |
÷ |
|||
(t) = ç |
÷ |
(t) = ç |
÷ |
|||
|
è |
- sint ø |
|
ècost ø |
|
решений системы дифференциальных уравнений (4.23).
4.4. Системы линейных однородных дифференциальных уравнений с постоянными коэффициентами
Пусть дана система |
|
|
|
|
dY |
= A×Y , |
(4.25) |
|
dt |
||
|
|
|
где A = ||ai j ||, i, j =1,K,n , aij - действительные числа.
Так как постоянные aij непрерывные на любом отрезке [a, b] функции, то применима сформулированная выше теорема 4.5.
Частные решения системы (4.25), составляющие фундаментальную систему решений, ищем в виде
æ a1 ö
= ç ÷ kt . (4.26)
Y (t) ç L÷e
çèan ÷ø
Подставляем (4.26) в систему уравнений (4.25)
ì ka1ekt = a11a1ekt + a12a2ekt +K+ a1nanekt |
|||||||||
ï |
|
ekt = a a ekt + a |
|
a |
ekt +K+ a |
|
a |
ekt |
|
ïka |
22 |
2n |
|||||||
í |
2 |
21 1 |
2 |
|
n |
|
|||
ï |
|
|
|
... |
|
|
|
|
|
ïka |
ekt = a |
a ekt + a |
n2 |
a |
ekt +K+ a |
nn |
a |
ekt |
|
î |
n |
|
n1 1 |
2 |
|
n |
|
89
PDF created with pdfFactory Pro trial version www.pdffactory.com
или
ì(a - k)a + a a |
2 |
+K+ a |
a |
n |
= 0 |
|
||||||||
ï |
11 |
1 |
|
|
12 |
|
|
1n |
|
|
|
|||
ïa21a1 + (a22 - k2 )a2 +K+ a2nan = 0 |
(4.27) |
|||||||||||||
í |
|
|
|
|
... |
|
|
|
|
|
|
|
||
ï |
|
|
|
|
|
|
|
|
|
|
|
|
||
ïa |
a + a |
n2 |
a |
2 |
+K+ (a |
nn |
- k)a |
n |
= 0. |
|
||||
î |
n1 1 |
|
|
|
|
|
|
|
|
|
Система (4.27) - система n линейных однородных алгебраических уравнений относи-
тельно α1, …, αn. Она имеет нетривиальные решения, если
|
a11 - k |
a12 |
... |
a1n |
|
|
|
|
|
|
|||||
D = |
a21 |
a22 - k2 |
... |
a2n |
= 0 . |
(4.28) |
|
... |
... |
... |
... |
||||
|
|
|
|||||
|
an1 |
an2 |
... |
ann - kn |
|
|
Уравнение (4.28) называется характеристическим уравнением для системы дифференциальных уравнений (4.25), а его корни - характеристическими корнями системы
(4.25).
Случай 1. Характеристические корни действительны и различны, тогда их ровно n. Пусть это числа k1, …, kn. Подставляя ki, i = 1, …, n в (4.27), находим соответствующие
a(i) ,...,a(i) , i =1,...,n . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
1 |
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
æa(i) ö |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ç 1 |
÷ |
kit |
|
линейно независимы на любом отрезке [a, b]. |
|
|||||||||||
|
Векторы Yi = |
ç |
÷ |
, i =1,...,n |
|
|||||||||||||
|
ç ... |
÷e |
|
|
||||||||||||||
|
|
ça(i) ÷ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
è n |
ø |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Покажем это. Допустим, существуют β1, …, βn, такие, что b1Y1(t) + ... + bnYn (t) = 0 . Тогда |
|||||||||||||||||
выполняются следующие равенства: |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
ìb a(1)ek1t + b a(2)ek2t |
+...+ b a(n)eknt = 0 |
|
|||||||||||
|
|
|
|
|
ï |
1 1 |
|
|
2 |
1 |
|
|
n |
1 |
|
|
|
|
|
|
|
|
|
ïb a(1)ek1t + b |
a(2)ek2t |
+...+ b |
a(n)eknt = 0 |
(4.29) |
|||||||||
|
|
|
|
|
í |
1 2 |
|
|
2 |
2 |
|
... |
n |
2 |
|
|
|
|
|
|
|
|
|
ï |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ï |
(1) |
e |
k1t |
|
(2) |
e |
k2t |
|
(n) |
e |
knt |
= 0. |
|
|
|
|
|
|
îb1an |
|
+ b2an |
|
+ ...+ bnan |
|
|
Так как система функций ek1t ,ek2t ,...,eknt линейно независима на любом отрезке [a, b], из равенств (4.29) следует n2 равенств
ìb a(1) |
= 0, b |
a(2) |
= 0,...,b |
a(n) |
= 0, |
|
||||
ï |
1 |
1 |
|
2 |
1 |
|
n |
1 |
|
|
ïb a(1) |
= 0, b |
2 |
a(2) |
= 0,...,b |
n |
a(n) |
= 0, |
(4.30) |
||
í |
1 |
2 |
|
2 |
|
2 |
|
|||
ï |
|
|
|
|
... |
|
|
|
|
|
ï |
|
(1) |
|
|
(2) |
|
|
(n) |
= 0. |
|
îb1an |
= 0, b2an |
= 0,...,bnan |
|
90
PDF created with pdfFactory Pro trial version www.pdffactory.com
|
Так как Y1 |
≠ |
θ, первый столбец равенств в (4.30) означает, что |
β1 |
= 0 (среди a(i) ,...,a(i) |
хотя бы одно число ≠ 0), аналогично заключаем, что β2 = … = βn = |
|
|
1 |
n |
|
0. Тогда система вектор-функций Y1, ..., Yn линейно независима на любом отрезке [a, b].
|
|
æa(1) |
ö |
|
|
æa(1) |
ö |
|
|
|
|
ç 1 |
÷ |
k1t |
|
ç 1 |
÷ |
knt |
|
Вектор-функции |
Y1 = |
ç |
÷ |
, ..., Yn = |
ç |
÷ |
составляют фундаментальную систему |
||
ç ... |
÷e |
|
ç ... |
÷e |
|
||||
|
|
ça(1) |
÷ |
|
|
ça(1) |
÷ |
|
|
|
|
è n |
ø |
|
|
è n |
ø |
|
|
решений системы дифференциальных уравнений (4.25). Тогда по приведенной теореме 4.5 общее решение системы уравнений (4.25) имеет вид
æa(1) |
ö |
|
ç 1 |
÷ |
k1t |
ç |
÷ |
|
Yобщ = C1ç ... |
÷e |
|
ça(1) |
÷ |
|
è n |
ø |
|
|
æa(1) |
ö |
|
ç 1 |
÷ |
+...+ Cn |
ç ... |
÷eknt . |
|
ç |
÷ |
|
ça(1) |
÷ |
|
è n |
ø |
Пример 4.4. Проинтегрировать систему дифференциальных уравнений
ì |
x& = 5x + y |
. |
(4.31) |
í |
= 10x + 2y |
||
îy& |
|
|
Решение. Записываем соответствующее характеристическое уравнение
|
5 - k |
1 |
|
= 0 . |
|
|
|||
|
10 |
2 - k |
|
|
Или
k2 – 7k = 0.
Его корни: k1 = 0, k2 = 7. Частные решения ищем в виде
|
|
|
|
|
|
æ a |
ö |
0t |
æ a |
|
ö |
æ b |
ö |
|
7t |
|
|
|
|
|
|
|
Y1(t) = |
ç 1 |
÷ |
ç 1 |
|
÷ |
ç 1 |
÷ |
|
. |
|
|
|||
|
|
|
|
ç |
÷e |
|
= ç |
÷ , |
Y2 (t) = ç |
÷e |
|
|
|
|||||
|
|
|
|
|
|
èa2 |
ø |
|
èa2 |
ø |
èb2 |
ø |
|
|
|
|
|
|
Подставляем Y1 в систему (4.31) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
ì 5a1 + a2 = 0 |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
í |
|
|
|
|
= 0, |
|
|
|
|
|
|
|
|
|
|
|
|
|
î10a1 + 2a2 |
|
|
|
|
|
|
|||||
|
|
|
æ |
a |
ö |
|
|
|
|
|
|
|
æ |
|
1 |
ö |
|
|
откуда α2 |
= –5α1 |
и Y1 |
ç |
1 |
÷ |
|
|
|
|
|
|
|
ç |
|
|
÷ |
|
|
(t) = ç |
- 5a |
÷ |
. Пусть α1 = 1, тогда Y1 = ç |
- 5÷ . |
|
|||||||||||||
|
|
|
è |
1 |
ø |
|
|
|
|
|
|
|
è |
|
|
ø |
|
|
Аналогично находим Y2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
ì- 2b + b |
2 |
= 0 |
Þ |
|
|
|
|
æ b |
ö |
7t |
|
||||
|
|
|
|
1 |
|
b2 = 2b1 и Y2 (t) = |
ç |
|
1 |
÷ |
. |
|||||||
|
|
|
í10b - 5b |
2 |
= 0 |
|
|
ç2b |
÷e |
|
||||||||
|
|
|
î |
1 |
|
|
|
|
|
|
|
è |
|
1 |
ø |
|
|
91
PDF created with pdfFactory Pro trial version www.pdffactory.com
Полагаем β1 = 1, тогда Y2 |
æ |
1 |
ö |
7t |
. |
ç |
|
÷ |
|||
= ç |
2 |
÷e |
|
||
|
è |
ø |
|
|
Итак, общее решение системы (4.31)
æ |
1 ö |
æ |
1 |
ö |
7t |
. |
|
ç |
|
÷ |
ç |
|
÷ |
||
Yобщ = C1ç |
- 5 |
÷ |
+ C2ç |
2 |
÷e |
|
|
è |
ø |
è |
ø |
|
|
Случай 2. Среди характеристических корней есть комплексные. Так как коэффициенты в системе уравнений (4.25), а следовательно, и в уравнении (4.28) действительные числа, комплексные корни входят сопряженными парами.
Пусть k1 = α + iβ и k2 = α – iβ являются корнями характеристического уравнения (4.28).
æ a1 ö
В этом случае вместо комплекснозначных вектор-функций Y = ç... ÷e(α+iβ)t и
1 ç ÷ çèan ÷ø
|
|
æ b |
ö |
|
|
|
~ |
|
Y |
= |
ç 1 |
÷ |
e(α−iβ)t |
|
|
= ReY |
|
ç ... |
÷ |
в фундаментальную систему решений включаем вектор-функции Y |
||||||
2 |
|
|
|
|
1 |
1 |
||
|
|
ç |
÷ |
|
|
|
|
|
|
|
èbn |
ø |
|
|
|
|
|
и |
~ |
|
|
|
~ |
~ |
линейно независимы. |
|
Y2 |
= ImY1 (теорема 4.4). Самостоятельно докажите, что Y1 |
и Y2 |
Тогда по теореме 4.5 общее решение системы уравнений (4.25) имеет вид
|
|
|
æa(1) |
ö |
|
|
æa(n−2) |
~ |
~ |
+ C3 |
ç 1 |
÷ |
k |
t |
ç 1 |
Yобщ = C1Y1 |
+ C2Y2 |
ç ... |
÷e 3 |
|
+ ... + Cn ç ... |
||
|
|
|
ç |
÷ |
|
|
ç |
|
|
|
èa(n1) |
ø |
|
|
èa(nn−2) |
ö
÷
÷eknt .
÷
ø
Пример 4.5. Проинтегрировать систему дифференциальных уравнений
ìdx |
= x + y, |
||
ï |
|
|
|
|
|||
í dt |
(4.32) |
||
ï |
|
dy |
= -2x + 3y. |
î dt |
Решение. Соответствующее характеристическое уравнение:
|
1- k |
1 |
|
= 0 |
|
|
|||
|
- 2 |
3 - k |
|
|
или
k2 - 4k - 5 = 0 Þ k1,2 = 2 ± i .
Частные решения ищем в виде
æ a |
ö |
(2+i)t |
|
ç 1 |
÷ |
. |
|
Y(t) = ç |
÷e |
|
|
èa2 |
ø |
|
|
92
PDF created with pdfFactory Pro trial version www.pdffactory.com