Нелинейная динамика. 2012. Т. 8. № 1. С. 3-28. Полнотекстовая версия в свободном доступе http://nd.ics.org.ru
УДК: 517.9
М8С 2010: 37С70, 37С29, 37С25
К вопросу о сценариях возникновения хаоса у трехмерных отображений
А. С. Гонченко, С. В. Гонченко,
В настоящей работе изучаются вопросы хаотической динамики трехмерных гладких отображений (диффеоморфизмов). Показывается, что здесь существует два основных сценария развития хаоса от устойчивой неподвижной точки либо к спиральному аттрактору, либо к странному аттрактору лоренцевского или «восьмерочного» типа. Дается качественное описание этих аттракторов, приводятся условия, при которых они могут быть «настоящими» (псевдогиперболическими странными аттракторами). В работу также включены соответствующие результаты численного исследования аттракторов в трехмерных отображениях Эно.
Ключевые слова: странный аттрактор, хаотическая динамика, спиральный аттрактор, тор-хаос, гомоклиническая траектория, инвариантная кривая, трехмерное отображение Эно
1. Введение
Одной из основных проблем математической теории динамического хаоса и нелинейной динамики является задача описания сценариев перехода к странным аттракторам (СА) от простых притягивающих режимов — устойчивых состояний равновесия и периодических траекторий. В случае двумерных отображений и трехмерных потоков в этом направлении,
Получено 6 октября 2011 года После доработки 24 декабря 2011 года
Работа выполнена в рамках гранта Правительства Российской Федерации для государственной поддержки научных исследований, проводимых под руководством ведущих ученых в российских образовательных учреждениях высшего профессионального образования, договор №11.034.31.0039. Работа поддержана грантами РФФИ N0.10-01-00429, N0.11-01-00001 и РФФИ N0.11-01-97017-р-поволжье.
Гонченко Александр Сергеевич agonchenko@mail.ги Гонченко Сергей Владимирович [email protected] Шильников Леонид Павлович lpshilnikov@mail.ги
НИИ Прикладной математики и кибернетики 603005, Россия, Нижний Новгород, ул. Ульянова, 10
Л. П. Шильников
как хорошо известно, получено большое число весьма интересных и фундаментальных результатов. Здесь достаточно отметить такие из них, как создание математической теории спирального хаоса [1-3], аттракторов Лоренца [4-7], а также теории СА (квазиаттракторов) двумерных диффеоморфизмов [8, 9]. Естественно, все эти результаты могут быть использованы и при исследовании хаотической динамики многомерных систем (размерности ^ 3 для отображений и ^ 4 для потоков). Однако, как недавно выяснилось, такие системы могут обладать СА новых типов, так называемыми дикими гиперболическими аттракторами [10, 11]. Главной особенностью этих аттракторов является то, что они допускают гомокли-нические касания, но не содержат устойчивых периодических траекторий, которые не появляются также и при возмущениях.1 Соответственно, дикие гиперболические аттракторы нужно относить к «настоящим» СА, к которым, как известно, до недавнего времени можно было приписывать только лишь гиперболические и квазигиперболические (аттракторы Лоренца) СА.
В связи с этим возникает естественный интерес к проблемам хаотической динамики многомерных систем, связанный, в частности, с нахождением сценариев возникновения СА, в том числе и нового типа — диких гиперболических. В настоящей работе мы имеем дело как раз с этой задачей и предпринимаем попытку анализа хаотической динамики трехмерных гладких отображений с позиций качественной теории. Мы показываем, что, независимо от природы конкретных математических моделей, если только они допускают описание с помощью трехмерных диссипативных (с якобианом меньше единицы) отображений, существуют два основных сценария развития хаоса, приводящих либо к аттрактору Шильникова, который может быть спиральным или седловым, либо к СА лоренцевского типа или «восьмерочного» типа (см. ниже). Заметим, что сценарий со спиральным аттрактором был рассмотрен еще в работе [3] (см. также обзор [7]), а здесь мы предлагаем его модификацию, приспособленную именно к трехмерным отображениям.
Отметим, что указанный выше спиральный аттрактор (содержащий неподвижную точку типа седло-фокус) является квазиаттрактором [12, 15], то есть либо он сам содержит устойчивые периодические траектории, либо последние неизбежно появляются при возмущениях. Квазиаттракторы являются по сути «физическими СА», поскольку часто бывает, что устойчивые траектории внутри них имеют весьма большой период и малые области притяжения и, соответственно, не наблюдаются в эксперименте в силу неизбежного наличия шумов и/или ошибок счета (см. дискуссию на эту тему в [7, 16, 17]). Однако указанные СА лоренцевского или «восьмерочного» типа вполне могут претендовать на роль «настоящих» аттракторов (диких гиперболических). Косвенно это подтверждается численными исследованиями, проведенными в работе [18], в которой изучались, в частности, СА лоренцевского типа в трехмерных квадратичных отображениях Эно.
В настоящей работе будут представлены также и результаты численных экспериментов, иллюстрирующие соответствующие теоретические положения.2 В качестве конкретных
1 Системы с дикими гиперболическими аттракторами принадлежат областям Ньюхауса, то есть открытым (в С2-топологии) областям в пространстве динамических систем, в которых плотны системы с гомоклиническими касаниями. Однако эти касания таковы, что их бифуркации (в отличие от двумерного случая) не приводят к рождению устойчивых периодических траекторий [12, 13]. Сам термин «дикий» восходит к Ньюхаусу [14], который ввел понятие «дикого гиперболического множества» — равномерно гиперболического (базисного) множества, у которого касания устойчивых и неустойчивых многообразий сохраняются при малых С2-возмущениях.
2При этом многие компьютерные иллюстрации в статье были получены с помощью программного комплекса '^п8е1 [19, 20].
моделей мы рассматриваем как раз трехмерные отображения Эно — сокращенно 3БНМ (аббревиатура от «3Б Нёпоп Мар»). Напомним, что, по определению, это такие трехмерные квадратичные отображения, которые обладают следующими тремя свойствами: (1) они имеют постоянный якобиан, (и) отображения, обратные к ним, также является квадратичными, (ш) координаты «не отщепляются».3 Как известно [21, 22], любое 3БНМ с помощью аффинных замен координат может быть приведено к такому виду (нормальной форме)
х = у, у = г,
(1.1)
г = Ы\ + Ы^у + Ы%г + Вх + ау2 + Ьуг + ег2,
где всегда можно потребовать, что выполняется какое-нибудь условие нормировки (например, можно сделать равным ±1 один из коэффициентов квадратичной формы), которое достигается соответствующим масштабированием (рескейлингом) координат; также, если |а| + |е| = 0, то (сдвигом координат) всегда можно добиться, чтобы один из коэффициентов, М^ или Ы|, обращался тождественно в нуль. Заметим, что коэффициент В — это якобиан отображения (1.1). Таким образом, нормальная форма для трехмерных отображений Эно содержит не более пяти независимых параметров.
Выбор в качестве эталонных именно трехмерных отображений Эно объясняется многими причинами. Во-первых, они являются в определенном смысле простейшими по форме нелинейными трехмерными отображениями. Во-вторых, несмотря на эту простую форму, они могут демонстрировать весьма нетривиальную динамику, связанную, в частности, с существованием странных аттракторов разнообразных типов. В-третьих, равно как и хорошо известное двумерное отображение Эно, его трехмерные аналоги не являются чисто искусственными. В частности, два следующих вида 3БНМ: первый — с Ы\ = 0, Ь = 0, е = 0 и а = —1, то есть
х = у, у = г, г = Мг + М2у + Вх - г2, (1.2)
и второй — с Ы2 = 0,а = 0,Ь = 0 и е = —1, то есть
х = у, у = г, 1 = М\ + М22 + Вх — у2, (1.3)
представляют особый интерес для нелинейной динамики, поскольку они имеют «гомоклини-ческое происхождение», то есть являются нормальными формами отображений первого возвращения вблизи гомоклинических касаний некоторых типов (см., например, [12, 15, 22-25]).
Содержание работы. В § 2 формулируется постановка задачи и дается краткое описание двух основных сценариев возникновения хаоса у трехмерных отображений; §§ 3 и 4 составляют основную часть работы, в них приводится развернутое изложение предложенных сценариев, теоретические положения которых иллюстрируются результатами численного счета; в § 5 обсуждается один из главных механизмов возникновения хаотической динамики (имеющий место и в двумерном случае), который присутствует в обоих сценариях — это переход к хаосу в результате разрушения тора (замкнутой инвариантной кривой — в случае отображений).
3Последние два условия весьма существенны: они исключают некоторые типы квадратичных отображений с постоянным якобианом. Например, отображение, обратное к х = х, у = у + ж2, 1 = = г + у2, будет уже четвертой степени; условие (ш) исключает «тривиальные случаи», когда общая динамика отображения определяется только частью координат (как, например, у отображения х =
= X, у = г, £ = г2 - у).
2. Постановка задачи и описание основных сценариев возникновения хаоса у трехмерных диффеоморфизмов
Пусть Т^ — однопараметрическое семейство трехмерных ориентируемых гладких отображений. Предположим, что это семейство Т^ непрерывно зависит от параметра у и удовлетворяет при этом следующим условиям.
A1 Существуют такие уо и fi\ (положим для определенности, что уо < fi\), что отображение Тц имеет при всех у Е (уо,у1) асимптотически устойчивую неподвижную точку
Оц. При этом точка О^ является глобальным аттрактором в том смысле, что положительные итерации всех точек из некоторой достаточно большой области фазового пространства сходятся к Ор.
A2 Значение у = у1 параметра — бифуркационное, при котором точка О^ теряет устойчивость и становится седловой при у > Ц1.
A3 Существуют такие у2 и у3, где у3 > у2 > у1, что при у Е (у2,у3) точка Ор имеет гомоклинические к ней траектории, то есть устойчивое и неустойчивое инвариантные многообразия точки О^ имеют непустое (гомоклиническое) пересечение.
В этой работе мы будем считать, что
Но — это бифуркационное значение параметра у, при котором в семействе Т^ происходит невырожденная седло-узловая бифуркация: при у < уо у Т^ нет неподвижных точек (вообще говоря, в некоторой области фазового пространства), а при у > уо у у появляются две неподвижные точки: асимптотически устойчивая О^ и седловая О^, причем последняя точка имеет тип (2,1), то есть dim Ws(0р) = 2, dim Wи(О^) = 1;
у1 — это бифуркационное значение параметра у, при котором точка О^ мягко теря-
ет устойчивость: типично одним из двух основных способов — либо в результате (i) бифуркации Андронова-Хопфа, либо в результате (ii) бифуркации удвоения периода; при этом сама точка О^ становится седловой, а в ее окрестности рождается либо устойчивая замкнутая инвариантная кривая — в первом случае, либо устойчивый цикл периода 2 — во втором;
у2 — это некоторое значение параметра у, при котором у седла О^ появляется гомокли-ническая траектория: в этом случае W**(Оц) П Wи(О^) = 0:
у3 — это некоторое значение параметра у, при котором появляется гетероклиническая траектория между седлами Оц и О^: в этом случае Wи(О^) П Ws(0ц) = 0;4
— при у Е (у2,уз) существует странный гомоклинический аттрактор Л^, содержащий точку Ор вместе с ее неустойчивым инвариантным многообразием.
4Здесь Ws(0ц) всегда двумерно, но W“(0M) может быть либо одномерным (как в случае (i) бифуркации Андронова-Хопфа), либо двумерным (как в случае (ii) бифуркации удвоения периода). Соответственно, возникновение гетероклинических пересечений между W“(Ом) и Ws (Ом) — это бифуркационное явление коразмерности один, то есть оно может свободно наблюдаться в однопараметрических семействах общего положения.
Замечание 1. Мы считаем, что ро и р\ — это точные бифуркационные значения параметра р. Однако это нельзя, вообще говоря, утверждать относительно р2 и рз. Дело в том, что в трехмерном случае, в отличие от двумерного, не всегда возможно определить (даже теоретически) момент первого гомоклинического или гетероклинического пересечения (касания).5 Поэтому нужно считать, что р2 и рз — это некоторые значения параметра р, при которых обнаруживаются соответствующие гомоклиническое и гетероклиническое пересечения.
Замечание 2. В этой статье мы предполагаем, для определенности, что при р € (р2, Рз) точка 0^ принадлежит некоторому «гомоклиническому» странному аттрактору Лм. Однако понятно, что это не всегда так. Например, если Лм — квазиаттрактор, то на интервале (р2,рз) могут существовать так называемые «окна устойчивости», когда аттрактором является некоторая устойчивая периодическая точка или замкнутая инвариантная кривая. Более того, в реальных моделях могут наблюдаться и другие ситуации. Например, в известной механической системе «кельтский камень» [26, 27], хаос может возникать жестким образом (уже после того как образуются гомоклинические и гетероклинические структуры) в результате потери устойчивости предельного цикла (так называемого цикла Карапетяна). Здесь в критический момент этот цикл либо сливается с седловым того же периода и исчезает, либо в него «влипает» седловой цикл удвоенного периода (то есть происходит так называемая «обратная бифуркация удвоения»). В обоих случаях сразу после соответствующей бифуркации наблюдается хаотическая динамика.6
(а) (Ь) (с)
Рис. 1. Типы седло-узлов у ориентируемых 3D-отображений. На рисунке пунктирная окружность имеет радиус, равный модулю ведущего устойчивого мультипликатора.
Что касается седло-узловой бифуркации (при р = ро), то в случае трехмерных ориентируемых отображений существуют три основных «гладких» типа таких бифуркаций, в зависимости от расположения устойчивых мультипликаторов (внутри единичной окружности). Так как Т^ ориентируемо, то есть произведение всех мультипликаторов положительно, это отвечает трем основным типам седло-узлов, которые схематически представлены на рисунке 1: здесь либо оба устойчивых мультипликатора (а) действительны и положительны, либо
5Например, в случае неподвижной точки седло-фокус (2,1) ее неустойчивое многообразие может навиваться на самоё себя при приближении к моменту касания. Поэтому определить, когда в точности наступит этот момент, практически невозможно, поскольку любой наперед выбранный кусок многообразия будет окружен со всех сторон другими его же кусками.
6 Которую, однако, нельзя автоматически ассоциировать со странным аттрактором. Дело в том, что указанная система является обратимой — реверсивной (но не консервативной), и хаос может включать периодические траектории с якобианом как больше, так и меньше единицы. Поэтому здесь должна наблюдаться так называемая смешанная динамика [28-30], когда гиперболические периодические траектории всевозможных типов (устойчивые, вполне неустойчивые и седловые) сосуществуют, и типично, когда они в совокупности не отделимы друг от друга.
(Ь) действительны и отрицательны, либо (с) комплексно-сопряжены. Тогда седло в случае (а) будет иметь топологический тип (2+, 1+) (т. е. все мультипликаторы — два устойчивых и один неустойчивый — действительны и положительны), в случае (Ь) — топологический тип (2—, 1+) (устойчивые мультипликаторы отрицательные, неустойчивый — положительный) и в случае (с) точка — седло-фокус (2,1) (устойчивые мультипликаторы комплексно сопряжены, а неустойчивый — действительный положительный). (См. рис. 2.)
Рис. 2. Типы седло-узловых бифуркаций в соответствии с типами седло-узловых точек (а), (Ь) и (с) рисунка 1. Во всех случаях возникают седла типа (2,1). В случаях (а) и (Ь) — картинки одинаковые, но на ШЯ(О) имеется либо узел-плюс (а), либо узел-минус (Ь). В случае (с) рождается седло-фокус.
Отметим, что существование седловой точки Оц имеет весьма важное значение для хаотической динамики в целом, поскольку ее двумерное устойчивое многообразие образует естественную границу области притяжения как для самой точки О^ при ро < Р1, так и для всех тех аттракторов, которые она порождает (в том числе и странных) при р > р1. Поэтому естественно ожидать, что кризисы7 таких СА и их разрушения будут связаны с образованием гетероклинических пересечений между неустойчивыми инвариантными многообразиями точек на аттракторе и границей устойчивости Шя(0^). И здесь мы можем сделать следующее важное замечание.
Замечание. О кризисах СА, содержащих точку Ом. По предположению А3, при значениях р € (р2, рз), точка Оц принадлежит некоторому «гомоклиническому» аттрактору Л^, так как Ш“(ОМ) ПШЯ(ОМ) = 0. Момент р = р3, то есть когда Ш“(ОМ) начинает пересекать границу Шв(0^) глобальной области устойчивости, отвечает кризису этого аттрактора. Здесь аттрактор либо мгновенно разрушается8, либо, скачком, возникает больший СА, поглощающий Лм и содержащий точку
Ом. С другой стороны, если рз < р2, то есть если такое гетероклиническое пересечение образуется раньше, чем гомоклиническое пересечение инвариантных многообразий точки Ом, то странные аттракторы, содержащие точку Ом, вообще не появятся.
7Кризис аттрактора, в отличие от внутренних бифуркаций, резко меняет его структуру так, что данное предельное множество (в смысле топологического предела при стремлении значения параметра к критическому значению) — уже не аттрактор (см. точное определение в [32]): аттрактор здесь либо разрушается, либо занимает совсем другую область фазового пространства, либо поглощается большим аттрактором. Последнее, например, может происходить при жесткой потере устойчивости периодическим или квазипериодическим (тором) аттрактором, находящимся в поглощающей области.
8 Почти все его точки убегают либо на бесконечность, либо на какой-либо другой далекий аттрактор; в случае, когда Лм — квазиаттрактор, могут, однако, сохраняться некоторые устойчивые периодические траектории больших периодов вместе с их областями притяжения.
(аНЪ)
(с)
Таким образом, в системах общего типа седло-узловая бифуркация играет двоякую роль: в результате такой бифуркации, во-первых, возникает динамика и образуется тривиальный аттрактор (устойчивая точка)9 и, во-вторых, формируются границы глобальной области притяжения аттракторов (в том числе и странных), ассоциированных с устойчивой точкой (точкой О^ в нашем случае).
Что касается бифуркаций потери устойчивости точки О^ (при р = р1), то мы будем считать, что они во всяком случае невырожденны. Тогда это либо бифуркация Андронова-Хопфа, либо бифуркация удвоения периода. В первом случае при р = р1 на единичную окружность выходит пара мультипликаторов в±г1р (при условии отсутствия сильных резонансов, т. е. при ф = 0,п,п/2, 2п/3), а во втором — один мультипликатор —1. Эти бифуркации настолько важны для дальнейшего, что по их типу мы различаем два основных сценария развития хаоса у трехмерных отображений. Эти сценарии представлены ниже.
Сценарий 1 (сценарий Ш^ильникова [3])
1) При р = р1 у точки О^ появляется пара комплексно-сопряженных мультипликаторов в±г1р, где ф = 0,п,п/2, 2п/3, и здесь происходит мягкая бифуркация Андронова -Хопфа, то есть при р > р1 точка О^ становится седловой — седло-фокус типа (1, 2), а в ее окрестности рождается асимптотически устойчивая замкнутая инвариантная кривая Ьц.
2) При дальнейшем изменении р кривая Ь^ разрушается одним из известных способов (см. § 5), в результате чего в ее окрестности возникает СА (вообще говоря, типа «тор-хаос»). Точка Оц при р < р2 лежит в стороне от этого СА.
3) При р > р2 образуется гомоклинический аттрактор — аттрактор Шильникова, содержащий точку Ом. Этот аттрактор может существовать вплоть до момента р = р3 (когда возникает гетероклиническое пересечение между Ши(Ои Шя(0^)).
Сценарий 2
1) При р = р1 у точки О^ появляется мультипликатор —1 и происходит мягкая бифуркация удвоения периода, то есть при р > р1 в окрестности точки О^ рождается устойчивый цикл периода 2 (состоящий из двух точек Р1 и Р2, таких, что Т^(р{) = Р2 и ТДР2) = Р1), а сама точка становится седловой неподвижной с мультипликаторами Аъ А2, А3, такими, что \1 < —1, 0 < \2 < 1, —1 < \3 < 0.
2) При дальнейшем изменении р с циклом происходит серия бифуркаций, приводящая к возникновению СА (например, типа тор-хаос, когда отображение Т^ не является сильно-диссипативным). При этом (седловая) точка О^ лежит в стороне от СА при р < р2.
3) При р € (р2,рз) точка О^ принадлежит гомоклиническому СА Л^, тип которого зависит главным образом от мультипликаторов А2 и Аз седла О^:
я) если |А21 > |Аз|, то аттрактор «квазилоренцевский».
Ъ) если |А21 < |Аз|, то аттрактор «восьмерочный».
9В случае, когда седло-узел не имеет гомоклинической к нему траектории; в противном случае хаос может возникать даже при его исчезновении [8, 31].
В обоих случаях СА пунктов 3 момент р = р3 является кризисным, то есть сразу после того как образуется гетероклиническое пересечение инвариантных многообразий Ши(0^) и Шв(0^), аттрактор либо разрушается, либо поглощается большим аттрактором, содержащим точку Оц.
В §§ 3 и 4 мы изучим подробнее предложенные сценарии и покажем, как они реализуются в ряде конкретных трехмерных отображений (в квадратичных трехмерных отображениях Эно).
3. О сценарии Шильникова (сценарий 1)
На первом шаге в критический момент р = р1 на единичную окружность выходит пара комплексно-сопряженных мультипликаторов А12 = е±ш точки Ом. И мы предполагаем здесь, что в семействе Тц при р = р1 происходит мягкая невырожденная бифуркация Андронова-Хопфа. Это условие общего положения для однопараметрических семейств, и оно означает (в трехмерном случае), что неподвижная точка О^ при р > р1 становится седло-фокусом типа (1, 2), то есть с двумерным неустойчивым многообразием, а в ее окрестности рождается гладкая замкнутая инвариантная кривая Ь^. Эта кривая асимптотически устойчивая, и она служит границей многообразия Ши(О) — по крайней мере для всех значений р, близких к р1 (см. рис. 3).
бифуркация Андронова - Хопфа
Рис. 3. (а) При р < р\ аттрактором является неподвижная точка Ом типа фокус. При этом все траектории, кроме траекторий точек на неведущем (сильно устойчивом) инвариантном многообразии Ш88, входят в Ом, касаясь ведущей плоскости Р8. (Ь) Сразу при р > р\ аттрактором является замкнутая инвариантная кривая Ьм, а точка Ом становится седло-фокусом типа (1, 2).
Отметим, что условие невырожденности бифуркации Андронова-Хопфа в случае отображений включает также требование отсутствия сильных резонансов (резонансов 1:1, 1:2, 1:3 и 1 : 4), то есть и = 0,п,п/2, 2п/3. Когда в бифуркационный момент присутствует резонанс р/д, где д ^ 5, замкнутая гладкая инвариантная кривая все равно рождается, но на ней будут лежать чередующиеся седла и узлы периода д. Поэтому нужно считать, что для всех р > р1, близких к р1, (глобальным) аттрактором отображения Т^ является либо вся инвариантная кривая Ь^, либо устойчивая периодическая траектория на ней в резонансном случае.
При дальнейшем изменении параметра р кривая Ь^ может потерять устойчивость. Основные бифуркационные механизмы потери этой устойчивости в целом известны (см., например, [32, гл. 2, § 1.6 и гл. 3, § 8.5]. В принципиальном плане, можно выделить следующие:
• разрушение инвариантной кривой (по Афраймовичу-Шильникову [8]);
• «удвоение тора» (так называют серию бифуркаций, в результате которых исходная инвариантная кривая теряет устойчивость, передавая ее новой инвариантной кривой «двойной длины»), подробнее об этом см. в [37];
• «жесткая потеря устойчивости тора», когда, типично, либо (1) инвариантная кривая исчезает (в результате «седло-узловой бифуркации» по Шенсине (А. СЬепстег) -так обозначают некоторую серии бифуркаций, в результате которой устойчивая инвариантная кривая «сливается» с неустойчивой (седловой) и исчезает (может быть и в двумерном случае), либо (п) происходит «обратное удвоение тора», когда, в результате опять же некоторой серии элементарных бифуркаций, в устойчивую инвариантную кривую «влипает» седловая удвоенной длины: при этом исходная инвариантная кривая становится также седловой;
• «прохождение через резонанс», когда на (асимптотически устойчивой) замкнутой инвариантной кривой появляются, типично, две периодические траектории одного и того же периода, устойчивая и седловая: если трехмерное отображение сильно диссипативное (или период достаточно большой), то здесь, в силу нормальной гиперболичности, ситуация похожа на двумерную (см. [8, 32]); но когда якобиан отображения близок к единице и период относительно не очень большой, трехмерность играет существенную роль, и здесь могут возникать весьма необычные СА (см., например, рис. 7, на котором показаны этапы образования «суперспирального» аттрактора при разрушении инвариантной кривой вблизи резонанса 1:4).
В первых двух случаях возникающий СА называют обычно тор-хаосом, и в § 5 мы дадим описание основных бифуркационных явлений, сопровождающих его появление. В третьем случае тип СА зависит главным образом от специфики задачи. Случай «прохождения через резонанс» представляет особый интерес еще и потому, что здесь, как нам кажется, нет пока достаточно законченной теории.
Заметим, что потеря устойчивости кривой Ь^ и возникновение тор-хаоса естественно должно сопровождаться изменением геометрии неустойчивого многообразия точки Ом, а именно: можно предположить, что в случае нашего семейства Т^ существует такое р1, что при р1 < р < р\ многообразие Ши(0^) является диском с границей Ь^ (см. рис. 4а), а при р > р\ оно начинает «наматываться» на Ь^ (см. рис. 4Ь). Момент р = р1 не является по-настоящему бифуркационным: здесь, как говорят, происходит «гладкая бифуркация», то есть меняется «гладкий тип» инвариантной кривой Ь^. Именно, при р\ < р < р1 инвариантная кривая Ь^ имеет узловой тип, а сразу при р > р\ — фокусный (фокальный): в первом случае ляпуновские показатели Л1, Л2, Л3 на Ь^ все действительные и такие, что Лз < Л2 < 0 < Л1 = 0; во втором случае мы имеем, что И,е(Лз) = И,е(Л2) < 0 < Л1 = 0.
Геометрически это означает, что при р > р\ образуется так называемая воронка Шиль-никова, в которую «втягиваются» все траектории из поглощающей области (все, кроме седла О^ и траекторий точек одной из его устойчивых сепаратрис). При дальнейшем увеличении р это «втягивание» будет сохраняться, однако поведение инвариантных многообразий седла О^ может меняться. Более того, по предположению А3, наступит такой момент, р = р2, когда эти многообразия пересекутся. Тогда точка О^ будет принадлежать аттрактору вместе с Ши(0^), и здесь можно говорить об образовании аттрактора Шильникова. Если при этом точка О^ — седло-фокус, то такой аттрактор является «спиральным», и он относится к классу квазиаттракторов [12, 15]. Однако если О^ — седло, то есть имеет пару действительных и различных по модулю неустойчивых мультипликаторов, то такой аттрактор, в принципе, может быть уже «настоящим» (см. дискуссию в §4.1).
образование
«воронки»
V
(а) цх < ц < р\
Рис. 4. На рисунке (Ь) показана «воронка Шильникова», часть которой как бы срезана, чтобы лучше показать внутреннюю структуру многообразия Ш“(0М).
Подобный сценарий возникновения хаоса описал Шильников в статье [3]. Однако нужно заметить, что наше предположение А3 о том, что гомоклиническое пересечение многообразий Ши(0р) и Шя(0р) наступает раньше, чем гетероклиническое пересечение многообразий Ши(0р) и Шя(0р), то есть ^2 < рэ, является принципиально важным. В противном случае, если р2 > рэ, аттрактор Шильникова не возникает и развитие хаоса заканчивается на стадии, например, тор-хаоса, который разрушается (при р > рЭ > Рэ), как только возникает гетероклиническое пересечение между неустойчивыми многообразиями его точек
В качестве иллюстрации рассмотренного сценария мы приведем некоторые результаты численного счета для трехмерных отображений Эно. На рисунке 5 показаны стадии формирования аттрактора Шильникова в отображении (1.3) с фиксированными В = 0.7 и М\ = 0 при изменении М2 от 0 до 0.129 (при М2 ~ 0.13 аттрактор разрушается). На рисунке 6 показаны стадии образования (по Афраймовичу-Шильникову) тор-хаоса в 3БИМ (1.2): здесь мы зафиксировали В = 0.7 и М2 = —0.5 и меняли М\ от 0 до 0.3. В этом случае аттрактор разрушался при М2 ~ 0.31 (здесь имеет место случай рэ < Р2). На рисунке 7 показан пример образования необычного («суперспирального» или «спирально-винтового») СА при прохождении вблизи резонанса: в данном случае — резонанса 1 : 4.
4. Сценарий 2: возникновение лоренцевского или «восьмерочного» СА
На первом шаге в критический момент р = р\ у точки Ор появляется мультипликатор — 1. Мы предполагаем, что соответствующая бифуркация удвоения периода является мягкой и невырожденной. Тогда при р < р\ неподвижная точка Ор асимптотически устойчива (это точка типа «узел-минус», так как ее ведущий мультипликатор отрицательный), а при р > р\ она становится седлом типа (2,1), то есть с двумерным устойчивым инвариантным многообразием,10 и в окрестности точки Ор рождается устойчивый цикл периода два.
10Здесь мы как раз и видим одно из принципиальных отличий от сценария 1, в котором, напомним, при потере устойчивости неподвижная точка становится седловой с одномерным устойчивым инвариантным многообразием.
и Ш8(0р).
0.60
0.55
0.50
0.45
0.40
0.40 0.45 0.50 _ 0.55 0.60
Хх
0.8
0.7
0.6 . ■ .
Щ. "С /•'
0.5 ■
0.4
0.3
0.40 0.45 0.50 0.55 0.60 0.65 0.70 0.75 X,
1.2
1.0
0.8
Хо ^0.6
0.4
0.2
0.0
0.0 0.2 0.4 0.6 0.8 1.0 1.2
X,
(а)
(Ь)
*2
-Г"
// (/
1.2 1.0 0.8 -0.6 -0.4 0.2 -0.0 -
. '
■
.. ■
■■ ■■
-
■■
-0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4
0.0 0.2 0.4 0.6 0.8 1.0 1.2
0.0 0.2 0.4 0.6 0.8 1.0 1.2
XI
№
(е)
(!)
щ.
(8)
.; '■ •■■ ■ ■* ч' у ■ .
'
0.0 0.2 0.4 0.6 0.8 1.0 1.2 Ху
(Ь)
1.4
1.2 1.2
1.0
0.8 0.8
Ъ 0.6 0.4
0.4
0.2 0.0
0.0
-0.2 -0.4
-0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 -0.4 0.0 0.4 0.8 1.2
(1)
О)
Рис. 5. Эволюция аттракторов в 3DHM (1.3) при фиксированных В = 0.7 и Ы\ = 0: (а) М2 = 0.8 (неподвижная точка); (Ь) М2 = 0.875 (замкнутая инвариантная кривая); (с) М2 = 1.21 (инвариантная кривая «удвоилась»); (^ М2 = 1.24; (е) М2 = 1.245 (второе удвоение); (1) М2 = 1.25; (8) М2 = = 1.255; (Ь) М2 = 1.258 (развитие тор-хаоса); (1) М2 = 1.27 (тор-хаос); 0) М2 = 1.29 (аттрактор
Шильникова).
0.5
0.2 0.8
0.4
0.6
0.1 0.3 0.4 ' ' \.
0.2 !
0.0 0.1 ' 0.2 ■ N • . ■ ■
X2 *! \ ' о о ьГ . ■. ■ ■ V ( .
-0.1 \ -0.2
-0.1 ' ■
-0.4 ■ • ■ ■. ■ ■ ■ ■ ЦШ \
-0.2 -0.2 ■ V . ' .■■■ . ' ■■■■ 1 ■ . ..•••—' ■ ■ '
-0.6 ■
-0.3 ч-г ^ -----
-0.3 -0.8
-0.4
-0.20-0.15-0.10-0.01 0.00 0.05 0.10 0.15 0.20 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6
Ху Ху Ху
(а) (Ь) (с)
-0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 -0.8 -0.4 0.0 0.4 0.8 -1.2 -0.8 -0.4 0.0 0.4 0.8
Ху Ху Ху
№ (е) (Г)
Рис. 6. Эволюция аттракторов в 3DHM (1.2) при В = 0.7, М2 = —0.5: (а) Ы\ = 0 (неподвижная точка); (Ь) М\ = 0.05 (замкнутая инвариантная кривая); (с) М\ = 0.15 (на кривой виден резонанс типа 4 : 13); (^ М1 = 0.2 (кривая негладкая); (е) М1 = 0.25 (появление тор-хаоса); (Г) М1 = 0.3
(развитый хаос — аттрактор с условным названием «корабль зергов»).
Этот цикл состоит из двух точек р1 и Р2, таких, что Тр('р\) = Р2 и Тр(р2] = Р1, а его мультипликаторы, т. е. собственные числа матрицы линеаризации ОТ2\ ,1 = 1, 2, все по модулю
№ фг
меньше единицы.
В момент, близкий к ц = Ц1, (глобальным) аттрактором отображения Т^ является (по предположению А2) этот устойчивый цикл {р1,Р2} (см. рис. 8). При дальнейшем изменении параметра ц с этим циклом происходят бифуркации, соответственно тип аттрактора меняется, и он становится странным типично одним из известных способов, например, в результате таких цепочек бифуркаций:
(1) бесконечная последовательность бифуркаций удвоения ^ СА (фейгенбаумовский квазиаттрактор),
(п) бифуркация Андронова-Хопфа ^ инвариантный тор ^ тор-хаос,
(ш) конечная последовательность бифуркаций удвоения ^ инвариантный тор ^ тор-хаос.
Заметим, что последовательность (1) бифуркаций характерна в основном только для сильно-диссипативных систем11 или сохраняющих площадь отображений. Тор-хаос в случаях (п) и (ш) может возникать опять одним из известных способов, о которых мы говорили в предыдущем параграфе. (См. также §5, в котором мы более подробно рассмотрим сценарии «разрушение тора», «удвоение тора» и «седло-узловая бифуркация тора».)
пНапример, у отображения Эно (1.2) при В > ^ уже нет второго удвоения у устойчивой точки.
Ху
(а)
№
0.0 0.2 0.4 0.6 0.8 1.0
Х\
(Ь)
0.4 0.8
(е)
-0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2
Ху
(с)
>>
у
ш ^ '
/ 1 /У' ^ .
1. ' '">■ \ ;«\>\| \
-0.4 0.0 0.4 0.8 1.2
Ху
(Г)
Рис. 7. Эволюция аттракторов в 3DHM (1.3) при В = 0.7, М2 = 0.8: (а) М1 = 0 (неподвижная точка); (Ь) М1 = 0.1; (с) М1 = 0.26; (^ М1 = 0.335 (замкнутая инвариантная кривая); (е) М1 = = 0.34 (тор «рассыпался»); (Г) М1 = 0.35 («суперспиральный» аттрактор, вблизи резонанса 1 : 4, он содержит два седло-фокуса периода 4: один с одномерным, другой с двумерным неустойчивыми многообразиями).
Естественно, что СА в указанных цепочках могут быть совершенно разной структуры и разной природы. Однако при р = р2, то есть в тот момент, когда образуются гомокли-нические пересечения многообразий исходного седла Ор, происходит определенного рода унификация. В этом случае точка Ор (автоматически) будет принадлежать аттрактору вместе с Ши(Ор). Тогда (вообще говоря, при р ^ р2 и достаточно близких к Р2) можно выделить два основных геометрических типа таких глобальных гомоклинических аттракторов, в зависимости только от характера мультипликаторов А1, Х2, Аз самой точки Ор. Мы предполагаем здесь, что
А1 < —1, 0 < А2 < 1, —1 < А3 < 0
(т. е. седло Ор сохранило свой топологический тип, какой у него был после бифуркации удвоения периода, происходившей при р = Р1). Тогда СА в рассматриваемом случае (при р2 < < р < рз) можно разбить на два следующих типа.
• Квазилоренцевские аттракторы, если неведущим (сильно устойчивым) мультипликатором является отрицательный, то есть Аз; соответственно, здесь |Аз| < А2.
(а)
(Ь)
(с)
Рис. 8. Основные этапы развития гомоклинического аттрактора лоренцевского типа.
• Восьмерочные аттракторы, если неведущим (сильно устойчивым) мультипликатором является положительный, то есть А2; соответственно, здесь А2 < |Аз|.
В первом случае можно говорить об образовании аттрактора, очень похожего на аттрактор Лоренца. Сходство здесь может быть как чисто условным (см. рис. 8с), так и геометрическим, в смысле фазового портрета (см., например, рис. 12е, £). В принципе, можно указать следующие характерные признаки квазилоренцевских аттракторов:
• точка Ор является седлом, лежащим строго на границе аттрактора, а кусок многообразия Ши(Ор) с точкой Ор является ребром этого аттрактора, к которому «приклеены его ветви-листы» (см. рис. 9),
• сильно устойчивое многообразие Ш33 (О р) делит Ш3(Ор) (локально) на две половинки, и гомоклинические точки принадлежат только одной из них,
• аттрактор чаще всего имеет две «дырки», внутри которых типично лежит седловая точка периода два (однако иногда эти «дырки» могут «замазываться» за счет образования соответствующих гетероклинических пересечений, см. рис. 13 и пояснительный текст к нему), он также может иметь лакуны (здесь: трубчатые окрестности седловых замкнутых инвариантных кривых, не пересекающиеся с аттрактором, как на рис. 12е
12 Это условие является одним из необходимых для того, чтобы сам аттрактор или все близкие к нему не содержали бы устойчивых периодических траекторий (периодических аттракторов). Оно автоматически выполняется в случае классического аттрактора Лоренца (для отображения сдвига на единицу времени). Однако не ясно, насколько оно необходимо для геометрии аттрактора. В частности, кажется, что бифуркация «объединения» инвариантных кривых в одну, иллюстрируемая на рис. 14, происходит благодаря их «прохождению через гомоклиническую восьмерку», у которой а < 1, — когда а > 1, такое невозможно, так как кривые устойчивые. Аттрактор в это время является, конечно, квазиаттрактором.
и 14Ь),
• седловая величина а = |А1А21 больше единицы.12
(с) (а)
Рис. 9. Вид СА лоренцевского типа в случае отображения (1.2) (при В = 0.7, Ы\ = 0, М2 = 0.85) в различных ракурсах (повороты относительно вертикальной оси). Жирными точками отмечены неподвижные точки: Ом — на аттракторе и Ом — вне его.
«Восьмерочный аттрактор» не будет похож на квазилоренцевский по многим признакам. В частности, можно указать такие, как
• точка Ор является седлом, лежащим «в середине» аттрактора, а «ветви-листы» аттрактора «приклеены» к Ши(Ор) гладко,
• сильно устойчивое многообразие Ш33 (О р) делит Ш3(Ор) (локально) на две половинки, и гомоклинические точки принадлежат обязательно им обеим.
Это можно извлечь из простого наблюдения того, как могут быть устроены гомокли-нические пересечения многообразий седла Ор в случае Х\ < — 1, 0 < А 2 < 1, —1 < Аз < 0 при |А21 < |Аз | (см. рис. 10). Однако по-прежнему такой аттрактор имеет две «дырки», и условие а = |А1 Аз| > 1 является необходимым для того, чтобы сам аттрактор и все близкие к нему не содержали бы устойчивых периодических точек.
(а)
(Ь)
Рис. 10. Гомоклинические структуры для СА: (а) квазилоренцевского, (Ь) восьмерочного.
4.1. О псевдогиперболических (диких) странных аттракторах
Для обоих типов странных аттракторов, «квазилоренцевского» и «восьмерочного», можно поставить вопрос: при каких условиях эти аттракторы являются настоящими, то есть когда они сами и все близкие не содержат устойчивых периодических траекторий? Как вытекает из результатов Тураева и Шильникова [10, 11], они тогда должны обладать так называемой псевдогиперболической структурой: отображение на аттракторе должно быть устроено так, что всегда есть направление сильного сжатия, а на трансверсальном сечении отображение должно растягивать площади. В частности, это означает, что мультипликаторы А1, А2, Аз в седле Ор должны быть такими, что
• |А1А2Аз| < 1 (сжатие объемов необходимо для существования аттрактора),
• |А11 > 1 и |А1 > |А21 > |Аз | (все мультипликаторы должны быть простыми, в частности, это означает, что точка Ор у всех близких систем является седлом, а не седло-фокусом) — это условие обеспечивает существование как сильно неустойчивого, так и сильно устойчивого направлений;
• | А1А21 > 1 — растяжение площадей в ортогональной к сильно устойчивому направлению площадке.
Такие аттракторы допускают гомоклинические касания и, следовательно, могут принадлежать областям Ньюхауса. Однако при наличии псевдогиперболической структуры ни сам аттрактор, ни все близкие не будут содержать устойчивых периодических траекторий. Такие аттракторы были названы дикими гиперболическими аттракторами. Терминология восходит к Ньюхаусу, который ввел понятие дикого гиперболического множества, а Тураев и Шильников предложили использовать соответствующий термин для обозначения указанных аттракторов, математическую теорию которых они, собственно, и построили [10, 11]. Важно, что существование гомоклинических касаний здесь неизбежно. С другой стороны, известно (критерий Гонченко, Тураева и Шильникова, [12, 13]), что определенного типа гомоклинические касания могут порождать устойчивые периодические траектории. Это означает, например, что дикие аттракторы не могут иметь гомоклинических касаний
к точкам с |Л11 = |Л 21 > 1 или с |Аз| = |Л2| < 1, то есть, в частности, седло-фокусы запрещены.
-0.15 -0.1 -0.05 0 0.05 0.1
Рис. 11. Карта ляпуновских показателей (из [18]) для области Б\: {—0.15 < Мх < 0.1; 0.7 < М2 < 1} параметров при В = 0.7 в случае 3DHM (1.2). На цветной версии рисунка* зеленым, голубым, красным, синим и желтым показаны области значений параметров, отвечающие траекториям с разным набором ляпуновских показателей Лх, Л2 и Л3. На черно-белом рисунке эти цвета расположены слева-направо вдоль, например, линии М2 = 0.8.
Отметим, что в работе [18] были численно проверены условия псевдогиперболичности для аттракторов квазилоренцевского типа в случае отображения (1.2). Было показано, что эти условия выполняются для значений параметров из некоторых открытых областей. В частности, такие области были обнаружены при В = 0.7. На рисунке 11, взятом из [18] показано разбиение области 0\: {—0.15 < М\ < 0.1; 0.7 < М2 < 1} на участки, отвечающие разным типам асимптотического поведения траекторий на аттракторе. Другими словами, здесь изображена так называемая карта ляпуновских показателей. На этой карте разными цветами отмечены значения параметров, отвечающих траекториям с разным набором ляпуновских показателей Л\, Л2 и Л3. Во-первых, белые участки отвечают тем значениям параметров, при которых итерации начальной точки покидают некоторую заданную окрестность фазового пространства. Положим для определенности, что Л1 ^ Л2 ^ Л3. Тогда
1. зеленый цвет отвечает случаям, когда Л1 < 0, то есть периодическим аттракторам,
2. голубой цвет отвечает случаям, когда Л1 =0 > Л2 > Л3, то есть устойчивым замкнутым инвариантным кривым,
3. красный цвет отвечает случаям, когда Л1 > 0 = Л2 > Л3,
*Для читателя печатной версии: здесь и далее полноцветные версии рисунков см. в эл. версии статьи — http://nd.ics.org.rU/doc/r/pdf/1966/0
4. синий цвет отвечает случаям, когда Л1 > Л2 > 0 > Л3,
5. желтый цвет отвечает случаям, когда Л1 > 0 > Л2 > Л3.
Последние 3 случая отвечают странным аттракторам. Причем здесь знак «=» (при красном цвете) означает только, что отличие Л2 от 0 не больше, чем условленные границы точности вычисления ляпуновских показателей (в [18] эти границы выбирались между 10_5 и 10_6 в зависимости от типа вычислений).
Замечание 3. Для СА, представленного на рисунке 12^ была тщательно проверена зависимость старшего ляпуновского показателя Лх при В = 0.7,М2 = 0.85 от М\. Было численно установлено, что график Л^Мх) в окрестности точки М\ = 0 является гладкой кривой, одной и той же для разных начальных условий (даже когда исходный интервал по М\ длиной порядка 10~2 был разбит на 108 интервалов в эксперименте К.Симо). Обычная ситуация в случае квазиаттракторов — прямо противоположная: здесь (если захотеть) на любом участке значений параметра можно обнаружить фрактальную структуру соответствующего графика.
Используя такую карту, как на рисунке 11, можно легко проследить, при изменении параметров, эволюцию притягивающих множеств 3БИМ (1.2) с В = 0.7 от простого аттрактора к странному. На рисунках 12 и 13 показаны результаты при фиксированном М2 = 0.85, а на рисунке 14 — при фиксированном М2 = 0.75.
Некоторые результаты счета для 3БИМ (1.3) были показаны выше на рисунке 5. Здесь финальный аттрактор — это аттрактор Шильникова. Поскольку он спиральный (содержит седло-фокус Ор), то это квазиаттрактор по критерию из [12]. Однако промежуточные аттракторы типа тор-хаос (рис. 5£, g и Ь) вполне могут претендовать на роль настоящих СА, обладающих псевдогиперболической структурой.
5. О бифуркационных механизмах возникновения тор—хаоса
Как мы видим, бифуркации, связанные с возникновением тор-хаоса, в рассмотренных сценариях играют весьма важную роль, и, в принципе, если ^3 < Ц2, то как раз они могут приводить к финальным СА, не содержащим точку О р. Естественно, таких аттракторов, в отличие от финальных гомоклинических (содержащих Ор), может быть бесконечное множество различных типов. Вряд ли их возможно как-то разумно классифицировать, однако мы можем описать основные особенности бифуркационных механизмов возникновения тор-хаоса (см. §3). Два из них, «удвоение тора» и «жесткая потеря устойчивости тора», в главных чертах похожи на то, что происходит с циклами в (трехмерных) потоках (имеются в виду бифуркации удвоения периода циклов и бифуркации слияния двух циклов), и поэтому мы дадим ниже лишь их краткое описание. Другой механизм, «прохождение через резонанс», мы здесь не будем обсуждать, потому что он слишком специальный, да и возникающие здесь аттракторы могут иметь весьма слабую ассоциацию с обычным «тор-хаосом» (см., например, рис. 7£). Механизм же «разрушение тора» (по Афраймовичу-Шильникову) существенно отличается от того, что нам дает потоковая интуиция, и поэтому мы обсудим его сейчас более подробно.
5.1. О механизме Афраймовича —Шильникова развития тор—хаоса
Это один из типичных способов, когда возникает СА типа «тор-хаос» вследствие смены квазипериодического режима на хаотический в результате разрушения замкнутой инвариантной кривой (тора). Он может быть даже у двумерных отображений, собственно для которых и приведена иллюстрация (рис. 15). Здесь изображена окрестность (кольцо)
0.3
Ху
0.3 0.4 0.5 0.6
Ху
(с)
№
0.8 0.8
0.6 0.6
Ъ 0.4 ■ ■ 0.4 0.2
0.2 0.0
0.0
0.0 0.1 0.2 0.3 0.4 0.5
Ху
0.6 0.7 0.8
Ху
(е
(!)
Рис. 12. Эволюция аттракторов в 3DHM (1.2) с В = 0.7, М2 = 0.85 при изменении М\: (а) Ы\ = = —0.055 (неподвижная точка); (Ь) М1 = —0.045 (цикл периода 2); (с) М1 = —0.0342 (замкнутая инвариантная кривая периода 2); (^ М1 = —0.0341; (е) М1 = —0.045; (1) М1 = 0. В последних трех случаях наблюдается С А (лоренцевского типа).
резонансной инвариантном кривои, содержащей две неподвижные точки: устоичивыи узел N и седло 5. У последнего на рисунке показаны две половинки с устойчивым многообразием, которые на самом деле отождествлены. Три колонки на рисунке отвечают трем разным типичным способам разрушения инвариантной кривой.
(а)
3-2 0.4
-0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0
Ху -
(Ь)
Ху
(с
-0.8 -
0.0 0.4
Х1
№
Рис. 13. Эволюция аттракторов в 3DHM (1.2) с В = 0.7, М2 = 0.85 при изменении М1: (а) М1 = = 0.025; (Ь) М1 = 0.05; (с) М1 = 0.075; (^ М1 = 0.1. Во всех этих случаях наблюдается СА лоренцевского типа, который, однако, эволюционирует к «СА без дырок», показанному на рис. (^. Похоже, что здесь образуется нечто типа «быковских контуров» [33, 34], за счет появления пересечения между одномерными неустойчивыми многообразиями точек странного аттрактора и одномерным устойчивым многообразием седло-фокусного цикла {р1,р2} периода 2 (см. рис. 10), лежащего в «дырках».
В общем случае механизм Афраймовича-Шильникова начинается с того, что на инвариантной кривой Ьр появляются (резонансные) седла 5 и узлы N одного и того же периода. При этом поначалу инвариантная кривая (здесь это замыкание неустойчивых многообразий седел) является по-прежнему гладкой (но гладкость ее уже конечная и, вообще говоря, зависит от соотношений между мультипликаторами узлов13).
13Поясним это на примере линейного отображения х = ь>\х, у = 1/2у, где O<V2<V1<I. Точка 0(0, 0) является устойчивым узлом, у которого ось «у», х = 0, является сильно устойчивым инвариантным многообразием, а ось «х», у = 0, определяет ведущее направление. Нетрудно проверить,
что кривые вида у = Сха, где а = -7——^, являются инвариантными, то есть любая точка на такой
1п v1
кривой остается на ней при итерациях отображения. Это означает, что любая траектория (кроме траекторий тех точек, которые лежат на Ш88: {х = 0}) входит в О, касаясь ведущего направления у = 0. Соответственно, любая инвариантная кривая, входящая в О, будет лежать целиком между двумя какими-либо кривыми у = С1ха и у = С2ха, то есть если она Сг-гладкая, то ее гладкость не должна типично превышать целую часть числа а.
0.4575
0.4550
0.4525
0.4500
0.4475
”0.4450
0.4425
0.4400
0.4375
0.4450 0.4475 0.4500 0.4525 0.4550 Ху
0.8 0.7 0.6 0.5 0.4 0.3 ■ 0.2 -0.1
(а)
(Ф
0.4 0.5
XI
(8)
0.60
0.55
0.50
0.45 Хп 0.40
0.35
0.30
0.25
0.30 0.35 0.40 0.45 0.50 0.55 0.6 Ху
0.7 Г ~ \ \ 0.7
0.6 0.6
0.5 ~Г\ 0.5
Ь
0.4 \ \ \ \ 0.4
0.3 \ К 0.3
0.2 0.2
0.2 0.3 0.4 0.5 0.6 0.7 Ху
0.9
0.8
0.7
0.6
ю О О
0.3
0.2
0.1
0.0
(Ь)
0.2 0.3 0.4 0.5 0.1
Ху
(е)
I ’ ч чЛ. \
0.65
0.60
0.55
0.50
_ 0.45 Хп 0.40
0.35
0.30
0.25
.6 0.7
0.1 0.2 0.3
0.4 0.5 0.6 0.7 0.8
Ху
(Ь)
0.30 0.35 0.40 0.45 0.50 0.55 0.60 0.65 Ху
0.7 { 0.6 -0.5 0.4 ■ 0.3 ■ 0.2
(с)
0.2 0.3 0.4 0.5 0.6 0.7
Ху
(!)
Рис. 14. Эволюция аттракторов в 3DHM (1.2) при В = 0.7, М2 = 0.75: (а) М1 = —0.025 (неподвижная точка); (Ь) М1 = 0.025 (цикл периода 2); (с) М1 = 0.035 и (^ М1 = 0.039 (замкнутая инвариантная кривая периода 2); (е) М1 = 0.0395 — момент близкий к бифуркационному, когда замкнутая инвариантная кривая периода 2 «влипает в гомоклинику» и превращается в кривую периода 1; (!) М1 = 0.4 и (8) М1 = 0.05 (замкнутая инвариантная кривая периода 1); (Ь) М1 = 0.06 (СА — лоренцевского типа).
Однако, при изменении параметра, кривая Ьр может стать негладкой, поскольку инвариантное неустойчивое многообразие седла может идти в узел, не касаясь определенного направления. Это понятно во втором и третьем случаях (средняя и правая колонки рис. 15), когда Ши(5) пересекает Ш) и когда точка N становится фокусом. В первом случае (левая колонка рис. 15) требуется некоторое пояснение (см. рис. 16). На рисунке 16 пунктирными линиями изображены границы так называемых фундаментальных областей. Мы рассматриваем области, границы которых пересекают ровно по одной точке инвариантные кривые, которые гладко входят в узел N. Если инвариантное многообразие Ши(5) также пересекает каждую из этих границ в одной точке, то Ши(5) будет гладкой кривой
V5' ‘"■■к... - 5І
: т - - 1
(тор гладкий — устойчивое многообразие седла (<9) входит в узел (ЛГ) гладко)
Случай 1
(тор негладкий)
(бифуркационный момент — когда седло-узел(57У) исчезает, тор разрушается)
д
(хаос может возникать сразу после исчезновения седло-узла (если выполнено, например, условие «большой петли»[35]))
Случай 2
(тор негладкий)
(бифуркационный момент — образование гомоклиники седла 5, тора нет)
д
(СА появляется сразу после исчезновения седло-узла, и здесь наблюдается так называемое явление «перемежаемости»)
Случай 3 (тор негладкий)
(после бифуркации удвоения периода устойчивой точки тора нет)
д
(аттрактором является устойчивый цикл периода 2, СА появляется только после новых бифуркаций удвоения периода или бифуркаций, проходящих по первым двум сценариям)
Рис. 15. К механизму Афраймовича-Шильникова возникновения тор-хаоса.
(рис. 16а). Если же таких точек пересечения 3 (или больше), то, очевидно, кривая Ши(Б) теряет гладкость в точке N (рис. 16Ь и с).
В первом случае инвариантная кривая Ши(Б) существует (постепенно теряя гладкость) вплоть до того бифуркационного момента, когда седло сливается с узлом. Здесь сразу же после исчезновения седло-узла может возникать СА. Однако для этого требуется, чтобы выполнялись некоторые дополнительные условия, например, условие «большой петли» из [35] (см. рис. 16с). Заметим, что некоторые достаточные условия, при выполнении которых сразу после исчезновения седло-узла возникает СА, установлены, например, в [36].
Во втором случае инвариантная кривая Ши(Б) пересекает неведущее (сильно устойчивое) многообразие узла. Здесь сразу же тор Ши(Б) негладкий, и он разрушается, как только образуется гомоклиническое касание у седла Б (см. рис. 15, второй столбец). Однако для появления СА нужно, чтобы произошла седло-узловая бифуркация, как и в первом
случае. При этом здесь возникает сразу же развитый хаос и наблюдается динамическое
14
явление, называемое перемежаемостью.
14Иногда говорят о «перемежаемости по Помо-Манневилю», которые дали скорее физическую интерпретацию этого явления. Однако математическая его теория, в основе которой лежит бифуркация исчезновения седло-узловой периодической траектории, имеющей гомоклиническую орбиту,
(с) «большая петля»
Рис. 16. К вопросу о гладкости инвариантного тора.
В третьем случае узел сначала становится фокусом и тор автоматически теряет гладкость. При этом он разрушается, когда, например, происходит бифуркация удвоения периода устойчивой точки. Однако для появления СА нужно, чтобы произошла еще некоторая серия бифуркаций уже с устойчивой точкой периода 2. Например, здесь может иметь место бесконечная последовательность бифуркаций удвоения периода или бифуркации по типу случаев 1 и 2 механизма Афраймовича-Шильникова.
5.2. О возникновении тор—хаоса в результате удвоений инвариантных кривых
Предположим здесь, что существует такое ^\* > ц1, при котором опять меняется гладкий тип кривой Ьр — она становится при ц > ^1* снова узлового типа, но с прыгающими направлениями (для Т2 это обычная узлового типа инвариантная кривая). Здесь есть некоторая аналогия с периодической траекторией, имеющей отрицательные действитель-
была установлена Лукьяновым и Шильниковым [31]. Сам термин перемежаемость (Ш;егт1Мепсу) возник как отражение специфики характера временных реализаций хаотических колебаний, возникающих здесь: участки явно нерегулярных, случайных колебаний чередуются (перемежаются) с участками колебаний, более или менее похожих на регулярные. Это легко объясняется, так как в «нерегулярной фазе» траектория находится в области развитого хаоса, а в «регулярной фазе» — вблизи бывшего седло-узла.
ные мультипликаторы. Аналогия простирается даже дальше: такая замкнутая инвариантная кривая может потерять свою устойчивость в результате бифуркации удвоения — сама кривая становится седловой, а в ее окрестности появляется устойчивая замкнутая инвариантная кривая, дважды обвивающая исходную. Очевидно, такие бифуркации не могут быть у двумерных диффеоморфизмов, но в трехмерном случае они вполне возможны. Более того, с новой кривой может произойти опять бифуркация удвоения и т. д. Может ли этот процесс удвоений идти до бесконечности — неясно. Однако несколько таких бифуркаций можно наблюдать при численных экспериментах (см. рис. 5Ь-£).
5.3. О возникновении хаоса в результате жесткой потери устойчивости инвариантной кривой
Как известно (и это показывают, в частности, численные эксперименты), с инвариантными замкнутыми кривыми могут происходить бифуркационные явления, очень похожие на те, которые наблюдаются у предельных циклов. Конечно, теория бифуркаций замкнутых инвариантными кривых гораздо сложнее, однако если не вдаваться в определенные тонкие детали, то аналогия покажется весьма значительной. Так, известные бифуркация (Андронова-Хопфа) рождения цикла из сложного фокуса, бифуркация удвоения периода предельного цикла, бифуркации слияния и исчезновения двух предельных циклов, а также бифуркации предельных циклов при прохождении через петлю седла имеют свои аналоги и в случае замкнутых инвариантных кривых. Первые два типа бифуркаций мы обсуждали выше и, кроме того, приводили результаты численного счета, которые иллюстрируют эти бифуркации. В этой статье мы не останавливаемся подробно на описании бифуркаций последних двух типов. Однако такие «живые» бифуркации можно наблюдать на рисунке 14. на котором видно как две устойчивые замкнутые инвариантные кривые сливаются в одну, пройдя через «гомоклиническую восьмерку» седла Ор (рис. d-g), и как получившаяся инвариантная кривая (рис. g) жестко теряет устойчивость,15 так что сразу же проявляется хаос (образовавшийся ранее за счет гомоклиник седла Ор).
Авторы благодарят А. Борисова, И. Мамаева и Д. Тураева за весьма полезные замечания.
Список литературы
[1] Шильников Л. П. Об одном случае существования счетного множества периодических движений // Докл. АН СССР, 1965, т. 160, №3, с. 558-561.
[2] Шильников Л. П. К вопросу о структуре расширенной окрестности грубого состояния равновесия типа седло-фокус // Матем. сб., 1970, т. 81, №1, с. 92-103.
[3] Шильников Л. П. Теория бифуркаций и турбулентность // Методы качественной теории дифференциальных уравнений: Межвуз. сб. / Е. А. Леонтович-Андронова. Горький: ГГУ, 1986. С.150-163.
[4] Афраймович В. С., Быков В. В., Шильников Л. П. О возникновении и структуре аттракторов Лоренца // Докл. АН СССР, 1977, т. 234, №2, с. 336-339.
[5] Афраймович В. С., Быков В. В., Шильников Л. П. О притягивающих негрубых предельных множествах типа аттрактора Лоренца // Тр. ММО, 1982, т. 44, с. 150-212.
15Типично здесь возможны два варианта: «седло-узловая бифуркация» — устойчивая кривая сливается с такой же седловой и исчезает, и «обратное удвоение» — в устойчивую кривую влипает седловая инвариантная кривая двойной длины, передавая первой свою неустойчивость.
[6] Шильников Л. П. Теория бифуркаций и модель Лоренца // Бифуркация рождения цикла и ее приложения / Дж. Марсден, М. Мак-Кракен. М.: Мир, 1980. С. 317-335.
[7] Shilnikov L. P. Mathematical problems of nonlinear dynamics: A Tutoraial // J. Franklin Inst. B, 1997, vol. 334, nos. 5-6, pp. 793-864.
[8] Афраймович В. С., Шильников Л. П. Инвариантные торы, их разрушение и стохастичность // Методы качественной теории дифференциальных уравнений: Межвуз. тематич. сб. науч. тр. / Е. А. Леонтович-Андронова. Горький: ГГУ, 1983. С. 3-26.
[9] Benedicks M., Carleson L. The dynamics of the Henon map // Ann. of Math (2), 1991, vol. 133, pp.73-169.
[10] Тураев Д. В., Шильников Л. П. Пример дикого странного аттрактора // Матем. сб., 1998, т. 189, №2, c. 137-160.
[11] Тураев Д. В., Шильников Л. П. Псевдогиперболичность и задача о периодических возмущениях аттракторов лоренцевского типа // Докл. РАН, 2008, т. 418, №1, с. 23-27.
[12] Гонченко С. В., Тураев Д. В., Шильников Л. П. Динамические явления в многомерных системах с негрубой гомоклинической кривой Пуанкаре // Докл. РАН, 1993, т. 330, №2, с. 144-147.
[13] Gonchenko S.V., Shilnikov L.P., Turaev D. V. Dynamical phenomena in systems with structurally unstable Poincare homoclinic orbits // Chaos, 1996, vol. 6, no. 1, pp. 15-31.
[14] Newhouse S.E. The abundance of wild hyperbolic sets and nonsmooth stable sets for diffeomorphisms // Inst. Hautes Etudes Sci. Publ. Math., 1979, vol. 50, pp. 101-151.
[15] Гонченко С. В., Тураев Д. В., Шильников Л. П. О динамических свойствах диффеоморфизмов с гомоклиническими касаниями // Итоги науки и техн. Сер. Соврем. мат. и ее прил. Темат. обз., 2003, т. 7, с. 92-118.
[16] Newhouse S.E. Diffeomorphisms with infinitely many sinks // Topology, 1974, vol. 13, pp. 9-18.
[17] Aframovich V. S., Shilnikov L. P. Strange attractors and quasiattractors // Nonlinear Dynamics and Turbulence / G. I. Barenblatt, G. Iooss, D. D. Joseph. Boston: Pitmen, 1983. P. 1-34.
[18] Gonchenko S. V., Ovsyannikov 1.1., Simo C., Turaev D. V. Three-dimensional Henon-like maps and wild Lorenz-like attractors // Internat. J. Bifur. Chaos Appl. Sci. Engrg., 2005, vol. 15, pp. 34933508.
[19] Морозов А.Д., Драгунов Т.Н. Визуализация и анализ инвариантных множеств динамических систем. М.-Ижевск: НИЦ «Регулярная и хаотическая динамика», Инст. компьютерн. исслед., 2003. 304 с.
[20] Morozov A. D., Dragunov T. N., Boykova S. A., Malysheva O. V. Invariant sets for Windows. (World Sci. Ser. Nonlinear Sci. Ser. A Monogr. Treatises, vol. 37.) Singapore: World Sci., 1999. 259 p.
[21] Lomeli H. E., Meiss J. D. Quadratic volume preserving maps // Nonlinearity, 1998, vol. 11, pp. 557574.
[22] Gonchenko S. V., Meiss J.D., Ovsyannikov 1.1. Chaotic dynamics of three-dimensional Henon maps that originate from a homoclinic bifurcation // Regul. Chaotic Dyn., 2006, vol. 11, no. 2, pp. 191-212.
[23] Gonchenko S.V., Shilnikov L.P., Turaev D.V. On dynamical properties of multidimensional diffeomorphisms from Newhouse regions: 1 // Nonlinearity, 2008, vol. 21, pp. 923-972.
[24] Gonchenko S. V., Gonchenko V. S., Shilnikov L. P. On homoclinic origin of Henon-like maps // Regul. Chaotic Dyn., 2010, vol.15, nos. 4-5, pp. 462-481.
[25] Gonchenko S. V., Shilnikov L. P., Turaev D. V. On dynamical properties of multidimensional diffeomorphisms from Newhouse regions: 2 // Nonlinearity (to appear).
[26] Борисов А. В., Мамаев И. С. Странные аттракторы в динамике кельтских камней // УФН, 2003, т. 173, №4, с. 407-418.
[27] Борисов А. В., Мамаев И. С. Странные аттракторы в динамике кельтских камней // Неголоном-ные динамические системы: Интегрируемость, хаос, странные аттракторы: Сб. ст. / А. В. Борисов, И. С. Мамаев (ред.). М.-Ижевск: Инст. компьютерн. исслед., 2002. С. 293-316.
[28] Гонченко С. В., Тураев Д. В., Шильников Л. П. Об областях Ньюхауса двумерных диффеоморфизмов, близких к диффеоморфизму с негрубым гетероклиническим контуром // Тр. МИАН, 1997, т. 216, с. 76-125.
[29] Lamb J.S.W., Stenkin O.V. Newhouse regions for reversible systems with infinitely many stable, unstable and elliptic periodic orbits // Nonlinearity, 2004, vol. 17, no. 4, pp. 1217-1244.
[30] Delshams A., Gonchenko S.V., Gonchenko V. S., Lazaro J.T., Stenkin O.V. Abundance of attracting, repelling and elliptic periodic orbits in two-dimensional reversible maps // Nonlinearity (to appear)
[31] Лукьянов В. И., Шильников Л. П. О некоторых бифуркациях динамических систем с гомокли-ническими структурами // Докл. АН СССР, 1978, т. 243, №1, с. 26-29.
[32] Арнольд В. И., Афраймович В. С., Ильяшенко Ю. С., Шильников Л. П. Теория бифуркаций // Динамические системы — 5 / В. И. Арнольд. (Итоги науки и техн. Сер. Соврем. пробл. мат. Фундам. направления, т. 5.) М.: ВИНИТИ, 1986. C. 9-218.
[33] Быков В. В. О бифуркациях динамических систем, близких к системам с сепаратрисным контуром, содержащим седло-фокус // Методы качественной теории дифференциальных уравнений: Межвуз. темат. сб. научн. тр. Горький: ГГУ, 1980. С. 44-72.
[34] Bykov V. V. The bifurcations of separatrix contours and chaos // Phys. D, 1993, vol. 62, pp. 290-299.
[35] Афраймович В. С., Шильников Л. П. О некоторых глобальных бифуркациях, связанных с исчезновением неподвижной точки типа седло-узел // Докл. АН СССР, 1974, т. 219, №6, с. 12811285.
[36] Тураев Д. В., Шильников Л. П. Бифуркации квазиаттракторов тор-хаос // Математические механизмы турбулентности. Киев, 1986. С. 113-121.
[37] Анищенко В. С. Сложные колебания в простых системах. М.: Наука, 1990. 312 с.
Towards scenarios of chaos appearance in three-dimensional maps
Alexander S. Gonchenko1, Sergey V. Gonchenko2, Leonid P. Shilnikov 3
1,2,3Research Institute for Applied Mathematics and Cybernetics Ul’yanova st. 10, Nizhnii Novgorod, 603605, Russia
[email protected], [email protected], [email protected]
We study questions of chaotic dynamics of three-dimensional smooth maps (diffeomorphisms). We show that there exist two main scenarios of chaos developing from a stable fixed point to strange attractors of various types: a spiral attractor, a Lorenz-like strange attractor or a “figure-8” attractor. We give a qualitative description of these attractors and define certain condition when these attractors can be “genuine” ones (pseudo-hyperbolic strange attractors). We include also the corresponding results of numerical analysis of attractors in three-dimensional Henon maps.
MSC 2010: 37C70, 37C29, 37G25
Keywords: strange attractor, chaotic dynamics, spiral attractor, torus-chaos, homoclinic orbit, invariant curve, three-dimensional Henon map
Received October 6, 2011, accepted December 24, 2011
Citation: Rus. J. Nonlin. Dyn., 2012, vol. 8, no. 1, pp. 3-28 (Russian)
Leonid P. Shilnikov