CC BY
28
УДК 530.182
Д.И. Данилов, М.О. Журавлев, А.В. Иванов, А.А. Короновский,
М.К. Куровская, О.И. Москаленко, А.О. Сельский, С.А. Шурыгина
ИССЛЕДОВАНИЕ ХАРАКТЕРИСТИК ПЕРЕМЕЖАЮЩЕГОСЯ ПОВЕДЕНИЯ ВБЛИЗИ ГРАНИЦЫ СИНХРОННОГО РЕЖИМА ПРИ ИЗМЕНЕНИИ УПРАВЛЯЮЩЕГО ПАРАМЕТРА
Исследуются характеристики перемежающегося поведения, наблюдающегося вблизи границы синхронизации осциллятора Ван дер Поля под внешним гармоническим воздействием в случае, когда амплитуда внешнего сигнала изменяет свое значение случайным образом через некоторые интервалы времени.
Хаотическая синхронизация, перемежаемость, хаотический осциллятор
D.I. Danilov, M.O. Zhuravlev, A.V. Ivanov, A.A. Koronovsky,
M.K. Kurovskaya, O.I. Moskalenko, A.O. Selskii, S.A. Shurygina
RESEARCH INTO THE CHARACTERISTICS OF INTERMITTENT BEHAVIOR NEAR THE SYNCHRONIZATION BOUNDARY FOR THE VARYING CONTROL PARAMETER
Characteristics of intermittent behavior near the synchronization boundary are studied for the Van der Pole oscillator under the external harmonic signal where the amplitude is changed randomly at the given time intervals.
Chaotic synchronization, intermittency, chaotic oscillator
Явление перемежаемости широко встречается в природе и весьма активно изучается исследователями. Одной из важных особенностей перемежаемости является ее фундаментальный характер, приводящий к тому, что явление перемежающегося поведения наблюдается в системах самой различной природы. Известно, что перемежаемость является одним из классических сценариев перехода от периодических колебаний к хаотическим [1, 2]. На сегодняшний день существует определенная классификация перемежающегося поведения, в частности, выделяют перемежаемость типов I-III [1-3], on-off перемежаемость [4, 5], перемежаемость кольца [6], перемежаемость «игольного ушка» [7], тесно связанную с перемежаемостью типа I в присутствии шума [8].
Как правило, в подавляющем большинстве случаев перемежающееся поведение, наблюдающееся в той или иной системе, рассматривается при фиксированных во времени (стационарных) значениях управляющих параметров. Однако, вполне возможны ситуации, когда параметры системы, демонстрирующей перемежающееся поведение, будут изменяться с течением времени. В частности, подобная ситуация может быть характерна для живых систем. Например, известно, что перемежающееся поведение наблюдается вблизи границы возникновения режимов хаотической синхронизации [9]. Известно также, что режимы хаотической синхронизации (например, фазовой хаотической синхронизации) имеют место в живых системах, например, в кардиоваскулярной системе человека [10]. Очевидно, что на границе установления синхронного режима будет наблюдаться перемежающееся поведение, при этом на его характеристики будет оказывать существенное влияние изменение управляющих параметров системы с течением времени, присущее живым системам.
В настоящей работе рассматривается вопрос о характеристиках перемежающегося поведения осциллятора вблизи границы синхронного режима, если один из управляющих параметров системы будет менять свое значение случайным образом через некоторые промежутки времени. Пусть A -управляющий параметр, Ac - значение этого управляющего параметра, соответствующее границе установления синхронного режима. Будем считать, что при A > Ac в системе наблюдается синхронная динамика, а при A < Ac - асинхронное поведение, демонстрирующее черты перемежаемости. Предположим, что значение управляющего параметра находится вблизи критической точки Ac и подвержено случайным флуктуациям с течением времени, причем в разные
моменты времени величина A может быть как выше, так и ниже критического значения Ac .
В качестве исследуемой системы был выбран осциллятор Ван дер Поля, находящимся под внешним гармоническим воздействием, амплитуда которого A изменялась случайным образом:
x — (X - x )x + x = (Aq + D^)sin(wet), (1)
где A и - амплитуда и частота внешнего гармонического воздействия, соответственно, X -
параметр нелинейности, D - величина, характеризующая интенсивность флуктуаций, Ъ - случайная величина, характеризующаяся нормальным законом распределения
1 Г ъ21
Р(Ъ) = /г- exP — ^ , (2)
где а2 = 1. Значения указанных выше параметров были выбраны фиксированными: X = 0.1, =0.98, Aq =0.024, что больше граничного значения Ac =0.0238, соответствующего порогу синхронизации в отсутствии случайной добавки D^.
В качестве упрощающего предположения будем считать, что амплитуда внешнего воздействия изменяет свое значение случайным образом через фиксированные (равные) интервалы времени. Параметр, отвечающий за интенсивность флуктуаций D = 5 • 10-4, величина Ъ, а соответственно и величина A = Aq + D'b, отвечающая за амплитуду внешнего воздействия, принимала новое значение через интервал времени = 900. Очевидно, для плотности распределения вероятности p( A) величины A можно записать
p( A) = ^^exp л2па
2а
2
(3)
при этом а = Dа^ =5-10 4,ц = A0 = 0.024.
Рассмотрим, что происходит в системе между флуктуациями управляющего параметра A с точки зрения возникновения фазовых проскоков, соответствующих возникновению турбулентных участков поведения. Очевидно, что в случае, когда Л > Лс, система оказывается в синхронизованном состоянии, и, как следствие, в этот интервал времени длительностью т^ фазовых проскоков не
возникает и рассматриваемый интервал времени полностью принадлежит участку ламинарной (синхронной) динамики. Ситуация оказывается сложнее, если в результате флуктуации управляющий параметр Л принял значение, меньшее Лс. Хорошо известно, что в случае малой расстройки
параметров при постоянном значении параметра Л ниже границы синхронизации наблюдается поведение, демонстрирующее черты перемежаемости типа I, при этом зависимость средней длительности ламинарных участков поведения (участков синхронной динамики) для данного типа перемежающегося поведения имеет вид
Т(Л) = С1 - (Лс - Л)-1/2, (4)
где С - некоторая константа. В рассматриваемом случае для указанных значений управляющих параметров, С1 ~ 26. Очевидно, что если средняя длительность Т(Л) участка ламинарного поведения оказывается существенно меньше длительности интервала времени т^, в течение которого 28
значение управляющего параметра Л остается неизменным, то с вероятностью почти наверное в системе будет наблюдаться фазовый проскок, соответствующий окончанию участка ламинарного
(синхронного) поведения. Введем в рассмотрение значение Л*, которое соответствует такому
значению параметра Л, при котором средняя длительность ламинарной фазы Т(Л*) = т^. Из соотношения (4) нетрудно получить, что
Л* = Лс -(С^)2. (5)
Таким образом, если т^ = 900, то А ;0.023. Нетрудно заметить, что при всех Л < Л*, Т(Л) < т^, и, следовательно, можно считать, что до смены значения параметра Л в рассматриваемой системе с вероятностью, близкой к единице, произойдет фазовый проскок. В том же случае, когда значение управляющего параметра оказывается в диапазоне А* < А < Лс, вероятность фазового проскока может быть оценена как т^/Т(Л). Таким образом, для вероятности наблюдения фазового
проскока в течение интервала времени между флуктуациями значения управляющего параметра Л можно записать
1, Л < Л*,
Р( Л) =
, Л* < Л < Л, (6)
Т (Л) с
0 л > л .
На рис. 1 представлено распределение величины управляющего параметра Л для рассматриваемой системы (кривая 1), полученное в результате численного моделирования системы (1).
Одновременно на этом же рисунке показано распределение значений того же самого параметра Л , при котором наблюдался фазовый проскок (кривая 2). Оба распределения получены для одного и того же временного ряда по всем значениям, которые принимал параметр Л за время наблюдения. На рисунке также показаны стрелками значения Л* и Лс , играющие важную роль в формировании распределения плотности вероятности фазового проскока в рассматриваемой системе. Из рисунка отчетливо видно, что для значений параметра Л < Л* оба распределения совпадают друг с другом, что говорит о том, что для данных значений управляющего параметра в 100% случаев происходил фазовый проскок, то есть в системе начинался турбулентный участок поведения, и, соответственно, вероятность фазового проскока в данном случае равняется единице, в полном соответствии с соотношением (6). Точно так же наблюдается полное соответствие результатов численного моделирования для случаев, когда Л > Лс - в этой ситуации фазовые проскоки практически
отсутствуют, что говорит о том, что для значений параметра, принимающих в ходе флуктуаций значения больше порога установления синхронного режима, вероятность возникновения фазового проскока равняется нулю. Наконец, видно, что в диапазоне значений Л* < Л < Ас с ростом значения
Л количество фазовых проскоков уменьшается, что говорит об уменьшении вероятности Р. Для сопоставления поведения вероятности Р возникновения фазовых проскоков, полученной численно, с аналитической зависимостью (6) в диапазоне Л* < А < Ас, на рисунке приведена также зависимость
р( А)АЛР( А), показанная пунктирной линией (АЛ - диапазон агрегации представления данных).
Как видно из рисунка, и в этом диапазоне значений параметра Л соотношение (6) хорошо описывает процессы, происходящие в анализируемой системе.
С учетом всего вышеизложенного, можно говорить о том, что при рассмотрении динамики системы (1) в интервале времени т^ между флуктуациями значения параметра Л вероятность наблюдать уход фазы будет определяться соотношением
“ Л* Лс т р(Л)
р = \р(Л)Р(Л)йЛ = \р(Л)йЛ + \ йА. (7)
0 0 Л* Т (Л)
ЩА)
1000
600
200
0.022 0.023 0.024 0.025 А
Рис. 1: Плотность распределения вероятности р(А) значений параметра А для системы (1).
Распределение 1 построено по всем значениям, которые принимал параметр А за время наблюдения,
распределение 2 получено для значений параметра А, при которых наблюдался фазовый проскок.
Обе зависимости получены для одного и того же временного ряда по всем значениям,
которые принимал параметр А за время наблюдения. Штриховой линией показана зависимость р(А)ААР(А)
Тогда вероятность того, что в системе реализуется ламинарная фаза длительностью £ , соответствующей к интервалам флуктуаций параметра А (иными словами, £ = кт*), будет определяться как
Р2(к) = (1 - Р1)к-1 • Р2, (8)
где вероятность Р определяется соотношением (7). Соответственно, распределение длительностей ламинарных участков поведения, характеризующихся синхронной динамикой во времени, в этом случае будет подчиняться закону
*/Тс — 1 п
N(£) ~ (1 — Р) * • р2, (9)
и, соответственно, данное соотношение может быть переписано в виде
N(5) = С2 ехр(кх), (10)
где
к = (1/т* )1п(1 — Р), (11)
С2 - коэффициент пропорциональности, являющийся нормировочным множителем,
обеспечивающим выполнение условия нормировки. Таким образом, следует ожидать, что в случае медленных флуктуаций управляющего параметра А распределение длительностей ламинарных фаз в рассматриваемой системе будет подчиняться экспоненциальному закону (10). В случае, если распределение длительностей ламинарных участков поведения строить для длительностей ламинарных фаз, нормированных на интервал времени т*, показатель степени К в соотношении (10)
будет определяться исключительно величиной вероятности наблюдать уход фазы за интервал времени между флуктуациями управляющего параметра А:
К = 1п(1 — Р). (12)
На рис. 2 приведены распределения длительностей ламинарных фаз, полученные в результате численного моделирования рассматриваемой системы при указанных выше значениях параметров, для различных значений интервала т*, в течение которого управляющий параметр А считается неизменным. Для сопоставления с теоретическими данными на рисунке приведены также зависимости, соответствующие аналитическому соотношению (9), значения величины Р1 для которых определялись с помощью соотношения (7). Точки на графике соответствуют данным,
полученным численно, сплошными линиями показана зависимость (9) для каждого конкретного значения т^. Нетрудно видеть хорошее соответствие численных данных и теоретического
соотношения, независимо от величины интервала времени т^, по истечении которого происходило
изменение значения параметра А, моделирующее медленную флуктуацию данного параметра.
т
0.1
0.01
0.001
0 40 80 $/ц
Рис. 2: Распределения длительностей ламинарных фаз для системы (1), полученные при различных значениях т^. Точки на графике соответствуют численным данным, аналитическое
соотношение (9) показано сплошной линией. Значения интервала времени т^ и вероятности P , соответствующие каждому конкретному распределению следующие:
1 - т^ =1200, р =0.247 ; 2 - т^ =900,р =0.197 ; 3 - т^ = 600, р =0.133 ;
4 - т^ = 300, р = 0.066. Значения по оси ординат отложены в логарифмическом масштабе, чтобы показать экспоненциальный характер распределения длительностей участков ламинарного поведения
Таким образом, в настоящей работе исследованы характеристики перемежающегося поведения, наблюдающегося близи границы синхронизации осциллятора ван дер Поля, находящегося под внешним гармоническим воздействием в случае, когда амплитуда этого воздействия изменяла свое значение случайным образом через фиксированные интервалы времени. Был установлен аналитический закон для распределения длительностей ламинарных участков поведения и показано хорошее соответствие выявленной теоретической зависимости с данными, полученными в ходе численного моделирования системы при различных значениях интервала времени, через который происходила смена значения управляющего параметра.
Исследование выполнено при поддержке Министерства образования и науки Российской Федерации (соглашения 14.B37.21.0569 и 14.B37.21.0751) и Фонда некоммерческих программ «Династия».
ЛИТЕРАТУРА
1. Berge P. L'ordre dans le chaos / P. Berge, Y. Pomeau, and Ch. Vidal Hermann, Paris, 1988.
2. Берже П. Порядок в хаосе / П. Берже, И. Помо, К. Видаль. М.: Мир, 1991.
3. Dubois M. Experimental evidence of intermiasttencies associated with a subharmonic bifurcation /
M. Dubois, M. Rubio, P. Berge // Phys. Rev. Lett. 1983. V 51. P. 1446-1449.
4. Platt N. On-off intermittency: a mechanism for bursting / N. Platt, E. A. Spiegel, and C. Tresser // Phys. Rev. Lett. 1993. V. 70. P. 279-282.
5. Heagy J. F. Characterization of on-off intermittency / J.F. Heagy, N. Platt, S.M. Hammel // Phys. Rev. E. 1994. V. 49. P. 1140-1150.
6. Hramov A. E. Ring intermittency in coupled chaotic oscillators at the boundary of phase synchronization / A.E. Hramov, A.A. Koronovskii, M.K. Kurovskaya, S. Boccaletti // Phys. Rev. Lett. 2006. V. 97. P. 114101.
7. Pikovsky A.S. Attractor-repeller collision and eyelet intermittency at the transition to phase synchronization / A.S. Pikovsky, G.V. Osipov, M.G. Rosenblum, M. Zaks, J. Kurths // Phys. Rev. Lett. 1997. V. 79. P. 47-50.
8. Hramov A.E. Type-I intermittency with noise versus eyelet intermittency / A.E. Hramov, A.A. Koronovskii, M.K. Kurovskaya, O.I. Moskalenko // Phys. Lett. A. 2011. V. 375. P. 1646-1652.
9. Rosenblum M.G. From phase to lag synchronization in coupled chaotic oscillators / M.G. Rosenblum, A.S. Pikovsky, J. Kurths // Phys. Rev. Lett. 1997. V. 78. P. 4l93-4196.
10. Hramov A.E. Detecting synchronization of self-sustained oscillators by external driving with varying frequency / A.E. Hramov, A.A. Koronovskii, V.I. Ponomarenko, M.D. Prokhorov // Phys. Rev. E. 2006. V. 73. P. 026208.
Данилов Дмитрий Игоревич -
аспирант кафедры «Физика открытых систем» факультета нелинейных процессов Саратовского государственного университета имени Н.Г. Чернышевского
Журавлев Максим Олегович -
младший научный сотрудник научнообразовательного центра «Нелинейная динамика сложных систем»
Саратовского государственного технического университета имени Гагарина Ю.А.
Иванов Алексей Викторович -
кандидат геолого-минералогических наук, профессор, заведующий кафедрой «Г еоэкология и инженерная геоэкология» декан факультета экологии и сервиса Саратовский государственный технический университет имени Гагарина Ю.А.
Короновский Алексей Александрович -
доктор физико-математических наук, профессор, ведущий научный сотрудник научно-образовательного центра «Нелинейная динамика сложных систем» Саратовского государственного технического университета имени Гагарина Ю.А.
Куровская Мария Константиновна -
кандидат физико-математических наук, доцент кафедры «Нелинейная физика» Саратовского государственного университета имени Н.Г. Чернышевского
Москаленко Ольга Игоревна -
кандидат физико-математических наук, доцент, старший научный сотрудник научно-образовательного центра «Нелинейная динамика сложных систем» Саратовского государственного технического университета имени Гагарина Ю.А.
Сельский Антон Олегович -
младший научный сотрудник научнообразовательного центра «Нелинейная динамика сложных систем»
Саратовского государственного технического университета имени Гагарина Ю.А.
Шурыгина Светлана Андреевна -
младший научный сотрудник научнообразовательного центра «Нелинейная динамика сложных систем»
Саратовского государственного технического университета имени Гагарина Ю.А.
Dmitry I. Danilov -
Postgraduate
Department of Physics of Open Systems,
Gagarin Saratov State University
Maxim O. Zhuravlev -
Junior Researcher Scientific and Education Center «Nonlinear Dynamics of Complex Systems», Gagarin Saratov State Technical University
Alexey V. Ivanov -
Ph. D., Professor
Head: Department of Geoecology
and Engineering Geoecology,
Dean: Faculty of Ecology and Service,
Gagarin Saratov State Technical University
Alexey A. Koronovsky -
Dr. Sc., Professor, Leading Scientific Researcher: Scientific and Education Center «Nonlinear Dynamics of Complex Systems», Gagarin Saratov State Technical University
Maria K. Kurovskaya -
Ph. D., Associate Professor Department of Nonlinear Physics,
Faculty of Nonlinear Processes, Chernyshevsky Saratov State University
Olga I. Moskalenko -
Ph. D., Associate Professor,
Senior Researcher:
Scientific and Education Center «Nonlinear Dynamics of Complex Systems», Gagarin Saratov State Technical University
Anton O. Selsky -
Junior Researcher:
Scientific and Education Center «Nonlinear Dynamics of Complex Systems», Gagarin Saratov State Technical University
Svetlana A. Shurygina -
Junior Researcher:
Scientific and Education Center «Nonlinear Dynamics of Complex Systems», Gagarin Saratov State Technical University
Статья поступила в редакцию 12.12.12, принята к опубликованию 20.02.13
