Мультистабильность синхронных режимов в ансамблях неидентичных осцилляторов: два элемента Текст научной статьи по специальности «Физика»

Изв. вузов «ПНД», т. 17, № 2, 2009

УДК 621.391.01

МУЛЬТИСТАБИЛЬНОСТЬ СИНХРОННЫХ РЕЖИМОВ В АНСАМБЛЯХ НЕИДЕНТИЧНЫХ ОСЦИЛЛЯТОРОВ:

Два элемента

А.К. Крюков, Г.В. Осипов, А.В. Половинкин

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

Ключевые слова: Синхронизация, мультистабильность, синхронный режим, численные методы, моделирование, стохастическая система, ФитцХью-Нагумо, Бонхоффер-ван дер Поль.

Введение

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

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

16

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

1. Модель

Система Бонхоффера-ван дер Поля [13] (иногда называемая системой Фитц Хью-Нагумо) часто применяется для моделирования нейронных сетей и генераторов ритма сердца [14]. Динамика цепочки локально диффузионно связанных неидентичных осцилляторов Бонхоффера-ван дер Поля описывается уравнениями

f x3

dj j

£Хj = Xj - -3 - yj + $(Xj+i - 2Xj + Xj-1),

yj = (xj + aj ), (1)

, j =

Здесь в - параметр связи; e - малая величина; aj - параметр, отвечающий за период и симметрию колебаний; N - количество элементов.

2. Один элемент

Один элемент описывается системой уравнений второго порядка

X

eX = x - — - y = F(x,y),

y = (x + a) = G(x,y),

\ a \< 1.

(2)

Система (2) имеет неустойчивое состояние равновесия (фокус при a > \/1 - 2^/г или узел при a < л/1 - 2\/ё) и устойчивый предельный цикл. Малость параметра e обуславливает наличие быстрых и медленных движений в элементе.

Типичный фазовый портрет для одного элемента показан на рис. 1, а. h- и h+ - устойчивые участки кривой медленных движений, ho - неустойчивый участок. Для удобства описания динамики участок предельного цикла, лежащий в окрестности h-, будем называть пассивной фазой, участок в окрестности h+ - активной фазой. Переход из пассивной фазы в активную и обратно будем называть спайком2.

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

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

17

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

3. Два связанных элемента

Рассмотрим два связанных элемента Бонхоффера-ван дер Поля

ex'i = F(xi,yi) + в(х2 - xi),

yi = (х + ai),

ex' 2 = F(Х2,У2) + P(xi - X2),

У2 = (x + a2),

x3

^ F(xj, yj) = xj —у — Уз> 3 = 2-

В компьютерных экспериментах было показано, что в этой системе при некоторых значениях параметров существуют два устойчивых предельных цикла. Первый соответствует синфазному синхронному режиму (осциллограмма приведена на рис. 2, а), второй - противофазному (осциллограмма на рис. 2, б). Будем считать, что два элемента синхронизованы синфазно, если для межспайкового интервала Ti2 выполняется следующее неравенство: Ti2 = t2 — ti ^ T, где tj - время перехода j-го элемента из пассивной фазы в активную, T - период синхронного режима. Если же Ti2 = t2 — ti ~ T/2, то элементы синхронизованы в противофазе. Реализация в системе того или иного синхронного режима зависит от начальных условий.

Для аналитического доказательства сосуществования синфазного и противофазного синхронных режимов в цепочке модель (3) можно упростить. Для этого нужно положить e = 0 и заменить нелинейную функцию F(xj,yj) на ее кусочно-

18

Рис. 2. Осциллограммы элементов xi(t), x2(t) в синфазном (а) и противофазном (б) режимах линейную аппроксимацию

-4

—Xj - 2 - yj, для Xj ^ -1,

F = 4 oXj - yj’

для -1 < Xj < 1,

-4

—Xj + 2 - yj, для Xj ^ 1.

ч 3

Тогда получим кусочно-линейную систему второго порядка

(4)

0 = F(Xi,X2,yi,У2)+ |3(X2 - Xi),

dyi = ,

—г— = Xi + ai,

ax

(5)

0 = F (Xi, X2, yi ,y2) + |3(Xi - X2),

dy2

—— = X2 + a2. ax

Характерное поведение каждого из элементов показано на рис. 1, б. Динамику связанных элементов удобно рассматривать, отмечая их взаимное положение на ветвях h±. Тогда в данной системе возможны четыре случая взаимного расположения элементов на цикле:

1) «+—+» - оба элемента находятся в активной фазе (этот случай показан на рис. 1, б);

2) «+—» - первый элемент расположен в активной фазе, второй - в пассивной;

3) «—+» - первый элемент расположен в пассивной фазе, второй - в активной;

4) «--» - оба элемента находятся в пассивной фазе.

19

Для каждого из этих случаев система (5) имеет аналитическое решение. Например, для случая «—+» решение выглядит так:

Vr+(t)

1

2

(Vi(to) - V2(to)) + 2 - (ai - а2)в

exp(

t - to

2p

) +

+

1

2

(Vi(to) + V2(to))

2

3

(ai + a2)

( 3t - to) 2

exP(-----~) - 2

|3a2 + |3ai,

(6)

V- +(t) =

+

2(V2(to) - Vi (to)) - 2 + (ai - a2)e

exp(-t) + v 2p ;

2(Vi(to >+V2(to)) - 2(ai+a2)

3t to

exp(-----■——) + 2 - (3a2 + pai,

(7)

где Vi (to), V2(to) - начальные условия соответствующей кусочно-линейной аппроксимации, р = 2/3 + в и р = 4/3 + р. Сшивая решения с различных временных интервалов, можно получить состояние системы в произвольный момент времени по заданным начальным условиям.

С использованием (6) и (7) было построено отображение разности df (n + 1) = = V2(tn+i) - Vi(tn+i) в момент времени tn+i (момент начала (n + 1)-го прохождения предельного цикла) от такой же разности в момент tn (момент начала n-го прохождения предельного цикла). Наличие неподвижной точки такого отображения означает существование синхронного режима с соответствующей разностью df. Для различных значений параметра связи были численно построены данные отображения (рис. 3).

Как видно из рис. 3, при малой связи в системе существует два различных синхронных режима - синфазный (V2 - Vi ~ 0) и противофазный. При большой связи - только синфазный.

Для характеристики синхронного режима в исходной модели вычислялись частоты всех элементов. Частота элемента вычислялась как средняя частота генерации спайков в элементе.

Рис. 3. Отображение разности df для слабой связи р = 0.002 (а) и сильной связи р = 0.05 (б) в системе двух элементов

20

Зависимость средних частот элементов исходной системы связанных осцилляторов Бонхоффера-ван дер Поля от параметра связи при синфазных и противофазных начальных условиях представлена на рис. 4. При слабой (но достаточной для синхронизации) связи оба эти режима устойчивы. Частота синфазного режима близка к частоте самого быстрого (высокочастотного) элемента цепочки и практически не зависит от параметра связи р. Частота противофазного синхронного режима при относительно слабой связи близка к собственной частоте самого низкочастотного элемента и уменьшается с увеличением связи. Численно показано, что, если параметр связи больше некоторого порогового значения, противофазный режим становится неустойчивым и в системе реализуется только синфазный синхронный режим.

3.1. Возникновение и разрушение синхронных режимов. Рассмотрим детально механизмы возникновения и разрушения синхронных режимов в ансамбле двух связанных нелинейных осцилляторов типа (3) с увеличением параметра связи в в случае а1 = 0.995, е = 0.02 при Aa = а1 — а2 = 0.001 (случай относительно малых расстроек) и при Aa = 0.01 (случай относительно больших расстроек).

В обоих случаях существует критическое значение связи между элементами p^nin такое, что при в < р^п в системе наблюдаются биения - разность времен прохождения характерных значений (например, максимумов) координатами xi(t) и x2(t) или y1(t) и y2(t) непрерывно дрейфует, изменяясь от нуля до максимального значения.

Анализ отображений Пуанкаре и характеристических показателей Ляпунова, проведенный с использованием [9], позволяет заключить, что при изменении р в пределах от 0 до P1lnin в системе реализуются различные сложные режимы. При этом в случае Aa = 0.001 (рис. 5, а) с точностью до вычислительной погрешности по-

Рис. 4. Зависимость частот синхронизации от параметра связи при а1 = 0.995, a2 = 0.985, е = 0.02 (серые точки соответствуют синфазным начальным условиям, черные - противофазным)

Н|хго

-0.6768

-0.6769

-0.6770

-0.6771

-0.6772

П|хго

X

-0.675 - X N

-0.676

gауУ f frii,|iirш fitiimiiMmiiiM ■ -0.677 X , ^

-0.678

а

0.0 2.0 4.0 6.0 8.0 10.0 piO5 0.0 2.0 4.0 6.0 8.0 piO4

Рис. 5. Бифуркационная диаграмма отображения Пуанкаре (значение у1 при переходе элемента из пассивной фазы в активную при х1 = 0 в зависимости от параметра р) в системе связанных нелинейных осцилляторов ван дер Поля в области фазовых биений р < рП^ Aa = 0.001 (а) и 0.01 (б)

21

казателей Ляпунова можно сделать вывод о хаотическом характере наблюдаемых колебаний. В случае Да = 0.01 (рис. 5, б) эти режимы являются достоверно хаотическими и перемежаются окнами n : m синхронизации парциальных элементов.

C приближением в снизу к значению P^nin частота биений уменьшается, и при вПт = 1-23 • 10_4 (в случае Да = =0.001) или при Pnin=1-16 • 10_3 (в случае Да = 0.01) в системе возникает (c максимальным значением мультипликатора отображения Пуанкаре ц = +1) устойчивый предельный цикл, соответствующий 1 : 1 противофазному режиму синхронизации элементов (рис. 6).

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

При дальнейшем увеличении P вплоть до значений вП^ = 1-70 • 10_4 (в случае Да = 0.001) или при вПт = 1-95 • 10_3 (в случае Да = 0.01) в системе рождает-

ся (при возникновении также c максимальным значением мультипликатора ц = +1 отображения Пуанкаре) ещё один устойчивый предельный цикл, соответствующий уже синфазному режиму синхронизации и, таким образом, реализуется бистабильность синфазного и антифазного режимов синхронизации.

Отметим, что противофазный режим синхронизации возникает при меньших значениях параметра связи в по сравнению с синфазным. Данному факту можно

Рис. 6. Области: биений (1), моностабильности противофазного режима синхронизации (2), бистабильности противофазного и синфазного режимов (3), моностабильности синфазного режима синхронизации (4). Бифуркационные границы рД^Дя), Ртт(Дя), Ртах(Да)

Рис. 7. Проекции фазовых траекторий на плоскость (xi,x2) в случае Да = 0.001 и р < Piim на временных отрезках, соответствующих колебаниям, близким к противофазным - (а), и колебаниям, близким к синфазным - (б)

22

дать простую интерпретацию: за возникновение синхронизации «ответственно» линейное слагаемое в (xj — Xi),j = i В синфазном режиме типичные значения модуля разности Xj — Xi меньше, чем в противофазном, поэтому для заметного воздействия на систему требуются большие значения параметра р.

При р = pmax = 1.89 • 10_2 (в случае Да = 0.001) предельный цикл, соответствующий противофазному режиму синхронизации исчезает в результате (р = — 1)-бифуркации (вследствие слияния с седловым предельным циклом двукратного периода). В случае Да = 0.01 исчезновение предельного цикла, соответствующего противофазному режиму синхронизации, происходит при р = pmax = 2.17 • 10_2 в результате (р = +1)-бифуркации и, таким образом, при р > pmax в любом случае в системе реализуется только синфазный режим синхронизации. Некоторые отличия случая относительно больших расстроек (Да = 0.01) связаны также с тем, что цикл, соответствующий противофазному режиму синхронизации, при определенных значениях параметров может претерпевать бифуркацию удвоения периода, однако, с точки зрения существования мультистабильных синхронных режимов, данный факт не является существенным, так как фазовые траектории цикла удвоенного периода проходят в малой окрестности исходного противофазного цикла периода Т, не нарушая при этом факт существования бистабильности. Интересно отметить, что линии вт1п(Да) и рт;„(Да) на плоскости р, Да близки к прямо пропорциональной зависимости, в то время как зависимость ртах(Да) является, безусловно, нелинейной.

3.2. Смена динамических режимов в системе двух релаксационных осцилляторов в присутствии флуктуаций. Для анализа влияния флуктуаций (реально существующих в любых физических и биологических системах) на смену динамических режимов в ансамбле двух связанных релаксационных осцилляторов Бонхоффера-ван дер Поля при различных значениях параметра р численно исследовались статистические характеристики решений системы (8), которая отличается от (3) заменой времени и добавлением независимых источников белого гауссова шума ^i(t) и ^2(t) в правые части 1-го и 3-го уравнений,

X3

X1 = Xi — у — yi + р(х2 — Xi) + h(t),

yi = e(Xi + ai), (li(t)) = 0,

X3

X 2 = X2 — у — У2 + р^! — X2) + ^2 (t),

У2 = e(X2 + ai — Да), (ЫЩ (t + т)) = D bj Ь(т),

(8)

где D - интенсивность белого шума, bj - символ Кронекера. Для фиксированных параметров е, ai и Да использовались те же значения, что и в п. 3.1: е = 0.02, ai = 0.995, Да = 0.01.

Как показали результаты численного счета, при наличии флуктуаций биения в системе двух связанных релаксационных осцилляторов и, соответственно, несовпадение средних периодов колебаний 1-го и 2-го осцилляторов сохраняются и при

р > рт>п (рис. 8 a, б).

Скорость возрастания среднего периода биений (Ть) с ростом параметра р тем выше, чем меньше интенсивность шума. Так, при D = 10_8 величина (Ть) возраста-

23

Рис. 8. Зависимости от параметра связи P средних периодов колебаний {T1,2} первого и второго осцилляторов (а), среднего периода биений {Tb} (б), средней задержки {т12} момента возникновения спайка второго осциллятора относительно момента возникновения спайка первого осциллятора (в) при различных интенсивностях флуктуаций D: 0 - детерминированный случай (сплошная линия); 10~8 (ромбы); 10~6 (поперечные штрихи); 10~4 (штриховая линия). {Ts}, {Ta} - детерминированные периоды колебаний в синфазном и противофазном режиме синхронизации, соответственно. {x12s}, (ri2a} - детерминированные значения задержки в синфазном и противофазном режиме синхронизации, соответственно. В целях лучшей читаемости рисунков опущены зависимости {T1>2(P)} и {x12(P)} при D = 10~8. Видно, что в случае D = 10~6 при P > PLin зависимости {T1>2(P)} и {т12(в)} практически совпадают с соответствующими кривыми для детерминированного случая в синфазном режиме

ет до значений, превосходящих 105 при изменении P на практически неразличимом на рисунке интервале порядка 10_5 справа от P^m, при D = 10_6 этот интервал увеличивается до значений порядка Др = 0.001 и, наконец, при D = 10_4 занимает практически всю область от Pmin до Ртах (то есть включает диапазон значений Р, в пределах которого реализуется бистабильность синфазного - противофазного режимов колебаний в отсутствие флуктуаций). Интересно отметить, что при наличии относительно больших флуктуаций (D = 10_4) средние периоды колебаний как первого, так и второго осциллятора существенно возрастают (почти на 16% у осциллятора с параметром a = 0.995 при P < 10_3). Данный эффект связан с существенным изменением скорости движения изображающей точки осциллятора при флуктуационно-индуцированном отклонении от детерминированной траектории в область более медленного движения. Он аналогичен известному эффекту задержки распада неустойчивого состояния системы [10,11] под воздействием флуктуаций.

Происходящие под воздействием шума переключения между синфазным и противофазным режимами синхронизации проявляются не только в явлении биений, но и в изменении среднего временного интервала между моментами активации элементов (моментами перехода x\(t) и x2(t) через значение x = 0). На рис. 8, в изображена зависимость от силы связи P величины {т\2) - среднего значения разности т\2 момента перехода через 0 переменной xi(t) и предшествующего перехода через 0 переменной x2(t). В отсутствие флуктуаций четко различаются: линия, отражающая

24

однозначную зависимость (х12) = т12ь в области биений (0, Р^т), гДе значение х12 меняется с течением времени от 0 до T2, поэтому (Х12) — T2/2; кривая зависимости х12(|3) = т12а в области (P^in, Ртах), где реализуется либо только противофазная синхронизация (при pmin < в < Pmin), либо бистабильность синхронных режимов (при Pmin < P < Pmax); кривая зависимости Х12(P) = Ti2s в области (P > Pmin), где реализуется либо бистабильность синхронных режимов (при Pmin < P < Pmax), либо только синфазная синхронизация (при P > Pmax).

Как видно из рис. 8, а, в, переключения между сосуществующими режимами синхронизации, вызванные шумом малой интенсивности (D < 10_6 для исследуемых значений детерминированных параметров) приводят к срыву менее устойчивого по отношению к флуктуациям противофазного режима синхронизации. В результате подавляющую часть времени реализуется более устойчивый синфазный режим синхронизации, что можно трактовать как флуктуационно-индуцированное устранение бистабильности шумом малой интенсивности. Отсутствие в случае D = 10_6 уже при P = 0.0025 переключений между различными режимами синхронизации и, таким образом, высокую стабильность синхронного режима по отношению к слабым флуктуациям подтверждает гистограмма значений т12 (рис. 9, а). Наоборот, флуктуации относительно большой интенсивности (D > 10_6) приводят к относительно быстрым срывам не только противофазного, но и более устойчивого синфазного режима синхронизации, как бы восстанавливая бистабильные свойства системы. При этом усредненная задержка Х12 между спайками 1-го и 2-го элементов имеет величину, промежуточную между соответствующими значениями для синфазного и антифазного режимов синхронизации. Интересно отметить, что эта задержка может сохранять некоторое промежуточное значение и при значениях параметров, соответствующих отсутствию бистабильности синхронных режимов в детерминированном случае (в области Pmin < P < Pmin при D > 10“6, а также при P > Pmax, D > 10“4).

Рис. 9. Гистограммы значений задержки т12: при D = 10 6 и P = 0.0025 (а); при D = 10 6 и P = 0.0014 (б); при D = 10-4 и P = 0.0227 (в)

25

Чтобы разобраться в причинах данного явления, были построены гистограммы значений т12 при в = Pi = 0.0014, p^nin < Pi < P^in в случае D = 10-6

(рис. 9, б), и при p = р2 = 0.0227, р2 > Ртах в случае D = 10-4 (рис. 9, в).

Как видно из этих рисунков, на гистограммах четко выделяются два максимума. Один из максимумов соответствует задержке Т12, характерной для существующего устойчивого при данных значениях параметров в детерминированном случае режима синхронизации (противофазного при Р = р1 и синфазного при Р = p2). Другой максимум находится вблизи значения Т12, характеризующего другой режим синхронизации при бифуркационных значениях параметров, связанных с исчезновением этого режима синхронизации (синфазного при Р = pmin и противофазного при Р = pmax). Таким образом, при значениях параметров вне детерминированной области бистабильности синхронных режимов может наблюдаться индуцированная шумом бимодальность плотности вероятности распределения задержек Т12, свидетельствующая о флуктуационно-индуцированной перемежаемости режимов синхронизации. При Ргшт < Р < РПт флуктуационно-индуцированная перемежаемость сопровождается проскоками фaзы и приводит, как уже отмечалось, к биениям. Однако в случае относительно большого шума (D = 10-4) флуктуационно-индуцированная перемежаемость режимов синхронизации наблюдается и при Р > Ртах, когда биения отсутствуют, и средние частоты 1-го и 2-го осцилляторов практически совпадают (рис. 8, а, б). То есть в данном случае, в отличие от индуцированных шумом биений при малых значениях Р, происходят индуцированные шумом переходы между режимами без нарушения взаимной синхронизации осцилляторов.

Заключение

В ходе работы проанализирована система из двух связанных осцилляторов Бонхоффера-ван дер Поля. Экспериментально обнаружено, что при определенных значениях параметров в зависимости от начальных условий в системе могут устанавливаться синфазный и противофазный синхронные режимы. Аналитически доказано существование области параметров, соответствующих мультистабильности синхронных режимов. При сильной связи в системе реализуется только синфазный синхронный режим. При близких к нулю коэффициентах диффузионной связи Р между элементами в ансамбле двух осцилляторов наблюдаются биения. При малых значениях параметра связи в системе возможны как хаотические колебания, так и (п : т)-синхронизация (в узких диапазонах параметра Р). При Р = РП^ происходит рождение противофазного режима, а при Р = РП^ > РПт - также и синфазного режима 1 : 1 синхронизации элементов.

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

Данная работа выполнена при финансовой поддержке РФФИ, гранты 08-02-01259-а, 08-02-92004 и 08-02-97049.

26

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

1. Maurer J. and Libchaber A. // J.Phys (France) Lett. 1982. Vol. 41. P. 515.

2. Brun E., Derighette B., Meier D., Holzner R., and Raveni M. // J.Opt.Soc.Am. B.

1985. Vol. 2. P. 156.

3. DangoisseD., GlorieuxP., andHennequin D. //Phys. Rev. A. 1987. Vol. 36. P. 4775.

4. Thompson J.M.T. and Stewart H.B. Nonlinear Dynamics and Chaos. Wiley, Chichester,

1986.

5. Foss J., Longtin A., Mensour B., and Milton J. // Phys. Rev. Lett. 1996. Vol. 76.

P. 708.

6. Simonotto E., Riani M., Seife C., Roberts M., Twitty J., and Moss F. // Phys. Rev. Lett. 1997. Vol. 78. P. 1186.

7. Bonhoeffer K.F. // Naturwissenschaften. 1953. Vol. 40. P. 301.

8. Torre V // J. Theor. Biol. 1976. Vol. 61. P. 55.

9. Матросов В.В. Динамика нелинейных систем. Программный комплекс для исследования нелинейных динамических систем с непрерывным временем. Н. Новгород, ННГУ, 2002.

10. Mantegna R.N. and Spagnolo B. // Phys. Rev. Lett. 1996. Vol. 395, № 4. P. 563.

11. Malakhov A.N. and Agudov N.V // Phys. Rev. E. 1999. Vol. 60. P. 6333.

Поступила в редакцию 15.08.2008 После доработки 1.12.2008

VARIETY OF SYNCHRONOUS REGIMES IN ENSEMBLES OF NONIDENTICAL OSCILLATORS:

Two coupled elements

A.K. Kryukov, G.V Osipov, A.V. Polovinkin

We study synchronization of two coupled nonidentical Bonhoeffer-van der Pol oscillators. Coexistence of two different synchronous regimes is proved. Mechanisms of synchronous regimes origination and destruction are investigated. Fluctuations influence on syncronous regimes is considered. It is found that noise can cause: i) synchronization destruction and beating originations; ii) fluctuations-caused bistability destruction; iii) fluc-tuations-caused intermittency of synchronous regimes without synchronization destruction.

Keywords: Synchronization, multistability, synchronous regime, numerical methods, modeling, stochastic system, FitzHugh-Nagumo, Bonhoeffer-van der Pol.

27

Крюков Алексей Константинович - родился в 1984 году, окончил Нижегородский государственный университет им. Н.И. Лобачевского (2007). Аспирант кафедры теории управления и динамики машин факультета ВМК ННГУ. Область научных интересов: структуры синхронизации в нелинейных системах, моделирование динамики биологических сред. Стипендиат некоммерческого фонда «Династия» (2006). Имеет 12 публикаций (в сооавторстве) в том числе 2 журнальных. Лауреат стипендии имени академика Г.А. Разуваева (2008).

E-mail: [email protected]

Нижегородский государственный университет им. Н.И. Лобачевского 603950, Россия, г. Нижний Новгород, пр. Гагарина, 23

Осипов Григорий Владимирович - родился в Горьком (1960), окончил Горьковский государственный университет им. Н.И. Лобачевского (1982). В 19861988 годах учился в аспирантуре Научного совета по комплексной проблеме «Кибернетика» (НСК) АН СССР. Защитил диссертацию на соискание ученой степени кандидата физико-математических наук в НСК (1988) по математическому моделированию и доктора физико-математических наук в ННГУ (2004) по радиофизике. С 1988 года работает в ННГУ, с 2007 года в качестве заведующего кафедрой «Теории управления и динамики машин» факультета ВМК. Один из авторов монографий «Устойчивость, структуры и хаос в нелинейных сетях синхронизации» и «Synchronization in oscillatory networks». Опубликовал 120 научных статей (в том числе 2 обзора) по теории колебаний и волн и математическому моделированию.

E-mail: [email protected]

Нижегородский государственный университет им. Н.И. Лобачевского 603950, Россия, г. Нижний Новгород, пр. Гагарина, 23

Половинкин Андрей Владимирович - родился в 1959 году, окончил Нижегородский государственный университет в 1981 году. Докторант кафедры математики радиофизического факультета ННГУ, кандидат физико-математических наук (1987). Область научных интересов: флуктуационно-индуцированные процессы в нелинейных системах и неоднородных средах. Соавтор свыше 40 научных и 5 учебно-методических работ. Участник российских и международных научных проектов.

E-mail: [email protected]

Нижегородский государственный университет им. Н.И. Лобачевского 603950, Россия, г. Нижний Новгород, пр. Гагарина, 23

28