CC BY
41
Бифуркации ^^^^^^^^^^^^^^^^^
^ш в динамических системах
Изв. вузов «ПНД», т. 18, № 1, 2010 УДК 517.9
РОЖДЕНИЕ И РАЗРУШЕНИЕ МНОГОСЛОЙНОЙ ЗАМКНУТОЙ КРИВОЙ В НЕОБРАТИМЫХ ОТОБРАЖЕНИЯХ
Ж.Т. Жусубалиев, О.О. Яночкина
В статье обсуждается новый феномен, открытый авторами совсем недавно в двумерных эндоморфизмах, демонстрирующих переход к хаосу через разрушение замкнутой инвариантной кривой. Как известно, при рациональном числе вращения на этой кривой имеется четное число периодических орбит, половина из которых устойчивые, а половина - седловые, а сама кривая образована замыканием неустойчивых могогообразий сед-ловых циклов. Как оказалось, если отображение необратимо, то седловая периодическая орбита, лежащая на замнкутой инвариантной кривой, обычно претерпевает бифуркацию удвоения периода или бифуркацию вилки. В результате этого на инвариантной кривой мягко рождаются два новых «слоя», образованные неустойчивыми многообразиями сед-лового цикла удвоенного периода или неустойчивыми многообразиями двух седловых циклов того же периода в случае бифуркации вилки.
Ключевые слова: Двумерные эндоморфизмы, «многослойная» замкнутая инвариантная кривая, квазипериодическая динамика.
Введение
Возникновение хаотических колебаний через разрушение двумерного тора, на котором число вращения рационально, - широко распространенное явление. В.С. Афраймович и Л.П. Шильников доказали теорему о разрушении замкнутой инвариантной кривой в сечении тора с резонансной структурой и указали возможные пути хаотизации колебаний [1,2]. Общий характер выводов этой теоремы был подтвержден численно и экспериментально для широкого класса дискретных и потоковых систем (см., например, [3] и цитируемую там литературу).
В последние годы теория Афраймовича-Шильникова была существенно дополнена некоторыми модификациями, характерными для необратимых отображений (эндоморфизмов) [4,5] и для кусочно-гладких динамических систем [6,7]. Так например, в работах [4,5] показано, что разрушение инвариантной кривой может быть связано с особой нелокальной бифуркацией, приводящей к бесконечному числу самопересечений неустойчивых многообразий седловой периодической орбиты с образованием «петель» на этой кривой [8-10].
В недавней нашей работе [11] при анализе динамики трехмерной системы кусочно-гладких дифференциальных уравнений впервые была обнаружена «двухслойная» замкнутая инвариантная кривая в сечении резонансного тора. Такая кривая
является объединением неустойчивых многообразий двух седловых периодических орбит с устойчивым циклом. В последующих исследованиях [12-14] было установлено, что это явление - типичное свойство двумерных эндоморфизмов, причем число «слоев» замкнутой инвариантной кривой может быть и больше двух. Такую кривую мы назвали «многослойной».
В данной статье мы преследуем цель описать это явление на простейшем примере необратимого отображения.
1. Рождение и разрушение «многослойной» замкнутой кривой
Рассмотрим двумерное необратимое отображение [15,16]
/ : (!) - и У, ^ , (1)
у) \by — cx + x i '
принадлежащее к классу (Z0 — Z2).
Фазовая плоскость отображения (1) делится на две открытые области Zo и Z2 такие, что точка X = (x, y)T, принадлежащая Z2, имеет два разных прообраза f-l(X) ранга один, а точка X = (x,y)T из области Z0 - ни одного. Границей указанных областей является линия LC (ligne critique [15]), которая есть геометрическое место точек X, имеющих два совпадающих прообраза f-l(X), расположенных на так называемой кривой сливающихся прообразов LC-i (curve of merging pre-images [15]). Если отображение непрерывно и непрерывно дифференцируемо, то LC-i включена во множество, определяемое условием
\.1 (X )| =0,
где \,1 (X)\ - определитель матрицы Якоби.
Для рассматриваемого отображения определитель матрицы Якоби обращается в нуль вдоль вертикальной прямой
ЬС-1 = {(х,у) : х = с/2}.
Если отображение недифференцируемо, как в случае кусочно-гладких систем, то роль ЬС-1 обычно играет кривая недифференцируемости, на которой определитель матрицы Якоби меняет знак [15]. Кривую ЬС, удовлетворяющую /(ЬС-\) = = ЬС и /-1(ЬС) ^ ЬС-1, называют критической линией ранга один [15]. Для отображения (1) уравнение критической линии имеет вид
ЬС = /(ЬС-1) = {(х,у) : у = Ьх — с2/4}.
Образ ранга к кривой ЬС называется критической линией ранга к + 1: ЬСк = = /к(ЬС) [15]. Следует отметить, что отображения типа (^о — составляют класс необратимых отображений, обладающих простейшими свойствами. Мы не будем здесь обсуждать вопрос о роли критических линий при бифуркациях инвариантных кривых. Заинтересованных читателей отсылаем к источникам [9,10,15,16].
Отображение (1) имеет две неподвижные точки 0\ = (0,0)т и О2 = = (1 — Ь + с, 1 — Ь + с)т. На линии N = {(Ь,с) : с =1} плоскости параметров (Ь, с) неподвижная точка О1 = (0,0)т теряет устойчивость, когда комплексно-сопряженная пара мультипликаторов выходит на границу единичного круга. При этом из неподвижной точки рождается квазипериодический режим через суперкритическую бифуркацию Неймарка-Саккера.
На рис. 1 приведена карта динамических режимов, рассчитанная в плоскости параметров (Ь, с). Области периодичности, расположенные выше бифуркационной границы ^ф, соответствуют резонансным языкам (языкам Арнольда). Белая часть рисунка, обозначенная через , относится к области квазипериодичности и хаоса. На рис. 1 выделен резонансный язык 1: 5.
В классическом случае внутри каждого клюва на замкнутой инвариантной кривой, отвечающей рациональному значению числа вращения, устойчивый цикл сосуществует с седловым того же периода или сосуществуют несколько пар устойчивых и неустойчивых периодических орбит. Сама кривая является объединением неустойчивых многообразий седловых периодических орбит с устойчивыми циклами.
Оказывается, если отображение необратимое, то седловая периодическая орбита, лежащая на замкнутой инвариантной кривой, может претерпевать бифуркацию удвоения периода или бифуркацию вилки [12-14]. В результате на этой кривой появляются два новых «слоя», образованные неустойчивыми многообразиями седло-вой периодической орбиты удвоенного периода или неустойчивыми многообразиями двух седловых циклов того же периода в случае бифуркации вилки [12-14].
На рис. 2, а приведена бифуркационная диаграмма, иллюстрирующая основные стадии рождения «трехслойной» замкнутой кривой в отображении (1). В точке
1.56 ||
¡1 \\
ц" н
с
\ 7У+ 11
Пос \
•фи
1.06 1
0.36 Ъ 0.64
Рис. 1. Карта динамических режимов в плоскости управляющих параметров (Ь, с)
1.32
-0.80
т /
\ У 1С
V р* /
ьс
0.75
1.25
Рис. 2. а - бифуркационная диаграмма, иллюстрирующая рождение замкнутой инвариантной кривой с «трехслойной» структурой (Ь = 0.54). 1 - область квазипериодической динамики; 2 - область существования «однослойной» замкнутой кривой отвечающей рациональному значению числа вращения; 3 - область существования «трехслойной» инвариантной кривой. б - инвариантная кривая Т с иррациональным числом вращения, отвечающая квазипериодическому режиму (Ь = 0.54, с = 1.25). Здесь Е - неустойчивый фокус
/ / / Л' // 1 /з / / / / / \ / / / / ^ / \ / 1 т / / — \ т, \ Жи \ \ 1.4
Х\ г Л— \ V 1 / /—¿с ' \ я > / Г ( V \ / ¿Л" \ ЧХх ^ хРг ис У
X -0.8
-0.75
X
1.25
БУ / X
У'-..... /
/
// ' 1 У к
\ У'! /
Уг/ч /■+ /
V /ч / ^ ьс
-0.75
X
1.25
Рис. 3. а - «однослойная» замкнутая кривая (Ь = 0.54, с = 1.32); N, Б - соответственно устойчивый и седловой резонансные циклы; - устойчивые и неустойчивые многообразия седловой пери-
одической орбиты Б. б - «трехслойная» инвариантная кривая (Ь = 0.54, с = 1.375); ^о
- устойчивый резонансный цикл; Я - неустойчивый узловой 5-цикл (репеллер); Б* - седловой 10-цикл, возникающий через бифуркацию удвоения периода
1.4
с = еф, лежащей на линии ^, из неподвижной точки ()\ рождается инвариантная кривая Т с иррациональным числом вращения, отвечающая квазпериодическому режиму (рис. 2, б).
На рис. 3, а изображен фазовый портрет отображения (1) после седло-узловой бифуркации, при которой возникает «однослойная» замкнутая кривая. При дальнейшем увеличении параметра с резонансный цикл сначала переходит из устойчивого узла N в устойчивый фокус когда мультипликаторы становятся комплексно-сопряженными. При этом инвариантный тор теряет гладкость вблизи точек устойчивого цикла. Затем в точке с = С2, как это видно на рис. 2, а, мягко возникает неустойчивый 10-цикл через бифуркацию удвоения периода, а седловой 5-цикл становится
неустойчивым узлом Я (репеллером). В результате такого перехода из «однослойной» замкнутой инвариантной кривой мягко рождается «трехслойная» (рис. 3, б). Внешние два слоя являются притягивающими и образованы неустойчивыми многообразиями седло-вого 10-цикла Б *, а внутренний слой -многообразиями неустойчивого узла Я, соответствующими большему по модулю мультипликатору.
При изменении параметров «трехслойная» замкнутая кривая может перейти в «пятислойную» следующим образом (рис. 4). При движении по параметрам к левой границе N+ резонансного языка неустойчивый 5-цикл Я претерпевает бифуркацию вилки, в резуль-
-о.:
я2 / ьс,
\ ш ^^¡^
\ Ш АС!
щ/ р* [ ьс
-0.75
1.25
Рис. 4. «Пятислойная» замкнутая кривая (Ь = 0.4905, с = 1.379). Бо - седловой цикл периода 5; N - устойчивый резонансный цикл; Я\, Я2 - неустойчивые циклы периода 5 (репеллеры), возникающие через бифуркацию вилки неустойчивого узла Я
тате которой он становится седловым Б0, а в окрестности Я мягко рождается пара неустойчивых симметричных цикла Я1 и Я2. После этой бифуркации центральный слой вновь становится притягивающим, но возникают два новых «отталкивающих» слоя (рис. 4).
При дальнейшем изменении параметров наблюдается сложная последовательность бифуркаций, приводящая к разрушению инвариантной кривой и возникновению хаотических колебаний. Мы не будем детально обсуждать механизмы разрушения «многослойной» замкнутой инвариантной кривой, а ограничимся рассмотрением лишь двух сценариев, связанных с бифуркацией седло-узел и гомоклинической бифуркацией.
Рис. 5 и рис. 6 иллюстрируют случай, связанный с седло-узловой бифуркацией. Этот переход наблюдается при выходе из резонансного языка через левую границу В этом случае сначала «пятислойная» замкнутая кривая (см. рис. 4) переходит в «трехслойную», когда неустойчивые узлы Я1, Я2, лежащие на центральном
1.4
///5* - /// / / / г
1/........1......1..
/ / / У
6/ / /■' V + >< ьс
-0.8
-0.75
1.25
-0.72
1.22
Рис. 5. Переход от «пятислойной» к «трехслойной» инвариантной кривой. Средние два слоя (см. рис. 4) исчезают, когда неустойчивые узлы К\ и Я2 становятся неустойчивыми фокусами ,
слое, становятся неустойчивыми фокусами ^, Г2. Эта бифуркация проиллюстрирована на рис. 5. Затем седловой цикл £0, лежащий на центральном слое, «подтягивается» к устойчивому узловому 5-циклу N0, сливается с ним и исчезает на линии N+. В результате этого центральный слой «трехслойной» замкнутой кривой исчезает, а неустойчивые многообразия внешних слоев формируют «хаотическое множество», причем, связанные с ними седловая орбита £ * и неустойчивые фокусы , Г2 продолжают существовать (см. рис. 6).
При изменении параметров в направлении правой границы резонансного языка, «пятислойная» инвариантная кривая сначала переходит в «трехслойную» через бифуркацию вилки, а затем реализуются две гомоклинические бифуркации, в результате которых сначала разрушается внешний слой, а затем внутренний.
На рис. 7 изображен момент гомоклинического касания устойчивого и неустойчивого многообразий седловой периодической орбиты £*. Дальнейшие изменения параметров приводят к рождению гомоклинических точек трансверсального пересечения многообразий (рис. 8, а и рис. 8, б). Как можно видеть из рис. 8 и рис. 9, а, после этой бифуркации внешний слой превращается в довольно сложное множество с фрактальной структурой [4], однако центральный слой продолжает существовать. На рис. 8, в показан отдельно увеличенный фрагмент центрального слоя после го-моклинической бифуркации. Затем седловой 10-цикл £ * претерпевает «обратную» бифуркацию удвоения периода, что приводит к исчезновению внешнего слоя. На рис. 9, а изображен фазовый портрет перед бифуркацией удвоения периода. Рис. 9, б иллюстрирует фазовый портрет после «обратной» бифуркации удвоения периода. На рис. 9, в приведен увеличенный фрагмент «однослойной» замкнутой кривой.
Далее реализуется вторая гомоклиническая бифуркация, приводящая к разрушению «однослойной» инвариантной кривой. После бифуркации седло-узел на границе N+ резонансного языка движение становится хаотическим. На рис. 10, а показан фазовый портрет после гомоклинической бифуркации, а на рис. 10, б - увеличенный фрагмент фазового портрета.
Рис. 7. Разрушение внешнего слоя через гомоклиническую бифуркацию: а - фазовый портрет отображения вблизи гомоклинического касания; б - увеличенный фрагмент фазового портрета, выделенного прямоугольником
1.32
1.20
#Л> и/- \ \ ! 1 ! I
/Ж /// 4 } / у У
1 / х V 1 1 УУу^
-.....IX ■ 0.57
П *__\)|Й»\ 4
¿^¡¡Ш^ / \ N ) ^ \
\
/ \
/ ^
Рис. 8. а - фазовый портрет после разрушения внешнего слоя. б - увеличенный фрагмент фазового портрета, изображенного на рис. 8, а. в - центральный слой после гомоклинической бифуркации
Рис. 9. а - фазовый портрет перед «обратной» бифуркацией удвоения периода. б - фазовый портрет после бифуркации удвоения периода седлового 10-цикла. в - увеличенный фрагмент «однослойной» замкнутой кривой
1.20
0.57
-0.80
-0.75 х 1.25 ^ 0.2 х 1.0
а о
5* 1
ТУ * \
У Мн//
/ \
/ М. ^ ж
0.2 г 1.0
Рис. 10. Разрушение «однослойной» замкнутой кривой через гомоклиническую бифуркацию
Заключение
Представленная статья, как и предыдущие [12-14], посвящена обсуждению бифуркаций замкнутой инвариантной кривой с «многослойной» структурой в двумерных эндоморфизмах. В качестве конкретного примера приведены результаты бифуркационного анализа квадратичного необратимого отображения. В плоскости параметров выделены области с различным динамическим поведением. Показано, что при изменении параметров в пределах резонансного языка, седловая периодическая орбита, лежащая на замкнутой инвариантной кривой, претерпевает бифуркацию удвоения периода, становясь неустойчивым узлом (репеллером). В результате такой бифуркации на этой кривой мягко рождаются два новых «слоя», образованные неустойчивыми многообразиями седлового цикла удвоенного периода. Неустойчивый узловой цикл (репеллер), лежащий на центральном слое, может претерпевать бифуркацию вилки, при которой возникает «пятислойная» замкнутая кривая. При изменении параметров инвариантная кривая разрушается, и при выходе из резонансного языка через линию седло-узловой бифуркации, движение становится хаотическим. Детальное обсуждение механизмов разрушения «многослойной» замкнутой инвариантной кривой требует отдельной статьи. Поэтому, мы ограничились анализом лишь двух сценариев, связанных с гомоклинической бифуркацией и бифуркацией седло-узел.
Вопрос, связан ли описанный феномен в общем случае с необратимостью отображения, остается открытым. Во всяком случае, нам пока не удается доказать это строго и найти примеры диффеоморфизмов, демонстрирующих такое явление.
Библиографический список
1. Афраймович В. С., Шильников Л. П. Инвариантные двумерные торы, их разрушение и стохастичность //Методы качественной теории диффернециальных уравнений. Горький: Изд-во Горьк. ун-та, 1983. С. 3.
2. Арнольд В. И., Афраймович В. С., Ильяшенко Ю. С., Шильников Л. П. Теория бифуркаций //Современные проблемы математики: фундаментальные направления /Под ред. В. И. Арнольда. М.: ВИНИТИ, 1986.
3. Anishchenko V. S., Astakhov V. V., Neman A. B., Vadivasova T. E., Schimansky-Geier L. Nonlinear Dynamics of Chaotic and Stochastic Systems. Tutorial and Modern Development. Berlin: Springer, 2007.
4. Maistrenko V., Maistrenko Yu., Mosekilde E. Torus breakdown in noninvertible maps // Phys. Rev. E. 2003. Vol. 67. P. 046215.
5. Frouzakis C. E., Kevrekidis I. G., Peckhman B. B. A route to computational chaos revisited: noninvertibility and breakup of an invariant circle // Physica D. 2003. Vol. 177. P. 101.
6. Zhusubaliyev Zh. T., Mosekilde E. Bifurcations and Chaos in Piecewise- Smooth Dynamical Systems. Singapore: World Scientific, 2003.
7. Zhusubaliyev Zh. T., Mosekilde E., De S., Banerjee S. Transition from phase-locked dynamics to chaos in a piecewise-linear map // Phys. Rev. E. 2008. Vol. 77. P. 026206.
8. Lorenz E. N. Computational chaos - a prelude to computational instability // Physica D. 1989. Vol. 35. P. 299.
9. Krauskopf B., Osinga H. M., Peckham B. B. Unfolding the cusp-cusp bifurcation of planar endomorphisms // SIAM J. Applied Dynamical Systems. 2007. Vol. 6(2). P. 403.
10. England J. P., Krauskof B., Osinga H. M. Bifurcations of stable sets in noninvertible planar maps // Int. J. Bifurcat. Chaos. 2005. Vol. 15. P. 891.
11. Zhusubaliyev Zh. T., Mosekilde E. Birth of bilayered torus and torus breakdown in a piecewise-smooth dynamical system // Phys. Lett. A. 2006. Vol. 351. P. 167.
12. Zhusubaliyev Zh. T., Mosekilde E. Formation and destruction of multilayered tori in coupled map systems // Chaos. 2008. Vol. 18. P. 037124.
13. Zhusubaliyev Zh. T., Mosekilde E. Multilayered tori in a system of two coupled logistic maps // Phys. Lett. A. 2009. Vol. 373. P. 946.
14. Zhusubaliyev Zh. T., Mosekilde E. Novel routes to chaos through torus breakdown in noninvertible maps // Physica D. 2009. Vol. 283. P. 589.
15. Mira C., Gardini L., Barugols A., CathalaJ. C. Chaotic Dynamics in Two-Dimensional Noninvertible Maps. Singapore: World Scientific, 1996.
16. Frouzakis C. E., Gardini L., Kevrekidis I. G., Millerioux G., Mira C. On some properties of invariant sets of two-dimensional noninvertible maps // Int. J. Bifurcat. Chaos. 1997. Vol. 7. P. 1167.
Курский государственный Поступила в редакцию 15.03.2009
технический университет После доработки 04.06.2009
FORMATION AND BREAKDOWN OF A MULTILAYERED CLOSED CURVE IN NONINVERTIBLE MAPS
Z. T. Zhusubaliyev, O. O. Yanochkina
The paper describes the mechanism for the formation of closed invariant curves that are formed as layered structures of several sets of interlacing manifolds each with their associated stable or unstable resonance modes. Such invariant curves can arise, for instance, if the saddle cycle on a «simple resonance curves» undergoes period-doubling or pitchfork bifurcations transversely to the circumference of the closed curve.
Keywords: Two-dimensional endomorphisms, multilayered closed invariant curve, quasi-periodic dynamics.
Жусубалиев Жаныбай Турсунбаевич - доктор технических наук, профессор. Окончил Томский институт автоматизированных систем управления и радиоэлектроники (ТИАСУР) в 1982 г. После окончания института работал в НИИ автоматики и электромеханики при ТИАСУР. В 1989 г. защитил кандидатскую диссертацию, а в 2002 г. - диссертацию доктора технических наук в области исследования нелинейной динамики и хаоса в импульсных системах автоматического управления. Профессор кафедры вычислительной техники Курского государственного технического университета. Область научных интересов - теория бифуркаций и динамического хаоса; сложная динамика и хаос в кусочно-гладких динамических системах, импульсных системах управления и вибрационной механике. Опубликовал свыше 100 работ в отечественной и зарубежной научной печати, является соавтором 5 монографий и одного учебного пособия, включая монографию «Zhanybai T. Zhusubaliyev and Erik Mosekilde, Bifurcations and Chaos in Piecewise-Smooth Dynamical Systems. World Scientific, 2003». Член Международного общества физики и управления (IPACS). 305040 Курск, ул. 50 лет Октября, 94 Курский государственный технический университет E-mail: zhanybai@mail.kursk.ru, zhanybai@hotmail.com
Яночкина Ольга Олеговна - окончила машиностроительный факультет Курского государственного технического университета (КурскГТУ) в 2001 г. Инженер кафедры вычислительной техники КурскГТУ. Область научных интересов - нелинейная динамика, теория бифуркаций, математическое моделирование импульсных систем автоматического управления. Имеет 11 научных публикаций.
305040 Курск, ул. 50 лет Октября, 94
Курский государственный технический университет
E-mail: yanoolga@gmail.com
