Хаотическая динамика и классическая теория устойчивости движения
.pdf7.1. МЕТОД ТРИАНГУЛЯЦИИ ПЕРРОНА – ВИНОГРАДА |
71 |
7.1. Метод триангуляции Перрона – Винограда
Проведем процедуру ортогонализации Шмидта для решений zj (t), образующих матрицу (7.3).
v1(t) = z1(t) |
|
|
|
|
|
|
|
|
|
|
|
|
|
v2(t) = z2(t) − v1(t) z2(t) |
|
|
v1(t) |
|
|
|
|
||||||
|
|
v (t) 2 |
|
|
|
|
|
|
|
||||
|
| |
1 |
|
| |
|
|
|
|
|
|
|
||
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . |
|
|
|
|
|
(7.5) |
|||||||
|
|
|
|
v1(t) |
|
|
|||||||
vn(t) = zn(t) − v1(t) zn(t) |
|
|
|
|
|
||||||||
|
|
− . . .− |
|
|
|
||||||||
|
v (t) 2 |
|
|
|
|||||||||
|
| 1 |
|
| |
|
|
|
vn−1(t) |
|
|
||||
|
− vn−1(t) zn(t) |
|
|
. |
|||||||||
|
|
v |
(t) |
2 |
|||||||||
|
|
|
|
|
|
| |
n−1 |
| |
|
||||
Из соотношений (7.5) сразу следуют равенства |
|
|
|
|
|||||||||
v |
(t) v (t) = 0, |
|
j = i, |
|
|
(7.6) |
|||||||
i |
j |
|
|
|
|
|
|
|
|||||
|
|vj (t)|2 = vj (t) zj (t). |
|
|
(7.7) |
|||||||||
Из равенства (7.7) сразу вытекает следующее утверждение. |
|||||||||||||
Лемма 7.1. Имеет место оценка |
|
|
|
|
|
|
|
|
|
||||
|vj (t)| |zj (t)|, |
t 0. |
|
|
(7.8) |
|||||||||
Насколько сильно могут уменьшиться вектор-функции vj (t) по сравнению с исходной системой векторов zj (t)? Ответ на этот вопрос дает следующее утверждение.
Лемма 7.2. Если для некоторого числа C выполнено неравенство
n |
t |
|
|
|
|
! |
Tr A(s) ds t 0, |
|
|zj (t)| C exp |
(7.9) |
|
j=1 |
0 |
|
|
|
|
то существует число r > 0, для которого имеет место оценка |
||
|zj (t)| r|vj (t)|, |
t 0, j = 1, . . . , n. |
(7.10) |
72 |
ГЛАВА 7 |
ДОКАЗАТЕЛЬСТВО.
Из формулы Остроградского – Лиувилля [12] и неравенства (7.9) имеем соотношения
"
"
"det
"
"
|
z1 |
(t) |
zn(t) |
|
" |
||
|
|
" |
|||||
|
|
|
, . . . , |
|
|
" = |
|
| |
| |
| |
| |
||||
|
" |
||||||
|
z1 |
(t) |
zn(t) |
|
" |
||
= (|z1(t)|, . . . , |zn(t)|)−1| det(z1(0), . . . , zn(0))|×
× exp |
t |
Tr A(s)ds C−1| det(z1(0), . . . , zn(0))|, t 0. |
|
0 |
|
Отсюда следует, что для линейного подпространства L(t), натянутого на вектора z1(t), . . . , zm (t) (m < n), найдется число ε (0, 1), такое, что выполнена оценка
|zm+1(t) e(t)| |
1 − ε, t 0, |
(7.11) |
|
|zm+1(t)| |
|||
|
|
для всех e(t) L(t), удовлетворяющих равенству |e(t)| = 1. Соотношения (7.5) можно переписать следующим образом
vj (t) |
j−1 |
|
vi(t)vi(t) |
|
zj (t) |
|
|
|
! |
|
|
|
|
|
|
|
= i=1 |
I − |
|
|
|
. |
(7.12) |
|zj (t)| |
|vi(t)|2 |
|zj (t)| |
|||||
Предположим, что утверждение леммы не имеет места. В этом случае существует последовательность tk → +∞, для которой
vj (tk )
lim = 0.
k→+∞ |zj (tk )|
Но тогда из равенства (7.12) получим, что существует число l < j, для которого
# zj (tk ) |
|
|
|
vl(tk ) $ |
|
|||
lim |
|
|
|
|
|
|
= 0. |
(7.13) |
|
zj (tk ) |
| | |
|
|||||
| |
| |
|
||||||
k→+∞ |
|
− |
vl (tk ) |
|
||||
7.1. МЕТОД ТРИАНГУЛЯЦИИ ПЕРРОНА – ВИНОГРАДА |
73 |
Поскольку vl(t) L(t) соотношения (7.11) и (7.13) противоречат друг
другу. Полученное противоречие доказывает оценку (7.10). |
|
Из приведенных здесь рассуждений легко видеть, что |
усло- |
вие (7.9) является необходимым и достаточным условием существования числа r > 0, для которого выполнена оценка (7.10).
Заметим, что условие (7.9) является необходимым и достаточным условием невырожденности при t → +∞ нормированной фундаментальной матрицы системы первого приближения (7.4):
"
"
lim "det
"
t→+∞ "
|
z1(t) |
|
|
zn(t) |
|
" |
|
|
, · · · , |
|
" |
||
|
z1(t) |
| |
zn(t) |
| |
" |
|
| |
|
| |
" |
|||
Перейдем теперь к описанию процедуры триангуляции Перрона – Винограда.
Введем в рассмотрение унитарную матрицу
U (t) = v1(t) , · · · , vn(t)
|v1(t)| |vn(t)|
и сделаем следующую замену в системе (7.4): z = U (t)w.
Из унитарности матрицы U (t) следует, что для столбцов w(t) матрицы W (t) = (w1(t), . . . , wn(t)) = U (t) Z(t) выполнены соотношения |wj (t)| = |zj (t)|.
Из соотношений (7.5)– (7.7) вытекает, что матрица W (t) имеет следующий треугольный вид:
|v1(t)| · · · |
. |
|
(7.14) |
|
W (t) = |
. . . |
.. |
. |
|
|
0 |
|vn(t)| |
|
|
|
|
|
||
Матрица W (t) является фундаментальной матрицей для системы
dw |
= B(t)w, |
|
(7.15) |
|
dt |
|
|
|
|
где |
|
|
U (t)−1U˙ (t). |
(7.16) |
B(t) = U (t)−1A(t)U (t) |
− |
|||
|
|
|
|
|
74 ГЛАВА 7
Из того факта, что W (t) — верхняя треугольная матрица, сле-
−1 ˙
дует, что и W (t) , W (t) — верхние треугольные матрицы. Поэтому матрица
B(t) = W˙ (t)W (t)−1 |
|
|
||
является верхней треугольной матрицей вида |
|
|
||
(ln |v1(t)|)• |
· · · |
. |
|
(7.17) |
B(t) = |
. . . |
.. |
. |
|
|
|
|
|
|
0 |
|
(ln |vn(t)|)• |
|
|
Покажем, что если A(t) ограничена при t 0, то и матри- |
||||
˙ |
|
|
0. Ограниченность |
|
цы B(t), U (t) и U (t) также ограничены при t |
||||
матрицы U (t) имеет место всегда и очевидна. Поэтому ограничена также матрица U (t)−1A(t) U (t) = U (t) A(t)U (t).
Из унитарности U (t) следует тождество
% |
1 |
˙ |
& |
|
U (t)− |
U (t) |
|
= |
− |
U (t)−1U˙ |
(t). |
(7.18) |
|
|
|
|
Отсюда, из соотношений (7.16) и (7.17) следует, что модуль
−1 ˙
элемента матрицы U (t) U (t) совпадает с модулем некоторого эле-
−1 −1 ˙
мента матрицы U (t) A(t)U (t). Таким образом, матрица U (t) U (t) ограничена при t 0. Отсюда и из равенства (7.16) следует огра-
˙ |
= |
ниченность B(t). Из ограниченности B(t) и равенства U (t) |
|
˙ |
|
= AU (t) − U (t)B(t) следует ограниченность U (t). |
|
Таким образом, здесь доказан следующий результат. |
|
Теорема 7.1 (Перрона о триангуляции линейной системы).
[12]. Систему (7.4) с помощью унитарного преобразования z = = U (t)w можно привести к системе (7.15) с матрицей B(t) вида (7.17).
Если матрица A(t) ограничена при t 0, то и матрицы B(t),
˙
U (t), U (t) также ограничены при t 0.
Докажем теперь еще одну полезную оценку для вектор-функ- ции vn(t).
|
|
7.2. ТЕОРЕМЫ О НЕУСТОЙЧИВОСТИ |
|
|
75 |
|||||||||
Лемма 7.3. Имеет место оценка |
|
|
|
|
|
|
|
|
||||||
|
| |
vn(t) |
|
t |
n |
− |
1 |
vj (τ ) |
| |
|
|
|||
|
|
|
|
|
|
|
||||||||
|
| |
exp Tr A(s) ds ! |
| |
|
. |
(7.19) |
||||||||
|vn(τ )| |
|
|
|
|||||||||||
|
|
|
τ |
j=1 |zj (t)| |
|
|||||||||
ДОКАЗАТЕЛЬСТВО. |
|
|
|
|
|
|
|
|
|
|
|
|||
Из соотношения (7.14) следует равенство |
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
n−1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
det W (t) |
' vj (τ ) |
|
|
|
||||
|
|
|
|vn(t)| |
|
|
| |
|
|
| |
|
|
|
|
|
|
|
|
= |
|
j=1 |
|
|
|
. |
|
|
|
||
|
|
|
vn(τ ) |
det W (τ ) |
n−1 |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
| |
| |
|
' vj (t) |
|
|
|
||||||
|
|
|
|
|
|
|
j=1 | |
|
|
| |
|
|
|
|
Из формулы Остроградского–Лиувилля, (7.15), (7.16) и (7.18) имеем равенства
t |
t |
|
|
det W (t) = det W (τ ) Tr B(s) ds = det W (τ )exp |
Tr A(s) ds. |
τ |
τ |
Из этих равенств и оценки (7.8) сразу следует утверждение лем- |
|
мы. |
|
7.2. Теоремы о неустойчивости
Теорема 7.2. |
[31]. Если выполнено неравенство |
|
|||
|
|
|
t |
|
|
sup |
lim |
1 |
|
Tr A(s)ds − ln |zj (t)| > 0, |
(7.20) |
t |
|
||||
k |
t→+∞ |
|
0 |
j=k |
|
|
|
|
|
|
|
то решение x(t) ≡ 0 системы (7.1) неустойчиво по Красовскому.
76 ГЛАВА 7
ДОКАЗАТЕЛЬСТВО.
Не умаляя общности, можно считать, что в (7.20) супремум по k достигается при k = n.
Сделаем замену x = U (t)y в системе (7.1): |
|
||
|
dy |
= B(t)y + g(t, y), |
(7.21) |
|
dt |
||
|
|
|
|
где B(t) определяется по формуле (7.16) и
g(t, y) = U (t)−1f (t, U (t)y).
Таким образом, последнее уравнение системы (7.21) примет вид
y˙n = (ln |vn(t)|)•yn + gn(t, y). |
(7.22) |
Здесь yn и gn — n-е компоненты векторов y и g.
Предположим теперь, что решение x(t) ≡ 0 устойчиво по Красовскому. Это означает существование для некоторой окрестно-
сти x = 0 числа R > 0 такого, что выполнена оценка |
|
|x(t, x0)| R|x0|, t 0. |
(7.23) |
Здесь x(0, x0 ) = x0. |
|
Из условий (7.2) и (7.23) следует оценка |
|
|g(t, y(t))| κRν |y(0)|ν . |
(7.24) |
Из условия (7.20) по лемме 7.3 получим существование такого числа μ > 0, что при достаточно больших значениях t выполнена оценка
ln |vn(t)| μt. |
(7.25) |
Заметим, что решение yn(t) уравнения (7.22) можно представить в форме
|
|vn(t)| |
|
t |
|vn(0)| |
|
|
|
|
|
|
|
(7.26) |
|||
yn(t) = |vn(0)| |
yn(0) + |
|vn(s)|g(s, y(s))ds . |
|||||
|
|
||||||
0
7.2. ТЕОРЕМЫ О НЕУСТОЙЧИВОСТИ |
77 |
Из оценки (7.25) следует существование числа ρ > 0 такого, что выполнено неравенство
|
t |
|
|
|||
|
|vn(0)| |
|
|
(7.27) |
||
|vn(s)| ds ρ, |
t 0. |
|||||
0 |
|
|
||||
Возьмем теперь начальные условия |
x0 = U (0)y(0) так, чтобы |
|||||
yn(0) = |y(0)| = δ, где число δ удовлетворяет неравенству |
|
|||||
δ> ρκRν δν .
Вэтом случае из соотношений (7.24)–(7.27) получим при достаточно больших t 0 оценку
yn(t) exp(μt)(δ − ρκRν δν ).
Отсюда сразу следует соотношение
lim yn(t) = +∞.
t→+∞
Последнее противоречит предположению об устойчивости по Красовскому тривиального решения системы (7.1). Теорема доказана.
Замечание относительно метода доказательства теоремы 7.2.
Если бы мы предположили наличие устойчивости по Ляпунову и попытались вывести к этому противоречие, так же, как это было сделано в доказательстве по отношению к устойчивости по Красовскому, то в этом случае необходимо было бы доказать неравенство
+∞ |
|vn(0)| |
|
|
|
|
|
(7.28) |
||
|vn(s)|g(s, y(s)) ds = 0. |
||||
yn(0) + |
|
|||
0
78 ГЛАВА 7
Это неравенство легко обеспечить, когда речь идет об устойчивости по Красовскому и является труднопреодолимым препятствием при рассмотрении устойчивости по Ляпунову.
Аналогичная схема сведения задачи к одному скалярному уравнению типа (7.22) была использована Н. Г. Четаевым [65, 66] при получении критериев неустойчивости. И аналогичная трудность в доказательстве неравенства (7.28) имеется в схеме, предложенной Н. Г. Четаевым. Поэтому метод Н. Г. Четаева в настоящее время дает возможность получения критериев неустойчивости по Красовскому.
Методика получения критериев неустойчивости по Ляпунову требует дальнейшего развития. Такое развитие с некоторыми допол-
нительными ограничениями будет дано в теореме 7.4. |
|
Условие (7.20) теоремы 7.2 выполнено, если имеет место нера- |
|
венство |
|
Λ − Γ > 0. |
(7.29) |
Здесь Λ — максимальный характеристический показатель, Γ — коэффициент неправильности.
Условие неустойчивости по Красовскому (7.29) было получено Н. Г. Четаевым [65, 66] при дополнительном требовании аналитичности f (t, x).
Напомним здесь условие устойчивости (6.5) теоремы 6.3, которое, согласно теореме 3.3, может быть записано следующим образом:
(ν − 1)Λ + Γ < 0. |
(7.30) |
Поскольку теоремы 6.1–6.3 являются также критериями устойчивости по Красовскому, то можно сформулировать следующий результат:
Теорема 7.3. Если
Λ < − Γ , (ν − 1)
то решение x(t) ≡ 0 устойчиво по Красовскому, если
Λ > Γ,
то решение x(t) ≡ 0 неустойчиво по Красовскому.
7.2. ТЕОРЕМЫ О НЕУСТОЙЧИВОСТИ |
79 |
Для правильных систем, когда Γ = 0, теорема 7.3 полностью решает задачу об устойчивости по Красовскому в некритическом случае (Λ = 0).
Заметим, что для системы (4.6) имеет место равенство Γ = = Λ + 2a + 1. Поэтому условие (7.29) для системы (4.6) не выполняется.
Теорема 7.4. [23]. Будем предполагать, что для некоторых чисел C > 0, β > 0, αj < β (j = 1, . . . , n − 1) выполнены следующие условия:
1) при n > 2
|
|
|
|
n |
t |
|
|
|
|
|
|
|
|
|
|
|
|
! |
Tr A(s) ds, t 0 |
|
|
|
|
|
|zj (t)| C exp |
(7.31) |
|
|
|
|
|
j=1 |
0 |
|
|
|
|
|
|
|
|
2) |
|
|
|
|zj (t)| C exp(αj (t − τ ))|zj (τ )|, |
|
|
|
|
|
|
(7.32) |
||
|
|
|
|
t τ 0, |
j = 1, . . . , n − 1 |
|
3) |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
1 |
|
t |
n−1 |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
Tr A(s) ds > β + αj , t τ 0. |
(7.33) |
|
|
(t |
− |
τ ) |
|||
|
|
|
τ |
j=1 |
|
|
Тогда нулевое решение системы (7.1) неустойчиво по Ляпунову.
ДОКАЗАТЕЛЬСТВО.
Напомним, что фундаментальная матрица W (t) системы (7.15) имеет вид (7.14). Вектор–столбцы этой матрицы wj (t) имеют вид:
|
w11(t) |
|
0 |
w1 |
(t) = . |
|
.. |
|
|
|
0 |
|
|
|
|
wn−1,1(t) |
|
||
|
|
|
|
|
. |
|
|
|
, wn |
|
.. |
|
|||
, |
· · · |
1 |
(t) = |
1,n |
1(t) |
. |
|
|
− |
|
wn |
|
|||
|
|
|
|
− |
− |
|
|
|
|
|
|
|
0 |
|
|
80 |
ГЛАВА 7 |
|
|
Поэтому матрица |
|
|
|
w11(t) · · · |
wn−1,1(t) |
|
|
|
0 · · · |
|
|
W (t) = . |
. |
|
|
( |
. |
. |
|
. |
. |
||
|
|
|
|
0· · · wn−1,n−1(t)
является фундаментальной матрицей системы
w˙ |
= B(t)w, w |
|
Rn−1 |
|
||
|
|
|
|
|
|
|
с матрицей |
|
|
|
|
|
|
(ln |v1(t)|)• |
· · · |
|
. |
|
||
B(t) = |
|
|
. . . |
|
.. |
. |
|
|
|
|
|
|
|
|
0 |
|
|
|
(ln |vn−1(t)|)• |
|
|
|
|
|
|
||
Из условия (7.31) и тождеств |wj (t)| ≡ |zj (t)|, которые вытекают из унитарности матрицы U (t), следуют оценки
|wj (t)| C exp(αj (t − τ ))|wj (τ )|,
(7.34)
t τ 0, j = 1, . . . , n − 1.
Кроме того, из условия (7.31) по лемме 7.2 следуют оценки (7.10), а из условий (7.32) и (7.33) по лемме 7.3 вытекает оценка
|
|
|vn(t)| |
|
n−1 |
|
|||
|
|
exp(β(t − τ ))C1−n ! |
|vj (τ )| |
|
(7.35) |
|||
|
|
|vn(τ )| |
|
|||||
|
|
|
j=1 |zj (τ )| |
|
||||
|
|
|
|
t τ 0. |
|
|
|
|
Из (7.10) и (7.35) получим неравенство |
|
|
|
|
||||
|
|vn(t)| |
exp(β(t − τ ))(Cr)1−n, |
t τ 0. |
(7.36) |
||||
|
|
|||||||
|
|vn(τ )| |
|
|
|
|
|
||
