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

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

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

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

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) x2 + 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:

2xb .

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) и

 

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 + 4

 

 

 

 

 

 

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

Матрица Якоби правой части системы (9.17) имеет вид:

 

−σ

σ

0

 

J =

 

 

z

 

1

 

x

 

r

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

x

−b

 

 

 

 

 

Введем матрицу следующего вида:

 

 

−a1

0

0

 

S =

 

σ1(b

 

1) 1

0

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

1

 

где a =

 

 

σ

 

 

 

. В этом случае

 

*

 

 

 

 

 

 

 

+ (b

1)(σ

b)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b − σ − 1

σ/a

 

 

 

 

 

 

 

σ

 

 

 

 

 

 

 

SJ S1

 

 

 

az

b

 

 

 

 

 

a

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a(b − 1)

 

 

 

 

 

 

 

ay +

x

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ

 

 

0

−x

.

−b

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

125

Поэтому характеристический полином матрицы

 

 

 

 

 

1

((SJ S1) + (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 S1) + (SJ S1)]

 

 

 

 

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

 

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 + 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

 

21σ + 2γ2

 

 

 

σ

 

 

x

 

22y

 

 

 

 

 

(b

1)2

+ γ3

 

4

22bz2 + 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

 

+ 21σ + 2γ2

 

 

σ

 

+

 

 

 

σ2

 

 

x

+

 

+ [a2 22]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

 

 

+ 21σ + 2γ2

 

σ

 

 

+

 

σ2

x

 

.

Легко доказать, что из условия (9.19) при определенных выше параметрах γ1, γ2, γ3 имеем неравенство

 

a(b − 1)2

 

a2(b − 1)2

21σ + 2γ2

 

+

 

0.

σ

σ2

Таким образом, при всех x, y, z

ϕ(x, y, z) *4+ (σ − 1)2 .

128

 

 

 

 

ГЛАВА 9

Отсюда следует, что для любого числа

 

 

= *

 

 

2b − σ − 1

s < s

0

4+ (σ − 1)2

 

 

 

*

 

+ σ + 1

 

 

 

4+ (σ − 1)2

существует число ε > 0 такое, что при всех x, y, z выполнена оценка

λ1(x, y, z) + λ2(x, y, z) + 3(x, y, z) + v˙ (x, y, z) < −ε.

Устремив s к s0 справа и используя теорему 9.2, получим оценку

dimL K 3

2(σ + b + 1)

 

.

(9.21)

σ + 1 + *

 

 

 

 

(σ

1)2

+ 4

 

 

 

 

 

 

 

Используя лемму 9.3, легко видеть, что локальная ляпуновская размерность стационарной точки x = y = z = 0 системы (9.17) равна

3

2(σ + b + 1)

 

.

(9.22)

σ + 1 + *

 

 

 

 

(σ

1)2

+ 4

 

 

 

 

 

 

 

Из соотношений (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

лученный частотный критерий дает условия неустойчивости на аттракторах.

Для дискретных систем обычно нет той наглядности представлений в фазовом пространстве, которая характерна для непрерывных систем с гладкими траекториями, заполняющими это пространство. Однако этот недостаток может быть частично ликвидирован, если вместо наблюдения за изменением в дискретном времени расстояния между двумя решениями рассматривать в начальный момент времени отрезок прямой, связывающий начальные данные, и дальнейшие его итерации — последовательность непрерывных кривых. Уменьшение или увеличение длин этих кривых дает информацию об устойчивости или неустойчивости и о других свойствах решений дискретных систем. Рассмотрение длин таких непрерывных кривых может быть положено и в основу определений устойчивости и неустойчивости, что особенно важно для неинъективных отображений, определяющих соответствующие дискретные системы (а именно они являются

внастоящее время широко известными генераторами хаоса). Дело

втом, что в некоторых случаях две траектории могут «слипнуться» за конечное число итераций. Однако итерируемый отрезок, связывающий начальные данные, может иметь все большую длину даже с учетом такого слипания.

Оказалось, что механизм растягивания в ряде случаев является достаточно сильным и его можно обнаружить и оценить, применяя весь богатый арсенал прямого метода Ляпунова, который можно соединить с эффективным алгебраическим аппаратом: лем-

мой Калмана – Сеге и леммой Шепелявого. Последние дают оценки растяжения длин итерируемых кривых в терминах частотных неравенств, оперирующих с передаточными функциями линейных блоков.

Эффект «слипания», который является менее сильным, чем «растягивание» в малых окрестностях рассматриваемых траекторий, усиливается по мере увеличения размеров окрестности (а это, как правило, неизбежно при превалировании эффекта растягивания). В какой-то момент окрестность становится такой большой, что «сли-

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