Синхронизация и формирование структур во взаимодействующих системах с локальными связями

Синхронизация колебаний взаимодействующих систем — одно из фундаментальных свойств природы, имеющее широкое применение в различных областях науки и техники. Впервые научное описание этого явления было сделано в 17 веке X. Гюйгенсом [1], который исследовал взаимную подстройку хода двух маятниковых часов, висящих на общей балке. В радиофизике работы по синхронизации начинаются в первой половине 20 века, когда было обнаружено свойство захвата частоты колебаний триодного генератора периодическим сигналом [2,3]. Подробный теоретический анализ данного явления был сделан в работах A.A. Андронова [4,5]. Существенную роль в исследовании задач синхронизации взаимодействующих периодических осцилляторов при наличии исчезающе малых связей сыграли фундаментальные работы А. Пуанкаре [6] и A.M. Ляпунова [7], развитие которых приведено в монографии И. Г. Малкина [8]. Фундаментальные результаты по синхронизации периодических колебаний с точки зрения качественной теории дифференциальных уравнений и теории бифуркаций были сделаны В. И. Арнольдом [9]. В дальнейшем, синхронизация периодических автоколебаний была описана во множестве научных работ, включая ряд монографий [10−12]. Обзор основных результатов, а также библиографию по этому вопросу можно найти в [13−15].

Явление синхронизации имеет множество разнообразных проявлений в природе, технике, экономике и обществе, поэтому ему трудно дать достаточно строгое и полное определение. Наиболее удачным представляется определение, данное в монографии И.И. Блех-мана [10]: «синхронизацию можно определить как свойство материальных объектов самой различной природы вырабатывать единый ритм совместного существования, несмотря на различие индивидуальных ритмов и на подчас крайне слабые взаимные связи». Выработка единого ритма заключается в том, что, при синхронизации:

1) происходит «захват собственных частот» автоколебаний, когда система N осцилляторов, каждый из которых имеет свою индивидуальную частоту tUi (г = 1,2 при наложении связей начинает колебаться с некоторой единой для всех частотой и, Inf {wi, tu2, Wjv} < LO < Sup{uJl, LU2, ., UJNy,.

2) устанавливаются определенные стационарные значения разностей текущих фаз между колебаниями осцилляторов: ipi (t) —.

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

При синхронизации периодических осцилляторов может существовать несколько значений стационарных разностей фаз, соответствующих разным устойчивым синхронным состояниям (в этом случае и устанавливаемые общие частоты и> будут, как правило, также разными). Выбор между сосуществующими синхронными состояниями определяется начальными условиями. Данное явление получило название «фазовой мультистабиль-ности» [17]. Если взаимодействует два осциллятора, то при форме колебаний близкой к гармонической, обычно наблюдается одно синхронное состояние — синфазное или противофазное (хотя и здесь возможны исключения, см, например, [18]). Однако, если форма колебаний начинает усложняться, то число сосуществующих режимов увеличивается. Классическим примером этого является развитие фазовой мультистабильности в системе двух взаимодействующих осцилляторов при субгармоническом каскаде [17,19−26]. Так, при исследовании динамики двух идентичных систем с симметричной связью, каждая из которых демонстрирует переход к хаосу через последовательность бифуркаций удвоения периода, было обнаружено, что в ходе каскада бифуркаций происходит последовательное увеличение числа предельных циклов в два раза: в фазовом пространстве существует один устойчивый однооборотный предельный цикл (1С0) с некоторым исходным периодом Т (одна периодическая орбита периода один в случае дискретных отображений) — два устойчивых двухоборотных предельных цикла (2С0 и 2С1) с периодом 2 Т (две периодические орбиты периода два в случае дискретных отображений), четыре цикла (4С°, 4С1, 4С2 и 4С3) с периодом 4 Т и т. д. На границе перехода к хаосу число сосуществующих синхронных режимов стремится к бесконечности. Для любого из указанных режимов 2NCM (N = 0,1,2,., М = 0,1,., 2N — 1) форма колебаний в каждом из осцилляторов.

1 Здесь и далее будет использована система обозначений мультистабильных состояний, впервые примененная в работе [20] и затем ставшая традиционной при описании фазовой мультистабильности: nRl первый числовой индекс п означает «тактность» режима, то есть число петель аттрактора, или число связных областей в сечении ПуанкареR- тип колебаний (С — цикл, Т — тор, Ахаотический аттрактор) — верхний числовой индекс г — задержку колебаний одного осциллятора относительно колебаний второго в периодах однооборотного цикла. идентична колебаниям одиночной системы, а их временные реализации отличаются на некоторый временной интервал т: хг (£) = Хх (£ — т). Причем, как показали исследования, величина временного сдвига между колебаниями подсистем равна целому числу периодов исходного однооборотного цикла: т = МТ. Подробный анализ этого свойства фазовой мультистабильности и его объяснение с точки зрения захвата Фурье-фаз в спектрах колебаний осцилляторов приведен в работе [25]. Закономерности формирования фазовой мультистабильности, бифуркационный сценарий этого явления, как будет показано далее, имеют непосредственную связь с механизмами формирования и разрушения синхронизации сложных непериодических колебаний в системах с удвоениями периода. Выявлению этой взаимосвязи посвящена первая и, частично, — вторая глава диссертационной работы.

В ансамблях из большего чем два числа осцилляторов сосуществование множества синхронных режимов с разными фазовыми сдвигами наблюдается уже для гармонических колебаний. Одними из характерных режимов, которые могут возникать в цепочках локально связанных осцилляторов с периодическими граничными условиями, являются автоволны, бегущие вдоль ансамбля с постоянной фазовой скоростью [27−30]. При этом, колебания в соседних осцилляторах имеют равную амплитуду и отличающуюся на постоянное значение фазу колебаний. В цепочках конечной длины число возможных режимов бегущих волн конечно. Аттракторы, соответствующие режимам с разной длиной волны, могут сосуществовать в фазовом пространстве при одних и тех же значениях параметров, демонстрируя тем самым явление мультистабильности.

Синхронизация колебаний в ансамблях осцилляторов с регулярным поведением — одна из традиционных областей исследований для радиофизики. Первые работы в этом направлении известны с середины прошлого века [31] и рассматривали, как правило, задачу частотной синхронизации в цепочке осцилляторов с гармоническим поведением [32−35]. В работе [34] было обращено внимание, что в подобных системах возможны режимы с разными фазовыми сдвигами между осцилляторами — то есть сосуществуют разные пространственные моды. Исследование мультистабильных состояний в ансамбле идентичных замкнутых в кольцо осцилляторов подробно описано в монографии П. С. Ланда [36]. Условия существования и устойчивости разных пространственных мод в ансамблях осцилляторов, колебания в которых возникают через бифуркацию Андронова-Хопфа, были получены для разных типов связей в цикле работ Эрментроута [37−39]. Детальное описание динамики пространственно-однородных и неоднородных волн, возникновение пространственно — разупорядоченных колебательных режимов и особенности переходов между режимами с разными длинами волн для ансамбля автогенераторов с жестким возбуждением можно найти в работе [40]. Рассмотрение процессов образования частотных кластеров при различных видах расстройки собственных параметров гармонических осцилляторов было проведено в работах [41−43]. В работах [44−46] было обнаружено явление подавления автоколебаний, когда в цепочке осцилляторов появляются области с пренебрежимо малой амплитудой колебаний. Это явление получило название «вымирание автоколебаний» или «амплитудная смерть» .

Большинство указанных работ используют в качестве модели либо фазовые уравнения, либо укороченные уравнения для амплитуд и фаз, решениями которых являются гармонические колебания. В то же время, представляется интересным вопрос, как ведут себя ансамбли «реальных» регулярных осцилляторов, в которых присутствуют такие явления, как ангармоничность и неизохронность. Одной из попыток ответить на данный вопрос можно считать работу Даидо [47], в которой он моделирует ангармоничность введением дополнительного слагаемого в уравнение для фазового осциллятора. При этом, при большой ангармоничности в кольце осцилляторов наблюдается появление пространственной разупорядоченности (т.н. «странные волны»). Моделирование ангармоничности в работе Даидо носит достаточно искусственный характер. Кроме того, обнаруженные эффекты наблюдаются в его модели при очень больших значениях ангармоничности, не наблюдаемых в более реалистичных моделях. Будут ли подобные режимы возникать в реальных ангармонических осцилляторах? Какие еще особенности поведения там возможны? Как постепенный переход от гармонических колебаний к релаксационным влияет на формирование мультистабильности в цепочках осцилляторов? Пространственно — периодические режимы не исчерпывают всего разнообразия структур, которые могут возникать даже в небольших цепочках. Возникает вопрос — насколько типичными являются эти режимы. Будут ли они возникать при случайных начальных условиях из некоторой окрестности однородного состояния? Какова структура их бассейнов притяжения? Все эти вопросы оставались открытыми на момент начала работы над данной диссертацией. Их рассмотрение проводится в третьей главе диссертации.

Усложнение формы колебаний в ансамбле осцилляторов может быть связано не только с переходом от гармонических колебаний к релаксационным, но и с уже отмеченным ранее каскадом бифуркаций удвоения периода в каждой из подсистем, который ведет к появлению хаотических колебаний в каждом из осцилляторов, то есть к хаотическим бегущим волнам в ансамбле. Возможность существования подобных колебательных режимов рассматривалась в ряде работ [48−51] (в работах [48, 50] они были названы как «rotating waves» — «вращающиеся волны»). Тем не менее, детальный анализ эволюции таких режимов с изменением параметров системы не был сделан. Представляется интересным, каким образом усложнение временной динамики влияет на пространственную упорядоченность колебаний в ансамбле. Сохраняется ли пространственная периодичность вдоль всего каскада субгармонических бифуркаций? На каком этапе при переходе к временному хаосу и как именно происходит разрушение пространственной периодичности и переход к пространственному хаосу? Решение этих вопросов также рассматривается в третьей главе диссертационной работы.

После открытия динамического хаоса, когда было обнаружено, что детерминированные динамические системы могут демонстрировать сложное непредсказуемое поведение, внимание исследователей было обращено на возможность синхронизации взаимодействующих хаотических осцилляторов. Такую синхронизацию стали называть «синхронизацией хаоса» или «хаотической синхронизацией». На сегодняшний день нет единого подхода к определению хаотической синхронизации. Этот термин объединяет несколько разных способов взаимной подстройки колебаний хаотических осцилляторов: полная синхронизация хаоса (ПС) [52−55], частотная синхронизация хаоса (ЧС) [56−58], фазовая синхронизация (ФС) [59,60], обобщенная синхронизация (ОС) [61−63], синхронизация с задержкой по времени (lag-synchronization в англоязычной литературе) [64]. В последнее время появилась попытка рассмотреть разные виды синхронизации с единых позиций, как постепенный взаимозахват фаз вейвлет-образов в разных частотных диапазонах [65−67]. Данный метод анализа получил название «синхронизация временных масштабов» .

Одним из видов хаотической синхронизации является «полная синхронизация хаоса» [52−55], при которой в ансамбле идентичных осцилляторов в результате действия связей колебания всех подсистем становятся одинаковыми: Xi (i) = X2(i). = Этот вид синхронизации подразумевает наличие в фазовом пространстве системы инвариантного симметричного подпространства U: Xi = Х2 =. = хдг, что, в свою очередь, является следствием инвариантности системы уравнений, описывающей ансамбль, к взаимным перестановкам динамических переменных Xj Xj, i, j = 1,2,., N. Прилагательное «полная» отражает тот факт, что при данном виде синхронизации происходит максимально возможное согласование поведения всех осцилляторов ансамбля. Несмотря на достаточно жесткое условие, полная синхронизация — типичное явление, наблюдаемое не только в специально созданных математических моделях, но и в реальных природных, технических и социальных системах [68−75].

Первыми работами по исследованию режимов полной синхронизации хаоса были статьи Ямады и Фуджисаки [52], А. С. Пиковского [53], С. П. Кузнецова [76]. Было установлено, что полная хаотическая синхронизация наблюдается в системе двух идентичных осцилляторов при диффузионной связи. Был определен критерий устойчивости синхронных колебаний по отношению к возмущениям, трансверсальным к симметричному подпространству ij через трансверсальный показатель Ляпунова Aj. Как оказалось, полная синхронизация хаоса реализуется при достаточно сильной связи (коэффициент диффузионной связи 7 должен быть больше некоторого критического значения 7с), где транс-версальный показатель Ляпунова отрицателен. При переходе через критическое значение связи 7с трансверсальный показатель Ляпунова становится положительным, при этом синхронные колебания теряют трансверсальную устойчивость и перестают наблюдаться в эксперименте (данный переход называется бифуркацией прорыва [77]). Граница области синхронизации на плоскости управляющих параметров имеет достаточно сложный характер. В ее окрестности наблюдается перемежаемость Ямада-Фуджисака [78,79] (также называемая «переключающейся синхронизацией» или «оп-о:И ^егтИ^епсу» в англоязычной литературе): интервалы почти синхронного поведения чередуются с короткими «всплесками» несинхронных колебаний.

При более детальном исследовании разрушения полной синхронизации хаоса было установлено [77,80−91], что рассинхронизация колебаний возможна и при отрицательном значении трансверсального показателя Ляпунова. Это связано с наличием в хаотическом аттракторе трансверсально неустойчивых траекторий, с положительными трансверсаль-ными показателями Ляпунова (данное свойство получило название «локальная изрешечен-ность хаотического аттрактора»), из окрестности которых фазовая точка, при наличии сколь угодно малых несимметричных возмущений, выбрасывается в сторону от симметричного подпространства. При наличии локальной изрешеченности хаотический аттрактор перестает быть аттрактором в традиционном понимании и становится т.н. аттрактором Милнора [92,93]2 Локальная изрешеченность ведет к двум возможным последствиям в наблюдаемой динамике системы:

• если, области трансверсальной неустойчивости аттрактора направлены в сторону седловых предельных траекторий, то фазовая точка, уйдя от симметричного подпространства вдоль этих областей, через некоторое время вследствие глобального притяжения к аттрактору вернется в его окрестность. Поскольку хаотический аттрактор в среднем остается трансверсально притягивающим, то после каждого ухода она будет все ближе и ближе подходить к подпространству /?, при этом следующие «выбросы» будут происходить все реже и реже. В итоге, режим хаотической синхронизации оказывается асимптотически орбитально устойчивым. Однако, добавление в систему малейшего шума или малейшая неидентичность взаимодействующих осцилляторов делают режим перемежаемости перманентным. Таким образом, полная синхронизация продолжает существовать, но перестает быть грубой. Такой характер синхронизации получил название «слабой синхронизации», а поведение системы.

2Аттрактор в смысле Милнора — это наименьшее замкнутое множество, имеющее область притяжения ненулевой меры. В отличие от классического определения, здесь не требуется, чтобы любая точка прикосновения аттрактора входила в его область притяжения. в этом режиме — «пузырящимся поведением» (перекрестная проекция фазового аттрактора как бы «вскипает» под действием малого шума).

• Если области трансверсальной неустойчивости направлены в сторону других аттракторов, то фазовая точка, уйдя по ним от симметричного подпространства, притягивается к другому аттрактору. При этом, режим хаотической синхронизации заменяется другим колебательным режимом. Поскольку области неустойчивости исходят из самого аттрактора, в его непосредственной окрестности оказываются области из бассейнов других аттракторов. Данное явление называется «изрешеченностью бассейна притяжения хаотического аттрактора» .

Как было показано в ряде работ [86,94], потеря синхронизации непосредственно связана с бифуркациями седловых циклов, встроенных в хаотический аттрактор. Как известно, в хаотический аттрактор может быть встроено бесконечное число седловых циклов. В связи с этим, представляется интересным вопрос: какие именно циклы являются определяющими при потере хаотической синхронизации. Этот вопрос рассматривался в работе [95], где авторы обнаружили, что потеря трансверсальной устойчивости начинается, как правило, с бифуркаций седловых циклов малого периода. Однако, детального исследования бифуркационного сценария процессов разрушения полной синхронизации сделано не было. Каковы типичные бифуркационные переходы, ведущие к постепенной потере хаотической синхронизации? Как они связаны с бифуркациями, ответственными за формирование синхронного хаоса и с бифуркациями, ответственными за формирование фазовой муль-тистабильности? Как меняется бифуркационный сценарий и наблюдаемые явления (если меняются) при переходе к неидентичным системам? Существуют ли различия в полной синхронизации для обратимых и необратимых систем? Эти вопросы рассматриваются в первой главе диссертационной работы на примере связанных осцилляторов, переход к хаосу в которых происходит через каскад бифуркаций удвоения периода.

Естественным обобщением понятия полной синхронизации стала «обобщенная синхронизация» («generalized synchronization»), при которой между сигналами, генерируемыми подсистемами, существует некоторая детерминированная взаимосвязь: x2(i) = F (Xl (i)).

Первое определение обобщенной синхронизации было дано в работе Афраймовича [54]. Открытая в середине 90-х годов в работах Рулькова и Кокарева [61−63, 96], она затем была обнаружена во множестве различных динамических систем [97−103], включая синхронизацию распределенных сред [67,104] и систем с задержкой по времени [106]. Детерминированная функция^Р (ж), связывающая состояния подсистем, может иметь достаточно сложный характер и даже быть недифференцируемой [107], что подразумевает нетривиальные методы диагностирования данного вида синхронизации. Обобщенная синхронизация наблюдается, как правило, в связанных неидентичных осцилляторах. Однако, как будет показано в третьей главе, режимы обобщенной синхронизации могут существовать и в ансамблях идентичных осцилляторов с симметричной диффузионной связью.

Полная синхронизация, как это видно из ее определения, является частным случаем обобщенной синхронизации, если в качестве детерминированной функции взять тривиальную зависимость: Г (х) = х. Другим частным случаем обобщенной синхронизации, получившим, однако, свое собственное название, является «противофазная синхронизация хаоса» [108−112], при которой временные реализации, генерируемые подсистемами, удовлетворяют соотношению XI (?) = —х2(£). Соответственно, синхронные траектории лежат в симметричном инвариантном подпространстве 1а: = — х23.

Противофазная синхронизация — типичное явление для регулярных осцилляторов. Достаточно вспомнить, что именно это явление было впервые описано X. Гюйгенсом, как пример синхронизации. Однако, для хаотических осцилляторов противофазная синхронизация гораздо более редкое явление, механизмы которого оставались не изученными на момент создания данной работы. Как происходит формирование хаотического предельного множества, соответствующего режиму противофазной синхронизации? Каков сценарий и бифуркационный механизм потери синхронизации? Сходен ли он с соответствующим сценарием для полной синхронизации хаоса или имеет свои особенности?

Как показывают исследования, противофазная синхронизация хаоса не реализуется в системах с диффузионным типом связи, поскольку соответствующие предельные множества трансверсально неустойчивы. Однако, для реализации таких режимов можно использовать методы принудительной синхронизации (или методы синхронизации посредством управления). Задача принудительной синхронизации колебаний — есть задача стабилизации траекторий в симметричном подпространстве по отношению к трансверсальным к этому подпространству возмущениям. Метод решения этой задачи — добавление в систему цепи дополнительной обратной связи, либо внешнего, явным образом зависящего от времени воздействия, которые меняют характер устойчивости синхронных колебаний. При этом управляющее воздействие не должно оказывать влияния на динамику системы внутри симметричного подпространства, то есть не должно менять форму синхронных.

3Существует и другое определение для противофазной хаотической синхронизации (см., например, [113]), которое ближе к понятию фазовой синхронизации: противофазной синхронизацией считается такой хаотический режим, при котором текущие фазы оказываются захвачены около значения Д<�р = ж, или, что тоже самое, колебания в одном из осцилляторов отстают от колебаний во втором на половину среднего периода. В настоящей работе указанный вид синхронизации рассматриваться не будет. колебаний. В противном случае мы имеем дело не с синхронизацией (в том смысле, в котором она определяется в работе), а с индуцированными в системе новыми колебательными режимами при воздействии на нее. Известно, что устойчивые и грубые режимы полной синхронизации хаоса во взаимодействующих системах могут быть реализованы только при определенных типах связи выше некоторого порогового значения. Часто независимо от величины коэффициента связи в системе существуют синхронные хаотические движения, которые являются неустойчивыми к несимметричным возмущениям. В этих случаях в системе можно осуществить переход из режима несинхронных хаотических колебаний к режиму синхронизации, используя методы управления хаосом [114−124]. Под управлением хаосом обычно понимают целенаправленное воздействие на систему, с помощью которого различные непритягивающие предельные множества можно превратить в устойчивые по определенным собственным направлениям. В большинстве работ, посвященных синхронизации посредством управления хаосом, рассматривается простейший случай синхронизации — синфазная синхронизация [125−135]. Возможно также распространить методы принудительной синхронизации хаоса для случая противофазной синхронизации. Синхронизация и принудительная синхронизация противофазных колебаний рассматривается во второй главе диссертационной работы.

В качестве метода управления хаосом можно рассматривать предложенный в ряде работ [136−139] подход, предусматривающий периодическую модуляцию одного из параметров системы. Аналитически и методами численного и физического экспериментов в этих работах было показано, что периодическая модуляция параметра может привести к подавлению хаоса и к переходу на периодический режим, на базе которого возник исходный хаотический аттрактор. Указанный эффект оказывается возможен, если частота модуляции кратна частоте предельного цикла, то есть данное явление носит резонансный характер. В другой работе для подавления хаоса и перехода на регулярный режим использовалась высокочастотная модуляция параметра системы, когда частота модуляции много больше собственной частоты осциллятора [140]. Было показано, что движение системы с высокочастотным воздействием может быть представлено как сумма «медленного» движения с характерной частотой системы без модуляции и «быстрого» движения с характерной частотой параметрического воздействия. Уравнение для полной системы разделяется на уравнение для «быстрых» и для «медленных» переменных, причем параметры уравнения для «быстрых» переменных оказываются зависящими от амплитуды и частоты высокочастотного воздействия. Таким образом, можно говорить о том, что высокочастотная модуляция параметра может менять средние значения параметров системы и таким образом индуцировать переход к другим колебательным режимам. Можно ли этот способ управления использовать для принудительной синхронизации хаотических колебаний? Идея использования параметрического воздействия для синхронизации хаотических систем связана с хорошо известной классической задачей о маятнике с вибрирующей точкой подвеса [141,142]. В ней было показано, что высокочастотная модуляция параметра может изменить устойчивость состояния равновесия. Следуя этой идее в работах [143,144] было показано, что периодическая модуляция параметра связи приводит к синхронизации двух диффузионно связанных осцилляторов Дуффинга в хаотическом режиме. Возможно ли этот метод использовать для синхронизации в больших ансамблях осцилляторов? Есть ли ограничения на число взаимодействующих осцилляторов в ансамбле? Данные вопросы рассматриваются в четвертой главе диссертации.

Как известно, в ансамблях из большего чем два числа элементов, наряду с полной синхронизацией может наблюдаться кластерная синхронизация, при которой существуют наборы осцилляторов (кластеры), работающие в режиме полной синхронизации, между которыми полная синхронизация отсутствует. Если число осцилляторов в ансамбле невелико, то вместо термина «кластерная синхронизация» принято использовать термин «частичная синхронизация» [145] (термин «частичная синхронизация» представляется не совсем удачным, поскольку, между «несинхронными» осцилляторами может существовать другой тип синхронизации). Теоретические исследования режимов кластерной синфазной и кластерной противофазной синхронизации колебаний были проведены в работе [109], где были установлены конфигурации инвариантных множеств, соответствующих разным случаям частичной синхронизации при разном числе осцилляторов в системе. В работе [146] методом Ляпунова были найдены условия асимптотической устойчивости для частично синхронных режимов. Простейшая модель для систем, в которых возможно наблюдать явление частичной синхронизации — кольцо из трех осцилляторов. На ее примере проще всего рассмотреть сосуществование режимов полной и частичной синхронизации, закономерности их появления и разрушения и соответствующие бифуркационные механизмы. «Сильная» и «слабая» (то есть грубая и негрубая) частичная синхронизация хаоса в системах с различными топологиями рассматривалась в работе [147]. В работе [148] исследовался режим частичной синхронизации в системе из трех и четырех связанных осцилляторов Рёсслера. Тот же подход был применен к системе из большего числа осцилляторов, в результате была выяснена зависимость числа элементов системы, совершающих синхронные колебания, от величины связи между элементами.

Явление частичной синхронизации в системе, состоящей из трех связанных логистических отображений с несимметричной связью, исследовалась в работе [149], в которой был проведен детальный бифуркационный анализ механизмов разрушения полной и формирования частичной хаотической синхронизации. Однако для систем с симметричной диффузионной связью такая задача решена не была. Возникновения и исчезновение режимов частичной синхронизации регулярных и хаотических колебаний в кольце из трех отображений с симметричной диффузионной связью рассматривалось в работе [150]. В ней в качестве основного метода анализа использовался расчет трансверсальных показателей Ляпунова. Подобный подход дает возможность оценить область устойчивости для некоторых синхронных режимов на плоскости управляющих параметров, но не позволяет выявить бифуркационные механизмы их формирования, а также не может дать достаточно полной картины устройства пространства параметров при исследовании систем с развитой мультистабильностью. Чтобы ответить на эти вопросы необходимо проведение детального бифуркационного анализа седловых периодических орбит, формирующих скелет синхронного хаотического аттрактора. Проведение такого анализа, включая сопоставление механизмов формирования и разрушения режимов частичной и полной синхронизации, представляется интересной задачей, важной для понимания процессов синхронизации хаоса в больших ансамблях локально связанных элементов.

Наличие разных видов синхронного поведения хаотических осцилляторов, а также тот факт, что выход из режима синхронизации хаоса происходит постепенно, через частично синхронное состояние к полностью несинхронному, делает актуальным задачу диагностики и количественного анализа хаотической синхронизации. Для каждого из видов синхронизации вводились свои количественные меры. Легче всего ввести меру полной синхронизации или синхронизации с задержкой, как степень «схожести» временных реализаций взаимодействующих осцилляторов. В работе [64] для этих целей предлагалось использовать минимальное значение функции подобия S®.

В случае полной синхронизации минимум функции 5(т) достигается при г = 0, при синхронизации с задержкой — при некотором положительном т. В обоих случаях свидетельством полностью синхронного поведения является равенство нулю индекса синхронизации Постепенное разрушение синхронизации, например, через процесс перемежающейся синхронизации, будет выражаться в монотонном росте ?1. Функция подобия естественным образом может быть использована и для количественного анализа обобщенной синхронизации. Однако, в этом случае необходимо знать в явном виде функцию детерминированный связи между колебаниями подсистем.

В случае фазовой синхронизации, когда «амплитуды» колебаний оказываются некорц = min (5®) где a-2(i + r)-a-i (i))2) релированными, а «фазы» — захвачены, можно ввести формальное определение в виде неравенства:

Ит М) ~ Mt) < М (1) t—"00 где мгновенные фазы первого и второго осцилляторов, а М некоторая положительная константа. Формула (1) означает, что разность фаз двух осцилляторов остается ограниченной во времени, то есть, что средние «периоды» колебаний одинаковы. Слабым местом указанного определения является само понятие мгновенной фазы ip (t) в применении к хаотическим колебаниям. Если в случае периодических колебаний, когда аттрактор одномерен, фаза однозначно характеризует установившееся состояние системы (что собственно и является основанием для использования термина «фаза» как синонима термина «состояние»), то уже при квазипериодических колебаниях можно говорить о нескольких независимых фазах (например, фаза несущей и фаза огибающей). Что касается хаотических колебаний, то хотя формальное введение мгновенной фазы тем или иным способом возможно, ее физическая интерпретация остается неясной. Поэтому, на сегодняшний день концепция мгновенной фазы хорошо работает лишь для «слабого» то есть когерентного хаоса, который близок к периодическим колебаниям и не работает для развитого хаоса, при котором спектр не содержит выраженных пиков.

Естественной количественной мерой для фазовой синхронизации может служить коэффициент эффективной диффузии разности мгновенных фаз [151]:

Deff = ~ [<М*) — Мф — - Vi (*))>2].

Эта величина показывает насколько быстро происходит диффузия разности фаз взаимодействующих осцилляторов и связана, таким образом, со средним временем, в течение которого имеет место захват фаз. Для полной фазовой синхронизации Д.// = 0. При разрушении фазовой синхронизации, когда длительные интервалы синхронного поведения прерываются внезапными «проскальзываниями» фаз, Deff начинает монотонно расти по мере ухода от режима синхронизации. В работе [152] было предложено использовать для анализа фазовой синхронизации «параметр порядка» Я, введенный Курамото [153] для определения степени упорядоченности ансамбля фазовых осцилляторов, связанных через среднее поле:

TV—1 ехрО’Д^(Ш)) i=о где Aip (iAt) представляет собой разность текущих фаз колебаний подсистем в момент времени iAt.

Указанные количественные характеристики хаотической синхронизации являются част.

R = N ными, поскольку применимы лишь для определенного вида синхронизации. Кроме того, получаемые из них величины имеют относительный характер: их конкретные значения ничего не скажут исследователю, насколько сильной или слабой является синхронизация в каждом рассматриваемом случае. Имеет смысл лишь сопоставление значений одной и той же величины на разных этапах перехода от синхронного состояния к несинхронному. Хотелось бы иметь универсальную количественную меру синхронизации, позволяющую сравнивать степень согласованности движений для разных проявлений взаимоподстройки осцилляторов. Разумно потребовать, чтобы эта количественная характеристика удовлетворяла следующим критериям:

1) Универсальность, то есть применимость к разным видам синхронного поведения.

2) Ясный физический смысл, позволяющий легко интерпретировать получаемые результаты.

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

4) Устойчивость по отношению к действию малых возмущающих факторов, таких как внешний шум или малые искажения генерируемых сигналов.

5) Кроме того, количественная характеристика должна представлять собой вещественное число, принимающее значение в заданном интервале, например, от нуля до единицы. Причем выбранный интервал определения должен быть единым для всех рассматриваемых динамических систем и видов синхронного поведения.

Первыми работами, рассматривающими задачу количественной меры хаотической синхронизации, можно, по-видимому, считать работы П. С. Ланды и М. Розенблюма [154,155], в которых для этой цели предлагалось использовать размерность аттрактора (речь шла о корреляционной размерности, однако, предполагалось, что можно применять и другие определения размерности хаотического аттрактора). В частности, полностью синхронными полагались такие колебания, для которых размерность аттрактора в полном фазовом пространстве системы совпадала с «частичной размерностью», то есть с размерностью, подсчитанной в проекции фазового подпространства на один осциллятор. В дальнейшем появились и другие подходы к построению количественной характеристики хаотической синхронизации, которые базировались на измерении взаимной информации между подсистемами [156], измерении «близости» фазовых точек в синхронизируемой подсистеме, отвечающих соседним фазовым точкам в синхронизирующем осцилляторе (так называемый метод «нелинейных взаимозависимостей») [156−160], а также на анализе отображений возврата («recurrent plots approach») [161]. Все указанные методы основывались на непосредственном анализе временных реализаций сигналов, генерируемых взаимодействующими подсистемами. С другой стороны, в радиофизике часто используется спектральный анализ сигналов, то есть переход от временного представления сигналов к частотным. Соответственно, появились алгоритмы построения количественной меры синхронизации, основанные на анализе спектров сигналов. Так в работе [162] для количественного измерения степени синхронности колебаний применялось построение распределений разностей Фурье-фаз.

Наличие большого числа разнообразных алгоритмов, дающих подчас противоречивые результаты свидетельствует о том, что задача нахождения «хорошей» универсальной количественной меры хаотической синхронизации пока не решена. В диссертационной работе проводятся количественные исследования полной и фазовой хаотической синхронизации, при использовании нормированной функции взаимной информации. Кроме того, в работе предлагается собственный метод подсчета степени согласования колебаний в подсистемах, использующий усредненную взвешенную функцию когерентности.

Методы нелинейной динамики носят универсальный характер и применяются для исследования колебательных и волновых процессов в системах самой различной природы. Традиционной областью применений для нелинейной динамики является химия [153,163, 164]. Химические процессы могут демонстрировать колебательный характер, когда концентрации веществ, участвующих в реакциях, регулярно или хаотически меняются во времени. Колебательный характер многих химических реакций неоднократно наблюдался в экспериментах, начиная со знаменитой реакции Белоусова-Жаботинского [165,166]. Кроме того, при определенных условиях, возможно построение динамических моделей таких реакций, которые позволяют выявлять их закономерности аналитически или методами численного эксперимента. Один из подходов к моделированию реакций — метод среднего поля, при котором вероятность реакции в каждой точке пространства одинакова и зависит от средних концентраций реагирующих веществ. Использование данного метода оправдано, если благодаря тем или иным процессам (например, сильная диффузия, конвекция или перемешивание) в пространстве реакций поддерживается все время однородная концентрация всех компонентов. Использование метода среднего поля позволяет построить модель реакций в виде системы обыкновенных дифференциальных уравнений, переменными которых являются концентрации химических реагентов. Если в ходе реакций существенными является диффузия веществ, то учет процессов диффузии приводит к моделям в виде системы уравнений в частных производных (уравнение «реакция — диффузия»). Существует, однако, класс реакций, в которых перемешивания компонент не происходит. При этом в пространстве реакций формируется сильно неоднородная среда из компонентов реакций. В этом случае методы среднего поля перестают адекватно описывать процессы превращений, поскольку теперь необходимо учитывать их локальный характер. Возможным подходам для построения модели становится использование различных клеточных автоматов [167−169], одним из которых является метод Монте-Карло (МК) [170,171]. В МК пространство реакций моделируется дискретной решеткой определенных размеров и размерностей. Если реакции идут на плоской поверхности, то размерность решетки выбирается равной двум, если в пространстве — то трем. Существуют модели с одномерными решетками [172,173], а также с фрактальными решетками, обладающими нецелыми размерностями [174]. Каждая «клетка» решетки заполняется той или иной «частицей» или же остается «вакантной». В ходе реакций на случайно выбранных клетках происходят взаимопревращения частиц, их «рождения» или «уничтожение», в соответствии с схемой реакций, при этом учитывается состояние соседних клеток. Сами превращения носят случайный характер, то есть они либо происходят, либо не происходят в соответствии с выбранными вероятностями. Модели МК достаточно полно учитывают случайный и локальный характер химических превращений и дают более реалистическую картину происходящих процессов. Кроме того, они позволяют получить не только значения средних концентраций, но их локальное распределение по всему пространству реакций. Детальный анализ образования пространственно-временных структур в таких моделях, изучение процессов синхронизации в них, а также исследование возможности сопоставления результатов МК моделирования с результатами более простых моделей среднего поля, — все это является актуальной и интересной задачей для нелинейной динамики.

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

В ходе проведения исследований, изложенных в диссертации, решались следующие научные задачи:

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

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

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

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

5) Разработка универсальных методов диагностирования и количественного анализа хаотической синхронизации на основе использования функции взаимной информации и функции когерентности. Тестирование этих методов на примерах взаимодействующих осцилляторов с разными типами синхронного поведения.

6) Исследование процессов синхронизации и образования пространственных структур в системах химических превращений, моделируемых с помощью метода Монте-Карло. Сопоставление результатов, полученных этим методом с результатами использования метода среднего поля.

Научная новизна.

В диссертационной работе впервые.

• показано, что различные системы с дискретным временем, объединенные общим сценарием формирования хаоса через каскад бифуркаций удвоения периода и общим типом диффузионной связи, обладают также общим сценарием формирования и разрушения режима полной синхронизации хаоса, который практически не зависит от индивидуальных особенностей подсистем, таких, как, например мультистабиль-ность/ моностабильность и обратимость/ необратимостьпроведен детальный бифуркационный анализ процесса разрушения синхронизации и формирования мультистабильности в связанных системах с дискретным временем, переход к хаосу в которых происходит через последовательность бифуркаций удвоения периодавыявлено влияние слабой неидентичности подсистем на бифуркации, ведущие к разрушению полной синхронизациипоказано, что в диффузионно связанных отображениях с бифуркациями удвоения периода наблюдаются только регулярные противофазные колебания с любым периодом, поскольку противофазные хаотические колебания оказываются трансверсально неустойчивымиисследован механизм формирования фазовой мультистабильности на базе противофазных колебанийпродемонстрирована возможность управляемой противофазной синхронизации хаоса в двух диффузионно связанных отображениях, исследован бифуркационный механизм управляемой противофазной синхронизациипроведен детальный бифуркационный анализ формирования и разрушения режимов полной и частичной синхронизации в кольце из трех дискретных отображений с симметричной диффузионной связьюпостроена структура областей существования пространственно — периодических режимов с разными длинами волн на плоскости параметров «связь — возбуждение» для ангармонических осцилляторов и выявлено качественное различие в устройстве пространства параметров для гармонических и негармонических автоколебанийобнаружены и описано новое явление в ансамблях ангармонических осцилляторов, заключающееся в появлении на фоне волнового режима движущихся локальных сбоев фазыпроведен анализ закономерностей перехода к пространственному беспорядку от пространственной периодичности при усложнении формы колебаний в ходе каскада бифуркаций удвоения периода в кольце осцилляторов, показано, что разрушение пространственной упорядоченности связано с разрушением фазовой синхронизации хаоса между соседними осцилляторамипоказана возможность принудительной синхронизации в цепочках конечной длины посредством высокочастотной модуляции параметра связи;

• предложен метод диагностики и количественного анализа хаотической синхронизации, основанный на усредненной по диапазону частот функции когерентности, показана его применимость для разных видов синхронного поведения и слабая чувствительность к искажениям экспериментальных данных;

• показана принципиальная разница в хаотической синхронизации систем с непрерывным и дискретным временем: в последнем случае в режиме когерентного хаоса отсутствует порог синхронизации по параметру связи;

• показана возможность диагностирования фазовой синхронизации посредством функции когерентности;

• предложен метод расчета степени синхронизации хаотических колебаний, основанный на нормированной функции взаимной информации. Показано, что этот метод дает возможность анализировать разные виды синхронного поведения: полную и фазовую синхронизацию хаоса, включая синхронизацию на гармониках.

Основные положения и результаты, выносимые на защиту.

1) Диффузионно связанные системы с удвоениями периода обладают общим сценарием разрушения полной синхронизации хаоса, который практически не зависит от особенностей индивидуальной динамики осцилляторов. Разрушение полной синхронизации происходит через два промежуточных этапа: «пузырящееся поведение» и «изрешечивание бассейна притяжения синхронного аттрактора». Бифуркационный сценарий этих явлений совпадает с механизмом формирования в данных системах фазовой мультистабильности.

2) В диффузионно связанных отображениях режим противофазной синхронизации наблюдается лишь для регулярных колебаний. На плоскости параметров «связь — возбуждение» области устойчивости периодических противофазных колебаний сужаются по мере увеличения периода колебаний, и на пороге перехода к хаосу их ширина становится равной нулю.

3) Для полной и противофазной синхронизации слабая неидентичность подсистем, а также малая асимметрия связи не меняет сценарий разрушения синхронизации и структуру фазового пространства в окрестности синхронного аттрактора, но меняет его бифуркационный механизм в части бифуркаций, обусловленных симметрией системы.

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

5) Функция когерентности может служить средством диагностики фазовой синхронизации хаоса как для идентичных автоколебательных систем, так и систем с расстройкой по параметрам. Данный метод позволяет уйти от использования расчета мгновенных фаз хаотических колебаний, которые нельзя корректно определить в случае развитого хаоса.

Научно-практическая значимость результатов.

Проведенные исследования относятся к фундаментальным проблемам радиофизики и нелинейной динамики. Научно-практическая значимость результатов заключается в том, что.

• проведен детальный бифуркационный анализ и построены карты динамических режимов для ряда колебательных систем, представляющие интерес в радиофизике;

• предложен метод принудительной синхронизации противофазных хаотических колебаний и определены области значений параметров, где эта синхронизация реализуется;

• предложен и апробирован метод принудительной синхронизации хаотических колебаний в ансамбле осцилляторов посредством параметрического воздействия;

• разработан простой метод диагностирования и количественного анализа синхронизации хаоса, применимый для разных типов синхронизации;

• проведено сопоставление методов среднего поля и Монте-Карло для исследования динамики химических реакций на поверхности каталитической решетки, показаны границы применимости метода среднего поля.

Результаты работы используются в учебном процессе в Саратовском государственном университете при чтении общих курсов («Теория колебаний»), а также спецкурсов по специальности «радиофизика» .

Структура и объем работы.

Работа состоит из введения, шести глав, заключения и списка литературы. Она содержит 416 страниц, включая 155 рисунков. Библиография содержит 216 ссылок на литературные источники.

Апробация работы и публикации.

Основные результаты работы были представлены на следующих научных конференциях:

1) 5-th international school on chaotic oscillations and pattern formation CHAOS-98, Saratov, Russia, October 6−10, 1998.

2) 2nd International Conference «Control of Oscillations and Chaos», vol. 3, July 5−7, St. Petersburg, Russia, 2000.

3) Международная межвузовская конференция «Современные проблемы электроники и радиофизики СВЧ», Саратов, Российская Федерация, 20−24 марта 2001 г.

4) 6-th International School on Chaotic Oscillations and Pattern Formation «CHAOS'Ol», Saratov, Russia, October 2−7, 2001.

5) International conference «Synchronization of chaotic and stochastic oscillations», Saratov, Russia, September 22−28, 2002.

6) International conference «Physics and Control», S. Petersburg, Russia, Aug. 20−22, 2003.

7) International conference «Physics and Control (PhysCon2005)», S. Petersburg, Russia, Aug. 20−22, 2005.

8) CEC AM Workshop on Dynamics at the Mesoscale: Theory, Modelling and Experiments September 8−11, Lyon, France, 2004.

9) VII международная школа «Хаотические автоколебания и образование структур», 1−6 октября 2004 г. Саратов.

10) Научная школа-конференция «Нелинейные дни в Саратове для молодых», 2005.

11) Международная конференция молодых учёных по фундаментальным наукам «Ломоносов 2006». Москва.

Результаты работы неоднократно обсуждались на научных семинарах:

• кафедры радиофизики и нелинейной динамики СГУ;

• центра нелинейной динамики и биофизики при СГУ;

• института физической химии национального научного центра республики Греция «Demokritos» (Греция, Афины);

• лаборатории нелинейной динамики Потсдамского университета (Германия, Потсдам);

• кафедры нелинейной динамики технического университета г. Лодзь (Польша);

• кафедры нелинейных систем университета г. Лозанны (Швейцария).

По материалам работы опубликовано 46 работ: 33 статьи в реферируемых научных журналах (из них 27 статей в журналах, рекомендованных ВАК для опубликования результатов докторских диссертаций), 4 статьи в сборниках международных научных конференций и 9 тезисов докладов. Результаты работы использованы при выполнении госбюджетных тем «Радиофизика», «Радиофизика-2», «Радиофизика-3», при выполнении грантов РФФИ, совместного гранта министерства Образования и науки РФ и американского фонда развития гражданских исследований (СШЭР).

Работа выполнена в докторантуре на кафедре радиофизики и нелинейной динамики Саратовского государственного университета (научный консультант — доктор физ.-мат. наук профессор Астахов В.В.).

Краткое содержание работы.

В первой главе исследуются бифуркации, ведущие к разрушению полной синхронизации хаоса в двух связанных системах с удвоениями периода. В качестве моделей используются дискретные отображения, демонстрирующие при изменении управляющего параметра переход к хаосу через каскад бифуркаций удвоения периода: логистические отображения, кубические отображения, отображения Эно. Проводится двухпараметриче-ский бифуркационный анализ на плоскости параметров «коэффициент связи — параметр возбуждения». Результаты анализа сопоставляются с наблюдаемыми явлениями, которые сопровождают разрушение синхронного состояния: пузырение хаотического аттрактора и его изрешеченность бассейна притяжения. Исследуется влияние неидентичности на бифуркационный сценарий и наблюдаемые явления. Сопоставляется процесс разрушения синхронизации в обратимых и необратимых системах.

Вторая глава посвящена исследованию противофазной синхронизации в отображениях с дискретным временем и потоковых системах. Проводится бифуркационный анализ формирования противофазных синхронных режимов и формирования фазовой мульти-стабильности в окрестности инвариантного симметричного подпространства, содержащего аттракторы противофазных колебаний. Проводится анализ устойчивости противофазных режимов для диффузионно связанных отображений. Строятся области существования для устойчивых периодических орбит раличных периодов и исследуются бифуркации, в результате которых, они теряют трансверсальную устойчивость. Предлагается алгоритм управления, стабилизирующий трансверсально неустойчивые хаотические предельные множества и тем самым обеспечивающий принудительную синхронизацию противофазных хаотических колебаний. Исследуется бифуркационный механизм разрушения управляемой противофазной синхронизации. На примере связанных осцилляторов Чуа с симметричной емкостной связью исследуются процессы противофазной синхронизации регулярных и хаотических колебаний в осцилляторах с непрерывным временем. Исследуется сценарий выхода из режима хаотической синхронизации.

Третья глава посвящена вопросам синхронизации в небольших ансамблях отображений с дискретным временем и осцилляторов. На примере трех логистических отображений исследуются механизмы формирования и разрушения режимов полной и частичной синхронизации хаоса. Строятся двупараметрические бифуркационные диаграммы для семейств режимов, удовлетвояряющих условиям полной или частичной синхронизации. Исследуются режимы бегущих волн в ансамбле связанных периодических осцилляторов при переходе от гармонических к релаксационным колебаниям. Строятся области устойчивости для мультистабильных состояний с разными пространственными периодами. Рассматриваются новые колебательные режимы — волны с движущимися фазовыми дефектами. Определяются области притяжения для волн с разными пространственными периодами. Исследуется переход к хаосу через каскад бифуркаций удвоения периода в цепочках генераторов Чуа.

В четвертой главе изучается возможность полной синхронизации хаоса в ансамблях осцилляторов при периодическом изменениии коэффициента связи. Аналитически решается задача постоения области устойчивости синфазных колебаний в зависимости от амплитуды и частоты параметрической модуляции. Исследуется синхронизация переключений в связанных осцилляторах Дуффинга при параметрическом воздействиии. Определяется максимальная длина цепочки осцилляторов, при которой возможна принудительная синхронизация при заданной амплитуде воздействия.

Пятая глава посвящена задаче количественного анализа хаотической синхронизации. Предлагается использовать для этой цели усредненную функцию когерентности, а также, нормированную функцию взаимной информации. Проводится тестирование предлагаемых алгоритмов на связанных осцилляторах с дискретным и непрерывным временем. Рассматриваются явления полной и фазовой синхронизации. Исследуется влияние внешнего шума, а также искажений обрабатываемых данных в канале связи. Предлагается новый метод диагностирования фазовой синхронизации, основанный на использовании функции когерентности.

В шестой главе проводится моделирование сложной пространственно — временной динамики химических реакций на поверхности катализатора при помощи метода Монте.

— 30.

Карло. Исследуется процесс образования гомогенных кластеров и определяется связь этого процесса с колебаниями средних концентраций реагирующих веществ. Проводится статистический анализ пространственного распределения веществ на решетке. Исследуется влияние малого перемешивания на динамику системы. Проводится сопоставление полученных результатов с результатами моделирования методом среднего поля.

Основные результаты диссертационной работы можно сформулировать следующим образом:

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

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

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

4) Усложнение формы колебаний в ходе каскада бифуркаций удвоения периода противофазных орбит приводит к развитию фазовой мультистабильности, закономерности которой схожи с закономерностями развития фазовой мультистабильности на базе синфазных колебаний. В этом случае также происходит последовательное увеличение в два раза числа сосуществующих периодических решений после каждого удвоения периода. Для сосуществующих орбит, колебания в одном из осцилляторов отличаются от колебаний во втором на целое число периодов исходного цикла с точностью до множителя —1.

5) Режим хаотической противофазной синхронизации можно обеспечить с помощью добавления дополнительной управляющей связи. Выход за область значений управляющего параметра, соответствующих трансверсально устойчивому противофазному хаосу, приводит к разрушению синхронизации, сценарий которой зависит от степени асимметрии управляющей связи. При симметричной и слабо несимметричной связи, разрушение противофазной синхронизации происходит через этап негрубой синхронизации. за которым следует бифуркация прорыва. При сильной асимметрии связи, включая однонаправленную связь, разрушение синхронизации осуществляется через этапы негрубой синхронизации и изрешечивания бассейна притяжения хаотического аттрактора.

6) В кольце из трех отображений с симметричной диффузионной связью кроме полной синхронизации наблюдаются режимы частичной синхронизации. Регулярные и хаотические режимы частичной синхронизации появляются на базе орбиты периода два через каскад тангенциальных бифуркаций удвоения периода. Каждый из таких режимов ограничен по параметру связи линией седло-узловой (для регулярных режимов) или седло-седловой бифуркации (для хаотических режимов). Каждая из орбит, участвующих в каскаде тангенциальных бифуркаций удвоения периода бифуцирует еще и в трансверсальном направлении, что приводит как к локальной изрешеченно-сти хаотического аттрактора, соответствующего режиму частичной синхронизации хаоса, так и к появлению новых седловых режимов вне подпространства частичной симметрии. Вновь появившиеся в результате трансверсальных бифуркаций удвоения периода несинхронные периодические орбиты впоследствии становятся устойчивыми, что приводит к развитию мультистабильности.

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

8) Режимы бегущих волн наблюдаются в ансамбле осцилляторов с бифуркациями удвоения периода на всем протяжении перехода к хаосу. Зоны устойчивости этих режимов на плоскости параметров «связь — параметр возбуждения» представляют собой ограниченные области, расположенные таким образом, что области для более коротковолновых режимов оказываются включенными в области более длинноволновых. После перехода к хаосу строгая пространственная периодичность сменяется периодичностью «в среднем», которая постепенно разрушается при развитии временного хаоса или при уменьшении связи. Разрушение пространственной периодичности колебаний связано с разрушением фазовой синхронизации хаотических колебаний.

9) Возможна принудительная синхронизация хаотических колебаний посредством высокочастотной модуляции коэффициента связи. Данный эффект наблюдается в ограниченной области управляющих параметров на плоскости частота модуляции — амплитуда модуляции. Применение данного метода к бистабильным системам позволяет осуществлять синхронизацию процессов переключений между состояниями. Внутри это области бистабильная система демонстрирует высоко когерентный процесс переключений между состояниями системы. Найденный эффект не зависит от природы колебаний осциллятора. Он наблюдается и при динамической и при шумовой природе переключений. Параметрическое воздействие позволяет также осуществлять синхронизацию хаоса в цепочках осцилляторов, но только при конечной длине цепочки.

10) Функция когерентности является мощным средством для анализа явления хаотической синхронизации. Индекс синхронизации, построенный на ее основе может применяться для измерения и диагностики различных типов синхронного поведения: полной, обобщенной или фазовой синхронизации. В последнем случае функция когерентности дает возможность диагностировать синхронизацию и при отсутствии расстройки между осцилляторами. Проведенный анализ показал принципиальное отличие синхронизации в системах с непрерывным и дискретным временем. Если в первом случае индекс синхронизации всегда равен нулю при нулевой связи, то в случае дискретных отображений он может оставаться на на некотором положительном уровне даже при полном отсутствии связи.

Заключение

.