Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Хаотическая динамика и классическая теория устойчивости движения

.pdf
Скачиваний:
0
Добавлен:
15.08.2026
Размер:
1 Мб
Скачать

7.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 C1| 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 )

 

" > 0.

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(τ )|

 

 

 

 

 

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]