Хаотическая динамика и классическая теория устойчивости движения
.pdfУСТОЙЧИВОСТЬ ПО ЖУКОВСКОМУ |
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 )) dτ 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 ))dτ 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, 2−1 |
, 3−1, . . . |
} |
||||
|
|
|
{ |
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. Здесь через обозначена операция транспонирования (в вещественном случае) или эрмитова сопряжения.
