Хаотическая динамика и классическая теория устойчивости движения
.pdfФУНКЦИИ ЛЯПУНОВА В ОЦЕНКАХ РАЗМЕРНОСТИ |
121 |
Из этих равенств ясно, что αk (x)2 являются нулями полинома
λ2 − (4x2 + 2b)λ + b2,
который является характеристическим полиномом матрицы
−2x |
√ |
|
−2x |
√ |
|
|
||||
b |
|
b . |
||||||||
√ |
|
|
0 |
√ |
|
|
0 |
|
||
|
b |
|
b |
|
||||||
Таким образом, доказаны формулы (9.11).
Из теоремы 9.1 следует, что если существуют число s [0, 1) и непрерывно дифференцируемая положительная на K функция p(ξ), такие, что выполнено неравенство
|
|
p(F (ξ)) |
1+s |
|
||
sup α1(x)α2(x)s |
|
< 1, |
(9.12) |
|||
p(ξ) |
||||||
ξ K |
|
|
|
|||
|
|
|
|
|||
то |
|
|
|
|
||
|
|
dimL K 1 + s. |
|
|
||
Положим |
|
|
|
|
||
|
p(ξ)1+s = eγ(1−s)(x+by) , |
|
||||
где γ — положительный варьируемый параметр. |
|
|||||
Легко видеть, что выполнено соотношение |
|
|||||
p(F (ξ)) |
1+s = eγ(1−s)(a+(b−1)x−x2 ). |
|
||||
|
p(ξ) |
|
||||
|
|
|
|
|
||
Отсюда следует, что после логарифмирования условие (9.12) принимает вид:
sup , ln α1(x) + s ln α2(x) + γ(1 − s)(a + (b − 1)x − x2)- =
K
= sup ,(1 − s) ln α1(x) + s ln b + γ(1 − s)(a + (b − 1)x − x2)- < 0.
K
122 |
|
|
|
ГЛАВА 9 |
|
|
Это неравенство будет выполнено, если |
|
|||||
s ln b + (1 − s)ϕ(x) < 0, x (−∞, +∞), |
|
|||||
где |
|
|
|
|
|
|
.* |
|
|
|
/ |
2 |
). |
x |
2 |
|
||||
ϕ(x) = ln |
|
+ b + |x| |
+ γ(a + (b − 1)x − x |
|||
Из неравенств γ > 0, b − 1 < 0 следует оценка
ϕ(−|x|) ϕ(|x|).
Поэтому достаточно рассмотреть экстремальные точки функции ϕ(x) при x (−∞, 0]. Ясно, что на этом множестве
ϕ (x) = |
√ |
− 1 |
+ γ,(b |
− |
1) |
− |
2x-, ϕ (x) < 0. |
||||
|
x |
2 |
+ b |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|||
Выбрав |
|
|
|
|
|
|
1 |
|
|
|
|
|
γ = |
|
|
|
|
, |
|||||
|
|
|
|
|
|
|
|||||
*
(b − 1 − 2x−) x−2 + b
получим, что ϕ (x−) = 0 и, следовательно, при таком выборе γ выполнено неравенство
% |
* |
|
|
& |
2 |
|
|||
ϕ(x) ln |
|
x− + b |
+ |x−| |
= ln α1(x−). |
Таким образом, неравенство (9.12) выполнено для всех s, удовлетворяющих условию
s > |
ln α1(x−) |
. |
(9.13) |
|
ln α1(x−) − ln b |
|
|
Следовательно, выполнена оценка
dimL K 1 + s
для всех s, удовлетворяющих неравенству (9.13). Переходя к пределу, получим неравенство
dimL K 1 + |
1 |
. |
(9.14) |
|
1 − ln b/ ln α1(x−) |
||||
|
|
|
ФУНКЦИИ ЛЯПУНОВА В ОЦЕНКАХ РАЗМЕРНОСТИ |
123 |
Заметим, что точка x = x−, y = x− является неподвижной для отображения F . При этом
α |
(x |
− |
)α (x |
− |
)s = 1, |
(9.15) |
|
1 |
|
2 |
|
|
|
||
где |
|
|
1 |
|
|
|
|
s = |
|
|
|
. |
|
||
1 − ln b/ ln α1(x−) |
|
||||||
Легко видеть, что величины α1(x−) и α2(x−) являются собственными значениями матрицы Якоби Tξ F отображения F в неподвиж-
ной точке y = x = x−: |
−2x− b . |
|
Tξ F = |
||
|
1 |
0 |
Отсюда и из равенства (9.15) по лемме 9.2 следует, что локальная ляпуновская размерность последовательности отображений F i в этой стационарной точке равна
1 + |
1 |
. |
(9.16) |
1 − ln b/ ln α1(x−)
Отсюда и из неравенства (9.14) следует утверждение теоремы 9.3.
Заметим, что при a = 1, 4, b = 0, 3 из теоремы 9.3 имеем равен-
ство
dimL K = 1, 49532 . . . |
|
Рассмотрим систему Лоренца |
|
x˙ = −σx + σy, |
|
y˙ = rx − y − xz, |
(9.17) |
z˙ = −bz + xy, |
|
где r, b, σ — положительные числа. |
|
Предположим, что выполнены неравенства r > 1, |
|
σ + 1 b 2. |
(9.18) |
Рассмотрим оператор сдвига по траектория системы (9.17) GT , где T — произвольное положительное число. Пусть K — инвариантное множество относительно этого оператора GT . Будем предполагать, что множество K содержит стационарную точку x = y = z = 0. Такое множество изображено на рисунках 2.3 и 2.4.
124 ГЛАВА 9
Дадим здесь формулу для ляпуновской размерности dimL K множества K относительно последовательности отображений (GT )i.
Теорема 9.4. Пусть выполнены неравенства (9.18) и |
|
||||||
rσ2(4 − b) + 2σ(b − 1)(2σ − 3b) > b(b − 1)2. |
(9.19) |
||||||
Тогда |
2(σ + b + 1) |
|
|||||
dimL K = 3 − |
(9.20) |
||||||
|
|
|
|
. |
|||
σ + 1 + * |
|
|
|
||||
(σ |
− |
1)2 + 4rσ |
|||||
|
|
|
|
|
|
||
ДОКАЗАТЕЛЬСТВО.
Матрица Якоби правой части системы (9.17) имеет вид:
|
−σ |
σ |
0 |
|
||||
J = |
|
|
z |
|
1 |
|
x |
|
r |
− |
− |
− |
. |
||||
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
y |
|
x |
−b |
|
||
|
|
|
|
|||||
Введем матрицу следующего вида:
|
|
−a−1 |
0 |
0 |
|
|
S = |
|
σ−1(b |
|
1) 1 |
0 |
|
|
− |
, |
||||
|
− |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
0 |
1 |
|
|
где a = |
|
|
σ |
|
|
|
. В этом случае |
|
|||||
* |
|
|
|
|
|
|
|||||||
|
rσ + (b |
− |
1)(σ |
− |
b) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
b − σ − 1 |
σ/a |
||||||
|
|
|
|
|
|
|
σ |
|
|
|
|
|
|
|
|
SJ S−1 |
|
|
|
− |
az |
− |
b |
||||
|
|
|
|
|
a |
||||||||
|
|
= |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a(b − 1) |
|
|
|
||
|
|
|
|
ay + |
x |
x |
|||||||
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
σ |
|
|
|
0
−x
.
−b
ФУНКЦИИ ЛЯПУНОВА В ОЦЕНКАХ РАЗМЕРНОСТИ |
125 |
Поэтому характеристический полином матрицы
|
|
|
|
|
1 |
((SJ S−1) + (SJ S |
||||
|
|
|
|
|
2 |
|
|
|
|
|
|
b − σ − 1 |
σ |
− |
az |
1 |
|||||
|
|
a |
2 |
2 |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
σ |
|
|
az |
|
|
|
|
|
|
|
− |
|
|
− |
b |
|
|||
|
|
a |
2 |
|
|
|
||||
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
a(b − 1) |
|
|
|
|
|||
|
ay + |
|
x |
0 |
|
|
||||
|
2 |
|
|
σ |
|
|
|
|||
имеет вид:
#
(λ + b) λ2 + (σ + 1)λ + b(σ + 1 − b)
−1)) =
a(b − 1) ay + σ x
0
−b
|
|
|
|
$ |
− |
σ |
− |
az |
− |
a |
2 |
|||
|
# |
|
|
|
− |
a(b − 1) |
x + |
ay |
2σ |
2 |
Отсюда следует, что собственными значениями матрицы
$2
.
|
|
1 |
[(SJ S−1) + (SJ S−1)] |
|
|
||||||
|
|
2 |
|
|
|
|
|
|
|
|
|
являются величины |
|
|
λ2 = −b, |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||
|
σ + 1 |
|
1 |
# |
2 |
|
2σ |
|
2 |
||
λ1,3 = − |
|
± |
2 |
(σ + 1 − |
2b) + |
|
|
− az + |
|||
2 |
a |
||||||||||
|
|
|
|
|
a(b − 1) |
|
2 1/2 |
||||
|
|
|
|
+ |
|
|
|
|
x + ay |
. |
|
|
|
|
|
|
|
σ |
|
||||
Из соотношений (9.18) легко видеть, что λ1 λ2 λ3.
126 |
ГЛАВА 9 |
Введем в рассмотрение функцию ляпуновского типа
|
|
1 |
|
|
|
|
|
|
|
|
|
(b − 1)2 |
|
|
|
||
v(x, y, z) = |
2 |
(1 −s) γ1x |
2 |
+ γ2 |
|
|
2 |
+ z |
2 |
− |
|
2 |
+ γ3z , |
||||
2 aθ |
|
y |
|
|
|
x |
|
||||||||||
|
|
|
σ2 |
|
|||||||||||||
где s (0, 1), |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
) |
|
|
|
|
|
|
|
−1 |
|
|
|
|||
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|||||
θ |
2 |
= 2 (σ + 1 − 2b)2 + |
2σ |
|
, |
|
|
|
|||||||||
|
|
|
a |
|
|
|
|||||||||||
γ3 = − |
4σ |
, γ2 |
= |
a |
, |
ab |
2 |
||||
|
1 |
# |
rσ − (b − 1)2 |
||
γ1 = − |
2σ |
2γ2 |
|
|
σ |
+ γ3 |
+ 2 |
a(b − 1) |
$. |
|
σ |
||||
|
|
|
Рассмотрим выражение
, |
- |
2λ1 + λ2 + sλ3 + v˙ =
=−(σ + 1 + 2b) − s(σ + 1) + (1 − s)ϕ(x, y, z),
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
2σ |
|
|
|
2 |
|
|
|
|
|
||||||
ϕ(x, y, z) = (σ + 1 − 2b) |
+ |
|
|
|
|
− az + |
|
|
|
|||||||||
a |
|
|
|
|||||||||||||||
|
|
+ |
a(b − 1) |
x + ay |
2 |
|
1/2+ |
|
|
|||||||||
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
σ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 3 |
|
|
|
|
a(b − 1)2 |
|
|
2 |
|
2 |
|
||||||
+ th |
|
−2aγ1σ + 2γ2 |
|
|
|
σ |
|
|
x |
|
− 2aγ2y |
− |
||||||
|
|
|
|
|
rσ |
− |
(b |
− |
1)2 |
+ γ3 |
|
4 |
||||||
− 2aγ2bz2 + a |
2σγ1 + 2γ2 |
|
|
|
|
|
xy − γ3abz . |
|||||||||||
|
|
|
σ |
|
|
|
||||||||||||
ФУНКЦИИ ЛЯПУНОВА В ОЦЕНКАХ РАЗМЕРНОСТИ |
127 |
Используя очевидное неравенство |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
√ |
|
|
1 |
+ th2 u, |
|
|
|
|
|
|
||||||||
|
|
u |
|
|
|
|
|
|
||||||||||||
|
|
4 th2 |
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
получим оценку |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
23 |
|
|
|
|
2 |
|
2σ |
2 |
|
|
|
|||||||
ϕ(x, y, z) |
|
+ th |
|
(σ + 1 − 2b) |
+ |
|
|
a |
+ |
|
|
|||||||||
4 th2 |
|
|
|
|||||||||||||||||
# |
|
|
|
|
a(b − 1)2 |
|
a2(b − 1)2 |
$ |
2 |
|
||||||||||
+ −2aγ1σ + 2γ2 |
|
|
σ |
|
+ |
|
|
|
σ2 |
|
|
x |
+ |
|||||||
|
+ [a2 − 2aγ2]y2 + [a2 − 2γ2ab]z2 + |
|
|
|
||||||||||||||||
+ a |
2σγ1 + 2γ2 |
rσ − (b − 1)2 |
+ γ3 |
+ |
|
|||||||||||||||
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
σ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+ 2a2 |
b − 1 |
xy |
− |
[γ3ab + 4σ]z4. |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
σ |
|
|
|
|
|
||
Заметим, что параметры γ1, γ2 и γ3 выбраны так, чтобы выполнялось соотношение
ϕ(x, y, z) |
1 |
|
23 |
2 |
|
|
2σ |
2 |
|
|
|
||
|
+ th |
(σ + 1 − 2b) |
+ |
|
a |
+ |
|
|
|
||||
4 th2 |
|
|
|
||||||||||
|
|
|
|
a(b − 1)2 |
|
a2(b − 1)2 |
2 |
4 |
|||||
|
|
+ −2aγ1σ + 2γ2 |
|
σ |
|
|
+ |
|
σ2 |
x |
|
. |
|
Легко доказать, что из условия (9.19) при определенных выше параметрах γ1, γ2, γ3 имеем неравенство
|
a(b − 1)2 |
|
a2(b − 1)2 |
|
−2aγ1σ + 2γ2 |
|
+ |
|
0. |
σ |
σ2 |
|||
Таким образом, при всех x, y, z
ϕ(x, y, z) *4rσ + (σ − 1)2 .
128 |
|
|
|
|
ГЛАВА 9 |
||
Отсюда следует, что для любого числа |
|||||||
|
|
= * |
|
|
− 2b − σ − 1 |
||
s < s |
0 |
4rσ + (σ − 1)2 |
|||||
|
|
|
* |
|
+ σ + 1 |
||
|
|
|
4rσ + (σ − 1)2 |
||||
существует число ε > 0 такое, что при всех x, y, z выполнена оценка
λ1(x, y, z) + λ2(x, y, z) + sλ3(x, y, z) + v˙ (x, y, z) < −ε.
Устремив s к s0 справа и используя теорему 9.2, получим оценку
dimL K 3 − |
2(σ + b + 1) |
|
. |
(9.21) |
|||
σ + 1 + * |
|
|
|
|
|||
(σ |
− |
1)2 |
+ 4rσ |
||||
|
|
|
|
|
|
|
|
Используя лемму 9.3, легко видеть, что локальная ляпуновская размерность стационарной точки x = y = z = 0 системы (9.17) равна
3 − |
2(σ + b + 1) |
|
. |
(9.22) |
|||
σ + 1 + * |
|
|
|
|
|||
(σ |
− |
1)2 |
+ 4rσ |
||||
|
|
|
|
|
|
|
|
Из соотношений (9.21) и (9.22) следует формула (9.20). |
|
||||||
Использование аналогичного подхода с построением функций Ляпунова позволили получить формулы ляпуновской размерности для аттракторов диссипативного отображения Чирикова [29].
ГЛАВА 10
Частотный критерий cлабой экспоненциальной неустойчивости на аттракторах дискретных систем
В настоящей главе предлагается подход, связанный с нижними частотными оценками ляпуновских экспонент [36–39, 97]. Этот подход развивается в духе классической теории устойчивости и современной нелинейной теории управления: основные результаты формулируются в терминах передаточных функций и частотных характеристик линейных частей рассматриваемых систем. На характеристики нелинейных блоков ϕ(σ) накладываются условия, близкие к традиционным «неравенствам секторного типа» [1, 6, 11, 27, 51, 92, 94], характерным для теории абсолютной устойчивости нелинейных систем управления. Однако для неустойчивости на аттракторах оказалось необходимым «вывернуть наизнанку» такие секторные условия. Вместо привычных и традиционных условий ϕ(σ)/σ [k1, k2] или ϕ [k1, k2] для всех σ R1 здесь для почти всех σ требуется выполнение условий ϕ [k1, k2]. При этом на некоторых множествах ненулевой меры выполнены как неравенство ϕ < k1, так и неравенство ϕ > k2. Для непрерывных систем такие условия и соответствующие им частотные неравенства для передаточных функций линейной части выделяют в пространстве параметров пустые множества. Однако для дискретных систем такое рассмотрение позволяет получить частотные критерии неустойчивости на аттракторах. Здесь приведены примеры конкретных дискретных динамических систем первого и второго порядка, для которых по-
130 |
ГЛАВА 10 |
лученный частотный критерий дает условия неустойчивости на аттракторах.
Для дискретных систем обычно нет той наглядности представлений в фазовом пространстве, которая характерна для непрерывных систем с гладкими траекториями, заполняющими это пространство. Однако этот недостаток может быть частично ликвидирован, если вместо наблюдения за изменением в дискретном времени расстояния между двумя решениями рассматривать в начальный момент времени отрезок прямой, связывающий начальные данные, и дальнейшие его итерации — последовательность непрерывных кривых. Уменьшение или увеличение длин этих кривых дает информацию об устойчивости или неустойчивости и о других свойствах решений дискретных систем. Рассмотрение длин таких непрерывных кривых может быть положено и в основу определений устойчивости и неустойчивости, что особенно важно для неинъективных отображений, определяющих соответствующие дискретные системы (а именно они являются
внастоящее время широко известными генераторами хаоса). Дело
втом, что в некоторых случаях две траектории могут «слипнуться» за конечное число итераций. Однако итерируемый отрезок, связывающий начальные данные, может иметь все большую длину даже с учетом такого слипания.
Оказалось, что механизм растягивания в ряде случаев является достаточно сильным и его можно обнаружить и оценить, применяя весь богатый арсенал прямого метода Ляпунова, который можно соединить с эффективным алгебраическим аппаратом: лем-
мой Калмана – Сеге и леммой Шепелявого. Последние дают оценки растяжения длин итерируемых кривых в терминах частотных неравенств, оперирующих с передаточными функциями линейных блоков.
Эффект «слипания», который является менее сильным, чем «растягивание» в малых окрестностях рассматриваемых траекторий, усиливается по мере увеличения размеров окрестности (а это, как правило, неизбежно при превалировании эффекта растягивания). В какой-то момент окрестность становится такой большой, что «сли-
