АВТОВОЛНОВЫЕ СТРУКТУРЫ В ДВУМЕРНЫХ РЕШЕТКАХ НЕЛОКАЛЬНО СВЯЗАННЫХ ОСЦИЛЛЯТОРОВ

УДК 53.043 https://doi.org/10.18500/0869-6632-2020-28-3-299-323

Автоволновые структуры в двумерных решетках нелокально связанных осцилляторов

А. В. Бух, Е. В. Рыбалова, В. С. Анищенко

Саратовский национальный исследовательский государственный университет имени Н.Г. Чернышевского Россия, 410012 Саратов, Астраханская, 83 E-mail: buh.andrey@yandex.ru, rybalovaev@gmail.com, wadim@info.sgu.ru Автор для переписки Вадим Семенович Анищенко, wadim@info.sgu.ru Поступила в редакцию 16.02.2020, принята к публикации 15.04.2020, опубликована 30.06.2020

Цель настоящего исследования заключается в сравнении динамики спиральных и концентрических структур, в том числе химерных, в ансамблях с различающимися парциальными элементами. Для этого выполняется численное моделирование автоволновых структур в двумерных ансамблях связанных генераторов ван дер Поля и отображений Рулькова. Рассмотрены случаи локальной и нелокальной связи между индивидуальными осцилляторами ансамблей. Методы. Эволюция динамики решетки отображений Рулькова в численном эксперименте напрямую определяется соответствующими рекуррентными соотношениями. Ансамбль осцилляторов ван дер Поля интегрируется методом Гюна. По результатам расчетов в обоих случаях строятся мгновенные значения амплитуд для всех элементов решетки, пространственно-временные диаграммы для её сечения, фазовые портреты и временные реализации для отдельных элементов решетки. Кроме того, вычисляются средние значения частот для всех парциальных элементов и зависимости мгновенных частот от времени для выбранных отображений или осцилляторов. Результаты сравниваются. Результаты. Показано, что при локальной связи в ансамблях реализуются режимы классических спиральных и концентрических волн. С введением нелокальной связи при вариации радиуса связи получены более сложные структуры, включая химерные. Описаны спирально-волновые химерные структуры с одним и более ядрами некогерентности. Показана возможность реализации нового типа химерной структуры на основе концентрических волн - концентрических химер. Представлены результаты анализа свойств ядер некогерентности для спиральных и концентрических химерных структур. Обсуждаются результаты реакции динамики ансамблей на изменения силы нелокальной связи и внешнее шумовое воздействие. Заключение. В обоих ансамблях для концентрических химер средние значения частот почти одинаковы для всех элементов. В случае спиральных химер средние значения частот в ядре отличаются от средних значений в когерентной области. В решетке генераторов ван дер Поля возможен переход от концентрических волн к спиральным при воздействии шумом достаточно большой интенсивности. В ансамбле отображений Рулькова реализуется индуцированный шумом переход от спиральной волны к концентрической волне.

Образец цитирования: Бух А.В., Рыбалова Е.В., Анищенко В.С. Автоволновые структуры в двумерных решетках нелокально связанных осцилляторов//Известия вузов. ПНД. 2020. T. 28, № 3. С. 299-323. https://doi.org/10.18500/0869-6632-2020-28-3-299-323

Статья опубликована на условиях лицензии Creative Commons Attribution License (CC-BY 4.0).

Финансовая поддержка. Исследование выполнено при финансовой поддержке Немецкого Физического Общества (проект SFB 910). Е.В.Р. и В.С.А. благодарят за финансовую поддержку РФФИ и DFG в рамках исследовательского проекта №20-52-12004 (получение результатов, написание текста), А.В.Б. благодарит за финансовую поддержку РФФИ в рамках научного проекта № 19-32-90005 (анализ результатов).

https://doi.org/10.18500/0869-6632-2020-28-3-299-323

Autowave structures in two-dimensional lattices of nonlocally coupled oscillators

A. V. Bukh, E. V. Rybalova, VS. Anishchenko

Saratov State University 83 Astrakhanskaya str., Saratov 410012, Russsia E-mail: buh.andrey@yandex.ru, rybalovaev@gmail.com, wadim@info.sgu.ru Correspondence should be addressed to Vadim S. Anishchenko, wadim@info.sgu.ru Received 16.02.2020, accepted 15.04.2020, published 30.06.2020

Objective. The aim of the research was to compare the dynamics of spiral and target structures including the dynamics of chimera states in ensembles with different nodes. Numeric simulations for autowave structures in two-dimensional ensembles of coupled van der Pol oscillators and Rulkov&s maps were performed. Cases of local and nonlocal coupling between ensemble nodes were considered. Methods. The evolution dynamics of Rulkov&s map lattice is strictly defined with corresponding recurrent formulae. The ensemble of van der Pol oscillators was integrated using the Heun&s method. Snapshots of amplitude values, spatio-temporal diagrams for corresponding sections, phase portraits and time series for single elements were constructed. Moreover, mean values of frequencies of all nodes, and dependencies of instantaneous frequencies on time for selected maps and oscillators. Results were compared. Results. It is shown that classical spiral and target wave regimes were realized at local coupling. More complex structures including chimera states were obtained when a nonlocal coupling was included. Spiral wave chimera structures with single and several incoherent cores were described. A new chimera structure based on target waves (target wave chimera) was shown to be possible. The analysis of features of incoherent cores for both spiral wave chimera and target wave chimera were presented. The results of coupling parameter effect and external noise influence on the dynamics of the ensembles were discussed. Conclusion. Mean values of frequencies were almost the same for all the elements in both ensembles at target wave chimera regimes. The mean values in the incoherent core differed from the ones in coherent area at spiral wave chimera regimes. A transition from target waves to spiral waves was possible when a noise with sufficiently large intensity was introduced. Noise-induced transition from spiral waves to target waves took place in the ensemble of Rulkov&s maps.

Reference: Bukh A.V., Rybalova E.V., Anishchenko V.S. Autowave structures in two-dimensional lattices of nonlocally coupled oscillators. Izvestiya VUZ. Applied Nonlinear Dynamics, 2020, vol. 28, no. 3, pp. 299-323. https://doi.org/10.18500/0869-6632-2020-28-3-299-323

This is an open access article distributed under the terms of Creative Commons Attribution License (CC-BY 4.0).

Acknowledgements. This work was funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) -Projektnummer. 163436311-SFB 910. E.V.R. and V.S.A. thank for the financial support provided by RFBR and DFG according to the research project #20-52-12004 (obtaining results, writing the manuscript), A.V.B. thanks for the financial support provided by RFBR, project number 19-32-90005 (analysis the results).

Введение

Одним из актуальных направлений исследований в нелинейной физике является анализ автоволновых структур в нелинейных средах и их моделях. Под автоволнами понимают самоподдерживающийся волновой процесс в неравновесной среде, остающийся неизменным при достаточно малых изменениях как начальных, так и граничных условий. Математическим аппаратом для описания автоволн чаще всего служат уравнения диффузионного типа с активной нелинейностью [1-11]. Понятие «автоволны» впервые было введено в научную лексику академиком Р.В. Хохловым и затем активно использовалось в исследованиях ученых Института Биофизики АН СССР (г. Пущино) применительно к задачам кардиологии. Оно явилось обобщением определения автоколебаний на случай пространственно распределенных активных нелинейных сред и систем [12,13]. Исследования автоволновых структур были начаты во второй половине

двадцатого века и активно продолжаются в настоящее время [1-42]. Интерес к исследованиям автоволновых процессов в различных моделях активных сред резко возрос в начале двухтысячных годов. Во многом это связано с введением в рассмотрение так называемых химерных структур. Анализ волновых структур проводился, как правило, с использованием математических моделей нелинейной среды в виде уравнений типа реакция - диффузия. Такой подход соответствует моделированию среды в виде ансамблей взаимодействующих активных осцилляторов с локальной связью. Каждый осциллятор связан только с ближайшими соседями в ансамбле. При анализе химерных структур использовалась более сложная, нелокальная связь. При нелокальной топологии связи каждый индивидуальный осциллятор симметрично связывается с различным числом соседних осцилляторов, которое определяется радиусом связи [43]. Использование нелокальной связи позволило найти ряд более сложных пространственно-временных структур, которые не наблюдались при локальной связи [10,32,44-48]. Так, в пионерских работах [24,26] при численном моделировании динамики двумерного ансамбля нелокально связанных фазовых осцилляторов впервые была обнаружена структура, впоследствии названная спирально-волновой химерой (spiral wave chimera, SWC). Важно отметить, что этот результат был получен раньше, чем в работе [46] было дано определение химерной структуры. В настоящей работе предметом изучения выбраны автоволновые структуры в виде спиральных и концентрических волн. Названные структуры были получены методом численного моделирования динамики двумерных решеток нелокально связанных осцилляторов, в качестве которых использовались осцилляторы ван дер Поля и отображения Рулькова [49, 50]. Подробно анализируются случаи локальной и нелокальной связи между осцилляторами решеток. Показано, что при локальной связи в решетках при вариации параметров реализуются только «классические» спиральные или концентрические волны. Увеличение радиуса нелокальной связи порождает более сложные автоволновые структуры, включая химерные. Представлены результаты по SWC с одним и несколькими ядрами некогерентности, а также новый тип химерной структуры, которая рождается на основе концентрической волны при увеличении радиуса нелокальной связи, концентрическая химера (target wave chimera, TWC) [49]. Детально анализируются свойства некогерентных кластеров найденных типов химерных структур (SWC и TWC). Рассматриваются свойства реализуемых автоволновых структур при вариации силы нелокальной связи и воздействии внешнего шума.

Рассмотрим динамику двумерных решеток связанных динамических подсистем с N х N элементами, которые описываются следующими уравнениями:

г к=1 3 1=1

где Вг х В^ - количество связей в матрицах о^, о^д для данного элемента; и^ - переменные, отражающие динамику отдельных осцилляторов; г,] = 1, 2 ..., N - индексы этих осцилляторов в квадратной решетке по соответсвующим направлениям; N х N - количество осцилляторов в решетке; й^ = 1, Ь = 1, 2,.. .Т в случае отображения, и = и^, £ е [0;Т] в случае, когда осциллятор задается системой ОДУ; Т - время, на котором рассматривается поведение системы; Е(и^) определяет динамику отдельного элемента и^; С (и*к и*^ определяет связь между

элементами с индексами (к, I) и (г, ]), С 1, и* ^ = 1 — и* ^ в случае диффузионной

связи и С (и^ , 1, и^^ = , 1) - Е(и^)) в случае инерционной связи. Для нелокальной связи матрицы связи со свободными граничными условиями определяются следующими уравнениями:

Огк =

где R - параметр нелокальности связи, определяющий количество связей отдельного элемента В{ х В^ в решетке. Графическое представление связей в случае нелокальной связи приведено на рис. 1.

В рамках работы рассматриваются решетки с разными типами парциальных элементов. Первая решетка представляет собой систему N х N = = 100 х 100 нелокально связанных осцилляторов ван дер Поля и задается следующими соотношениями в (1):

Рис. 1. Топологическая схема связей ансамбля (1) в нелокальном случае при R = 2

Fig. 1. Scheme of the coupling topology for the ensemble (1) in a nonlocal case with R = 2

иг A = U

fK , ) = (f-(x •> " )W

V fy (xi,j , У г,j )) V

te [0; T], G (ul, u{3 ) =Н — ХгА .

£(1 — Х2 j) • Уг,] — Ы2Хг,3

где £ и ш - параметры нелинейности и частоты отдельных осцилляторов, соответственно. Вторая решетка состоит из N хN = 200x200 отображений Рулькова, описывающих динамику нейронов. В этом случае соотношения в (1) определяются следующим образом:

f X (xi,i , У i.i )

t = 1, 2....T.

F(U )=ff*(xh,yhЛ

) U, ) &

(xi,j, у i,j)

У (xi.i, у i.i )

a/(1 — х\\3 ) + ух, , x\\j < 0

a + у* „•, 0 < х < a + у

id, уы) — fy (хгл > уи)

— 1, xj, > a + yfG (uU, ) = (xi

\\J У (xk,l, УкЛ

fy (x{j, y{j) = y\\,j — + 1) +

Функция /х (х1-^, у\\ ^) описывает процесс генерации спайка, а уравнение /у (х| ^, ) отвечает за относительно медленное изменение режима. Увеличение параметров у и а ведет к увеличению частоты колебаний в отдельном отображении Рулькова.

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

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

Автоволновые структуры, наблюдаемые при локальной связи. Как показали результаты численного моделирования, в зависимости от величины управляющих параметров е и ю при постоянном значении радиуса связи К = 1, в исследуемом ансамбле возможна реализация как спиральных, так и концентрических волновых структур. Исследование решетки (3) в случае локальной связи проводилось с использованием случайных начальных условий для ж0 ^, у0 ^ € [-0.9 : 0.9]. Карта режимов (рис. 2, a) иллюстрирует, что для каждого режима можно

Рис. 2. Карта режимов поведения системы (3) в случае локальной связи (R = 1, а = 0.7) (a). На карте отмечены области синхронных колебаний (Sync); области существования концентрических волн (TW, мгновенный снимок при ю = 0.9, е = 2.15 приведен на фрагменте (с)) и спиральных волн (SW, приведены мгновенные снимки при ю = 0.9, е = 0.5 (d) и £ = 2.15 (е)); точки, в которых возможно наблюдение обоих типов волн в одной структуре (STW, мгновенный снимок при ю = 1.65, £ = 2.1 приведен на фрагменте (b))

Fig. 2. Regime map of system (3) behaviour for local coupling case (R = 1, а = 0.7) (a). Synchronous oscillations (Sync); target waves (TW, snapshot is shown in (c) at ю = 0.9, £ = 2.15); spiral waves (SW, snapshots are shown in (d) ю = 0.9, £ = 0.5 and in (е) at ю = 0.9, £ = 2.15); and structures consisting of both spiral and target waves (STW, snapshot is shown in (b) at ю = 1.65, £ = 2.1) are depicted on the map

выделить конечную область его существования. При малых значениях параметра ш возможно наблюдение только спиральных волн (рис. 2, d, e) и некогерентных режимов. Когда параметр ш увеличивается, возникают концентрические волны (рис. 2, в) при достаточно больших е. Изучение поведения системы на плоскости параметров (е, ш) показывает, что существует возможность наблюдения спиральных и концентрических волн в одной пространственно-временной структуре. Один из таких режимов приведен на рис. 2, Ь.

Рассмотрение режимов, наблюдаемых при фиксированном значении параметра ш = 0.9 и изменении , позволяет обнаружить, что параметр управляет длиной волны как в случае спиральных, так и в случае концентрических волн. На рис. 2, d, e приведены мгновенные снимки спирально-волновых структур, рассчитанных для различных значений параметра возбуждения генератора. Как видно из рис. 2, d, e, увеличение значения параметра е приводит к уменьшению длины спиральной волны. Дальнейшее увеличение параметра е ^ 2.2 приводит к разрушению спиральной волны и переходу к режиму некогерентности. При уменьшении параметра е < 1.65 спиральные волны в решетке сохраняются, и их длина волны продолжает увеличиваться. Пример концентрической волны при е = 2.15 приведен на рис. 2, c. Если уменьшать параметр е для этой структуры, её длина волны будет увеличиваться, так же, как длина спиральной волны при уменьшении параметра е. Затем концентрические волны переходят в профиль, отвечающий режиму синхронизации, при е ^ 1.6. При увеличении параметра е концентрические волны разрушаются при е ^ 2.4. Таким образом, параметр е позволяет управлять длиной волны обеих структур: длины как спиральной, так и концентрической волн обратно пропорциональны значению этого параметра.

Исследования показали наличие мультистабильности в системе (3). Этот эффект проявляется в том, что при вариации начальных условий или смене направления движения по параметру можно для одних и тех же фиксированных значений управляющих параметров получать различные пространственно-временные структуры. В рамках исследования сначала были получены спиральная и концентрическая структуры при = 2, а затем параметр изменялся с сохранением начальных условий. В результате было выяснено, что в диапазоне е € [1.65; 2.2] в решетке наблюдаются как спиральные, так и концентрические волны, то есть мультистабильность сохраняется в некотором конечном интервале по параметру е. На рис. 2 показано, что при одних и тех же значениях параметров о = 0.7, ш = 0.9, R = 1, е = 2.15 реализуются как спиральная (рис. 2, e), так и концентрическая (рис. 2, в) волны. Таким образом, в исследуемом ансамбле (3) в силу мультистабильности могут быть реализованы как спиральные, так и концентрические волновые структуры. Отметим, что приведенные результаты получены для случая локальной связи ( R = 1) и для этого случая были описаны в работе [49].

Появление химерных автоволновых структур при введении нелокальности. Как показали исследования многих авторов, введение нелокальной связи в двумерных ансамблях из различных нелинейных осцилляторов приводит к усложнению структур и появлению спирально-волновых химерных состояний [24,26,29,30,36,42,43,49,51-53]. Рассмотрим усложнение пространственно-временных структур, которое имеет место с введением нелокальной связи и зависит от величины радиуса связи в ансамбле (3). На рис. 3, a приведена бифуркационная диаграмма для этой системы на плоскости параметров нелокальности связи R и параметра возбуждения е. В зависимости от величины параметра е, при увеличении радиуса нелокальной связи R, могут возникать или спирально-волновые химеры, или концентрические химеры.

На рис. 3, d, e представлены результаты расчетов влияния радиуса связи на спирально-волновую структуру в системе (3). С ростом R > 1 вначале увеличивается длина волны, а затем, при R = 16, рождается спирально-волновая химерная структура (появляется ядро некогерентности в окрестности центра возбуждения волны (рис. 3, e)). Качественно сходные с рис. 3, d, e

■ ■■■

. Д A A .

А Д A A A A AA A A A AA A A AAAAAAAAAAAAAA ДДДЛДДййЛДДДДДДйДйй4ДДйДДДйДЛД

АДДДАДДДДАДДДДАДДДДАДДДДАДДДДА

Л А Л

А д А Д

А А А Д А А А А А А А А АД АД А А А А А ЖЖ*ЖЖЖЖЖ

.AAA .AAA

AAAAAAAAAA^i

АддддАджжжзйо

;жжжжжжжжжж;

i (em swc

Рис. 3. Бифуркационная диаграмма режимов поведения системы (3) при введении нелокальной (R ^ 1) связи при ю = 0.9, а = 0.7 (a). На карте отмечены области существования: синхронных колебаний (Sync); концентрических и спиральных волн (TW, SW, мгновенные профили приведены на фрагментах (b) при е = 2, R = 1 и (d) при е = 1, R = 1, соответственно); концентрических и спирально-волновых химерных состояний (TWC, SWC, мгновенные профили приведены на фрагментах (с) при е = 2, R = 8 и (e) при е = 1, R = 16, соответственно)

Fig. 3. Bifurcations diagram of behaviour regimes of system (3) when the nonlocal coupling is introduced (R ^ 1) at ю = 0.9, а = 0.7 (a). The areas depicted on the diagram: synchronous oscillations (Sync); target and spiral waves (TW, SW snapshots are shown in (b) at e = 2, R = 1 and in (d) at e = 1, R = 1 correspondingly); target and spiral wave chimera states (TWC, SWC, snapshots are shown in (c) at e = 2, R = 8 and in (e) at e = 1, R = 16 correspondingly)

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

Рассмотрим влияние радиуса нелокальной связи на режим концентрических волн. Можно ожидать, что с введением нелокальной связи в решетке в режиме концентрических волн структуры будут усложняться и, возможно, возникнет новая химерная структура на основе концентрической волны. Эксперименты это подтвердили, о чем свидетельствуют результаты, представленные на рис. 3, Ь, с. Эти результаты иллюстрируют эволюцию концентрической волны при увеличении радиуса связи К. Уже при К = 4 в центре концентрической волны возникает ядро некогерентРис. 4. Переход в режим спирально-волновой химеры для (3) с увеличением радиуса связи R. Мгновенные снимки амплитуд решетки Xij при R = 1 (a); R = 5 (b) и R = 10 (с). Остальные параметры: а = 0.7, ю = 0.9, е = 0.8

Fig. 4. Transition to the spiral wave chimera for (3) when nonlocal coupling range R increases. Snapshots Xij at R =1 (a); R = 5 (b); and R =10 (c). Other parameters: а = 0.7, ю = 0.9, e = 0.8

Рис. 5. Переход от спиральной волны к химере на основе концентрической волны. Мгновенные снимки амплитуд решетки Xij при R = 1 (a); R = 2 (b); R = 8 (с) и R =16 (d). Остальные параметры: а = 0.7, ю = 0.9, е = 2

Fig. 5. Transition from the spiral wave to the target wave chimera. Snapshots Xij at R = 1 (a); R = 2 (b); R = 8 (с); and Д =16 (d). Other parameters: а = 0.7, ю = 0.9, e = 2

ности, которое увеличивается вместе со значением параметра R. На рис. 3, с приведена новая химерная структура при значении параметра R = 8, которую мы назвали концентрической химерой [49]. Как и в случае спирально-волновой химеры, в окрестности источника концентрической волны (рис. 3, b-е) формируется кластер осцилляторов (ядро) с некогерентной динамикой. Концентрические химерные структуры также могут быть получены на основе спиральной волновой структуры (рис. 5). При этом область некогерентности топологически отличается от области некогерентности TWC, полученной из концентрической волны.

Влияние параметра силы связи на пространственно-временные структуры. Исследовалось также влияние параметра силы связи на динамическое поведение решетки (3). С этой целью выбиралась волновая структура при а = 0.7. Далее плавно изменялось значение параметра а в сторону увеличения или уменьшения. Расчеты показали, что, в основном, влияние силы связи ведет к изменениям длин волн как спиральных, так и концентрических. При этом для тех и других волновых структур длина волны возрастала с увеличением параметра связи. Расчеты показали, что области существования этих режимов по параметру силы связи а различаются. Спиральная волна существует при а е [0.65; 1], а концентрическая волна - при а е [0.59; 0.83]. В расчетах значения величины параметра связи а изменялись в пределах от 0 до 1.

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

•<s>

a 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1

...........}..................[............j.....EHEt-...........Ф.......■■■-Et-BSE

Ьооо ■■■■фввс -••--ф>вф--тиЦзввЕ •Ьоооф ■■Мввс

........«¡ЭОев^Э-.....И*ВВЕ

•шооофо на

........•фгеефео.....МЕЕ

• оосроо i ■■

........•»»эефеое-КииЕ

»•с&офосоа ■[ ...........[-••-феофэ........I

••фОООфО

...........^---мвееоаэээээ.....•

OXOXQ -чЭ.....1-2©О \\J=

С=1 С=4 С=9 ring

SBEBEB EBBBEDI

пп -ёвв

Рис. 6. Спирально-волновые химеры в сети нелокально связанных отображений Рулькова (а). С=9 - девятиядерная спирально-волновая химера при а = 0.65, R = 19 (b); С=4 - четырехядерная спирально-волновая химера при а = 0.70, R = 23 (с); С=1 - одноядерная спирально-волновая химера при а = 0.65, R = 53 (d); ring - спирально-волновая химера с ядром в виде кольца при а = 0.85, R = 33 (е). Собственные параметры отображений: а = 4.0, = 0.001, у = 0.2. Время установления режимов t = 20000 + 1000

Fig. 6. Spiral wave chimera states in the network of nonlocally coupled Rulkov maps (a). C=9 is a 9-core spiral wave chimera at а = 0.65, R = 19 (b); C=4 is a 4-core spiral wave chimera at а = 0.7, R = 23 (с); C=1 is a single-core spiral wave chimera at а = 0.65, R = 53 (d); ring is a spiral wave chimera with core in a ring shape at а = 0.85, R = 33 (е). Local map parameters are а = 4, |j, = 0.001, у = 0.2. Transition time is t = 20000 + 1000

рис. 6, а, спирально-волновые химеры появляются только в случаях нелокального и глобального характера связи.

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

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

Переход от спиральных химер к концентрическим волнам. При вариации параметров и начальных условий в области, находящейся слева от границы существования режимов 9-ядерной спирально-волновой химеры (см. рис. 6, а), были обнаружены режимы концентрических волн. Наблюдаются как классические (регулярные) концентрические волны, так и концентрические химеры, содержащие одно некогерентное ядро. Выделить область параметров существования концентрических волн и TWC затруднительно, так как их установление сильно зависит от начальных условий, и они наблюдаются в узкой области параметров. Пример реализации классического режима концентрической волны приведен на рис. 7, который иллюстрирует переключение системы с режима спирально-волновой химеры на режим концентрических волн при протягивании радиуса нелокальной связи от К = 19 до К = 11 с шагом АК = -1, первоначальная

Рис. 7. Переход от спирально-волновых химер к концентрическим волнам при уменьшении радиуса связи в сети нелокально связанных отображений Рулькова (4). Первоначальная структура в слое при R = 19 (a). Динамика новой установившейся структуры при R =11 и t = 7 (b), t = 19 (c), t = 28 (d). Остальные параметры: a = 4, у = 0.2, И = 0.001, a = 0.65

Fig. 7. Transition from spiral wave chimeras to target waves when nonlocal coupling range decreases in the network of coupled Rulkov maps (4). Initial structure in the layer taken with R =19 (a). The dynamics of newly given structure at R =11 and t = 7 (b), t = 19 (c), t = 28 (d). Other parameters: a = 4, у = 0.2, и = 0.001, a = 0.65

Рис. 8. Концентрическая волновая химера (а) и увеличенный фрагмент некогерентного кластера (b) в решетке нелокально связанных отображений Рулькова. Мгновенный пространственный профиль амплитуд решетки. Параметры: a = 4, и = 0.001, у = 0.2, a = 0.45, R = 15

Fig. 8. Target wave chimera (а) and the enlarged fragment of incoherent cluster (b) in the lattice of nonlocally coupled Rulkov maps. Snapshot of amplitudes in the lattice. Parameters: a = 4, и = 0.001, у = 0.2, a = 0.45, R = 15

структура в сети представлена на рис. 7, а. На рис. 8 представлен пример химерной структуры на основе концентрической волны. Кластер осцилляторов с некогерентной динамикой находится в области 155 < i < 185, 145 < j < 180.

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

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

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

Рассмотрим взаимное поведение отдельно взятых соседних осцилляторов из когерентного кластера и из некогерентного кластера. На рис. 9, 10 для каждой выбранной пары осцилляторов были рассчитаны фазовые портреты (фрагменты е, /), мгновенные значения частот в зависимости от времени (фрагменты g, И) и временные реализации (фрагменты 1,]). Мгновенные значения частот вычислялись как величина, обратная длительности одного прошедшего колебания. Соседние осцилляторы, взятые из когерентного кластера, ведут себя одинаково как в случае решетки осцилляторов ван дер Поля (см. рис. 9, е, g, г), так и в случае решетки отображений Рулькова (см. рис. 10, е, g, г). При этом их поведение почти эквивалентно поведению соответствующего осциллятора, взятого отдельно - без ансамбля. В поведении соседних осцилляторов из ядра химерных

Рис. 9. Описание динамики SWC в случае R =16 в решетке ( ). Мгновенный снимок амплитуд решетки xij (a); средние значения частот осцилляторов (b); сечение мгновенного снимка амплитуд j = 49 (c); пространственно-временная диаграмма этого сечения (d); фазовые портреты (e, f), зависимости мгновенных частот от времени (g, h) и временные реализации (i, j) соседних осцилляторов из когерентной области (e, g, i) и некогерентной области f h, j). Остальные параметры: a = 0.7, ш = 0.9, е = 1

Fig. 9. Description of the SWC dynamics for R = 16 in lattice ( ). Snapshot xij of the lattice (a); mean values of oscillators frequencies (b); cut of the snapshot at j = 49 (c); spatio-temporal diagram of the cut (d); phase portraits (e,f), dependencies of instantaneous frequencies on time (g, h), and time series (i, j) of neighbouring oscillators from the coherent region (e, g, i) and from the incoherent region f, h, j). Other parameters: a = 0.7, ш = 0.9, e = 1

состояний наблюдаются отличия: в обоих случаях изменяется фазовый портрет (см. рис. 9, f 10,/); временные реализации осцилляторов ядра химеры в решетке (3) представляют собой смещенные друг относительно друга амплитудно-модулированные колебания (см. рис. 9,у), временные зависимости отображений ядра некогерентности химеры в решетке (4) представляют собой не отличающиеся по форме от колебаний отдельного осциллятора некогерентные друг с другом спайковые колебания (см. рис. 10, у); зависимость мгновенных значений частот от времени в решетке (3) свидетельствует также о наличии частотной модуляции колебаний внутри кластера некогерентности (см. рис. 9, к), мгновенные частоты соседних осцилляторов в ядре некогерентности решетки (4) меняются для каждого отображения нерегулярно во времени и несвязанно между собой (см. рис. 10, к).

Рис. 10. Описание динамики SWC в случае R =53 в решетке (4). Мгновенный снимок амплитуд решетки Xij (a); средние значения частот осцилляторов (b); сечение мгновенного снимка амплитуд j = 100 (c); пространственно-временная диаграмма этого сечения (d); фазовые портреты (e,f), зависимости мгновенных частот от времени (g, h) и временные реализации (i, j) соседних осцилляторов из когерентной области (e, g, i) и некогерентной области f, h, j). Остальные параметры: а = 0.65, a = 4, ц = 0.001, у = 0.2

Fig. 10. Description of the SWC dynamics for R =53 in lattice (4). Snapshot Xij of the lattice (a); mean values of oscillators frequencies (b); cut of the snapshot at j = 100 (c); spatio-temporal diagram of the cut (d); phase portraits (e,f), dependencies of instantaneous frequencies on time (g, h), and time series (i, j) of neighbouring oscillators from the coherent region (e, g, i) and from the incoherent region f, h, j). Other parameters: а = 0.65, a = 4, ц = 0.001, у = 0.2

стрирует рис. 11, а в решетке связанных отображений Рулькова - рис. 12. На фрагментах а, с, d рис. 11 и 12 приведены результаты расчета мгновенных амплитуд, сечение мгновенных амплитуд по ядру некогерентности и пространственно-временная диаграмма этого сечения, соответственно. Эти графики позволяют обнаружить кластер некогерентности в середине концентрической волны, характеризующийся некоторым негладким распределением амплитуд соседних генераторов как в случае связанных осцилляторов ван дер Поля, так и в случае отображений Рулькова. Расчет средних значений частот генераторов решеток в обоих случаях говорит об отсутствии заметных различий с точки зрения распределения величин средней частоты (см. рис. 9, Ь, 10, Ь).

Различия в поведении в выбранном режиме становятся заметными при рассмотрении взаимного поведения отдельно взятых соседних осцилляторов из когерентного кластера и из некогерентного кластера. На рис. 11, 12 для каждой выбранной пары осцилляторов были рассчитаны фазовые портреты (фрагменты е, /), мгновенные значения частот в зависимости от времени (фрагменты g, Н) и временные реализации (фрагменты ¡, /). Мгновенные значения частот вычислялись как величина, обратная длительности одного прошедшего колебания. Соседние осцилляторы, взятые из когерентного кластера, ведут себя одинаково как в случае решетки осцилляторов ван дер Поля (рис. 11, е, g, г), так и в случае решетки отображений Рулькова (рис. 12, е, g, г). Стоит отметить, что поведение отборажений Рулькова в когерентном кластере почти идентично поведению отдельного автономного отображения, но осциллятор ван дер Поля демонстрирует несвойственное ему поведение даже в когерентной области: его колебания становятся модулированными как по амплитуде (рис. 11, /), так и по частоте (рис. 11, g). В поведении соседних осцилляторов из ядра химерных состояний также наблюдаются различия: в обоих случаях изменяется фазовый портрет (рис. 11, /, 12, Временные реализации элементов из ядра химеры в решетке (3) представляют собой смещенные друг относительно друга амплитудно-модулированные, но почти гармонические и противофазные колебания (рис. 11,]), временные зависимости элементов из ядра некогерентности химеры в решетке (4) представляют собой колебания, не отличающиеся по форме от колебаний отдельного осциллятора. При этом колебания соседних элементов некогерентны (рис. 12, ]). Зависимость мгновенных значений частот от времени в решетке (3) свидетельствует также об уменьшении коэффициента частотной модуляции колебаний внутри

Рис. 11. Описание динамики TWC в случае R = 8 в решетке ( ). Мгновенный снимок амплитуд решетки Xij (a); средние значения частот осцилляторов (b); сечение мгновенного снимка амплитуд j = 49 (с); пространственно-временная диаграмма этого сечения (d); фазовые портреты (e, f), зависимости мгновенных частот от времени (g, h) и временные реализации (i, j) соседних осцилляторов из когерентной области (e, g, i) и некогерентной области f, h, j). Остальные параметры: а = 0.7, ш = 0.9, е = 2

Fig. 11. Description of the TWC dynamics for R = 8 in lattice ( ). Snapshot Xij of the lattice (a); mean values of oscillators frequencies (b); cut of the snapshot at j = 49 (с); spatio-temporal diagram of the cut (d); phase portraits (e,f), dependencies of instantaneous frequencies on time (g, h), and time series (i, j) of neighbouring oscillators from the coherent region (e, g, i) and from the incoherent region f h, j). Other parameters: а = 0.7, ш = 0.9, e = 2

Рис. 12. Описание динамики TWC в случае R =15 в решетке (4). Мгновенный снимок амплитуд решетки Xij (а); средние значения частот осцилляторов (b); сечение мгновенного снимка амплитуд j = 160 (c); пространственно-временная диаграмма этого сечения (d); фазовые портреты (e,f), зависимости мгновенных частот от времени (g, h) и временные реализации (i, j) соседних осцилляторов из когерентной области (e, g, i) и некогерентной области f, h, j). Остальные параметры: a = 0.45, a = 4, и = 0.001, у = 0.2

Fig. 12. Description of the TWC dynamics for R =15 in lattice (4). Snapshot Xij of the lattice (a); mean values of oscillators frequencies (b); cut of the snapshot at j = 160 (c); spatio-temporal diagram of the cut (d); phase portraits (e,f), dependencies of instantaneous frequencies on time (g, h), and time series (i, j) of neighbouring oscillators from the coherent region (e, g, i) and from the incoherent region f h,j). Other parameters: a = 0.45, a = 4, и = 0.001, у = 0.2

кластера некогерентности (рис. 11, к). Мгновенные частоты соседних отображений в ядре некогерентности решетки (4) меняются для каждого отображения нерегулярно во времени и несвязанно между собой, довольно значительно отклоняясь от среднего значения частоты (рис. 12, к).

некогерентности как для SWC, так и для TWC (сравните рис. 9, h и 11, h). Есть особенность, свойственная только для TWC в (3), заключающаяся в том, что колебания даже в когерентном кластере становятся модулированными по частоте (рис. 11, g). Распределение средних частот для SWC в решетке (3) отличается от случая SWC в решетке (4) тем, что её значения к центру некогерентного кластера уменьшаются (рис. 9, b), а не увеличиваются (рис. 10, b). Распределение средних частот для TWC в решетке (3) по характеру идентично распределению частот для TWC в решетке (4) и заметно не изменяется в некогерентном кластере (сравните рис. 11, b и 12, b). Из вышесказанного можно заключить, что колебания в некогерентном кластере TWC характеризуются большей регулярностью, чем колебания в некогерентном кластере SWC вне зависимости от выбора парциального элемента.

Рассматривается влияние внешнего шумового воздействия на различные автоволновые структуры, наблюдаемые в исследуемых ансамблях, и проводится сравнение результатов этого воздействия. Как показали проведенные исследования, в обеих решетках возможна реализация двух типов волновых структур: спиральных и концентрических волн. Эти структуры важны не только с чисто научной точки зрения, но и в приложениях. Известно, например, что режим концентрических волн в сердечной мышце отвечает нормальному функционированию сердца, а спиральные волны отвечают возникновению патологии, связанной с фибрилляцией, аритмией и тахикардией [56,57]. В связи с этим возникает вопрос о возможности подавления режима спиральных волн с переходом к концентрическим волнам с целью избежать патологий в работе сердца. Подобный вопрос исследовался ранее в работах [58-60], где авторы реализовали указанный переход путем воздействия периодического сигнала на одиночный элемент ансамбля. В нашей работе внешнее воздействие вводится в виде аддитивного белого Гауссова шума с интенсивностью D в уравнение fx в ( ) и ( ), то есть добавляется слагаемое \\/2D&%tij, где ^ - источник шума. Наши исследования показали, что переход от спиральной волны к концентрической в решетке локально связанных отображений Рулькова может быть реализован при воздействии шумового возмущения на конечный кластер осцилляторов ансамбля. В решетке связанных осцилляторов ван дер Поля как спиральные, так и концентрические волны оказываются устойчивыми к небольшому шумовому воздействию, но при увеличении шумового воздействия спиральная волна с одним источником перестраивается в спиральную волну с несколькими источниками, а концентрическая переходит в спиральную волну.

ван дер Поля. Рассмотрим влияние внешнего шумового воздействия одинаковой интенсивности D на все элементы решетки одновременно. С этой целью мы устанавливаем интенсивность шума на некотором значении D = Do и интегрируем уравнения ( ) 10000 единиц дискретного времени. Затем при нулевом значении D = 0 проводим вычисления в течение такого же времени. Мгновенный снимок амплитуд Xij осцилляторов после 20000 итераций рассматриваем как результат. На рис. 13 приведены данные расчетов влияния шумового воздействия для спиральной (рис. 13, a) и концентрической (рис. 13, d) волн, соответственно. При Dq < 0.01 обе структуры устойчивы к шуму (рис. 13, a, d). При Dq ^ 0.01 спиральная волна разрушается и осуществляется переход к спиральным волнам с множеством ведущих центров (рис. 13, b, c). Иначе ведет себя решетка в режиме концентрической волны: под действием шума интенсивности Dq ^ 0.01 реализуется индуцированный шумом переход к режиму спиральной волны с множеством ведущих центров (рис. 13, e,f).

J 1 25 50 75 J 1 25 50 75 J

Рис. 13. Сравнение результатов воздействия шума на ансамбль в режимах спиральных (a-c) и концентрических (d-f) волн. Мгновенные снимки амплитуд Xij решетки (3) при Do = 0.001 (a, d); Do = 0.01 (b, e); Do = 0.1 (c, f). Остальные параметры: R =1, ю = 0.9, e = 2, a = 0.7

Fig. 13. Comparison of noise influence results on the ensemble in the spiral wave regime (a-c) and in the target wave regime (d-f). Snapshots Xij of lattice ( ) at Do = 0.001 (a, d); Do = 0.01 (b, e); Do = 0.1 (c, f). Other parameters: R = 1, ю = 0.9, e = 2, a = 0.7

J 1 50 100 150 J 1 50 100 150 J

Рис. 14. Воздействие аддитивного белого Гауссова шума на спиральные волны. а - Спиральная волна, устанавливающаяся в системе ( ). b - Установившаяся структура t = 199970 после воздействие шума на квадратную область 2х2 (69 < г < 70, 120 < j < 121). c - Индуцированная концентрическая волна при воздействие шума на квадратную область 11х11 (65 < г < 75, 115 < j < 125), t = 199950. Параметры: а = 4, ц = 0.001, у = 0.2, а = 0.8, R = 1, D = 0.5

Fig. 14. Effect of additive white Gaussian noise on spiral waves. a - Initial spiral wave regime in the lattice ( ). b - Steady-state structure at t = 199970 after the noise effect on a square of 2 x 2 elements (69 < г < 70, 120 < j < 121), c - induced target wave under the influence of noise on a square region of 11 x 11 elements (65 < г < 75, 115 < j < 125) at time t = 199950. Parameters: а = 4, ц = 0.001, у = 0.2, а = 0.8, R = 1, D = 0.5

Рис. 15. Установление концентрических волн во времени при воздействие шумом на спиральную волну в решетке Рулькова (4) (см. рис. 14, a) для различных моментов времени: t = 9940 (a), t = 99960 (b), t = 199920 (c). d - Пространственно временная диаграмма сечения j = 120 после полного установления концентрических волн в системе. Параметры, как на рис. 14

Fig. 15. Target wave formation with the noise presence in the Rulkov lattice (4) (fig. 14, a) for différent time: t = 9940 (a), t = 99960 (b), t = 199920 (c). d - Spatio-temporal diagram for the cross section j = 120 for the steady-state regime of target wave. Parameters are as in Fig. 14

воздействии шума на малый квадрат размером 2х2 элемента (69 ^ г ^ 70, 120 ^ j ^ 121). Малой области воздействия шума не достаточно, чтобы разрушить спиральную волну, но это воздействие перемещает ядро спиральной волны (сравн. рис. 14, a и 14, b). Увеличение области воздействия шума приводит к разрушению спиральной волны и установлению концентрической волны (см. рис. 14, b). В данном случае область воздействия шума представляет собой квадрат размером 11х11 элементов (65 ^ i ^ 75, 115 ^ j ^ 125).

Отметим, что концентрические волны наблюдаются в системе только при воздействии шума. Если убрать шумовое воздействие (Dx,y = 0), то система возвращается к спиральным волнам или переходит в когерентный режим. Рис. 15 демонстрирует становление в системе концентрических волн при воздействии шума на спиральную волну. Как можно видеть, спиральная волна разрушается постепенно, и со временем вокруг области воздействия устанавливаются новые кольца концентрической волны, пока они полностью не вытеснят спирально-волновую структуру из среды. На рис. 1 , d показана динамика во времени сечения j = 120 после полного установления концентрических волн в системе. Данная концентрическая волна распространяется в направлении к ядру (области, на которую мы воздействуем шумом).

Заключение

Обнаружение нового типа автоволновой структуры - химерного состояния на основе концентрической волны - породило интерес к более детальному исследованию этого явления и сравнению его с широко известной спирально-волновой химерной структурой. В работе представлены результаты моделирования динамики двумерного ансамбля в случаях нелокально связанных осцилляторов ван дер Поля и отображений Рулькова. Индивидуальные осцилляторы ван дер Поля рассматривались в режиме релаксационных колебаний, а индивидуальные отображения Рулькова - в режиме генерации спайков. Проведено исследование зависимости разнообразия наблюдаемых структур от выбора парциального элемента. Показано, что в обоих ансамблях возможно наблюдение как классических спиральных и концентрических волн, так и химерных структур на их основе. Последние обнаруживаются при введении нелокальной связи при достаточно больших значениях параметра нелокальности. Для обоих ансамблей представлены соответствующие бифуркационные диаграммы на плоскости управляющих параметров. Спиральные и концентрические химеры были детально исследованы с точки зрения сравнения усредненных и мгновенных частотных характеристик отдельных элементов ансамблей, реализующих эти режимы. Результаты численных экспериментов показали, что для обоих ансамблей в случае TWC средние значения частот в ансамбле распределяются равномерно, тогда как в случае SWC средние значения частот в некогерентном кластере отличаются от значений в когерентных областях. Отсюда следует, что для двух ансамблей с различными динамическими элементами колебания в некогерентном кластере TWC характеризуются большей регулярностью, чем колебания в некогерентном кластере SWC. Результаты по влиянию внешнего шумового воздействия на классические спиральные и концентрические структуры показали, что в ансамбле осцилляторов ван дер Поля становится возможным индуцированный шумом переход от концентрических волн к спиральным. В ансамбле связанных отображений Рулькова была получена концентрическая структура, которая существует только при наличии постоянного шумового воздействия. Учитывая, что спиральные и концентрические волны широко используются при анализе динамики сердечной мышцы [7, 18,20,22,25], результаты настоящей работы могут быть полезны в приложениях к задачам кардиологии.

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

References

Бух Андрей Владимирович - родился в Палласовке, Волгоградской области (1993). Окончил с отличием бакалавриат на физическом факультете Саратовского государственного университета по направлению «Радиофизика» (2015). Затем окончил с отличием магистратуру на физическом факультете Саратовского государственного университета по направлению «Радиофизика» (2017). С 2015 года работает на кафедре радиофизики и нелинейной динамики Саратовского государственного технического университета в должности инженера. С 2017 года обучается в аспирантуре Саратовского государственного университета по направлению «Физика и астрономия», направленности «Радиофизика». Научные интересы -нелинейная динамика и теория колебаний, бифуркационный анализ, синхронизация, теория устойчивости, ансамбли связанных осцилляторов, химерные состояния. Опубликовал свыше 10 научных статей по указанным направлениям.

Россия, 410012 Саратов, Астраханская, 83

Саратовский государственный университет имени Чернышевского Н.Г. E-mail: buh.andrey@yandex.ru

Рыбалова Елена Владиславовна - родилась в Энгельсе, Саратовской области (1996). Окончила с отличием бакалавриат на физическом факультете Саратовского государственного университета по направлению «Радиофизика» (2018). Затем поступила в магистратуру на физическом факультете Саратовского государственного университета по направлению «Радиофизика» (2018). С 2018 года работает на кафедре радиофизики и нелинейной динамики Саратовского государственного технического университета в должности инженера. Научные интересы - нелинейная динамика и теория колебаний, синхронизация, влияние флуктуаций, ансамбли связанных осцилляторов, химерные состояния. Опубликовала свыше 10 научных статей по указанным направлениям.

Россия, 410012 Саратов, Астраханская, 83

Саратовский государственный университет имени Чернышевского Н.Г. E-mail: rybalovaev@gmail.com

Анищенко Вадим Семенович - родился в Саратове (1943). Окончил с отличием физический факультет Саратовского государственного университета по специальности «Радиофизика и электроника» (1966). Защитил кандидатскую диссертацию на тему «Шумовые свойства электронных потоков СВЧ-усилителей О-типа» под руководством Штырова А.И. (1970, СГУ). Защитил докторскую диссертацию на тему «Механизмы развития и свойства хаотических колебаний в радиофизических системах с конечным числом степеней свободы» (1986, СГУ). С 1966 года работает инженером в НИИМФе СГУ. В 1970 году стал ассистентом на кафедре радиофизики СГУ, в 1972 году занял должность старшего преподавателя, в 1976 стал доцентом, в 1987 - профессором, а с 1998 года является заведующим кафедрой радиофизики СГУ. Научные интересы - теория нелинейных колебаний, статистическая радиофизика, нелинейная динамика сложных систем, динамический хаос и влияние флук-туаций, синхронизация и управление. Опубликовал свыше 430 научных статей и более 20 книг по указанным направлениям.

Россия, 410012 Саратов, Астраханская, 83

Саратовский государственный университет имени Чернышевского Н.Г. E-mail: wadim@info.sgu.ru