Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Учебники / НЕЛИНЕЙНАЯ ДИНАМИКА СЛОЖНЫХ.pdf
Скачиваний:
612
Добавлен:
30.03.2022
Размер:
28.89 Mб
Скачать

 

http://profbeckman.narod.ru/

20.2 Отображение xn+1=а-xn2.

 

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

 

xn+1=a-xn2

(4)

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

Унимодальное отображение - гладкое отображение на себя при условии, что оно имеет единственную критическую точку (ноль производной).

Все унимодальные отображения демонстрируют похожее динамическое поведение. Основных типов динамических режимов два: 1) стационарные или периодические режимы, 2) хаотические режимы. В первом случае итерации в конечном итоге посещают только конечный набор разных значений, всегда повторяющиеся в фиксированном порядке. В случае (2) состояние системы никогда не повторяется точно и, по-видимому, эволюционирует неупорядоченным образом.

Бифуркационная диаграмма ограничена диапазоном a0<a<aR. Между a0=-1/4 и a1 =3/4 предельный набор состоит из одного значения. Это соответствует стационарному режиму, который есть периодическая орбита периода-1. При a=a1 происходит бифуркация, рождающая периодическую орбиту периода-2: итерации колеблются между двумя значениями – бифуркация удвоения периода. При a=a2=5/4 существует еще одна бифуркация удвоения периода, где орбита периода-2 уступает место орбите периода-4.

Бифуркации удвоения периода, возникающие при a=a1 и a=a2, являются первыми двумя членами бесконечной серии – каскада удвоения периода, в котором для каждого целого n создается орбита периода 2n. Бифуркация при a=a3, приводящая к орбите периода-8, легко видна на бифуркационной диаграмме на рис. 13, та, что при a=a4 едва заметна, а следующие совершенно неразличимы для невооруженного глаза. Это связано с тем, что значения параметра a, при которых создается орбита периода 2n, геометрически сходятся к точке накопления a=1.401155189... с коэффициентом сходимости,

существенно превышающим 1: lim an an 1 4,6692016091...

n an 1 an

Постоянная δ, обладает замечательным свойством: периоды удвоения каскадов, наблюдаемые в чрезвычайно большом классе систем, имеют скорость сходимости, определяемую δ. В точке накопления aпериод решения становится бесконечным. На рис. 13 справа от этой точки система переходит в хаотические режимы, что видно по обилию тёмных областей в этой части бифуркационной диаграммы, которые указывают на то, что система посещает много разных состояний. Каскад удвоения периода - один из наиболее известных путей хаоса и может наблюдаться во многих низкоразмерных системах.

Рис. 13: Бифуркационная диаграмма отображения (4): рассчитана для нескольких значений параметров от a=- 0.25 до a=2.0; проведено 50

последовательных итераций отображения с целью вымираня переходных процессов. Вертикальные линии отмечают создание (i) орбиты периода-2; (ii) орбиты периода-4; (iii) орбиты периода-8 и (iv) точка накопления каскада удвоения периода; (v) начальная точка окна периода-3.

http://profbeckman.narod.ru/

Рис. 14: Увеличенный вид хаотической зоны бифуркационной диаграммы рис. 13.

Внутри периодических окон периода до 8 вертикальные линии указывают значения параметров, где соответствующие орбиты наиболее стабильны (период указан над линией).

В хаотической области есть относительно большое периодическое окно, соответствующее области устойчивости орбиты периода-3, которое начинается с а=7/4, внутри хаотической зоны. Периодические окна и хаотические области произвольно часто чередуются; между любыми двумя периодическими окнами бесконечно много периодических окон. Периодические окна видны невооруженным глазом только у очень низких периодов. У более высоких периодов (1) периодическое окно слишком узкое по сравнению с масштабом графика, и (2) количество выборок достаточно велико, чтобы окно не могло отличаться от хаотических режимов.

Конечно, было бы полезно определить для каждого периодического решения диапазон значений параметров, по которым он устойчив. Увы! Эта задача недостижима, но всё же простые топологические методы позволяют классифицировать разные периодические режимы и ответить на такие вопросы, как продолжается ли последовательность различных динамических режимов, возникающих при увеличении параметра a, возможно диагностика сложности хаоса и др.

Поиск неподвижных точек просто сводится к решению квадратичного уравнения x=a-x2 (5)

который имеет два решения:

( ) = ; ( ) =

(6)

Эти выражения 2.4 имеют смысл только при a>a0=-1/4. Ниже этого значения все орбиты уходят в бесконечность. Бесконечно удаленную точку x, можно формально рассматривать как другую неподвижную точку

системы, хотя и нефизическую.

Рис. 14: Бассейн притяжения неподвижной точки x+ расположен между левой неподвижной точкой x- и ее прообразом, обозначенной двумя вертикальными линиями.

Орбиты, помеченные 1 и 2, находятся внутри бассейна и сходятся к х+. Орбиты, помеченные 3 и 4, находятся за пределами бассейна и уходят в бесконечность (т. е. сходятся к бесконечной точке x).

Важным качественным изменением, происходящим при a=a0, является повсеместное явление низкоразмерной нелинейной динамики,

http://profbeckman.narod.ru/

касательной или седлоузловой бифуркации: две неподвижные точки (6) становятся одновременно вещественными и вырождаются: x-(a0)=x+(a0)=-1/2. Два обозначения указывают на два разных (но связанных) свойства этой бифуркации. Эти две бифуркационные неподвижные точки имеют разные свойства устойчивости. Если значения а несколько выше a0, то орбиты, расположенные вблизи x+, сходятся к нему, тогда как те, которые начинаются в окрестности x, оставляют его либо сходящимся к x+, либо уходящим на бесконечность, в зависимости от того, расположены они справа или слева от х-. Неподвижная точка x+ (а также x) называется устойчивой, а x- неустойчивой. Они называются узлом и седлом соответственно. Так как траектории в окрестностях этих точек сходятся к ним, x+ и xявляются притягивающими множествами или аттракторами. Множества точек, орбиты которых сходятся к аттрактору системы, называются бассейном притяжения этой точки. Из рис. 14 видно, что неустойчивая неподвижная точка x- находится на границе между бассейнами притяжения двух устойчивых неподвижных точек x+ и x. Другая граничная точка – прообраз f-1(x-) значения x-. Устойчивость неподвижной точки зависит от производной отображения в неподвижной точке.

Как известно, неподвижная точка устойчива (соответственно, неустойчива), когда абсолютное значение ее мультипликатора меньше (соответственно больше), чем единица. Мультиплкаторы μ ± двух неподвижных точек отображения (4) задаются формулами

=

=

( )

( )

= −2

= 1 + √1 + 4

(7б)

= −2

= 1 − √1 + 4

(7а)

 

Ур.7 показывает, что x- безусловно неустойчиво по всей своей области существования и, следовательно, в целом не наблюдается как стационарный режим, тогда как x+ стабилен для параметров a чуть выше a0=-1/4. Поэтому на бифуркационной диаграмме, показанной на рис. 13, можно наблюдать только х+. Точнее, x+ устойчиво для a [a0, a1], где a1=3/4 таково, что μ+=-1. Это согласуется с бифуркационной диаграммой, изображенной на рис. 13. Заметим, что при a0= [a0, a1] множитель μ+=0 и, следовательно, возмущения затухают быстрее, чем по экспоненциальной зависимости: тогда неподвижная точка называется суперстабильной. При бифуркации седло-узел обе неподвижные точки вырождены, а их множитель равен +1. Это фундаментальное свойство связано с тем, что в точке бифуркации график отображения (4) касается диагонали (рис.

15), поэтому эта бифуркация также известна как касательная бифуркация. Касание двух гладких кривых (графика f и диагонали) является общим в точке множественного пересечения. Это пример структурно неустойчивой ситуации: произвольно малое возмущение f приводит к двум различным пересечениям или вообще не пересекает (альтернативно, двум действительным или двум комплексным корням).

Рис. 15: График отображения (2) на начальной стадии бифуркаций седло-узел.

Периода-1 орбита x+ существует для любого a>a0. Следовательно, она не исчезает при бифуркации удвоения периода, а просто становится неустойчивой. Она присутствует во всех динамических режимах, наблюдаемых после потери устойчивости, в том числе в хаотических режимах правой части бифуркационной диаграммы рис. 13. Это справедливо для всех периодических решений квадратичного отображения. Поэтому периодические орбиты играют важную роль в динамике даже после того, как они станут неустойчивыми.

Так как орбиту периода-2 можно рассматривать как неподвижную точку второй итерации отображения, то можно определить её диапазон стабильности. Две периодические точки {x1, x2} являются решениями квартирного уравнения

http://profbeckman.narod.ru/

 

 

=

( ) = −1( −x2

)

(8)

x+

x-

 

 

 

 

 

 

 

 

 

Для нахождения x

и

воспользуемся тем,

что фиксированные точки

 

и

являются решениями Ур.8. Решение квадратичного уравнения

 

 

( )

=

( )

= 1 − − + = 0

(9)

 

 

( )

 

 

 

есть

 

; =

 

 

 

 

 

=

 

 

(10)

 

 

Орбита периода 2 (x1, x2) появляется при a=a1=3/4 и существует для каждого a>a1. Используя правило цепочки для производных, получим множитель фиксированной точки x1 функции f2 как

,

=

( )

=

(

)

( )

= 4

= 4(1 − )

(11)

 

 

1

 

 

 

Заметим, что x

 

и x2, рассматриваемые как неподвижные точки f2, имеют один и тот

же множитель, который определяется как множитель орбиты (x1, x2). В точке бифуркации a=a1 имеем μ1,2=1, сигнатура двух периодических точек x1 и x2 вырождается при бифуркации удвоения периода. Однако структура этой бифуркации не полностью аналогична структуре касательной бифуркации, рассмотренной ранее. Действительно, две периодические точки x1 и x2 также вырождены с неподвижной точкой x+. Таким образом, бифуркация удвоения периода фиксированной точки x+ является ситуацией, когда вторая итерация f2 имеет три вырожденные неподвижные точки. Сигнатура этого тройного вырождения соответствует особой форме более высокого порядка чем сингулярность складчатости, упомянутой при обсуждении касательной бифуркации. Заметим, что x+ имеет мультипликатор -1 как фиксированную точку f при бифуркации и, следовательно, существует по обе стороны от бифуркации: он просто становится неустойчивым при a=a1. Напротив, x1 и x2 имеют мультипликатор 1 для первых итераций f, при которых они являются неподвижными точками, и, следовательно, существуют только с одной стороны бифуркации. Орбита периода-2 устойчива только в конечном диапазоне параметров. Другой конец области устойчивости имеет вид a=a2=5/4, где μ1,2=-1. При этом значении параметра происходит новая бифуркация удвоения периода, где орбита периода-2 теряет устойчивость и рождает орбиту периода-4. Как показано на рис. 13, удвоение периода происходит неоднократно, пока не будет создана орбита бесконечного периода.

20.3 Подобие окон периодической динамики

Общая особенность для хаотических динамических систем состоит в том, что при изменения параметров рождается устойчивая и неустойчивая орбиты с периодом n (касательная бифуркация). Далее происходит переход к хаосу через каскад удвоения периода. Все заканчивается кризисом, при котором неустойчивая орбита с периодом n рожденная в исходной касательной бифуркации сталкивается с хаотическим аттрактором. Область параметров между касательной бифуркацией и кризисом называется окном периодической динамики с периодом n. Заметим, что в этом случае центральная часть рисунка подобна всей бифуркационной диаграмме (рис. 16).

Так, при С=-1.75 в окне с периодом 3 в результате тангенциальной бифуркации отображения f o3 возникают устойчивая и неустойчивая орбиты с периодом 3 (в верхней части окна). Устойчивая орбита с периодом 3 показана ниже налево. Если задать N=3 (рис. 17б), то мы получим 8 пересечений (неподвижных точек отображения f o3). Два из них – неустойчивые неподвижные точки. Остальные шесть пересечений принадлежат устойчивой и неустойчивой орбитам с периодом 3.

На рис. 17а периодическая орбита проходит через два "линейных" (по краям) и один квадратичный участок параболы (в центре). Поэтому в окрестности x=0 отображение f o3 с точностью до растяжений совпадет с квадратичным (это преобразование называется