Научная статья на тему 'Эффект "гибели колебаний" и квазипериодические бифуркации в низкоразмерном ансамбле осцилляторов ван дер Поля'

Эффект "гибели колебаний" и квазипериодические бифуркации в низкоразмерном ансамбле осцилляторов ван дер Поля Текст научной статьи по специальности «Физика»

CC BY
171
40
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
СИНХРОНИЗАЦИЯ / SYNCHRONIZATION / ЭФФЕКТ ГИБЕЛИ КОЛЕБАНИЙ / OSCILLATOR DEATH / КВАЗИПЕРИОДИЧЕСКИЕ БИФУРКАЦИИ / QUASI-PERIODIC BIFURCATIONS

Аннотация научной статьи по физике, автор научной работы — Кузнецов Александр Петрович, Тюрюкина Людмила Владимировна

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — Кузнецов Александр Петрович, Тюрюкина Людмила Владимировна

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

"Oscillator death" and quasiperiodic bifurcations in low-dimensional ensemble of van der Pol oscillators

The dynamics of the four dissipatively coupled van der Pol oscillator is considered. Lyapunov chart is presented in the parameter plane and its arrangement is discusses. The effect of increase of the threshold for the <> regime and the possibility of complete and partial broadband synchronization are revealed. We discuss the bifurcations of tori in the system at large frequency detuning of the oscillators, in particular, quasi-periodic saddle-node and Hopf bifurcations.

Текст научной работы на тему «Эффект "гибели колебаний" и квазипериодические бифуркации в низкоразмерном ансамбле осцилляторов ван дер Поля»

Бифуркации ^^^^^^^^^^^^^^^^^

в динамических системах

Изв. вузов «ПНД», т. 21, № 2, 2013

ЭФФЕКТ «ГИБЕЛИ КОЛЕБАНИЙ» И КВАЗИПЕРИОДИЧЕСКИЕ БИФУРКАЦИИ В НИЗКОРАЗМЕРНОМ АНСАМБЛЕ ОСЦИЛЛЯТОРОВ ВАН ДЕР ПОЛЯ

А. П. Кузнецов, Л. В. Тюрюкина

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

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

Введение

Синхронизация ансамблей осцилляторов является важной проблемой, представляющей как общетеоретический интерес, так и интерес с точки зрения приложений в радиофизике, лазерной физике, биофизики и т.д. [1-6]. Картина возможных типов режимов взаимодействующих осцилляторов достаточно хорошо изучена для малого их числа, или, наоборот, для больших массивов осцилляторов. Здесь мы обратимся к случаю ансамбля, содержащего четыре осциллятора ван дер Поля, и обсудим основные характерные типы режимов. Оказывается, что увеличение числа осцилляторов до четырех приводит к ряду новых моментов, которые и являются предметом настоящей работы. Будем рассматривать случай диссипативной связи, величина которой определяется относительной скоростью осцилляторов. Будем обсуждать традиционную для таких задач плоскость параметров частотная расстройка - величина связи осцилляторов. Как известно, для случая двух осцилляторов на этой плоскости выделяются следующие характерные области [1-6]. Во-первых, это основная область синхронизации с соотношением захваченных частот 1:1. Далее, область квазипериодических режимов со встроенной системой языков синхронизации с другими рациональными соотношениями частот. Еще одна характерная область - об-

ласть «гибели колебаний» (oscillator death). Она располагается при достаточно больших частотных расстройках при выполнении условия на величину связи ц и параметр отрицательного трения X < ц [1]. Интересный режим наблюдается для неидентичных по параметрам отрицательного трения осцилляторов, когда Xi = X2 [7-9]. В этом случае между областями гибели колебаний и квазипериодических режимов возникает область периодических режимов Xi > ц > X2, простирающаяся в сторону больших частотных расстроек. Причина ее появления состоит в том, что при выполнении условий Xi > ц > X2 первый осциллятор доминирует над вторым, который подавлен диссипативной связью. Такие режимы в [8,9] названы широкополосной синхронизацией.

Для трех идентичных осцилляторов в дополнение к указанной картине возникает возможность двух типов квазипериодических колебаний, отвечающих двум и трем несоизмеримым частотам [1,2,10]. На плоскости параметров они очень просто разделены при большом рассинхронизме осцилляторов: трехчастотной квазипериодичности отвечают меньшие значения связи. При малом рассинхронизме наблюдаются резонансы, и возникает картина языков двухчастотных торов, погруженных в область трехмерных торов [10,11]. Полная широкополосная синхронизация для идентичных по параметру возбуждения таких осцилляторов не наблюдается, она возможна только, если есть отличия в параметрах отрицательного трения. Однако в области X/2 < ц < X возникает широкополосный квазипериодический двухчастотный режим (подробности см. в [11] и ниже).

В настоящей работе, как отмечалось, обсуждается случай четырех осцилляторов. Будем использовать метод ляпуновских карт [10,11], в рамках которого выбирается точка на плоскости параметров и вычисляется спектр ляпуновских показателей системы Li. В зависимости от знаков и числа нулевых старших показателей1, определяется тип режима (аттрактора) в системе.

• P - предельный цикл, L1 = 0, L2 < 0, L3 < 0, L4 < 0, L5 < 0,

• T2 - двухчастотный тор, L1 = L2 = 0, L3 < 0, L4 < 0, L5 < 0,

• T3 - трехчастотный тор, L1 = L2 = L3 = 0, L4 < 0, L5 < 0,

• T4 - четырехчастотной тор, L1 = L2 = L3 = L4 = 0, L5 < 0,

• C - хаос, L1 > 0, L2 < 0, L3 < 0, L4 < 0, L5 < 0.

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

Следует отметить, что наличие торов разной размерности приводит к возможности квазипериодических бифуркаций. Мы обсудим основные типы таких бифуркаций в анализируемой системе. К сожалению, численные алгоритмы поиска квазипериодических бифуркаций очень трудоемки и находятся, во многом, в стадии разработки [12]. Однако поведение ляпуновских показателей в ряде случаев позволяет идентифицировать тип бифуркации, что мы и используем.

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

1. Случай трех осцилляторов

Для дальнейшего понадобятся более детальные сведения о случае трех дис-сипативно связанных осцилляторов [11]. Цепочка таких осцилляторов описывается уравнениями

X — (X — x2 )X + x + ц(Х — yy) = 0,

y — (X — y2)y + (1 + Ai )y + ц(у — X) + ц(у — z) = 0, (1)

z — (X — z2)z + (1 + A2 )z + Ф — y) = 0.

Здесь X - параметр возбуждения автономных осцилляторов, характеризующий отрицательное трение; Ai - частотная расстройка второго осциллятора относительно первого; A2 - третьего относительно первого; ц - коэффициент диссипативной связи.

Карта ляпуновских показателей системы (1) на плоскости (Ai, ц) для X = 0.1 и A2 = 0.05 показана на рис. 1. На ней располагается область полной синхронизации трех осцилляторов P, а также области двухчастотных торов T2 и трехчастотных торов T3. Хаотические режимы при данном наборе параметров не выявляются. Еще одна характерная область показана серым цветом и обозначена OD. Она отвечает режиму «гибели колебаний». Такой эффект имеет место при условии ц > X, когда трение между осцилляторами превышает отрицательное трение, свойственное автоколебательным элементам. Для удобства восприятия значение ц = X отмечено стрелкой на правой границе карты. Можно видеть, что оно является асимптотической границей области гибели колебаний при возрастании частотной расстройки Ai .

При уменьшении параметра связи в области X/2 < ц < X появляется широкая по частоте полоса области двухчастотной квазипериодичности, обозначенная на рис. 1 PBS. Этот режим в [11] назван частичной широкополосной синхронизацией. Его природа состоит в том, что связь по-разному возмущает осцилляторы в цепочке, даже если их параметры возбуждения одинаковы. Действительно, центральный осциллятор испытывает трение со стороны двух соседей, и эффективный коэффициент трения для него повышается в два раза. Поэтому при условии X > ц > > X/2 второй осциллятор подавлен соседями, которые все еще достаточно возбуждены [11].

0.15

■0.075 А, 0.125

Рис. 1. Ляпуновскиая карта системы трех связанных осцилляторов ван дер Поля (1), OD - область гибели колебаний, PBS - область частичной широкополосной синхронизации. Значения параметров X =0.1, A2 = 0.05

2. Случай четырех осцилляторов. Эффект повышения порога области гибели колебаний

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

(2)

стота же второго осциллятора

Д1

X — (X — х2)Х + х + ц(Х — у) = 0,

У — (X — у2)у + (1 + Д1 )у + ц(у — X) + ц(у — г) = 0,

г — (X — г2)г + (1 + Д2)-г + ц(г — у) + ц(г — ад) = 0,

г — (X — ад2 )г + (1 + Д3)г + — г) = 0.

Карта ляпуновских показателей для X = 0.1 на плоскости (Д1, ц) в разных масштабах показана на рис. 2. Остальные параметры выбраны так, чтобы частоты первого, третьего и четвертого осцилляторов были достаточно близки: Д2 = 0.03, Д3 = 0.1. Ча-

является варьируемым параметром. При небольших уровнях связи на рис. 2, а в окрестности нулевой расстройки Д1 можно видеть проявление резонансных эффектов для торов разной размерности. Так, наблюдаются два характерных языка трехмерных торов, погруженных в область четырехмерных. Эти языки имеют основания на оси абсцисс, которые отвечают двум резонансным условиям Д1 = 0 и Д1 = Д2. Физически, согласно (1), они соответствуют резонансам, когда частота второго осциллятора совпадает с частотой первого (Д1 = 0) или с частотой третьего (Д1 = Д2). Резонансные значения отмечены стрелочками на рис. 2, а. В свою очередь, области двухчастотных торов также имеют вид языков, однако, все они имеют порог по величине связи. Основания некоторых языков разрушены с возникновением областей хаоса С. Отметим, что в целом картина на рис. 2, а при небольшой связи качественно близка к результатам фазового приближения [10].

Обратимся теперь к рис. 2, б, отвечающему более широким диапазонам изменения параметров. Теперь в «поле

Рис. 2. Карта динамических режимов четырех диссипативно связанных осцилляторов ван дер Поля (2) в разных масштабах, X = 0.1, Д2 = 0.03, Дз = 0.1. Вертикальными стрелками на рис. а обозначены значения резонансных условий Дх = 0 и Дх = Д2. Горизонтальные стрелки на рис. б отвечают различным бифуркациям (см. пояснения в тексте)

зрения» попадает область «гибели колебаний» ОЭ. Однако нижняя граница этой области при большой частотной расстройке Д1 не стремится к значению ц = X, как это было в случае двух и трех осцилляторов, а довольно значительно превышает этот порог. Таким образом, для четырех осцилляторов наблюдается эффект повышения порога гибели колебаний. Из рис. 2, б можно видеть, что это превышение значительное: граница области гибели колебаний (точка Н на рисунке) отвечает примерно ц ^ = 0.26.

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

Проиллюстрируем сказанное некоторыми оценками. Запишем укороченные уравнения (уравнения Ландау-Стюарта)2:

2а = Xа — |а|2 а — ц(а — Ь),

2Ь = XЬ — |Ь|2 Ь + гД1Ь — ц(2Ь — а — с),

2 (3)

2с = Xc — |с| с + гД2с — ц(2с — Ь — й).

2й = Xd — |й|2 й + гД3й — ц(й — с).

Здесь а, Ь, с, й - медленные комплексные амплитуды. Для определения области гибели колебаний, как известно, достаточно линейного приближения. Кроме того, полагаем Д1 ^ то, а собственные частоты первого, третьего и четвертого осцилляторов малыми (по сравнению с Д1). Для оценок положим, что они просто совпадают, и Д2 = Д3 = 0. Тогда

2а = Xa — ц(а — Ь),

2Ь = Xb + гД1Ь — ц(2Ь — а — с),

(4)

2с = Xc — ц(2с — Ь — й).

2й = Xd — ц(й — с).

2Их вывод является стандартным для метода медленно меняющихся амплитуд [1,2] и поэтому здесь

не приводится.

Ищем решение этой линейной системы в виде экспоненциальной подстановки ехр(а£/2). Тогда получаем

(5)

(а — X + ц)а = цЬ, (а — X + 2ц — гД1)Ь = ц(а + с), (а — X + 2ц)с = ц(Ь + й), (а — X + ц)й = цс.

Условие рассинхронизма второго осциллятора Д1 ^ то в соответствии со вторым уравнением (5) требует, чтобы а/Ь ^ то, либо с/Ь ^ то. Физически это отвечает тому, что либо первый, либо третий осциллятор доминируют над вторым. (Напомним, что переменные а, Ь, с отвечают комплексным амплитудам соответствующих осцилляторов.)

Оценим первое отношение. Из первого уравнения (5) получаем

а/Ь =-^-. (6)

' а — X + ц 4 у

Отсюда видим, что а/Ь ^ то, если а = X — ц. Условие гашения колебаний определяется отрицательным декрементом затухания а < 0, так что получаем ц > X, то есть точно такое условие, как для двух или трех осцилляторов.

Оценим теперь второе отношение, связанное с доминированием третьего осциллятора. Из второго и третьего уравнений (5) получаем

с/Ь =-^-д-. (7)

Условие с/Ь ^ то выполняется при обращение в нуль знаменателя, что приводит к квадратному уравнению относительно декремента затухания а, которое легко решается

г

а = к---—(х. (8)

Для возникновения режима гашения колебаний должно быть а < 0 для обоих корней, что приводит к условию

ц > —и 2.618 • А. (9)

3 — V 5

Таким образом, ситуация доминирования третьего осциллятора при условии его эффективного взаимодействия с четвертым значительно повышает порог эффекта гибели колебаний - более чем в два раза.

Заметим, что если Д2 = 0, Д3 = 0 можно провести аналогичное рассмотрение, только получаются громоздкие выражения. При этом оценка (9) довольно эффективна и для случая, представленного на рис. 2, б. Действительно, при X = 0.1 она дает ц ~ 0.26, что очень близко к асимптотической границе области гибели колебаний на рис. 2, б. Это не удивительно, так как соответствующие частотные параметры малы: Д2 = 0.03, Д3 = 0.1.

На рис. 3 показаны фазовые портреты всех четырех осцилляторов для разного уровня связи при достаточно большой частотной расстройке второго осциллятора Ai = 2.0. (Обращаем внимание на разные масштабы по осям координат на фазовых портретах различных фрагментов.) Рис. 3, а отвечает случаю, когда ^ < X, но близко к этому значению. На карте это соответствует области частичной широкополосной синхронизации PBS. Видно, что второй осциллятор, выведенный из условий синхронизации, достаточно сильно подавлен связью. При этом его траектория сильно возмущена квазипериодическим образом по сравнению с автономным случаем. Все остальные осцилляторы возбуждены примерно равноценно. Причем это относится и к третьему осциллятору, испытывающему двойное трение со стороны соседей и для которого порог ^ = X/2 для выбранного значения параметра связи превышен.

Рис. 3, б относится к области периодических режимов. При этом выбрано значение ^ = 0.2, когда заметно превышен соответствующий порог ^ = X, отвечающий компенсации диссипативной связью отрицательного трения отдельного осциллятора. Можно видеть, что квазипериодический режим сменился периодическим. Теперь сильно подавлены связью как второй, так и первый осцилляторы. Несмотря на сильную диссипацию ^ > X, третий и четвертый осциллятор достаточно сильно возбуждены. Их орбиты примерно одинакового (по порядку) радиуса и не отличаются существенно от случая, приведенного на рис. 3, а. Таким образом, действительно возникла заметная синхронизация третьего и четвертого осцилляторов, которая и «затягивает» возникновение режима гибели колебаний.

Описанный механизм приводит также еще к одному эффекту, а именно к возможности полной широкополосной синхронизации CBS четырех осцилляторов для случая идентичных управляющих параметров X. Действительно, на рис. 2, а можно видеть характерную полосу CBS в диапазоне уровня связи 0.1 < ^ < 0.27 (между точками H и NS на карте). В случае трех идентичных осцилляторов такой режим не возникает (ср. рис. 2, б и рис. 1). Для его реализации для трех (как и для двух) осцилляторов необходима неидентичность по управляющим параметрам [11].

Рис. 3. Фазовые портреты четырех осцилляторов: а - режим частичной широкополосной синхронизации, = 0.08; б - режим полной широкополосной синхронизации, = 0.2. Цифры отвечают номеру осциллятора. Значения параметров X = 0.1, Ах = 2.0, А2 = 0.03, Аз = 0.1

3. Квазипериодические бифуркации в системе четырех осцилляторов

Продолжим обсуждение рис. 2, б. При больших частотных расстройках Ai с уменьшением параметра связи наблюдается каскад бифуркаций рождения инвариантных торов все более высокой размерности. Обсудим картину бифуркаций, сопровождающих наблюдаемый эффект. Для этого обратимся к графику зависимости пяти старших ляпуновских показателей Л^ от величины связи (рис. 4). Параметр частоты второго осциллятора выберем достаточно большим (Ai = 2.0), чтобы наблюдалась существенная рассинхронизация этого осциллятора с остальными.

При большой связи реализуется эффект гибели колебаний, и все показатели на рис. 4 отрицательны. В точке ^ ^ 0.27 первый показатель обращается в ноль, и происходит бифуркация Андронова-Хопфа H рождения устойчивого (многомерного) предельного цикла.

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

Затем в точке QNS в ноль обращается и третий показатель Л3. При этом наблюдается характерное поведение: до точки бифуркации третий и четвертый показатель равны друг другу (Л3 = Л4). За точкой бифуркации показатель Л3 обращается в ноль, а четвертый показатель Л4, коснувшись горизонтальной оси, снова уходит в отрицательную область. Это характерная картина квазипериодической бифуркации, приводящей к мягкому возникновению трехчастотного тора из двухчастотного [12]. В [12] такую бифуркацию называют квазипериодической бифуркацией Хопфа (quasi-periodic Hopf bifurcation), хотя более логичным было бы говорить о квазипериодической бифуркации Неймарка-Сакера QNS3.

В точке SNT размерность тора вновь повышается. Однако это бифуркация другого типа - теперь четвертый показатель просто обращается в ноль, а пятый показатель все время остается отрицательным. Это означает, что имеет место седло-

узловая бифуркация торов [12]: устойчивый трехмерный тор сливается с сед-ловым, исчезает, и возникает четырехмерный устойчивый тор. Отметим, что аналогичная бифуркация реализуется и в фазовом приближении [10]. Остальные бифуркации в фазовой модели невозможны.

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

3Поскольку именно бифуркация Неймарка-Сакера отвечает за рождение тора.

Рис. 4. Графики пяти старших ляпуновских показателей системы четырех связанных осцилляторов ван дер Поля, I = 0.1, Ах = 2, А2 = 0.03, Аз = 0.1. Здесь Н - точка бифуркации Андронова-Хопфа, N8 - точка бифуркации Неймарка-Сакера, - точка квазипериодической бифуркации Неймарка-Сакера двухчастотного тора, $N7 - точка седло-узловой бифуркации трехчастотного тора

Заключение

В системе четырех диссипативно связанных осцилляторов ван дер Поля наблюдаются определенные особенности в устройстве характерных областей на плоскости параметров. Так, режим «гибели колебаний» возникает при значительном превышении величины связи значения параметра отрицательного трения, что объясняется взаимодействием подсистемы из двух осцилляторов. При большой частотной расстройке наблюдаются последовательно бифуркация Андронова-Хопфа, Неймарка-Сакера, мягкая квазипериодическая бифуркация Неймарка-Сакера и жесткая седло-узловая квазипериодическая бифуркация. При этом их линии являются асимптотическими границами областей широкополосной синхронизации (полной и частичной) для квазипериодических режимов разной размерности.

Авторы выражают благодарность И.Р. Сатаеву за полезное обсуждение.

Работа выполнена при поддержке РФФИ, проект 12-02-00541-а.

Библиографический список

1. Пиковский А., Розенблюм М., Куртс Ю. Синхронизация. Фундаментальное нелинейное явление. М.: Техносфера, 2003. 494 с.

2. Ланда П.С. Автоколебания в системах с конечным числом степеней свободы. М.: Наука, 1980. 360 с.

3. Анищенко В.С., Астахов В.В., Вадивасова Т.Е. Регулярные и хаотические автоколебания. Синхронизация и влияние флуктуаций: Учебник-монография. Долгопрудный: Издательский Дом «Интеллект», 2009. 312 с.

4. Блехман И.И. Синхронизация в природе и технике. М.: Наука, 1981. 351 с.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

5. Glass L. and MacKey M.C. From Clocks to Chaos. Princeton University Press, 1988. 272 p.

6. Balanov A.G., Janson N.B., Postnov D.E., Sosnovtseva O. Synchronization: From simple to complex. Springer, 2009. 426 p.

7. Ivanchenko M., Osipov G., Shalfeev V., Kurths J. Synchronization of two non-scalar-coupled limit-cycle oscillators // Physica D. 2004. Vol. 189. № 1-2. P. 8.

8. Кузнецов А.П., Паксютов В.И., Роман Ю.П. Особенности синхронизации в системе неидентичных связанных осцилляторов ван дер Поля и ван дер Поля-Дуффинга. Широкополосная синхронизация // Известия вузов. Прикладная нелинейная динамика. 2007. Т. 15, № 4. С. 3.

9. Kuznetsov A.P., Roman Ju.P. Properties of synchronization in the systems of non-identical coupled van der Pol and van der Pol-Duffing oscillators. Broadband synchronization // Physica D. 2009. Vol. 238, № 16. P. 1499.

10. Кузнецов А.П., Сатаев И.Р., Тюрюкина Л.В. Синхронизация и многочастотные колебания в цепочке фазовых осцилляторов // Нелинейная динамика. 2010. Т. 6, № 4. С. 693.

11. Емельянова Ю.П., Кузнецов А.П., Тюрюкина Л.В. Динамика трех неидентичных по управляющим параметрам связанных осцилляторов ван дер Поля // Известия вузов. Прикладная нелинейная динамика. 2011. Т. 19, № 5. С. 76.

12. Broer H., Simo C., Vitolo R. Quasi-periodic bifurcations of invariant circles in low-dimensional dissipative dynamical systems // Regular and Chaotic Dynamics. 2011. Vol. 16, № 1-2. P. 154.

Саратовский государственный Поступила в редакцию 21.09.2012

университет им. Н.Г. Чернышевского После доработки 19.12.2012

«OSCILLATOR DEATH» AND QUASIPERIODIC BIFURCATIONS IN LOW-DIMENSIONAL ENSEMBLE OF VAN DER POL OSCILLATORS

A. P. Kuznetsov, L. V. Turukina

The dynamics of the four dissipatively coupled van der Pol oscillator is considered. Lyapunov chart is presented in the parameter plane and its arrangement is discusses. The effect of increase of the threshold for the «oscillator death» regime and the possibility of complete and partial broadband synchronization are revealed. We discuss the bifurcations of tori in the system at large frequency detuning of the oscillators, in particular, quasi-periodic saddle-node and Hopf bifurcations.

Keyword: Synchronization, oscillator death, quasi-periodic bifurcations.

Кузнецов Александр Петрович - родился в 1957 году. Доктор физико-математических наук, ведущий научный сотрудник Саратовского филиала Института радиотехники и электроники РАН, профессор Саратовского госуниверситета, заведующий базовой кафедрой динамических систем СГУ в СФ ИРЭ РАН. Специалист по нелинейной динамике, теории динамического хаоса и теории критических явлений. Занимается использованием идей теории катастроф и теории бифуркаций, а также развитием концепции сценариев перехода к хаосу применительно к многопараметрическим модельным и физическим нелинейным системам. Лауреат гранта Президента Российской Федерации для молодых российских ученых - докторов наук. Лауреат персонального гранта Фонда содействия отечественной науки. Соросовский профессор (2000, 2001), научный руководитель студенческой лаборатории «Теоретическая нелинейная динамика». Опубликовал более 100 научных работ. Автор нескольких оригинальных учебных курсов для факультета нелинейных процессов и Лицея прикладных наук СГУ, 10 учебных пособий и монографии «Нелинейные колебания» (совместно с С.П. Кузнецовым и Н.М. Рыскиным. М.: Физматлит, 2002).

410012 Саратов, ул. Астраханская, 83

Саратовский государственный университет им. Н.Г. Чернышевского E-mail: [email protected]

Тюрюкина Людмила Владимировна - родилась в 1977 году. Окончила факультет нелинейных процессов в Саратовском госуниверситете (2000). Кандидат физико-математических наук (2003, СГУ), имеет звание доцента по специальности радиофизика (2009). В настоящее время старший научный сотрудник Саратовского филиала Института радиотехники и электроники РАН, доцент базовой кафедры динамических систем СГУ в СФ ИРЭ РАН. Область научных интересов - новые аспекты явления синхронизации в системах различной физической природы (радиофизические системы, модели турбулентности, модели биофизических систем и др.); контроль (управление) неустойчивыми режимами; динамический хаос; физические системы с гиперболическими аттракторами. Автор более 80 публикаций, в том числе 30 статей в российских и зарубежных журналах и 3 учебно-методических пособий.

410019 Саратов, ул. Зеленая, д. 38 Саратовский филиал Института радиотехники и электроники им. В.А. Котельникова РАН E-mail: [email protected]

i Надоели баннеры? Вы всегда можете отключить рекламу.