
ГЛАВА 3. ПЕРВЫЙ МЕТОД ЛЯПУНОВА_______________________________________________________135
Подставляя найденное для χ[W(t)] выражение в (**), придём к неравенству Ляпунова ■
З а м е ч а н и е 2. В случае комплекснозначной матрицы A(t) имеем SpA = ReSpA +i Im SpA . А так как согласно формуле Эйлера мнимая часть следа матрицы ограничена, то неравенство Ляпунова принимает вид
|
t |
|
|
|
1 |
|
t |
|
SX ≥ χ exp |
Re SpA(τ)dτ |
= |
lim |
|
∫ |
Re SpA(τ)dτ. |
||
|
∫ |
|
|
t→+∞ t |
|
|
||
|
t0 |
|
|
|
|
|
t0 |
|
20. Равенство Ляпунова.
Следствие 1. Если для ф.с.р. X(t) линейной системы (1) с непрерывной, ограниченной вещественной матрицей A(t) имеет место равенство Ляпуно-
ва
|
|
t |
|
|
|
1 |
|
t |
|
SX |
= χ exp |
SpA(τ)dτ |
= |
lim |
|
∫ |
SpA(τ)dτ, |
||
|
|
∫ |
|
|
t→+∞ t |
|
|
||
|
|
t0 |
|
|
|
|
|
t0 |
|
то такая ф.с.р. является нормальной.
□ Для доказательства, не ограничивая общности, предположим, что в
фундаментальной совокупности |
X(t) ={x(1) ,..., x(n) } |
решения записаны в по- |
|||
рядке возрастания их характеристических показателей. |
|
||||
Если |
ф.с.р. |
X(t) , напротив, не является |
нормальной, то |
найдётся |
|
k {2,..., n} |
и линейная комбинация решений |
|
|
||
|
|
|
k |
|
|
|
|
z = ∑cil x(il ) (t) , |
|
(*) |
|
|
|
|
l=1 |
|
|
{i1,...,ik } {1,2,..., n}, i1 < ... < ik , |
с отличными от нуля коэффициентами (в ча- |
||||
стности, cik |
≠ 0), |
для которой |
|
|
|
|
|
χ[z] < max χ[x(il ) (t)] = χ[x(ik ) (t)] . |
(**) |
||
|
|
il =1,...,k |
|
|
Рассмотрим совокупность решений
Z ={x(1) ,...,x(ik −1) ,z,x(ik +1) ,...,x(n)}.
Она является ф.с.р., так как если предположить обратное, т.е.
∑a j x(j) (t) +aik z(t) ≡ 0 |
n |
при некоторых aj: ∑| a j |≠ 0, |
|
j≠ik |
j=1 |
ГЛАВА 3. ПЕРВЫЙ МЕТОД ЛЯПУНОВА_______________________________________________________136
то в силу линейной независимости системы |
{x(j) (t)} |
|
верно a |
≠ 0, а значит с |
||||||||
учетом (*) получаем |
|
|
|
|
|
|
j≠ik |
|
ik |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∑(a j +aik c j ) x(j) (t) +aik cik |
x(ik ) (t) + ∑a j x(j) (t) ≡ 0 . |
|
|
|||||||||
j<ik |
|
|
|
|
|
|
j>ik |
|
|
|
|
|
Отсюда, благодаря линейной независимости системы векторов |
{x(j) (t)}n |
|
, сле- |
|||||||||
|
|
|
|
|
|
|
|
|
|
j=1 |
|
|
дует равенство aik cik = 0, которое невозможно в силу cik ≠ 0, |
aik ≠ 0 . |
|
|
|||||||||
Полученное противоречие говорит о том, что |
Z является ф.с.р., |
причём |
||||||||||
благодаря (**) выполнено |
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
t |
|
|
|
|
|
SZ < SX = χ exp |
∫ |
SpA(τ)dτ |
= lim 1 |
∫ |
SpA(τ)dτ. |
|
|
|||||
|
|
|
|
|
t→+∞ t |
|
|
|
|
|||
|
|
0 |
|
|
|
|
|
0 |
|
|
|
|
Это неравенство не совместимо с неравенством Ляпунова для ф.с.р. Z ■ Пример 1. Приведем пример ф.с.р. линейной однородной системы урав-
нений
|
& |
x = y [sin(lnt) +cos(lnt)], |
|
|
(1 ≤ t < +∞), |
|
|
|
& |
y = x [cos(lnt) +sin(lnt)], |
для которой равенство Ляпунова не выполняется.
После сложения и вычитания уравнений данной системы, находим
x + y = 2C et sin(lnt) , |
|
|
x − y = 2C |
2 |
e−t sin(lnt) , |
||||||
1 |
|
|
|
|
|
|
|
|
|
|
|
откуда получаем её общее решение |
|
|
|
|
|
|
|
|
|
||
x = C et sin(lnt) +C |
|
e−t sin(lnt) , |
|
|
|
||||||
|
1 |
|
|
2 |
|
|
|
|
|
|
|
|
|
t sin(lnt) |
|
|
|
−t sin(lnt) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
y = C1 e |
−C |
2 e |
. |
|
|
|
|||||
|
|
|
|
|
|
||||||
Нетрудно вычислить, что |
χ[x] = χ[y] =1 (при C2 |
+C2 |
≠ 0). Поэтому ф.с.р. |
||||||||
|
|
|
|
|
|
1 |
|
|
|
2 |
|
et sin(lnt)
X = et sin(lnt)
e−t sin(lnt) |
|
|
|
|
является нормальной и S |
X |
= 2 . Однако для матрицы |
|
|
|
|
−e−t sin(lnt) |
|
|
|
0 |
sin(lnt) +cos(lnt) |
|
A(t) = |
|
|
|
0 |
|
cos(lnt) +sin(lnt) |
|
ГЛАВА 3. ПЕРВЫЙ МЕТОД ЛЯПУНОВА_______________________________________________________137
данной системы верно SpA(t) ≡ 0, а значит равенство Ляпунова нарушается:
2 = SX > |
|
1 |
t |
SpA(τ)dτ = 0 . |
lim |
||||
|
t→+∞ t |
∫ |
|
|
|
|
|
1 |
|
30. Правильные системы.
Существует известная связь между верхним и нижним пределами одной и той же функции. Благодаря этой связи и неравенству Ляпунова, получаем
|
|
1 |
t |
|
|
t |
|
S ≥ |
|
∫SpA(τ)dτ ≥ lim 1 |
∫SpA(τ)dτ, |
|
|||
lim |
|
||||||
t→+∞ t |
0 |
|
t→+∞ t |
0 |
|
||
а значит, |
|
|
|
||||
|
|
t |
|
|
|||
|
|
|
|
lim 1 |
|
|
|
|
|
|
S ≥ |
∫SpA(τ)dτ. |
|
||
|
|
|
t |
→+∞ t |
0 |
|
|
|
|
|
|
|
|
|
|
Линейная система, для которой в последнем неравенстве имеет место ра- |
|||||||
венство, носит специальное название. |
|
|
|
||||
О п р е д е л е н и е 2. Линейную систему (1) с непрерывной, ограничен- |
|||||||
ной, вещественной матрицей A(t) и спектром α1,...,αm (m ≤ n) |
называют пра- |
||||||
вильной по Ляпунову (или просто правильной), если |
|
||||||
|
|
|
|
lim 1 |
t |
|
|
|
|
|
S = |
∫SpA(τ)dτ, |
(4) |
t→+∞ t 0
где
m
S = ∑ns αs
s=1
есть сумма характеристических показателей (с учётом кратностей ns ) некоторой н.ф.с.р. системы (1).
З а м е ч а н и е 3. В аналогичном определении для системы (1) с комплекснозначной матрицей A(t) равенство (4) следует заменить на
S = lim 1 |
t |
∫ReSpA(τ)dτ. |
|
t→+∞ t |
0 |
|
|
Пример 2 (неправильная система). Рассмотрим линейную систему второ- |
|
го порядка |
|
|
& |
x = x cos(lnt) + y sin(lnt), |
|
|
(t ≥ 1). |
|
|
|
& |
y = x sin(lnt) + y cos(lnt),

ГЛАВА 3. ПЕРВЫЙ МЕТОД ЛЯПУНОВА_______________________________________________________138
При помощи подстановки нетрудно убедиться, что вектор-функции
|
t sin(lnt) |
|
|
|
t sin(lnt) |
|
x(t) = e |
|
, |
y(t) = e |
|
|
|
|
|
|
|
|
|
|
et cos(lnt) |
|
|
−et cos(lnt) |
являются решениями данной системы. Более того, они образуют линейно независимую систему и, тем самым, представляют собой ф.с.р. Кроме того, в силу
χ[C1x +C2y] = max{χ[x],χ[y]} =1,
это будет н.ф.с.р. Для такой н.ф.с.р. верно равенство S = 2. С другой стороны, имеем
lim 1 |
t |
lim 1 |
(τ (sin(ln τ) +cos(ln τ)))ττ==1t |
|
∫2cos(lnτ)dτ = |
= |
|||
t→+∞ t |
1 |
t→+∞ t |
|
|
|
|
|
|
= lim (sin(ln t) +cos(ln t)) = − 2 ≠ 2 = S.
t→+∞
Полученное в итоге неравенство свидетельствует о том, что данная система неправильная.
40. Приводимые системы.
О п р е д е л е н и е 3. Пусть x = (x |
,x |
2 |
,...,x |
n |
)T , |
y = (y ,y |
2 |
,...,y |
n |
)T . Ли- |
1 |
|
|
|
1 |
|
|
||||
нейное преобразование |
|
|
|
|
|
|
|
|
|
|
y = L(t) x |
|
|
|
|
|
|
|
|
|
(5) |
называют преобразованием Ляпунова, если непрерывно дифференцируемая при
t ≥ 0 |
матрица n-го порядка L(t) удовлетворяет следующим условиям: |
1) |
sup L(t) < +∞; |
t≥0
2)sup L&(t) < +∞;
t≥0
3)K > 0, t ≥ 0 : det(L(t)) ≥ K .
Важное свойство преобразования Ляпунова раскрывает следующая
Лемма 1. Преобразование Ляпунова сохраняет значение характеристического показателя, т.е.
y = L(t) x χ[y] = χ[x] .
□ Из (5) на основе условия 3) получаем x = L−1(t) y . Поэтому
ГЛАВА 3. ПЕРВЫЙ МЕТОД ЛЯПУНОВА_______________________________________________________139
|| y ||≤|| L(t) || || x || |
|
χ[y] = |
|
|
|
|
|
−1 |
|
|
|
|
|
|
|
|| x ||≤|| L |
(t) || || y || |
χ[x] = |
Следовательно, χ[y] = χ[x] ■
χ[|| y ||] ≤ χ[|| L(t) ||]+ χ[|| x ||] = χ[|| x ||] = χ[x]
χ[|| x ||] ≤ χ[|| L−1 (t) ||]+ χ[|| y ||] = χ[|| y ||] = χ[y].
О п р е д е л е н и е 4. Линейная система (1) называется приводимой, если при помощи некоторого преобразования Ляпунова (5) (с вообще говоря комплексной матрицей L ) её можно привести к системе
dy |
= B y |
(6) |
|
dt |
|||
|
|
с постоянной матрицей B.
Как показывает следующая теорема, класс приводимых систем содержится в классе правильных систем.
Теорема 2. Всякая приводимая линейная система (1) является правиль-
ной.
□ Пусть система (1) является приводимой и X(t) – некоторая её н.ф.м.р. По определению приводимой системы найдётся матрица L(t) :
X(t) = L(t) Y(t), |
(*) |
где Y(t) – ф.м.р. линейной системы (6) с постоянной матрицей |
B. Из (*) сле- |
дует |
|
det(X(t)) = det(L(t)) det(Y(t)). |
|
Отсюда с использованием формулы Остроградского-Лиувилля (см. п.10 в §6) получаем
t |
|
|
t SpB |
|
|
|
|
|
= det(L(t)) det(Y(0)) e |
, |
|
det(X(0)) exp |
∫SpA(τ)dτ |
|
|||
|
0 |
|
|
|
|
откуда |
|
|
|
|
|
t |
|
|
t SpB |
|
|
|
|
|
= c(0) | det(L(t)) | e |
, |
|
exp |
∫SpA(τ)dτ |
|
|||
|
0 |
|
|
|
|
где
c(0) =| det(Y(0) X−1 (0))|.
Поэтому
ГЛАВА 3. ПЕРВЫЙ МЕТОД ЛЯПУНОВА_______________________________________________________140
1 |
t |
|
|
|
= 1 ln[c(0) | det(L(t)) |]+1 SpB = |
|
||||
∫SpA(τ)dτ |
|
|||||||||
t |
0 |
|
|
|
t |
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
const |
|
1 |
|
огр. |
|
1 |
|
|
|
678 |
|
6447448 |
SpB → SpB. |
|
|||||
= |
t |
lnc(0) |
+ |
t |
ln | det(L(t)) | + |
t |
|
|||
|
|
|
|
|
|
t→+∞ |
|
|||
Таким образом, |
|
|
|
|
|
|
|
|
|
|
|
1 |
t |
|
|
|
lim 1 |
t |
|
|
|
lim |
∫SpA(τ)dτ = |
∫SpA(τ)dτ = SpB. |
(**) |
|||||||
t→+∞ t |
0 |
|
|
|
t→+∞ t |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Обозначим через SX и SY |
суммы характеристических показателей ф.с.р. |
|||||||||
X и Y соответственно. Согласно лемме 1 справедливо равенство |
|
|||||||||
|
|
|
|
|
SX = SY, |
|
|
|
(***) |
причём, так как по условию X – н.ф.с.р., то и Y – н.ф.с.р. Но для Y характеристическими показателями являются вещественные части корней характеристического уравнения det(B – λE) = 0, где каждый корень считается столько раз, какова его кратность. Поэтому
SY = ∑Reλk = Re ∑λk = ReSpB = SpB .
k |
k |
Отсюда вместе с учетом (**)–(***) следует
(***) |
(**) |
lim 1 |
t |
SX = |
SY = SpB = |
∫SpA(τ)dτ, |
|
|
|
t→+∞ t |
0 |
|
|
|
т.е. система (1) действительно правильная ■ Из доказанной теоремы очевидным образом вытекает
Следствие 2. Любая линейная однородная дифференциальная система с постоянной матрицей является правильной.
Можно доказать (см. [3]), что в случае, когда матрица A(t) линейной системы (1) является периодической, эта система также будет правильной.
Упражнение
1. Доказать, что линейная система (1) второго порядка с матрицей

ГЛАВА 3. ПЕРВЫЙ МЕТОД ЛЯПУНОВА_______________________________________________________141
0 |
1 |
|
|
|
2 |
|
(при t ≥ 1) |
A(t) = |
|
||
|
|
0 |
|
|
|
||
t2 |
|
|
неприводима.
§ 7. Теорема Перрона
10. Взаимно сопряжённые системы. Рассмотрим линейную однородную систему дифференциальных уравнений
|
|
dx |
= A(t) x |
(t ≥ 0) |
(1) |
|
|
|
dt |
||||
|
|
|
|
|
||
матрица A(t) которой в общем случае может быть комплекснозначной. |
|
|||||
О п р е д е л е н и е 1. Линейную однородную систему |
|
|||||
|
dy |
|
= −A* (t) y |
(t ≥ 0), |
(2) |
|
|
dt |
|||||
|
|
|
|
где A*(t) = AT (t) есть сопряжённая матрица для A(t) , называют сопряжённой системой для системы (1).
Если A(t) – вещественная матрица, то A*(t) = AT (t) и сопряжённая система для системы (1) принимает вид
ddty = −AT (t) y .
Нетрудно проверить, что система (1), в свою очередь, является сопряжённой для системы (2). На этом основании о системах (1) и (2) говорят как о вза-
имно сопряжённых системах.
Лемма 1. 1) Для любых решений x и y взаимно сопряжённых систем (1) и (2) выполняется
y* x = x, y = C − const , |
(3) |
где y* – |
вектор, сопряжённый для y. |
|
2) |
Для любых ф.м.р. X=X(t) и Y=Y (t) указанных систем имеет место |
|
равенство |
|
|
|
Y* X = C − постоянная матрица, |
(4) |
ГЛАВА 3. ПЕРВЫЙ МЕТОД ЛЯПУНОВА_______________________________________________________142
где Y* – сопряжённая для Y матрица. Обратно, из выполнения равенства (4)
с неособой матрицей |
|
C и |
ф.м.р. |
|
|
|
X системы (1) следует, что Y – есть |
|||||||||||||||||||||||||||||
ф.м.р. сопряженной системы (2). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
)T |
|
|
)T |
||||||||||||
□ 1) Зафиксируем произвольные решения |
x =( x |
,...,x |
n |
и y =( y ,...,y |
n |
|||||||||||||||||||||||||||||||
взаимно сопряжённых систем (1) и (2). |
|
|
|
|
|
|
|
|
|
1 |
|
|
|
1 |
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
Поскольку |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
& |
* |
= (−A |
* |
|
y) |
* |
|
|
&* |
= −y |
* |
A , |
|
|
|
|
|
|
|||||||||||||
|
|
|
|
(y) |
|
|
|
|
|
|
y |
|
|
|
|
|
|
|
||||||||||||||||||
строка y* является решением системы |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
dy |
* |
|
= −y |
* |
A(t) . |
|
|
|
|
|
|
|
|
(5) |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Из равенств (1) и (5) следует |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
y |
* |
|
dx |
= |
y |
* |
A(t) x , |
|
|
|
|
|
dy |
* |
|
|
* |
A(t) x . |
|
|
|
|||||||||||||||
|
dt |
|
|
|
|
|
|
|
dt |
x = −y |
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Складывая почленно последние два равенства, получим |
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
y |
* |
|
|
dx |
+ |
dy* |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
dt |
|
dt |
x = 0 , |
|
|
|
|
|
|
|
|
|
|
||||||||||||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
d |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
(y* x) = 0. |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Отсюда вытекает равенство (3). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X и |
Y устанавлива- |
||||||||||||
2) Первая половина этого утверждения для матриц |
||||||||||||||||||||||||||||||||||||
ется так же, как и для векторов x |
|
|
и y. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
Для доказательства обратного предположим, что имеет место равенство |
||||||||||||||||||||||||||||||||||||
(4). Из него следует |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Y = (C X−1 )* = (X* )−1 C* , |
|
|
|
|
|
|
|
|
(6) |
|||||||||||||||||||||
где матрица X* такова, что |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
& * |
= |
X |
* |
|
* |
(t) . |
|
|
|
|
|
|
|
|
|
|
|
(7) |
||||||||||
|
|
|
|
|
|
|
|
X |
|
|
|
A |
|
|
|
|
|
|
|
|
|
|
|
|
Дифференцируя (6) по правилу дифференцирования обратной матрицы, с учётом (7) получаем
ГЛАВА 3. ПЕРВЫЙ МЕТОД ЛЯПУНОВА_______________________________________________________143
& |
* |
−1 |
& * |
* |
−1 |
* |
(7) |
* |
−1 |
* |
* |
* |
−1 |
* |
(6) |
* |
(t) Y , |
|
Y = −(X ) |
|
|
X |
(X ) |
|
C |
=−(X ) |
|
X |
A |
(t) (X ) |
|
C |
=− A |
||||
|
144424443 |
|
|
14243 |
|
|
|
|
|
|
|
|||||||
|
|
|
d |
(X* )−1 |
|
|
|
|
E |
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
причём
det(Y) = det((X* )−1 C* )= det(X* )−1 det(C* ) = det(X)−1 det(C) ≠ 0 .
Полученное означает, что Y является ф.м.р. системы (2) ■
20. Критерий правильности системы.
Лемма 2. Линейная однородная система (1) с непрерывной, ограниченной и в общем случае комплекснозначной матрицей A(t) является правильной тогда и только тогда, когда
lim 1 |
t |
|
∫ReSpA(τ)dτ = S , |
(8) |
|
t→+∞ t |
0 |
|
|
|
где S – сумма всех характеристических показателей данной системы.
□ Часть «достаточность» вытекает из определения правильной системы. Для доказательства необходимости введём обозначения
|
|
1 |
t |
σ = lim 1 |
t |
σ = |
|
∫ReSpA(τ)dτ, |
∫ReSpA(τ)dτ. |
||
lim |
|||||
|
t→+∞ t |
0 |
t→+∞ t |
0 |
|
|
|
|
|
Согласно неравенству Ляпунова S ≥ σ. Но из-за того, что система (1) правильная, верно равенство S = σ. Следовательно, σ ≥ σ, что вместе с неравенством σ ≤ σ, связывающим нижний и верхний пределы, влечёт равенство σ = σ. Отсюда с учётом S = σ приходим к (8) ■
30. Теорема Перрона.
Теорема 1 (О. Перрон). Линейная система (1) с непрерывной, ограниченной и в общем случае комплекснозначной матрицей A(t) является правильной тогда и только тогда, когда полный спектр
α1 ≤ α2... ≤ αn
этой системы и полный спектр
β1 ≥ β2 ≥... ≥ βn
сопряжённой системы (2) симметричны относительно нуля, т.е.
αs +βs = 0 , |
s =1,2,...,n . |
ГЛАВА 3. ПЕРВЫЙ МЕТОД ЛЯПУНОВА_______________________________________________________144
□ Необходимость. Пусть система (1) правильная и |
X(t) = (x jk )n×n |
– её |
||||||||||||||||||||||||||||
н.ф.м.р., состоящая из столбцов |
|
x(k) = (x1k (t),...,xnk (t))T , |
k =1,2,...,n , таких, |
|||||||||||||||||||||||||||
что χ[x(k) ] = αk , причём характеристические показатели |
αk |
упорядочены по |
||||||||||||||||||||||||||||
возрастанию (как в условиях теоремы). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
Согласно лемме 1 матрица |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
Y(t) = [X−1 (t)]*= (y jk )n×n |
|
|
|
(9) |
||||||||||||||||||||
является ф.м.р. сопряжённой системы (2), так как Y* X = E. |
|
|
|
|||||||||||||||||||||||||||
Введём обозначение для столбцов матрицы Y: |
|
|
|
|
||||||||||||||||||||||||||
y(k) = (y |
(t),..., y |
nk |
(t))T , |
|
|
|
|
χ[y(k) ] = β |
k |
, k =1,2,...,n , |
|
|||||||||||||||||||
|
|
|
1k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
причём характеристические показатели |
|
|
βk |
|
упорядочены в порядке убывания |
|||||||||||||||||||||||||
(как в условиях теоремы). |
|
|
|
|
|
(s)* |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
Равенство |
* |
|
|
|
|
|
|
|
y |
x |
(s) |
|
=1, s =1,2,...,n . |
Переходя в этом |
||||||||||||||||
Y X = E влечёт |
|
|
|
|
|
|
||||||||||||||||||||||||
равенстве к характеристическим показателям, получим |
|
|
|
|
||||||||||||||||||||||||||
0 = χ[y(s)* x(s) ] ≤ χ[y(s)* ] + χ[x(s) ] = αs +βs |
|
|
αs |
+βs ≥ 0 . |
(*) |
|||||||||||||||||||||||||
С другой стороны, обозначая через X jk (t) |
|
алгебраическое дополнение |
||||||||||||||||||||||||||||
элемента x jk (t) |
матрицы X , согласно правилу обращения матрицы получаем |
|||||||||||||||||||||||||||||
|
|
|
|
|
|
Xsj(t) |
|
|
* |
|
|
|
|
|
js (t) |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
y js (t) = |
|
= |
|
|
X |
|
|
|
, |
|
j =1,2,..., n , |
|
|||||||||||||||
|
|
|
det(X(t)) |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
det(X(t)) |
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
где в соответствии с формулой Остроградского-Лиувилля |
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
det(X(t)) = det(X(0)) exp |
∫ |
SpA(τ)dτ ≠ 0. |
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
14243 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
≠0 |
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
||
Поэтому |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
1 |
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
χ[y js (t)] ≤ χ |
|
|
|
|
+ χ exp |
− ∫SpA(τ)dτ + χ[X js (t)], |
|
j =1,2,..., n . |
(**) |
|||||||||||||||||||||
det(X(0)) |
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Первое слагаемое в правой части (**) равно нулю. Рассмотрим подробнее второе и третье слагаемые. Прежде всего заметим, что система (1) правильная, а значит, в силу леммы 2