Линейные дифференциальные уравнения и системы линейных уравнений. Учебное пособие
.pdf
3. Линейныесистемы с постоянными коэффициентами
Решение. Будем отыскивать решение системы в виде
x 1ekt . y 2ekt , z 3ekt .
Составляем и решаем характеристическое уравнение
1 k |
1 |
2 |
|
1 k 2 1 k 8 4 |
1 k 4 |
1 k 0, |
|||
|
|||||||||
4 |
1 k |
0 |
0, |
||||||
k1 |
1, k2 |
k3 1. |
|
|
|||||
2 |
1 |
1 k |
|
|
|
||||
1) Строим собственный вектор, соответствующий простому корню k1 1. Решим систему (3.7), подставляя k 1 :
2 1 2 2 3 0,4 1 0,2 1 2 2 3 0.
Из последней системы находим 1 |
0 , 2 2 3 . Тогда собствен- |
|||||||
ный вектор |
|
1 и первое частное решение X , соответствующие |
||||||
A |
||||||||
|
|
|
|
|
|
|
|
1 |
корню k1 1, запишутся в виде |
|
|
|
|||||
|
|
|
|
0 |
|
|
0 |
|
|
|
|
|
1 |
2 |
, X |
|
2 et . |
|
|
A |
||||||
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
1 |
||
2) Построим два линейно независимых частных решения, соот-
ветствующих характеристическому числу |
k 1, кратность кото- |
|||
рого 2 . При k 1 получаем матрицу системы (3.7): |
||||
0 |
1 |
2 |
||
|
4 2 |
0 |
|
|
A kE |
. |
|||
|
2 |
1 |
0 |
|
|
|
|||
Её размерность n 3 , |
ранг |
r 2 . Число линейно независимых |
||||||
собственных векторов |
равно |
m n r 1 . |
Поскольку m , то |
|||||
решение следует искать в виде |
|
|
|
|
||||
x a |
a |
t e t , y b |
b t e t , z c |
c t e t . |
(3.11) |
|||
1 |
2 |
|
|
1 |
2 |
1 |
2 |
|
71
Нежельская Л. А. Линейные дифференциальныеуравнения
Для нахождения коэффициентов ai , bi , ci , i 1, 2 , подставим (3.11) в исходную систему и получим
a e t a |
a t e t |
b b t e t 2 |
c |
c t e t |
a |
a |
t e t , |
||||||||||
|
2 |
1 |
2 |
|
|
|
1 |
2 |
|
|
|
1 |
2 |
|
1 |
2 |
|
b e t b |
b t e t 4 |
a |
a t |
e t |
|
b |
b t e t , |
|
|
|
|
||||||
|
2 |
1 |
2 |
|
|
|
1 |
2 |
|
|
|
1 |
2 |
|
|
|
|
c e t c |
c t e t 2 |
a |
a t e t |
|
b b t e t c c t e t , |
||||||||||||
|
|
1 |
2 |
|
|
|
1 |
2 |
|
|
|
1 |
2 |
|
1 |
2 |
|
2 |
|
|
|
|
|
|
|
|
|||||||||
или, в результате сокращения на множитель e t 0 |
и приведения |
||||||||||||||||
подобных членов, запишем |
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
a |
2 |
b b t |
2c |
|
2c t 0, |
|
|
|
|
|||||
|
|
|
|
|
|
1 |
2 |
|
1 |
|
2 |
|
|
|
|
|
|
|
|
|
4a1 4a2t 2b1 2b2t b2 0, |
|
|
|
|
||||||||||
|
|
|
2a |
|
|
2a |
t b |
b t c |
0. |
|
|
|
|
||||
|
|
|
|
1 |
|
2 |
|
1 |
2 |
|
2 |
|
|
|
|
|
|
Приравнивая коэффициенты при t0 и t1 в левой и правой частях каждого уравнения, приходим к системе
|
|
t0 : a |
2 |
b |
|
2c |
0, t1 : |
b |
2c 0, |
|
|
|
|
|||
|
|
|
1 |
|
1 |
|
|
2 |
2 |
|
|
|
|
|||
|
|
4a1 2b1 b2 0, |
4a2 2b2 0, |
|
|
|
||||||||||
|
|
2a1 b1 c2 0. |
|
2a2 b2 0. |
|
|
|
|
||||||||
Поскольку кратность корня |
k 1 |
равна 2 , |
то произволь- |
|||||||||||||
ными будем считать два коэффициента a1 |
и a2 в решении (3.11). |
|||||||||||||||
Коэффициенты b1 , b2 , |
c1 , |
c2 |
выразим через a1 и a2 . Из последней |
|||||||||||||
системы находим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
b1 2a1 a2 , b2 2a2 , c1 a1 a2 , c2 a2 . |
|
||||||||||||||
Тогда решение (3.11) принимает вид |
|
|
|
|
|
|
||||||||||
x a a |
t e t , y |
2a a |
2 |
2a t e t , |
z a a |
2 |
a |
t e t . |
||||||||
1 |
2 |
|
|
|
|
1 |
|
2 |
|
1 |
|
2 |
|
|||
В качестве линейно независимых частных решений, соответствующих характеристическим числам k2 k3 1, можно взять
72
3. Линейныесистемы с постоянными коэффициентами
|
1 |
|
|
|
|
|
te t |
|
|
|
|
|
|
0 |
|
|
1 |
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
t |
. |
X2 |
2 e |
|
, X3 1 2t e |
|
|
|
1 |
t |
2 |
e |
|
||||||||||
|
|
|
|
|
|
1 t e |
t |
|
|
|
|
|
|
|
|
|
|
||||
|
1 |
|
|
|
|
|
|
|
|
|
1 |
|
|
1 |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3) Общее решение системы , согласно (3.5), имеет вид |
|
|
|||||||||||||||||||
x |
0 |
|
|
|
1 |
|
|
|
|
|
0 |
1 |
|
|
|
||||||
|
|
t |
|
|
|
t |
C3 |
|
|
|
|
|
|
t |
. |
|
|||||
y C1 2 e |
|
C2 |
2 e |
|
|
1 t |
2 e |
|
|
||||||||||||
z |
1 |
|
|
|
1 |
|
|
|
|
|
|
1 |
|
1 |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Пример 3.2. Найти общее решение системы |
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
x |
2x y z, |
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
x 2 y z, |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
y |
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
x y 2z. |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
z |
|
|
|
|
|
|
|
|
||||||||
Решение. Характеристическое уравнение имеет вид |
|
|
|
||||||||||||||||||
|
|
|
2 k |
|
1 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
1 |
|
2 k |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 , |
|
|
|
|
|
|
||||||
|
|
|
|
1 |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 k |
|
|
|
|
|
|
|
|
|||||||
2 k 3 1 1 1 2 k 1 2 k 1 2 k 0 , k3 6k 2 9k 4 0 .
Подбором находим корень k1 1; тогда характеристическое уравнение принимает вид k 1 k2 5k 4 0 или k 1 2 k 4 0 .
Таким образом, |
характеристические |
числа равны k1 k2 1, |
k3 4 . |
|
|
1) Для корня |
k 1 кратности 2 |
определим число линейно |
независимых собственных векторов. При k 1 система (3.7) принимает вид
73
Нежельская Л. А. Линейные дифференциальныеуравнения
|
|
2 |
|
3 |
0, |
1 |
|
|
|
||
1 2 3 0, |
|||||
|
2 |
3 0. |
|||
1 |
|||||
Размерность матрицы последней системы n 3 , ранг r 1. Число линейно независимых собственных векторов m n r 2 и совпадает с кратностью корня. Так как 1 2 3 , то в качестве
собственных векторов A 1 и A 2 можно взять
|
|
1 |
|
|
1 |
||
|
|
1 1 |
, |
|
2 |
0 . |
|
A |
A |
||||||
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
1 |
Два линейно независимых частных решения, соответствующих корню k 1 кратности 2 , запишутся в виде
|
|
|
|
|
|
|
1 |
|
|
|
|
|
1 |
|
||||
|
|
|
|
X |
|
|
1 et , X |
2 |
|
|
0 et . |
|||||||
1 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
1 |
|
|
|
|
2) Для простого корня k3 4 |
|
находим собственный вектор |
||||||||||||||
|
|
3 , решая систему |
|
|
|
|
|
|
|
|
|
|
|
|
||||
A |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
2 |
2 |
|
3 |
0, |
|
||||||||
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
||||
|
|
|
|
|
1 |
2 2 3 0, |
|
|
||||||||||
|
|
|
|
|
|
|
2 2 3 0. |
|
|
|||||||||
|
|
|
|
|
1 |
|
|
|||||||||||
Из последней системы находим |
|
1 3 |
и |
2 3 . Значит, соб- |
||||||||||||||
ственный вектор |
|
3 |
и частное решение |
X3 , соответствующие |
||||||||||||||
A |
||||||||||||||||||
простому корню k3 4 , принимают вид |
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
1 |
|
|
||
|
|
|
|
|
|
|
3 |
|
|
X3 |
|
|
4t |
|
||||
|
|
|
|
A |
|
. |
||||||||||||
|
|
|
|
|
1 , |
1 e |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
1 |
|
|
||
74
3.Линейныесистемы с постоянными коэффициентами
3)Общее решение исходной системы, согласно (3.5), запишется
ввиде
x |
|
1 |
|
|
1 |
|
1 |
|
|
||
|
|
|
|
|
t |
|
|
t |
|
4t |
. |
y |
C1 |
1 e |
|
C2 |
0 e |
|
C3 1 e |
|
|||
|
|
|
0 |
|
|
|
|
|
|
|
|
z |
|
|
|
|
1 |
|
1 |
|
|
||
Пример 3.3. Найти общее решение системы
x 3x y 2z,y 5y 3z,
z x 3y.
Решение. Составим и решим характеристическое уравнение
3 k |
1 |
2 |
|
|
k 3 |
k 5 k 3 2 5 k 9 3 k 0, |
||||||
|
0 |
5 k |
3 |
|
0 , |
|||||||
|
|
k |
3 |
|
8k |
2 |
22k 20 |
0. |
||||
|
1 |
3 |
k |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|||
Подбором находим корень |
k1 2 ; при этом характеристическое |
|
уравнение принимает |
вид |
k 2 k 2 6k 10 0 и уравнение |
k2 6k 10 0 имеет |
пару простых комплексных сопряжённых |
|
корней k2,3 3 i .
Таким образом, все корни характеристического уравнения различные. Построим фундаментальную систему решений.
1) Строим собственный вектор, соответствующий корню k1 2 . Решая систему
|
2 |
|
2 |
3 |
0, |
1 |
|
|
|
||
3 2 3 3 |
0, |
|
|||
|
|
|
2 3 0, |
||
1 3 2 |
|||||
находим 1 3 , 2 3 . Тогда собственный вектор A 1 и частное решение X1 , соответствующие простому корню k1 2 , имеют вид
75
Нежельская Л. А. Линейные дифференциальныеуравнения
|
|
|
|
|
|
|
|
1 |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 1 |
, X1 1 e2t . |
|
|
|
|
|
|
||||
|
|
|
|
|
A |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
1 |
|
|
|
|
|
|
|
|
2) Подставляя корень k2 3 i |
в систему (3.7), находим |
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
2 3 0, |
|
|
|
|
|
|
|
|
|
|
|
|
|
i 1 2 |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
2 i 2 3 3 0, |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
3 i 3 0. |
|
|
|
|
|
|
|||
|
|
|
|
1 3 2 |
|
|
|
|
|
|
||||||||
Приведём матрицу системы к ступенчатому виду |
|
|
|
|
||||||||||||||
1 3 |
3 i |
|
|
1 |
3 |
3 i |
1 |
3 |
3 i |
|||||||||
|
|
|
|
|
|
|
3 |
|
||||||||||
|
0 |
2 i |
3 |
|
|
0 2 |
i |
3 |
|
|
0 1 |
|
|
. |
||||
2 |
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
||
|
i |
1 |
2 |
|
|
|
|
0 1 3i |
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3i |
0 |
0 |
|
0 |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Выпишем связь между элементами собственного вектора:
2 23 i 3 , 1 3 2 3 i 3 .
Тогда собственный вектор можно задать в виде
|
|
|
|
|
|
2 i |
||
|
|
|
2 |
|
|
|
|
|
A |
3 |
|||||||
|
|
|
. |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 i |
||
Комплексное решение, соответствующее корню k2 3 i , есть
|
|
2 i |
|
|
2 i |
|
|
|
|
|
2 i |
|
|
|
|
|||||
|
|
|
|
3 |
|
|
3 i t |
|
|
|
3t |
e |
it |
|
3 |
|
3t |
cost i sin t |
||
|
|
|
|
|||||||||||||||||
X |
|
e |
|
|
|
3 e |
|
|
|
e |
|
|||||||||
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
2 i |
|
|
|
|
|
|
|
|
i |
|
|
2 |
i |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
2cost sin t |
|
|
|
cost 2sin t |
|
|
|
|||||||
|
|
|
|
|
|
3cost |
|
|
|
|
|
3sin t |
|
|
3t |
. |
||||
|
|
|
|
|
|
|
i |
|
|
e |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
2cost sin t |
|
|
cost 2sin t |
|
|
|||||||||
76
3. Линейныесистемы с постоянными коэффициентами
Взяв действительную и мнимую части от комплексного реше-
ния X , получим два линейно независимых действительных решения
|
|
|
2cost sin t |
|
|
|
|
|
cost 2sin t |
|
|
|
|
|
|
|
3cost |
|
3t |
, |
|
|
|
3sin t |
|
3t |
. |
|
|
|
|||||||||||
X2 Re X |
e |
|
X3 Im X |
e |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2cost sin t |
|
|
|
|
|
cost 2sin t |
|
|
||
3) Общее решение системы – линейная комбинация векторфункций Xi , i 1,3 :
x |
1 |
|
|
2cost sin t |
|
|
cos t 2sin t |
|
|
|
||
|
|
|
2t |
|
3cost |
|
3t |
|
3sin t |
|
3t |
. |
y |
C1 1 e |
|
C2 |
e |
|
C3 |
e |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
z |
1 |
|
|
2cost sin t |
|
|
cost 2sin t |
|
|
|||
3.2.2. Матричный метод
Линейную однородную систему всегда можно проинтегрировать матричным методом. Это значит, что всегда можно найти сразу всю фундаментальную матрицу системы.
Пусть дана система
dxi |
aji x j t , |
i 1, n , |
(3.12) |
||
dt |
n |
|
|
|
|
j 1 |
|
|
|
|
|
где все a ji – постоянные действительные числа. Введём в рассмотрение две матрицы
|
x11 t |
x12 t ... x1n t |
|
a |
|
a |
|
...a |
|
|
|||||
|
x |
|
|
|
|
t |
11 |
12 |
|
1n |
|
|
|||
X |
t |
x |
t ... x |
|
a |
|
a |
|
|
...a |
|
, |
|||
21 |
|
22 |
|
2n |
|
|
, A |
21 |
|
22 |
|
2n |
|||
|
................................. |
|
.................. |
|
|
||||||||||
|
|
t |
|
|
|
|
|
|
|
a |
|
|
...a |
|
|
|
|
xn2 |
|
|
|
|
a |
n1 |
n2 |
|
|
||||
|
xn1 |
t ... xnn t |
|
|
nn |
|
|||||||||
где X – фундаментальная матрица системы (3.12), A – транспонированная матрица коэффициентов системы (3.12).
77
Нежельская Л. А. Линейные дифференциальныеуравнения
Тогда систему (3.12) можно записать в виде одного равносильного матричного уравнения
|
|
|
dX |
XA . |
(3.13) |
|
|
|
|
dt |
|||
|
|
|
|
|
||
В (3.13) матрица |
dX |
называется производной от матрицы |
X и |
|||
dt |
||||||
|
|
|
|
|
||
находится дифференцированием всех элементов матрицы X . Всякая невырожденная матрица (с определителем, отличным
от нуля), обращающая уравнение (3.13) в тождество, называется
решением этого уравнения или интегральной матрицей.
Напомним, что показательной функцией eA |
матрицы A назы- |
||||||
вается сумма ряда |
|
|
|
|
|
|
|
eA E |
A |
|
A2 |
|
A3 |
... , |
(3.14) |
|
1! |
|
2! |
|
3! |
|
|
где E – единичная матрица той же размерности, что и матрица A . Ряд сходится для любой матрицы A .
Интегральной матрицей уравнения (3.13) является матрица ви-
да
X eAt , X 0 E . |
(3.15) |
Каждый элемент интегральной матрицы (3.15) представляет собой степенной ряд по степеням t , сходящийся при всех t . Все эти ряды можно просуммировать, так что интегральная матрица (3.15) будет состоять из элементарных функций.
Рассмотрим возможные случаи матрицы системы
1) Матрица A каноническая, при этом чисто диагональная
A k1, k2 ,...,kn ,
так что всем характеристическим числам k1,k2 ,...,kn матрицы A соответствуют простые элементарные делители
k k1, k k2 ,...,k kn ,
78
3. Линейныесистемы с постоянными коэффициентами
тогда как среди характеристических чисел k1,k2 ,...,kn могут быть и равные между собой. Пользуясь формулой (3.14), находим
eAt e k1,k2 ,...,kn t e k1t,k2t,...,knt ek1t ,ek2t ,...,eknt . |
||||
|
|
|
||
В этом случае интегральная матрица (3.15) принимает вид |
||||
X ek1t ,ek2t |
,...,eknt . |
(3.16) |
||
|
|
|
|
|
Матрица X чисто диагональная, |
с функциями |
ekit , i |
|
, на |
1, n |
||||
главной диагонали.
Замечание. Матрица A чисто диагональная, если:
а) все собственные значения ki простые, т.е. имеют кратность
i 1;
б) среди собственных значений ki имеются кратные, но каждо-
му кратному корню, кратность которого |
i 1, |
соответствует |
i |
||||||||||||||||||||||||
одномерных жордановых клеток. |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
2) Матрица A каноническая, но квазидиагональная |
|
|
|||||||||||||||||||||||||
A |
I |
1 |
k |
,I |
2 |
k |
2 |
,...,I |
s |
k |
s |
, |
|
2 |
... |
s |
n . |
|
|||||||||
|
|
|
1 |
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
||||||||||
Здесь I i |
ki |
|
– жордановая клетка, есть матрица размерности |
i , |
|||||||||||||||||||||||
имеющая вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
k |
|
|
0 |
|
0 ... |
0 |
0 |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
i |
ki |
|
0 ... |
0 |
0 |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
I i |
ki |
0 |
|
1 |
ki ... |
0 |
0 |
. |
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
.......................... |
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
0 |
0 ... |
1 k |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
Заметим, |
что |
|
I1 ki ki . |
|
Матрица |
|
I i |
ki |
|
соответствует элемен- |
|||||||||||||||||
тарному делителю |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
k ki i , i 1, s .
79
Нежельская Л. А. Линейные дифференциальныеуравнения
Среди характеристических чисел k1,k2 ,...,kn могут быть и равные
между собой.
В этом случае, пользуясь формулой (3.14), находим интегральную матрицу
X eAt |
I |
1 |
k |
,I |
2 |
k |
|
,...,I |
s |
k |
|
t |
|
|
I |
1 |
k |
t,I |
2 |
k |
|
t,...,I |
s |
k |
|
t |
|
||||||||||||
e |
1 |
|
|
|
2 |
|
|
|
|
s |
|
e |
1 |
|
|
|
|
2 |
|
|
|
|
s |
|
|||||||||||||||
|
|
|
|
|
|
|
I 1 k1 t |
,e |
I |
2 k2 t |
,...,e |
I s ks t |
, |
|
|
|
|
|
|
|
|
(3.17) |
|||||||||||||||||
|
|
|
e |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ekit |
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
... |
0 |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
tekit |
|
|
|
|
|
|
|
|
|
ekit |
|
|
|
|
|
|
... |
0 |
|
|
|
|
|
|
|
||||||
|
eI i ki t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
.................................................. |
. |
|
|
|
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t i 1 |
|
|
|
k t |
|
|
t i 2 |
|
|
|
|
k t |
|
|
|
|
|
k t |
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 !e |
|
|
|
|
2 |
!e |
|
|
... |
e |
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
3)Матрица A произвольного вида.
Вэтом случае её всегда можно привести к каноническому виду при помощи преобразования подобия, т.е. можно построить такую
невырожденную матрицу S det S 0 , чтобы матрица |
|
B SAS 1 |
(3.18) |
была канонической.
В матричном уравнении (3.13) выполним подстановку
X ZS , (3.19)
где Z – новая неизвестная интегральная матрица (решение нового дифференциального матричного уравнения). Получаем
ddtΖS ZSA , ddtZ ZSAS 1
или
dZ |
ZB . |
(3.20) |
|
dt |
|||
|
|
80
