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

книги / Нестационарные структуры и диффузионный хаос

..pdf
Скачиваний:
6
Добавлен:
12.11.2023
Размер:
14.86 Mб
Скачать

3.Отображение Т, действующее на множестве М и

имеющее инвариантную меру Р(х), называется перемешивающим,

если для любых f, g € Z«2

lim

S

f(Tn(x)) g(x) dP

=

f

fdP

S gdP.

 

 

 

(5.7)

n* a M

 

 

 

 

M

M

 

 

 

 

 

Смысл этого определения проясняется, если

вновь

 

рас­

смотреть

преобразование

тора

А (см. формулу

(5.2)).

Ин­

вариантнаямера

в

этом

случае имеет

постояннуюплотность.

Обозначим множество, на котором действует А

через

5,

две

произвольные области, принадлежащие S, через F и G. Обо­

значим меры

этих

множеств соответственно

через

 

mes

S,

mes

F и mes

G. Тогда из соотношения

(5.7)

следует, что

 

 

 

 

 

1 j mmes( ( AnF)

n

G)

_

mes

G

 

 

 

 

 

 

 

 

mes

F

 

 

mes

5

 

 

 

 

 

 

 

л-»oo

 

 

 

 

 

 

 

 

 

 

 

Это

означает,

что произвольное

множество

через

доста­

точно большое число итераций будет равномерно размазано по тору. «Черные» участки (см. рис. 5.1) окажутся перемешан­ ными с белыми. Доказательство соотношения (5.7) для отоб­ ражения (5.2) можно найти в книге [8]. Определение (5.7) естественно обобщается на случай непрерывного времени. Из свойства перемешивания следует эргодичность.

Проверка соотношения (5.7) в реальных системах оказы­ вается достаточно сложной. Поэтому, чтобы установить свой­ ство перемешивания, обычно, рассматриваются автокорреля­ ционные функции b(t)

m = Sf(x(t))f(x(0))dP(x) - Sfm)dP(x)Sf(x(0))dP(x).

В эргодических системах удобно заменить среднее по мере средним по времени для типичной траектории:

b(t) = \\mr'ff(x(t

T))f(x(T))dx -

Г-»а> Q

- ( l i m 7 ’- , J 7 ( * ( T ) ) d T ) 2 .

( 5 . 8 )

Г-» со

151

Если динамическая система обладает перемешиванием, то b(t) —* 0 при t —* оо, что свидетельствует о стохастичности

изучаемого процесса [175].

В физике это свойство известно

как свойство расцепления

временных корреляций [87]. В

частности, отображение (5.6) свойством перемешивания не обладает.

Закон убывания b(t) и спектральные характеристики ди­

намических систем тесно связаны. Представим b(t)

в виде

интеграла

 

b(t) = " ет а(ь>) du>.

(5.9)

-00

 

В случае эргодического движения без перемешивания функция

а(о>)

является

дискретной:

 

 

 

 

 

а(и>)

= £ ak 3(и> -

u>k),^

(5.10)

 

 

 

k

 

 

 

при

наличии

перемешивания -

непрерывной

[87].

 

4. Естественно

ожидать,

что

многие

свойства независи­

мых случайных величин будут характерными и для динамичес­ ких систем с хаотическим поведением. Одним из них является

центральная предельная теорема. Введем величину

т

т

(5.11)

. TU2 [ г 1 X f(x(t)) dt -

s f(x) dP(x)].

оо

Величина в скобках во многих случаях убывает как

7’-,/2 . Поэтому по аналогии с теорией вероятностей можно

ожидать,

что при Т —» со величина (5.11) должна подчиняться

гауссову

распределению

вероятностей. Физически

это

соот-

ветствует

тому,

что

 

Т

f(x(t))dt близок

к

сумме

интеграл X

слабо зависимых

 

 

0

 

 

 

слагаемых.

 

 

 

5. Одним из наиболее сильных стохастических свойств

является

экспоненциальное убывание корреляций:

 

 

b(t) —» ехр(-Л0*)

при

больших

t, hQ = 1/TQ.

 

(5.12)

152

В этом случае функция а(и) (см. формулу (5.9)) аналитична. Величину TQ часто называют временем расцепления корреляций, или радиусом корреляции системы.

При исследовании отображения (5.2) важным являлось

наличие сжимающегося и расширяющегося слоения. В

окрестности каждой траектории лежали бесконечно близкие

траектории,

одни

из

которых стремились

к ней при / —»<»,

другие -

при t

—>

-оо. Оказывается,

если таким образом

устроены все траектории динамической системы, то на основе представлений о локальных свойствах траекторий можно

исследовать глобальные свойства динамической системы. Такой подход развивается в гиперболической теории. Обратим внимание на некоторые ее результаты.

Свойство гиперболичности естественно обобщает свойст­

ва отображения (5.2) на многомерные отображения с более сложным фазовым пространством.

Рассмотрим динамическую систему с дискретным временем

вфазовом пространстве М (в отображении (5.2) это тор).

Обозначим некоторую

траекторию

через

(S*(x)},

а

касатель­

ное

пространство

в

точке

х

через

Т

(рис.

5.3).

Если

у

нас

есть

точка х

М

и бесконечно

X

 

ней

точки

1

близкие к

х ,

то

все

возможные

вектора

х1 -

х

и

составляют касательное

пространство Т .

 

 

 

 

 

 

 

 

 

 

153

Траекторию

{S*(x)}-

следуя

работе [158],

назовем

равномерно полно гиперболической, если касательное

про­

странство Т

может быть

разложено

на прямую сумму

под-

s\*)

 

 

 

 

 

пространств

 

 

 

 

 

Т= Es(S*(x)) + £ “(S'(x)).

S?(x)

причем эти подпространства являются инвариантными от­ носительно отображения Т (dSTEs(S?(x)) = Es(St+T(x)), dSTEu(S\x)) = Eu(St+T(x))) и существуют постоянные числа

О< Л < 1 < р

такие, что для всех t и т > О

\dST о|| s С Лт И ,

о € № '( * ) ) ;

|rf5T оI а С 'V х |о|. о € £“(S'(x));

 

 

 

 

 

 

 

3r(S*(x))

> const.

 

 

 

 

 

 

 

 

где

y(.S*(x))

-

угол

между

подпространствами

Es(S{(x))

и

£ “(S*(x))

(рис.

5.4).

 

Через

 

dSTEs(St(x))

обозначено

 

мно­

жество

векторов,

в

которое

через

время т

перейдут

векторы,

в момент

t

принадлежавшие

подпространству

Es{S\x)).

Под

о||

понимается

расстояние

между

траекторией,

начинаю­

||dS

щейся

в

точке

S°(x),

и

бесконечно

близкой

 

траекторией,

выходящей из

 

точки

S°(x)

+

v, в момент т. Чтобы найти

dSTv

 

в системах с

непрерывным

временем

х

=

Х(х),

решают

уравнение

в

вариациях v

=

Хх(х)

v на интервале 0

<

t

<

т.

Всистемах с дискретным временем вместе с хп+ = f(xn)

рассматривают

отображение »n+1 = fx(xn)v

(X = — ,

f =

df

х

Эх

х

= — , {, X, х -

векторы).

 

 

дх

 

 

 

Динамическая система называется У-системой, или диф­ феоморфизмом Аносова, если каждая ее траектория является равномерно полно гиперболической, а постоянные С, Л и д = = Л'1 можно выбрать одинаковыми для всех точек.

154

Диффеоморфизмы Аносова, так же как отображение (5.2), являются структурно устойчивыми. Число задач, где они

возникают

в чистом виде,

невелико.

К ним,

‘ например,

относится

задача о движении

по геодезическим на

поверхнос­

тях постоянной отрицательной кривизны [158].

 

Гиперболические траектории характерны для интересного

класса динамических систем,

называемых

биллиардами. Билли­

ардом называется система, отвечающая движению материальной

точки

по инерции внутри области с кусочно гладкой грани­

цей.

Предполагается выполненным условие упругого отражения

от границы. Если область ограничена гладкими выпуклыми внутрь кривыми, и кривизна положительна для каждой регу­ лярной точки границы, биллиард называется рассеивающим. Он представляет собой аналог гладких гиперболических систем. Биллиарды возникают в задачах акустики, радиофизики, кван­ товой механики. Хаотическое поведение в них во многих слу­ чаях удается детально исследовать [43].

_ Однако наибольший интерес при исследовании дис­ сипативных систем представляет собой использование резуль­ татов гиперболической теории для анализа странных аттрак­ торов.

Аттрактор динамической системы называется гиперболи­ ческим, если он состоит целиком из траекторий, удовлет­ воряющих условию равномерной полной гиперболичности с одними и теми же постоянными С и Л.

155

 

Гиперболические системы обладают хорошими стохасти­

ческими

свойствами.

В

частности,

отображение

Т является

перемешивающим

по

отношению к

распределению

вероятностей

Р,

если

аттрактор

удовлетворяет

условиям

гиперболичности

(при

некоторых

дополнительных

ограничениях

[224, 359]).

При этом для любого начального распределения PQ с плот­

ностью р0, сосредоточенного в окрестностиаттрактора W,

итерации

TnPQ сходятся

к Р. Однако ясно, что

с

помощью оп­

ределения гиперболичности можно исследовать свойства дос­

таточно узкого

класса динамических систем.

В работе

[175]

были приведены

условия гиперболичности,

которые в

ряде

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

 

Следуя работе [175], приведем

этот критерий

для сис­

темы

с

дискретным временем. Наряду с отображением Т в

точке

х

будем рассматривать его

матрицу Якоби

dT(x) =

=df/dx (отображение Т переводит УУ-мерный вектор х в

jV-мерный вектор /(х): xn+1

=

f(xn)). Вектор

е назовем рас­

тягивающимся, если |df(x)e|

>

||ej. Набор

всех растягиваю­

щихся векторов в данной точке х образует конус растягиваю­ щихся векторов. Изменив знак неравенства на противополож­ ный, получим определение сжимающегося вектора и конуса сжимающихся векторов. Для N = 2 этих определений было бы достаточно. В многомерном случае приходится определять еще расширяющееся L(u) и сжимающееся подпространства или

подпространства, состоящие соответственно из расширяющихся или сжимающихся векторов.

Допустим

теперь, что

для пары

целых чисел k, I, k+l =

= N и постоянной Л (1 < А < ») для

каждой

точки х некото­

рой области Q,

содержащей

аттрактор

(TQ с

Q), задано отк­

рытое множество &и\х) 6-мерных расширяющихся подпрост­ ранств lSu^ и открытое множество &s)(x) сжимающихся

/-мерных подпространств lSs\х) (рис. 5.5) и при этом

156

1)

dT(x)&u) (x)

с

&u)(Tx),

т. е. ДЛЯ любого lSa^€

&и)(х) Образ dTlSu)= L{u)

&и)(Тху,

более

того,

для

любого

вектора е

L ^ €

&и'(х)

справедливо

неравенство

 

 

 

 

№ ) е |

*

ЛЦе|;

 

 

 

 

 

2)

(dT(Tx))~*&s)(Tx)

с

&*\х),

т.

е.

для . любого

L{s)

&*>(Тх)

образ

 

 

=

L\s)

^ > (х );

более

того, для любого вектора е €

 

&s)(Tx)

справедливо

неравенство

 

 

 

 

 

 

 

 

 

 

 

 

 

IKdTrVnOefl

s Л«е|| .

 

 

 

 

Эти условия можно пояснить с помощью рис. 5.5 (N = 2).

 

Первое условие означает,

что

оператор

Т ‘ переводит

растягивающиеся векторы внутрь конуса, в котором лежат растягивающиеся векторы на следующей итерации. Второе условие получается, если заменить конус растягивающихся векторов на конус сжимающихся, а Г на Г-1..

В случае двумерных отображений (5.1) эти условия упрощаются и могут в явном виде быть выражены через произ­ водные функций / и g [44]

где ||h{x,y)\ = max \h(x,y)\. ( * у №

В таком виде критерий гиперболичности использовался, например, в работе [44] при исследовании модели Лоренца и

157

в работе [26] при изучении странного аттрактора в системе

дифференциальных уравнений, связанной с уравнением типа реакция - диффузия.

 

<?Си)(тх)

l f (u)(X)

dT<X)£(a)(x l

( s ) (r x )

d T - 1(Tx ) ( £ ( s ) (Tx } )

A— В

Рис.

5.6

 

Инвариантные множества

в

гиперболических системах

могут обладать сложной структурой

и необычными свойствами.

158

Примером может служить инвариантное множество отображения,

называемого

подковой Смейла.

Т участка

 

 

 

 

Представим

отображение

плоскости

в

себя,

при котором квадрат 5 с вершинами ABCD переходит в криво­

линейную фигуру

А‘ В' С' D' ,

изображенную на

рис.

5.6.

Для

наглядности

представим его

действие в

два

этапа.

Вначале

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

лении.

На

втором

этапе

полученная

полоска складывается в

виде

подковы. Понятно,

какой

будет

структура

множества

л

 

 

каждой

полоски возникают две меньшей ши­

Т (S) п S. Вместо

рины.

Они

заштрихованы

на рис.

5.7.

На

следующей

итерации

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

возникает 2" полос. Отсюда следует, что предельное мно­

жество 1i mTn(S) п S будет обладать структурой, повторяющей

п-»ео

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

ристики далее мы рассмотрим более подробно.

 

 

Для итераций разных точек квадрата можно предложить

простое символическое

описание. Обозначим две полоски, по­

являющиеся

на первом

шаге, V^ и У2На втором шаге

из пер­

вой

полосы

возникнут

полосы

и V12, из второй -

V21 и

V22

и т.

д.

Каждой

точке, остающейся в. пределах квадрата

5,

можно

сопоставить

бесконечную

последовательность

единиц

и двоек следующим способом. Если образ точки на первом ша­ ге попадает в полосу VJt то на первой позиции в этой по­

следовательности

будет

стоять

единица, если в

полосу V2,

то

двойка.

Пусть

на

&-ом шаге точка оказалась в

полосе

V .

,

. ,

......

Образ-

этой

полосы состоит

из

двух

частей

1

1

*2

и

lk

 

Если (£

+1)-я

итерация

точки

попадет в

-

левой

правой.

левую часть, то на (А>+1) месте в последовательности должна стоять единица, если в правую, то двойка. При этом каждая итерация отображения Т эквивалентна сдвигу на один элемент

вправо. В предельном множестве 1i mTn(S) n 5 все точки с п-><»

159

одной координатой у характеризуются одной и той же после­ довательностью единиц и двоек. Разные последовательности

соответствуют точкам

квадрата

с

различными

значениями

у

(см. рис. 5.7). В самом деле, в полосу V^, например, могут

попасть

только

точки

с

ординатой 0 - у

<

1/2,

в

полосу

V2 — с

ординатой

1/2

< у

-

1.

Поскольку

отображение Т

на

каждом

шаге растягивает прообраз вдоль оси

у,

то,

как

бы

ни

было

мало начальное расстояние по оси

у

между

точками

у ' -

у",

на некотором шаге

эти

точки попадут

в разные поло­

сы.

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 5.8

 

Описанное символическое представление итераций отоб­

ражения Т оказывается в точности таким же,

как

для случай­

ного процесса бросания монеты. Рассмотрим

бесконечную

се­

рию бросаний. Если в п

раз выпал орел, то

на л-м месте в

последовательности будем

ставить единицу,

если

решка,

то

двойку. Одно бросание эквивалентно сдвигу на один элемент

последовательности.

Таким образом, каждой «/-координате

точки в предельном

множестве lim Tn(S) n S однозначно c o -

 

n->00

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

ными объектами теории вероятностей. Такой подход, связан­

ный с сопоставлением различным траекториям символического описания, получивший название символической динамики, ока­ зался очень плодотворным. Методы символической динамики помогли исследовать хаотические режимы в ряде интересных динамических систем [152, 175, 289]. Чтобы построить инва­ риантное множество, поступим следующим образом.

160

Соседние файлы в папке книги