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

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

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

УСТОЙЧИВОСТЬ ПО ЖУКОВСКОМУ

101

Если для некоторых положительных чисел ρ и δ справедливы нера-

венства

tj+1

 

 

 

 

 

 

 

λ(x(t, x0)) dt ρ, j,

(8.29)

 

tj

 

 

z H(x(t, x0))z −δ|z|2 , t 0, z Rn,

(8.30)

то существуют положительные числа C и ε такие, что

 

 

|z(t)| Ceεt|z(0)|, t 0.

(8.31)

ДОКАЗАТЕЛЬСТВО.

Заметим сначала, что из ограниченности x(t, x0) и непрерывности λ(·) следует ограниченность при t 0 функции λ(x(t, x0 )).

Поэтому существует число α > 0, для которого

 

" t

"

 

 

 

 

 

 

 

 

"

"

 

 

 

 

 

 

 

 

 

"

 

 

 

 

 

 

 

 

"

α,

 

j,

 

t

 

[tj , tj+1] .

(8.32)

"

"

"

λ(x(t, x0)) dt"

 

 

 

"tj

"

 

 

 

 

 

 

 

 

"

"

 

 

 

 

 

 

 

 

Соотношение (8.25) можно переписать следующим образом:

 

 

 

t

 

 

 

 

 

 

 

λ(x(τ, x0 )) 0.

 

V (t) exp

 

 

0

 

 

 

 

 

 

 

 

Поэтому имеет место оценка

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

λ(x(τ, x0 )) dτ .

 

V (t) V (0) exp

(8.33)

0

Из неравенств (8.26) и (8.32) следует существование чисел ε > 0 и β > 0 таких, что

t

λ(x(τ, x0)) dτ −2εt + β, t 0.

(8.34)

0

Отсюда, из (8.33) и (8.27) следует оценка (8.28).

102

ГЛАВА 8

Если выполнены условия (8.29), то из (8.32), то существуют числа ε > 0 и β > 0 такие, что

t

λ(x(τ, x0 ))2εt − β, t 0.

0

Следовательно, неравенства (8.30), (8.33) и равномерная ограни-

ченность H(x, (t, x0)) влекут за собой оценку (8.31). Лемма доказана.

Ясно, что соотношение (8.28) соответствует асимптотической устойчивости по Жуковскому и (8.31) — неустойчивости по Жуковскому.

Заметим, что если (8.25) выполнено на решениях системы (8.8) лишь в некоторой окрестности Ω(0) точки z = 0, то утверждение леммы 8.3 о выполнении оценки (8.28) будет справедливым также для z(0) из некотрой достаточно малой окрестности Φ(0) Ω(0). Здесь, как и обычно это делают в прямом методе Ляпунова, предполагая, что z(0) дoстаточно мало, получим выполнение оценки (8.25) на достаточно большом промежутке времени [0, T ] (это следует, например, из того факта, что скорость изменения решения z˙(t) в окрестности Ω(0) мала). Тогда получим, что при t = T будет выполнена оценка (8.28). Из нее следует, что для больших T выполнено неравенство |z(T )| |z(0)|. Следовательно, при малых z(0) решение z(t) не может покинуть при некотoром t > 0 окрестность Ω(0). Следовательно, для этого z(t) оценка (8.25) выполнена при всех t 0. Но тогда и оценка (8.28) выполнена также при всех t 0.

Рассмотрим теперь симметричную матрицу M и ее собственные значения λ1 . . . λn.

Лемма 8.4. [86]. Для всех z Rn и y Rn имеют место неравенства

УСТОЙЧИВОСТЬ ПО ЖУКОВСКОМУ

103

(λn + λn−1)(|z|2 |y|2 (z y)2)

z M z|y|2 + |z|2y M y − 2z yz M y

(λ1 + λ2)(|z|2 |y|2 (z y)2)). (8.35)

ДОКАЗАТЕЛЬСТВО.

 

 

 

 

 

 

Поскольку ортогональное преобразование

 

,

,

S = S

1

 

 

 

 

z

= Sz

y = Sy

 

не изменяет скалярного произведения и может привести

M к диагональному виду:

 

 

 

 

 

 

 

 

 

λ1

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S M S =

. .

.

,

 

 

(8.36)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

λn

 

 

 

не умаляя общности, можно считать, что в доказываемых неравенствах (8.35) матрица M имеет вид (8.36). Но тогда эти неравенства примут следующий вид:

(λn + λn−1)(

 

zj2 yi2 (

zj yj )2)) =

 

 

 

 

j

i

j

 

 

 

 

 

1

 

n

n

 

 

 

 

(λn + λn−1)(

 

2

 

 

=

2

 

(zj yi − ziyj ) )

 

 

 

 

 

j=1 i=1

 

 

1

 

 

 

 

n n

 

 

 

 

2

 

2

2

(λi + λj )(zj yi − ziyj ) =

λj (zj yi − ziyj ) =

 

 

i,j

 

 

 

j=1 i=1

 

nn

=λj (zj zj yiyi + ziziyj yj 2zj ziyj yi) =

 

 

j=1 i=1

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

n

n

 

 

 

2

 

(λ1

+ λ2)(

 

2

=

2

(λi + λj )(zj yi − ziyj )

2

 

 

(zj yi − ziyj ) ) =

 

i,j

 

 

 

 

 

 

j=1 i=1

 

 

 

 

 

 

 

 

 

2

 

2

 

2

 

 

= (λ1 + λ2)(

zj

 

yi

(

zj yj ) ).

j

i

j

104

ГЛАВА 8

Поскольку приведенные выше неравенства очевидны, это доказывает лемму 8.4.

Возьмем теперь в качестве функции V (t) функцию

V (t) = |f (x(t, x0))|2|z(t)|2 .

Из соотношений (8.8)(8.10) следует, что для любого числа ν > 0 существует некоторая окрестность Ω(0) точки z = 0, в которой

V˙ (t)

|

f (x(t, x

 

)) 2z(t)

 

∂f

(x(t, x

 

))+

∂f

(x(t, x

 

))

 

z(t)+

 

 

 

 

 

 

 

 

 

 

 

 

 

0

|

 

 

 

∂x

0

 

 

 

 

 

∂x

 

 

 

0

 

 

 

 

 

z(t) 2f (x(t, x

 

))

∂f

(x(t, x

 

)) +

∂f

(x(t, x

 

 

))

 

f (x(t, x

))+

 

 

 

 

 

 

 

 

| |

 

 

0

 

 

 

∂x

 

 

0

 

∂x

 

 

 

 

0

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ ν|z(t)|2 .

Используя лемму 8.4 и полагая M =

 

∂f

 

+

 

∂f

 

, z = z(t), y =

 

∂x

 

∂x

= f (x(t, x0)), получим оценку

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

˙

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ν

 

 

 

 

 

 

 

 

λ1(x(t, x0)) + λ2(x(t, x0 )) +

 

 

 

 

 

 

 

 

 

 

V (t).

 

V (t)

 

|f (x(t, x0))|2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку для некоторого γ > 0 выполнена оценка |f (x(t, x0))| γ,t 0 и ν выбрано сколь угодно малым, то для выполнения условий леммы 8.3 достаточно, чтобы имела место оценка

tj+1

(λ1(x(t, x0)) + λ2(x(t, x0))) dt < −α, j

(8.37)

tj

 

 

 

для некоторого положительного α.

 

 

 

В самом деле, здесь

 

 

 

ρ = α −

ν

κ.

 

γ2

 

УСТОЙЧИВОСТЬ ПО ЖУКОВСКОМУ

105

Поэтому при малом ν выполнены условия (8.25)(8.27) в некоторой

окрестности Ω(0).

 

 

 

 

Таким образом, используя

характеристические значения

λ1(x) . . . λn(x) матрицы

 

 

 

 

 

∂f

 

+

∂f

 

,

 

∂x

 

 

 

∂x

 

cформулируем следующий результат.

Теорема 8.5 (Аналог признака Пуанкаре). Если для последовательности tj +выполнены неравенства tj+1 − tj κ и (8.37) с некоторыми положительными числами α и κ, то траектория x(t, x0) асимптотически устойчива по Жуковскому.

Аналогично, рассматривая функцию

V(t) = −|f (x(t, x0))|2|z(t)|2

ипримeняя леммы 8.3 и 8.4, получим следующий результат.

Теорема 8.6 (Аналог признака Пуанкаре). Если для последовательности tj +выполнены неравенства tj+1 − tj κ и

tj+1

(λn(x(t, x0)) + λn−1(x(t, x0 ))) dt α, j

(8.38)

tj

 

с некоторыми положительными числами α и κ, то траектория x(t, x0) неустойчива по Жуковскому.

Для n = 2 и T -периодического решения x(t, x0) из теорем 8.4 и 8.5 сразу следует широко известный [12, 92] признак Пуанкаре.

Теорема 8.7. Если для траектории x(t, x0) выполнено неравен-

ство

T

∂f

Tr ∂x (x(t, x0)) dt < 0,

0

106

ГЛАВА 8

то она асимптотически орбитально устойчива. Если

T

∂f

Tr ∂x (x(t, x0 )) dt > 0,

0

то она орбитально неустойчива.

Рассмотрим теперь функцию

V (t) = |z(t)|2 .

Из соотношений (8.8)–(8.10) следует, что для любого числа ν > 0 существует некоторая окрестность Ω(0) точки z = 0, в которой

V˙ (t)

 

2z(t)

∂f

(x(t, x

 

))z(t) + ν z(t) 2 .

∂x

 

 

 

 

0

| |

Таким образом, здесь будут выполнены все условия леммы 8.3 с δ = 1, κ = 1 и с малым положительным ρ, если имеет место неравенство

z(t)

∂f

(x(t, x0))z −α|z|2,

 

∂x

(8.39)

 

z

 

z f (x(t, x

)) = 0, z = 0

 

 

{

0

}

 

с некоторым положительным α. В самом деле, здесь достаточно положить λ(x(t, x0)) = α − ν, ρ = α − ν.

Таким образом, здесь получен следующий широко известный результат Г. Борга.

Теорема 8.5 (Борг [54, 77, 86]). Если выполнено условие (8.39), то траектория x(t, x0) асимптотически орбитально устойчива.

Дальнейшее развитие теории первого приближения для устойчивости и неустойчивости по Жуковскому имеется в [28, 91–93].

ГЛАВА 9

Функции Ляпунова в оценках размерности аттракторов

Числовыми характеристиками для гармонических колебаний являются амплитуда, период, частота, для периодических колебаний — их период. В результате многочисленных исследований стало ясно, что и более сложные колебания обладают числовыми характеристиками. Это — размерности аттракторов, соответствующих ансамблям таких колебаний.

Развитая в первой половине двадцатого века теория топологической размерности [88, 89] мало пригодна для того, чтобы дать шкалу размерностных характеристик аттракторов. Дело в том, что топологическая размерность может принимать лишь целочисленные значения, и составленная таким образом шкала размерностных характеристик оказывается слишком бедной.

Для исследования аттракторов более адекватной оказалась хаусдорфова размерность множеств. Эта размерностная характеристика может принимать любые неотрицательные значения и на таких привычных объектах в евклидовом пространстве как гладкая кривая, поверхность, счетное множество точек, совпадает с топологической размерностью.

Перейдем теперь к определению хаусдорфовой размерности множеств.

Рассмотрим компактное метрическое пространство X с метрикой ρ, подмножество E X и числа d 0, ε > 0.

Покроем E шарами с радиусами rj < ε и обозначим

μH (E, d, ε) = inf rjd,

j

где инфимум берется по всем таким ε-покрытиям E.

108

ГЛАВА 9

Очевидно, что μH (E, d, ε) не убывает при уменьшении ε. Поэтому существует предел (возможно бесконечный)

μH (E, d) = lim μH (E, d, ε).

ε→0

Определение 9.1. Функция μH (·, d) называется d-мерой Хаусдорфа.

При фиксированном d функция μH (E, d) обладает всеми свойствами внешней меры на X.

При фиксированном множестве E функция μH (E, ·) обладает следующим свойством. Существует такое dkp [0, ∞], что

μH (E, d) = ∞, d < dkp, μH (E, d) = 0, d > dkp.

Если X Rn, то dkp n. Здесь Rn — евклидово n-мерное пространство.

Положим

dimH E = dkp = inf

d μH (E, d) = 0 .

 

|

Определение 9.2. Величина dimH E называется размерностью Хаусдорфа множества E.

ПРИМЕР 9.1. Рассмотрим канторово множество из примера 1.3:

+

E = Ej ,

j=0

где E0 = [0, 1] и Ej состоит из 2j отрезков длины 3−j , получаемых из отрезков, входящих в Ej−1, удалением из них середин открытых интервалов длины 3−j (рис. 1.6).

В классической теории топологической размерности хорошо известно, что dimT E = 0.

ФУНКЦИИ ЛЯПУНОВА В ОЦЕНКАХ РАЗМЕРНОСТИ

109

Непосредственно из определения хаусдорфовой размерности легко получить, что μH (E, d) = 1 при d = log 2/ log 3 = 0,6309 . . . и, следовательно,

dimH E =

log 2

 

log 3 .

Топологическая размерность является инвариантом по отношению к гомеоморфизмам. Хаусдорфова размерность является инвариантом по отношению к диффеоморфизмам и нецелая хаусдорфова размерность не является инвариантом по отношению к гомеоморфизмам [88].

При исследовании аттракторов динамических систем в фазовом пространстве часто используются гладкие замены координат. Поэтому в таких рассмотрениях достаточно инвариантности по отношению к диффеоморфизмам.

Хорошо известно, что dimT E dimH E. Канторовское множество E дает пример, когда это неравенство является строгим.

Приведем теперь два эквивалентных определения фрактальной размерности.

Обозначим через (E) минимальное число шаров радиуса ε, необходимое для покрытия множества E X.

Введем в рассмотрение числа d 0, ε > 0 и положим

μF (E, d, ε) = N(E)εd ,

μF (E, d) = lim μF (E, d, ε).

ε→0

Определение 9.3. Фрактальной размерностью множества E называется величина

dimF E = inf{d| μF (E, d) = 0}.

Заметим, что это определение проведено по схеме, аналогичной определению хаусдорфовой размерности. Здесь только покрытие проводится шарами одного и того же радиуса ε.

110

ГЛАВА 9

 

 

 

 

Определение 9.4. Фрактальной размерностью множества E

называется величина

 

 

 

 

 

 

dimF E =

 

 

log N(E)

.

 

 

 

lim

 

 

 

 

 

 

 

 

ε→0 log(1)

 

 

 

 

Легко видеть, что

 

 

 

 

 

 

dimH E dimF E.

 

 

 

 

ПРИМЕР 9.2. Для X = [0, 1] и E = 0, 1, 21

, 31, . . .

}

 

 

 

{

1

 

 

dimH E = 0, dimF E =

 

 

 

2 .

 

Обобщения схем введения хаусдорфовой и фрактальной меры и размерности и определения различных метрических размерностных характеристик даны в [53]. Оказалось [19, 76, 81, 87, 92, 94], что верхней оценкой хаусдорфовой и фрактальной размерности инвариантных множеств является ляпуновская размерность. Дадим здесь определение ляпуновской размерности.

Рассмотрим непрерывно дифференцируемое отображение F открытого множества U Rn в Rn. Обозначим через TxF матрицу Якоби отображения F в точке x. Из непрерывной дифференцируемости F следует соотношение

F (x + h) − F (x) = (TxF )h + o(h).

В дальнейшем будем предполагать, что множество K U инвариантно относительно преобразования F : F (K) = K.

Введем в рассмотрение сингулярные числа n × n-матрицы A

α1(A) . . . αn(A).

Напомним, что сингулярным числом матрицы A называется квадратный корень из собственного значения матрицы A A. Здесь через обозначена операция транспонирования (в вещественном случае) или эрмитова сопряжения.

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