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

Динамический хаос (ИПИС, ФКС)

.pdf
Скачиваний:
64
Добавлен:
13.02.2015
Размер:
5.81 Mб
Скачать

 

 

kT

 

b exp

m

,

(162)

 

 

 

 

 

1

(1

b)P(x)

 

(x) (1 b)x k

 

 

 

 

.

(163)

Предположим,

что

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

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

импульсной силы таково, что

 

 

(x) 1 ax2 ,

(164)

где a – некоторая постоянная. Выбор функции (x) в виде (164) соответствует так называемой модели Эно.

Тогда преобразование (161) запишется в виде

 

2

 

 

xn 1

1 axn

b yn .

(165)

 

yn 1 xn

 

 

Это преобразование называется отображением Эно (Henon map). Оно впервые было рассмотрено французским астрофизиком Мишелем Эно в 1976 году. Как видим, отображение Эно относится к отображениям с дискретным временем. Помимо нашей частицы, отображением Эно описывается также динамика и других физических систем – например, диссипативный осциллятор и ротатор под импульсным периодическим воздействием (Heagy, 1992).

| b |.

Вычислим якобиан отображения Эно:

 

xn 1

xn 1

 

 

 

 

J

xn

yn

 

2axn

b

b

yn 1

yn 1

1

0

 

 

 

 

 

 

. (166)

 

xn

yn

 

 

 

 

Фазовое

пространство (x, y) рассматриваемой системы

двумерно, значит, его мерой является площадь. Рассмотрим область фазового пространства, обладающую изначально площадью S . Тогда равенство (166) означает, что эта область фазового пространства будет эволюционировать с течением времени так, что после каждой итерации её площадь будет умножаться на фактор

Так после первой итерации площадь станет равной S1 | b | S , после второй – S2 | b |2 S , после n-й – Sn | b |n S . В общем случае можно рассмотреть следующие варианты:

1)если | b | 1 , то с каждой итерацией, занимаемая системой площадь будет уменьшаться, стремясь к нулю. В этом случае отображение Эно является диссипативной системой, у которой должен присутствовать аттрактор;

2)если | b | 1 , то площадь, занимаемая системой, будет сохраняться. В этом случае отображение Эно является консервативной системой.

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

периодических толчков,

параметр b , определяемый формулой

(162),

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

неравенству: 0 b 1.

Значит,

отображение Эно, описывающее динамику нашей частицы,

является диссипативной системой. Причём при

b 0 (на-

пример, когда T ) отображение Эно (165) переходит в

логистическое отображение.

 

Эно в своей работе провёл исследование динамики системы (165) при различных значениях параметров a и b . Он обнаружил, что в зависимости от выбора этих параметров система может двигаться либо регулярно, либо хаотически. В области хаотической динамики Эно обнаружил также появление странного аттрактора. Так, на рис. 76 показано, как

выглядит

странный аттрактор при значениях параметров

a 1,4

и b 0,3 (такие значения параметров были рас-

смотрены в исходной работе Эно).

Рис. 76. Странный аттрактор отображения Эно (165) при значениях параметров a 1,4 и b 0,3 .

Как было сказано ранее, характер движения системы зависит от выбора параметров a и b . Изменяя параметры a и b , мы можем влиять на режим движения системы. Более наглядное представление о движении системы даёт карта динамических режимов – это диаграмма на плоскости, где по

осям отложены значения параметров (в нашем случае – это

a и b ), а каждая точка плоскости

(a,b)

закрашена

определённым цветом в зависимости от того, какой тип динамики реализуется в системе при данных значениях a и b.

Карту динамических режимов можно построить на компьютере следующим образом:

1)В каждой точке плоскости параметров (a,b) численно решается дифференциальное уравнение или итерируется отображение, задающее динамику исследуемой системы.

2)Полученные решения анализируются и определяется характер динамики системы при данных значениях a и b , т.е. определяется то, каким является движение – регулярным или хаотическим, периодическим или непериодическим, если периодическим, то с каким периодом и т.п..

3)Точке с координатами (a,b) присваивается цвет, зависящий от типа динамического режима, установленного в предыдущем пункте.

4)Каждая точка (a,b) отождествляется с пикселем графического изображения.

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

Из-за того, что нелинейным системам часто присуща мультистабильность (наличие нескольких аттракторов), карту динамических режимов таких систем, вообще говоря, надо представлять не как один лист, а как совокупность листов, перекрывающихся в тех областях параметров, где система имеет более одного аттрактора. Другими словами, в общем случае, каждой точке (a,b) может соответствовать не один, а несколько различных режимов движения, т.е. каждому пикселю карты динамических режимов может соответствовать не один, а несколько цветов.

Так, при определённых параметрах a и b может оказаться, что при выборе одних начальных условий система будет эволюционировать к одному аттрактору, а при выборе других начальных условий – к другому аттрактору. Следовательно, при данных a и b одновременно возможны сразу несколько режимов движения. Именно в таких случаях карту динамических режимов следует считать многолистной.

Рассмотрим карту динамических режимов отображения Эно. Она представлена на рис. 77.