ОХОТА НА ХИМЕР В ПОЛНОСВЯЗНЫХ СЕТЯХ НЕЛИНЕЙНЫХ ОСЦИЛЛЯТОРОВ Текст научной статьи по специальности «Математика»

Известия высших учебных заведений. Прикладная нелинейная динамика. 2022. Т. 30, № 2 Izvestiya Vysshikh Uchebnykh Zavedeniy. Applied Nonlinear Dynamics. 2022;30(2)

Научная статья УДК 517.926

DOI: 10.18500/0869-6632-2022-30-2-152-175

Охота на химер в полносвязных сетях нелинейных осцилляторов

Д. С. Глызин, С. Д. Глызин^, А.Ю. Колесов

Ярославский государственный университет имени П. Г. Демидова, Россия E-mail: glyzin@gmail.com, Elglyzm@uniyar.ac.ra, kolesov@uniyar.ac.ru Поступила в редакцию 21.12.2021, принята к публикации 16.02.2022, опубликована 31.03.2022

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

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

Благодарности. Исследования, изложенные в разделе 4, выполнены в рамках реализации программы развития РНОМЦ (ЯрГУ) при поддержке Минобрнауки России (соглашение № 075-02-2021-1397). Исследования, приведенные в разделах 1-3, выполнены при поддержке Российского научного фонда (грант № 21-71-30011).

Для цитирования: Глызин Д. С., Глызин С.Д., Колесов А.Ю. Охота на химер в полносвязных сетях нелинейных осцилляторов//Известия вузов. ПНД. 2022. T. 30, № 2. С. 152-175. DOI: 10.18500/0869-6632-2022-30-2-152-175

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

©Глызин Д. С., Глызин С. Д., Колесов А. Ю., 2022

Article

DOI: 10.18500/0869-6632-2022-30-2-152-175

Hunt for chimeras in fully coupled networks of nonlinear oscillators

D. S. Glyzin, S.D. Glyzin81, A. Yu. Kolesov

P. G. Demidov Yaroslavl State University, Russia E-mail: glyzin@gmail.com, Elglyzin@uniyar.ac.ru, kolesov@uniyar.ac.ru Received 21.12.2021, accepted 16.02.2022, published 31.03.2022

Abstract. The purpose of this work is to study the dynamic properties of solutions to special systems of ordinary differential equations, called fully connected networks of nonlinear oscillators. Methods. A new approach to obtain periodic regimes of the chimeric type in these systems is proposed, the essence of which is as follows. First, in the case of a symmetric network, a simpler problem is solved of the existence and stability of quasi-chimeric solutions — periodic regimes of two-cluster synchronization. For each of these modes, the set of oscillators falls into two disjoint classes. Within these classes, full synchronization of oscillations is observed, and every two oscillators from different classes oscillate asynchronously. Results. On the basis of the proposed methods, it is separately established that in the transition from a symmetric system to a general network, the periodic regimes of two-cluster synchronization can be transformed into chimeras. Conclusion. The main statements of the work concerning the emergence of chimeras were obtained analytically on the basis of an asymptotic study of a model example. For this example, the notion of a canonical chimera is introduced and the statement about the existence and stability of solutions of chimeric type in the case of asymmetry of the network is proved. All the results presented are extended to a continuous analogue of the corresponding system. The obtained results are illustrated numerically.

Keywords: fully coupled network of nonlinear oscillators, periodic modes of two-cluster synchronization, hunting for chimeras, stability, buffering.

Acknowledgements. Sections 1-3 of this work were supported by the Russian Science Foundation (project No. 21-71-30011). Section 4 was carried out within the framework of a development programme for the RSEMC of the P. G. Demidov Yaroslavl State University with financial support from the Ministry of Science and Higher Education of the Russian Federation (Agreement No. 075-02-2021-1397).

For citation: Glyzin DS, Glyzin SD, Kolesov AYu. Hunt for chimeras in fully coupled networks of nonlinear oscillators. Izvestiya VUZ. Applied Nonlinear Dynamics. 2022;30(2):152-175. DOI: 10.18500/0869-6632-2022-30-2-152-175

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

1. Общая стратегия охоты на химер

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

х^ = Р^ (Хз ,щ), ] = 1,2,...,т. (1)

Здесь т ^ 2, х^ = х^ (¿) е Мга, п ^ 2, точка — дифференцирование по I,

т

Щ = ^ GS(XS), (2)

8=1

«=3

а вектор-функции Р^ (х, и), С^ (х), ] = 1, 2,... ,т, со значениями в Мга бесконечно дифференцируемы по своим переменным (х,и) е Мга х Мга и х е Кга. Как правило, каждая из отвечающих сети (1), (2) парциальных систем

х = ^ (х, 0), ] = 1,2,...,т (3)

допускает экспоненциально орбитально устойчивый цикл, то есть представляет собой нелинейный осциллятор. Мы же рассматриваем ситуацию, когда т осцилляторов (3) взаимодействуют друг с другом по принципу «каждый со всеми».

В частном случае, когда

Рз (х, и) = Р^ (х) + Бз (х)и, ] = 1,2,...,т, где Dj(х) — квадратные матрицы размера пхп, система (1), (2) приобретает вид

т

%з = рз (X])+ Оз (Хз) ^ С3(х3), э = 1,2,...,т. (4)

8=1

2=3

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

Следует отметить, что цепочки и решетки связанных идентичных нелинейных осцилляторов используются в качестве математических моделей в различных областях естествознания: биофизике, экологии, оптике, химической кинетике, нейродинамике, генной инженерии и т. д. При этом в подавляющем большинстве работ, посвященных системам связанных нелинейных осцилляторов, изучаются специальные стационарные режимы упомянутых систем — так называемые химеры. Универсального математического определения химеры не существует. На эвристическом же уровне строгости одна из возможных ситуаций, когда мы имеем дело с химерой, состоит в следующем. Пусть т — количество осцилляторов в системе, а т1, т2 : т1 < т2 < т — некоторые фиксированные натуральные числа. Тогда в случае химерного стационарного режима имеются две зоны когерентности: для любых двух компонент Х3 с номерами 1 ^ ] ^ Ш1 или < ] ^ т наблюдается синхронизация (возможно, приближенная). Под термином «синхронизация» во избежание недоразумений мы понимаем выполнение равенств вида Ж1 = х2 = ... = хт1, хт2+1 = хт2+2 = ... = хт. Любые же два осциллятора с номерами .]1,.]2 : Ш1 < ,]1 < ,]2 ^ т2 колеблются асинхронно, то есть ф Хз2 (это так называемая зона некогерентности). Добавим еще, что, как правило, химеры характерны для случая т, т1, т2 ^ 1 и выявляются посредством численного анализа соответствующих математических моделей.

С момента своего открытия в 2002 г. (см. [1]) химеры стали активно изучаемым объектом исследования [2,3]. Было показано, что они могут возникать как в глобально связанных [4,5], так и в локально связанных [6] сетях нелинейных осцилляторов. Например, химеры удалось обнаружить в модели Стюарта-Ландау [7,8], в связанных системах ван дер Поля [9], ФитцХью-Нагумо [10], в моделях Ходжкина-Хаксли [11], Хиндмарша-Роуза [12] и т. д. Это, разумеется, далеко не полный библиографический список. Более подробную библиографию по химерам можно найти в книге [13].

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

На первом этапе мы рассматриваем аналогичную (1), (2) симметричную полносвязную сеть

т

±3 = Р (х^ ), Uj = ^ С(х3), ] = 1,2,...,т, (5)

8=1

2=3

где, как и выше, Р(х, и) е х Мга; Мга), С(х) е СМга). Далее, нас будут интересовать

периодические режимы системы (5), которые можно квалифицировать как квазихимеры. В описанной выше ситуации им соответствует случай Ш1 = Ш2, когда отсутствует зона некогерентного

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

С : Х1 = Х2 = ... = хк = Ь(к)(г), Хк+1 = Хк+2 = ... = хт = 'Ш(к)(г), (6)

где к : 1 ^ к ^ т — некоторое фиксированное натуральное число, а У(к)(Ь), -ы^^) — периодические вектор-функции. Предположим, что данный цикл в системе (5) существует и экспоненциально орбитально устойчив.

На втором этапе перейдем от (5) к несимметричной сети

т

= ^з (хз ,из ,£), из = С3(х3,е), ] = 1,2,...,т, (7)

8=1 э=3

(8)

где

^ (х,и,е) = Р (х,и) + (х,и), Су (х,е) = С(х) + еА2^ (х),

Аи (х,и) е х Мга; Мга), А2^ (х) е С ~(Мга; Мга), ] = 1, 2,...,т,

а е > 0 — малый параметр. Ясно, что система (7) допускает асимптотически близкий к (6) устойчивый цикл С (б), который при подходящих возмущениях А^-, из (8) может стать химерой.

Действительно, фиксируем натуральные т1, т2 : т1 < к < т2 < т и предположим, что

А1,1 = А1,2 = ... = А1,тх, А1,т2+1 = А1,т2+2 = ... = ^^ (9)

А2,1 = А2,2 = ... = А2,тг , А2,т2+1 = А2,т2+2 = ... = А2,m,

а при Ш1 + 1 ^ ] ^ Ш2 добавки А^-, А2^ нетривиально зависят от индекса ] (строгий смысл этому понятию будет придан ниже при рассмотрении конкретного примера сети (7)). Тогда, очевидно, система (7) допускает так называемые частично однородные решения (не обязательно периодические), для которых выполняются аналогичные (9) соотношения

Ж1 — Х2 — ... — , %т,2+1 — ^т.2+2 — ... — Хт. (10)

А отсюда, в свою очередь, следует, что равенства (10) справедливы и для цикла С (е) (поскольку он заведомо обладает частичной однородностью при е = 0). Что же касается компонент х^ этого цикла при Ш1 + 1 ^ ] ^ Ш2, то в общем случае они демонстрируют некогерентное поведение, то есть каждая пара из них колеблется асинхронно.

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

Фиксируем произвольно натуральное к : 1 ^ к ^ т — 1 и предположим, что совокупность индексов 1 ^ ] ^ т разбита на два непересекающихся множества А и В, состоящих из к и т — к элементов, соответственно, то есть

{1, 2,..., т} = А и В. (11)

Тогда, как нетрудно увидеть, система (5) допускает решения с компонентами

Xj = у^) при ] е А, Ху = ■ш^) при ] е В, (12)

где переменные V, ад удовлетворяют вспомогательной системе

V = Р(V, ик,1), й) = Р(ад, ик,2),

в которой

Uk,i = (к - l)G(v) + (т - k)G(w), ик>2 = к G(v) + (m - к - l)G(w).

(14)

А поскольку нас интересуют периодические решения системы (5) вида (12), то предполагаем выполненным следующее

Условие 1. Система (13), (14) имеет непостоянное периодическое решение

Ск : (у,ад) = (у{к)(1),ад{ф)) (15)

периода Т(к) > 0, удовлетворяющее требованию неоднородности

ад{к)(г). (16)

Сформулированное условие гарантирует существование у исходной системы (5) целого семейства циклов ик. Все циклы из этого семейства задаются равенствами (11), (12) при

V = v(k)(t), w = W(k)(t),

(17)

а их количество равно, очевидно, С^. Добавим еще, что в силу условия неоднородности (16) они представляют собой режимы двухкластерной синхронизации.

Для исследования устойчивости циклов (11), (12), (17) системы (5) нам потребуется серия Т(к)-периодических по £ матриц А3(Ь), В3(Ь), в = 1, 2, 3 размера пхп. Упомянутые матрицы зададим равенствами

dF

Ai(t) = Q.-(х,и)

dF

Л-2 (t) = — (х,и)

dF

A, (t) = — (х,и)

dF

Bl(t) = ~di(x,u)

dF

B2(t) = (x,u)

dF

Bä(t) = (x,u)

X=V(k)(t), U=Uk,l(t) X=V(k)(t), U=Uk,l(t) x=v(k)(t), u=uk)1(t) x=W(k)(t),u=uk2(t) x=W(k)(t),u=uk2(t) x=W(k)(t),u=uk2(t)

G'x(v(k)(t)),

Glx(w{k)(t)),

G'x(v(k)(t)),

Glx(w{k)(t)),

(18)

(19)

(20) (21) (22) (23)

где

uk,i(t) = (k - 1)G(v{k)(t)) + (m - k)G(w{k)(t)), uk,2(t) = kG(v{k)(t)) + (m - к - 1)G(w{k)(t)).

Далее, введем в рассмотрение линейные системы

¿1 = (Ai(t) + (к - l)A2(t))a + (т - k)A3(t)c2, , ч

(24)

¿2 = kB2(t)ci + (Bi(t) + (т - к - l)B3(t))c2, с =(Ai(t) - A2(t))c, c = (Bi(t) - B3(t))c, (25)

где c1; c2, с G Rra.

Из формул (18)—(23) непосредственно следует, что система (24) представляет собой линеаризацию вспомогательной системы (13), (14) на цикле (15). Тем самым, она заведомо имеет единичный мультипликатор. В дальнейшем считаем, что этот мультипликатор является простым. Более того, предположим, что выполнено очередное

Условие 2. Все мультипликаторы систем (24), (25) (за исключением простого единичного в случае системы (24)) по модулю строго меньше единицы.

Сформулированные ограничения позволяют доказать следующий базовый результат.

Теорема 1. При условиях 1, 2 система (5) допускает семейство U периодических режимов двухкластерной синхронизации, количество которых равно С^. Все эти циклы экспоненциально орбитально устойчивы.

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

Во-первых, несложная проверка показывает, что система (5) инвариантна относительно

замен

(Ж1,Ж2, . . .,Хт) ^ (Xj-1 ,Xj2 ,. . . , Xjm ), (26)

где (ji ,j2,..., jm) — произвольная перестановка набора индексов (1, 2,... ,т). Во-вторых, периодические режимы из семейства U допускают кодирование с помощью бинарных векторов

(0!, 02,..., 0m): 0J = 1 или 0, j = 1, 2,...,т. (27)

Точнее говоря, предполагаем, что j-я координата вектора (27) равна 1 или 0 при j G A или j G B, соответственно. В этом случае между векторами (27), содержащими к единиц и т, - к нулей, и циклами семейства U имеет место взаимно однозначное соответствие.

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

( , О,..^ ).

к m—к

Линеаризуя исходную систему (5) на указанном цикле, приходим к линейной системе вида

к m

hj = Ai(t)hj + ^ A2(t)hs + ^ A3(t)hs, 1 < j < к, (28)

s=i s=k+i

s=j

к m

hj = Bi(t)hj + B2(t)hs + Y^ B3(t)hs, к + 1 < j < m, (29)

«=i s=k+i

s=j

где As(t), Bs(t), s = 1,2,3 — матрицы (18)-(23), hj G Rn, j = 1,2,..., m.

Как оказывается, при подходящей замене переменных система (28), (29) принимает блочную структуру: распадается на одну 2п-мерную подсистему и т — 2 подсистем размерности п.

Действительно, введем в рассмотрение серию матриц размера пт х п, задающихся в виде матричных столбцов

elfl = ( I, I, О, О, у., О )*, е2,о = ( р,0,у.,0,,1,1,^..,1 )*, (30)

к т—к к т—к

ем = (I,0,...,0, —I,0,...,0)*, s = 1,...,к — 1, (31)

s + 1

e2,s = (0,...,0, I ,0,...,0, —I ,0,...,0)*, s = 1,...,т — к — 1. (32) к к+1 к+1+s

Здесь I и О — единичная и нулевая матрицы размера пхп, * — операция транспонирования, а нижние подписи обозначают номера позиций, которые занимают эти пхп-блоки в соответствующих матричных столбцах (30)-(32).

Опираясь на приведенные матрицы, выполним в системе (28), (29) замену

к—1 т—к—1

colon (hi,h2,..., hm) = ei,oCi,o + е2,ос2,о + ^ ei,sci,s + ^ e2,sc2,s, (33)

8=1 S=i

где ci}0, С2,о, ciyS, c2,s £ Rn — новые переменные. В результате приходим к выводу, что пара компонент (cit0, С2,о) = (cito(t), c2,0(t)) из (33) удовлетворяет линейной системе (24). Что же касается компонент ci}S = ci>s(t), s = 1,... ,к — 1, то все они являются решениями первой системы из (25). Аналогичным образом, группа координат c2,s = c2,s(t), s = 1,... ,т — к — 1, удовлетворяет второй системе из (25). А отсюда и из Условия 2 вытекает требуемая экспоненциальная орбитальная устойчивость всех циклов семейства U. Теорема 1 доказана. □

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

2. Анализ модельного примера

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

Для построения интересующей нас сети обратимся сначала к двумерной системе, имеющей в комплексной форме записи вид

z = z — d0lzl2z, z = х + iy, x,y £ R, d0 = 1 — ic0, c0 = const > 0. (34)

Несложно проверить, что эта система допускает экспоненциально орбитально устойчивый гармонический цикл Z0(t) = exp(ic0t), то есть является простейшим нелинейным осциллятором. Рассмотрим далее полносвязную сеть

т

2

Zj = zj — d>0\zj|2 Zj — —di z2 j = 1,2,...,fm

m

s=i

S=j

(35)

осцилляторов (34), где Zj = Xj + г yj, Xj, yj £ R, j = 1, 2,...,m, v = const > 0, d1 = 1 + i c1, c\ = const £ R. Как будет показано ниже, сеть (35) обладает интересующим нас свойством: все ее периодические режимы двухкластерной синхронизации, а также условия их устойчивости, выписываются в явном виде.

Анализ системы (35) начнем с получения условий ее диссипативности. Для этого нам потребуется функция

V (Z1, Z2,..., zm, Z1, Z2 • ••, Zm) = Y ' Zj I2.

(36)

3 = 1

Справедливо следующее утверждение. Лемма 1. При выполнении неравенства

v < m

(37)

имеет место оценка

±V ^ V - (1 - v) lv 2 \ m J m

2

(38)

где V — производная функции (36) в силу системы (35).

Доказательство. Для вычисления производной V заметим сначала, что в силу (35) при любом индексе ] : 1 ^ ] ^ т справедливо равенство

1 = - i1 - mb'1 - 2m'

('ill

d!Z2 ^

8=1

^+di ¿¡Y: ^2

8=1

Суммируя затем выражения (39) по j, приходим к выводу, что

(39)

1

(ч т

1 - m е

7 .7=1

o^ = V - 1 -- > | Zj |1 --

m

Е

3 = 1

<V - о-m) en

3 = 1

1

А отсюда с учетом очевидного свойства

(т ч 2 т

£ М2 ^ m Е М1

3=1 J 3=1

получаем требуемую оценку (38). Лемма 1 доказана. □

Установленная лемма позволяет без труда разобраться с интересующим нас вопросом о диссипативности системы (35). Действительно, пусть выполнено условие (37). Тогда в силу (38) при любом фиксированном К > т2/(т — у) все траектории нашей системы при увеличении £ втекают в шар

Z1, Z2 ,•••, Z,

г): Е 'ziI2

3 = 1

то есть имеет место требуемое свойство диссипативности. При выполнении же строго противоположного (37) неравенства диссипативность отсутствует.

Действительно, непосредственная проверка показывает, что система (35) допускает инвариантное многообразие вида

Zj = zехр(iyj), 3 = 1, 2,... ,т,

2

2

з

где вещественные постоянные у^ подчинены условию

^2ехр(2г ) = 0,

3 = 1

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

¿ = г-( (1о - — йЛ Ы2 г. (40)

\ т )

Что же касается уравнения (40), то в случае V > т любое его решение г(Ъ) = 0 определено на конечном полуинтервале вида [0, ¿0) и |г(Ь)1 ^ при £ ^ ¿0.

Перейдем теперь к вопросам о существовании и устойчивости периодических решений системы (35). Нетрудно проверить, что рассматриваемая сеть допускает однородный цикл

Х] = ехр( г ), 3 = 1,2,...,т, (41)

где

So = ^1/(1 + v(т — 1)/т), шо = (cq - S2,

(42)

а также периодические режимы двухкластернои синхронизации

Zj = S0 exp(iш0t), 1 ^ j ^ k, Zj = —S0 exp(гш0t), k + 1 ^ j ^ т (43)

при всех к : 1 ^ к ^ т — 1. Как уже отмечалось в разделе 1, каждый из циклов (43) порождает, в свою очередь, целое семейство U режимов двухкластерной синхронизации. Обозначим через U — объединение циклов из (J^Т^1 U с однородным циклом (41). Как оказывается, все эти периодические режимы устойчивы или нет одновременно.

Лемма 2. Любой цикл из U экспоненциально орбитально устойчив при выполнении неравенств

< 1, * (1 + <Ц) + С0С1 — 1 > 0 (44)

т т 1

и неустойчив при строгом нарушении хотя бы одного из них.

Доказательство. При обосновании леммы обратим внимание на два обстоятельства. Во-первых, как нетрудно увидеть, система (35) инвариантна по отношению к заменам переменных

(6i2i, Q2Z2, ..., Qrn Zm) ^ (Z1, Z2,..., Zm), (45)

где множители 6j, j = 1, 2,... ,т независимо друг от друга принимают значения 1 или —1. Во-вторых, любые два цикла семейства U переходят друг в друга при некоторой замене вида (45). А отсюда автоматически следует, что свойства устойчивости всех циклов из U одинаковы. Таким образом, достаточно разобраться с устойчивостью лишь одного из них — однородного цикла (41), (42).

Положим в (35) Zj = S0(1 + hj) exp(iш01), hj = h1,j + ih2,j, h1,j,h2,j e R и отбросим нелинейные по hj, hj слагаемые. В результате приходим к линейной системе

h3=— dotoh+h)+h3—¿1?о h—

2v 2 ™ (46)

—K j = 1,2,...,т.

S=1

Заметим, далее, что получившаяся система допускает инвариантное подпространство = Н2 = ... = Нт, на котором она записывается в виде

h = -¥0(do + А(h + h), h = -уhj,

т к ш ^—'

7 3=1

а также инвариантное подпространство (к\, Ъ2,..., Ьт) : ^^ = 0, на котором она приобре тает вид

h = -фо - ^^ di) h3 - Фо + h3, , = 1,

Bi = -но

Тем самым спектр устойчивости системы (46) совпадает с собственными значениями матрицы

_ Ут - 1- - Уш - х- (47)

, do +--й\ (1о +--й\ г

\ т т /

и т — 1 одинаковых матриц

(do — Ъ do + 22

л тп <т

(48)

U + d! do + v^-il ^

m m

do + ———— di do + ———— d\ I \ m m /

B2 = -H0

(do - di do + dl\

m m

, do +—(-) di do--(-) di I

\ m m /

Остается заметить, что матрица (47) имеет собственные значения Xi =0 и X2 = —2, а матрица (48) при условиях (44) является гурвицевой. При строгом же нарушении хотя бы одного из этих условий В2 допускает собственное значение в полуплоскости {X £ C : ReX > 0}. Лемма 2 доказана. □

Перейдем теперь от (35) к соответствующей несимметричной сети

т

z3 = (1 + ie^j )zj — d0\zj\2Zj--diZj^^ z2s, j = 1,2,...,m, (49)

s=i

S=j

где e £ R, \e\ ^ 1, И = const £ R, j = 1,2,...,m. Предполагаем существование таких натуральных mi, т2 : mi < т2 < т, что

И — И*2 — ... — Итог, Иго.2 — M-m^+i — ... — Иго. (50)

В случае же mi + 1 ^ j ^ т2 — 1 считаем зависимость величин и от индекса j нетривиальной. Последнее означает, что

Ил = И32 V ji,j2 £ [mi + 1,т2 — ^ h = 32. (51)

Как оказывается, для системы (49) приведенному в разделе 1 эвристическому описанию химеры можно придать строгий смысл. А именно, канонической химерой назовем цикл этой системы, компоненты Zj = Zj (t) которого удовлетворяют требованиям

¿1 (t) = Z2(t) = ... = Zmi (t), Zm2 (t) = Zm2 + i(t) = ... = zm(t), (52)

Z31(t) Ф zh(t) V h,32 £ [mi + 1,m2 — ^ h = ^ (53)

где натуральные т1, т2 заимствованы из (50), (51). Что же касается отрезков 1 ^ j ^ т1, т2 ^ 3 ^ т и т1 + 1 ^ ]^ т2 — 1, то, как обычно, их будем называть промежутками когерентного и некогерентного поведения осцилляторов.

При решении вопроса о существовании канонической химеры предположим, что параметры Со, С1, V системы (49) удовлетворяют требованиям

т

0 < V <-— , С0С1 > 1. (54)

т + 1

В этом случае справедливо следующее утверждение.

Теорема 2 (о существовании канонической химеры). При выполнении сформулированных выше ограничений и при всех е € [—е0, е0], где е0 > 0 достаточно мало, система (49) имеет экспоненциально орбитально устойчивый цикл

т

хз = Рз (е)ехР[ ъ Фз (е) + ъ 3 = 1,2,...,т, Е <фj (е) = 0, (55)

3 = 1

являющийся при е = 0 канонической химерой. Здесь вещественные функции Рj (е), фj (е), ю(е ) аналитически зависят от е и при е ^ 0 допускают асимптотику:

р^(е) = Р0 + £Р1^ + 0(е2), ф^-(е) = еф1^ + 0(е2), т(е) = Ю0 + + 0(е2), (56)

где

И,е х,- 1т х,- 1 ^л

Р0 = ^ р1,з = , Фи = -гГ~, Ю1 = ^

50 §0 т

8=1

vc1 — г(1 — v/т) , .

х,- = —----—-(ю1 — и,-), 1 = 1, 2,... ,т,

3 2v[cоCl — 1 + v(1 + с2)/тГ 1 * ' ' ' '

(57)

а 50, ю0 — величины (42).

Доказательство. В силу Условий (54) и Леммы 2 при е = 0 однородный цикл (41) системы (49) экспоненциально орбитально устойчив. Поэтому при малых по модулю значениях е он переходит в экспоненциально орбитально устойчивый цикл вида

Zj = Zj (е) ехр[г ю(е)£], 3 = 1,2,... ,т, (58)

где комплексные амплитуды Zj (е), Zj (0) = 50 и частота ю(е), ю(0) = Ю0 аналитически зависят от е. Далее подставим формулы (58) вместе с разложениями

Zj (е) = 50 + £%1,з + ..., т(£) = ®0 + е Ю1 + ... (59)

в систему (49), сократим результат на ехр)[ги приравняем коэффициенты при е в левой и правой частях получившихся соотношений. В результате для постоянных € С, Ю1 € М из (59) приходим к линейной неоднородной системе

• • лит , ь \ , v(т + 1) 1 г Ю1 = г и — а,0Ъ(Ъл + 51,]) +--т-015051,3 —

v(т — 1) * ьЕ 2v . 10

--а1$0$1Л--0,150} у51,з, 3 = 1,2,..., т.

т т

8=1

Последняя же, как нетрудно проверить, при значении Ш1, определенном соответствующей формулой из (57), имеет решение

х,-

5и = ^, 3 = 1, 2,...,т. (60)

5

0

Из установленных соотношений (59), (60) следует, что цикл (58) преобразуется к нужному виду (55)-(57). Необходимо только заметить, что выполнения фигурирующего в (55) равенства fi(£) = 0 всегда можно добиться посредством замены t на t + с, с £ R. Добавим, наконец, что при е = 0 цикл (55) является канонической химерой. Действительно, равенства вида (52) справедливы для него по тем же причинам, что и аналогичные равенства (10) для цикла С(е) системы (7) (см. соответствующее место в разделе 1). Что же касается требований

(53), то их справедливость при е = 0 гарантируют условие (51) и формулы (56), (57). Теорема 2 доказана. □

Обратим внимание на следующее обстоятельство. Поскольку при е = 0 система (49) по-прежнему инвариантна относительно замен (45), то каноническая химера (55) порождает целое семейство U, состоящее из 2т — 1 экспоненциально орбитально устойчивых циклов. Все эти циклы, получающиеся из одного цикла (55) посредством указанных замен переменных, будем называть химероподобными структурами. Ясно также, что если положить

Pj = v(s)ls=j/m, j = 1,2, ...,т, (61)

где ц(з) — некоторая непрерывная на отрезке 0 ^ s ^ 1 функция, то при выполнении аналогичных

(54) условий

0 < V < 1, coci > 1 (62)

и при согласованном стремлении е ^ 0, т ^ количество сосуществующих в системе (49) устойчивых циклов неограниченно растет. Иными словами, наблюдается известное явление буферности. Как показано в монографиях [16, 17], данное явление характерно для широкого класса динамических систем из различных областей естествознания.

Завершая обсуждение стратегии охоты на химер применительно к системе (49), отметим, что за счет подходящего выбора коэффициентов ^j мы можем гарантировать существование химеры с заранее запланированными свойствами. Действительно, предположим, что массив индексов 1 ^ j ^ т произвольным образом разбит на к, к ^ 2 отрезков, то есть

{1,2,...,т} = [1, ml] U [ml + 1,т2 ] U ... U [ms + 1, ms+l] U ... U [mk-l + 1,т]. (63)

Считаем далее, что на каждом из этих отрезков последовательность ^j либо постоянна, либо нетривиально зависит от j, и такого рода отрезки чередуются. В этом случае Теорема 2 остается в силе при условии, что мы надлежащим образом модифицируем определение канонической химеры. А именно в данном случае таковой уместно назвать химеру, у которой промежутки когерентного и некогерентного поведения осцилляторов совпадают с отрезками из (63) постоянства и непостоянства коэффициентов ^j.

3. Непрерывные химеры

В данном разделе установленные выше результаты для модельной системы (35) распространяются на ее непрерывный аналог. А именно исследуется эволюционное уравнение вида

i

z = z — d0lzl2z — vdlzj z2(t,s)ds, (64)

o

получающееся из (35) при m ^ Здесь, как обычно, точкой обозначена производная по t,

z = x(t, s) + iy(t, s), t ^ 0, s £ [0,1], x(t, s),y(t, s) £ R, d0 = 1 — icq,

(65)

c0 = const > 0, v = const > 0, dl = 1 + icl, cl = const £ R.

Остановимся сначала на вопросе о разрешимости задачи Коши для уравнения (64) с начальным условием из подходящего фазового пространства. В качестве такового возьмем вещественное банахово пространство Е, элементами которого являются векторы

5 = (x(s),y(s)): x(s), y(s) eL^(0,1). (66)

Что же касается нормы элемента (66), то ее зададим формулой

||5|| = ess sup УХ^^^Ту2^, (67)

где ess sup — существенная точная верхняя грань. Далее, введем в рассмотрение абстрактную функцию 5(t) = (x(t, s), y(t, s)) со значениями в Е, где x(t, s), y(t, s) — функции из (65). В результате уравнение (64) записывается в абстрактной форме

5 = F (5), (68)

где, как нетрудно увидеть, нелинейный оператор F : Е ^ Е есть сумма единичного оператора и слагаемого вида G(5,5,5), где G(■, ■, ) — непрерывная кубическая форма, действующая из

Е х Е х Е вЕ.

Элементы общей теории абстрактных уравнений вида (68) с ограниченной и гладкой по Фреше правой частью содержатся, например, в монографии [18]. Из этой теории вытекает, в частности, следующее утверждение.

Лемма 3. По любому фиксированному ограниченному множеству Q С Е можно указать такое t0 = t0(Q) > 0, что при V 5ü е Q решение 5 = 5(t) уравнения (68) с начальным условием

5|i=ü = 5о (69)

однозначно определено на отрезке 0 ^ t ^ t0.

В дальнейшем нам потребуется еще одно стандартное утверждение из теории абстрактных уравнений вида (68). Перед его формулировкой для решения 5 = 5(t) задачи Коши (68), (69) определим максимальный полуинтервал существования [0, tmax), где

¿max = sup {t0 : 5(t) существует на отрезке [0, t0]}. (70)

Согласно Лемме 3, множество значений to, по которому берется sup в (70), заведомо не пусто. Кроме того, любое решение 5(t) однозначно продолжается на свой максимальный полуинтервал. Если же оказалось, что tmax < ж, то справедлива следующая

Лемма 4. Предположим, что для некоторого решения 5(t) уравнения (68) конечна величина (70). Тогда имеет место предельное равенство

lim о ||5(i)|| = +ж, (71)

t—^tmax 0

где здесь и далее || ■ || — норма (67).

Доказательство сформулированной леммы проводится от противного. Действительно, предположим, что соотношение (71) не выполняется. В этом случае существуют такая последовательность моментов времени tn, tn ^ tmax при п ^ +ж, и такая постоянная М > 0, что ||5(in)|| ^ М. Далее, в условиях Леммы 3 в качестве ограниченного множества Q возьмем {5(tn), п ^ 1}. Тогда, очевидно, решение 5(t) будет определено на отрезках вида [tn, tn +t0], п ^ 1, где t0 = t0(Q) > 0. Ясно также, что при достаточно больших п выполняется включение imax е (tn, tn+to). Последнее же противоречит определению tmax (см. (70)). □

Перейдем теперь к вопросу о диссипативности уравнения (68). Как оказывается, в отличие от дискретного случая (35), где условие диссипативности имеет вид (37), непрерывная система (64) диссипативна при каждом v > 0. Точнее говоря, справедлива очередная

Лемма 5. Существует постоянная К0 > 0 со следующими свойствами. По любому ограниченному множеству □ С Е найдется такое > 0, что при каждом Но £ ^ решение Н^) соответствующей задачи Коши (68), (69) определено на полуоси £ ^ 0, а при всех £ ^ удовлетворяет неравенству

1Ш1 < Во. (72)

Доказательство. Фиксируем произвольно ограниченное множество □ из Е и элемент Но £ Далее, обозначим через Н^) решение задачи Коши (68), (69), определенное на своем максимальном полуинтервале [0,£тах), где Ьтах — величина (70). В последующем нам потребуется также представление

НСО = (х(г, 8),у&з)) : х(г, з),у(Ъ з) £ Ьх(0,1) по 8 при V £ £ [0,^ах) (73)

и связанная с ним функция = х(Ъ,,з) + гу{Ь,,з), удовлетворяющая при всех значениях

£ £ [0,1тах) уравнению (64) (почти при всех § £ [0,1]).

Обоснование Леммы разбивается на два этапа. На первом из них будет установлена оценка сверху для функции

1

V(I) = ! 8)12й8, 0 < 1< 1тах. (74)

0

Получение требуемой оценки проводится следующим образом. Сначала дополним уравнение (64) комплексно сопряженным уравнением. Умножим затем первое из этих уравнений на ~г, а второе — на г. В результате после сложения получившихся выражений и последующего интегрирования по в £ [0,1] приходим к равенству

(t) - j lz(t,s)l4ds - v J z2(t,s)ds j .

V = 2\ V(г) — у |z(t,s)|4ds — V о о

Далее, из соотношения (75) и из формулы (74) очевидным образом следует, что

4л„ 1 ^ о/Т/ЛЛ т/2/

(75)

(t) - j lz(t,s)l4dsS]j

где

F < 2( V(t) - j lz(t,s)l4ds) < 2(V(t) - V2(t)), V(0) < Vo, o

Vo= sup ||Ho||2 > 0 lo€Q

(случай V0 = 0 тривиален, поскольку тогда Q = {0} и t) = 0). А отсюда, в свою очередь имеем

V(t) < V*(t), 0 < t<tmax, (76)

где V*(t),t ^ 0 — решение задачи Коши

V = 2(V - V2), V|t=o = Vq. (77)

Добавим еще, что

V*(t) < Ri V t ^ 0, lim V* (t) = 1, (78)

где Ri = Ri(Q) — некоторая положительная постоянная.

Второй этап обоснования Леммы связан с получением при почти всех значениях € [0,1] оценки сверху на функцию

№(I, з) = |г(1, з)|2, 0 < ¿шах. (79)

Для этого привлечем аналогичное (75) равенство

№(г, в) =2(Ш (г, в) — № 2(г, в)) — 2vRe[dl а(£) г2(Ь, 8)], (80)

где, как обычно, точка — дифференцирование по Ь, а функция а(Ь ) имеет вид

2( ,

a(t) = j z2(t, s)ds, 0 ^ t < tr

Далее, объединяя формулы (79), (80) с очевидным фактом |a(i)| ^ V(t) и уже установленной оценкой (76), приходим к выводу, что (почти всюду по s)

W(t, S) < 2[(1 + V |d^V*(t))W(t, S) - W2(t, S)j, 0 < t < imax,

def 2 (81)

W(0, s) ^Wcd=f sup ||5c||2.

?o€Q

И наконец, из (81) в силу теоремы о дифференциальных неравенствах при почти всех s £ [0,1] имеем

W(t, S) ^W*(t), 0 < t < imax, (82)

где W*(t),t ^ 0 — решение аналогичной (77) задачи Коши

W = 2[(1 + V ldilV*(t))W - W2], W|i=0 = Wo.

Отметим еще, что в силу (78) решение W* (t) этой задачи обладает свойствами

W*(t) ^R2 Vt ^ 0, \imW*(t) = 1 + v|di|, (83)

где R2 = R2(Q) > 0.

Подводя итог, перейдем в неравенстве (82) к существенной точной верхней грани по s £ [0,1]. В результате, принимая во внимание способ задания нормы в пространстве Е (см. (67)), для функции (73) получаем оценку

|Ш||2 ^W*(t), 0 < imax. (84)

Далее, из (84) автоматически следует, что imax = (в противном случае имеет место предельное равенство (71), противоречащее оценке (84)). Остается отметить, что в силу (83), (84) все траектории системы (68) при увеличении t втекают в шар (5 £ Е : ||5|| ^ Ro} при любом фиксированном R0 > \J 1 + v |d1|. Тем самым, при подходящем выборе t* = t*(Q) > 0 требуемое неравенство (72) заведомо выполняется на полуоси t ^ t*. Лемма 5 доказана. □

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

z = 50 exp( i ю01), (85)

где

50 = 1/V1 + v, ®0 = (С0 - vci)50. (86)

Как оказывается, этот цикл порождает целое семейство циклов уравнения (64).

Действительно, опираясь на представление

[0,1] = А и В, (87)

где А, В — произвольные измеримые по Лебегу непересекающиеся подмножества отрезка [0,1] положительной меры, рассмотрим функцию

т ( \ \ 1 при 8 е А, (88)

I А, В00 = ^ (88) I —1 при § е В.

Заметим далее, что поскольку уравнение (64) инвариантно относительно замены

1а, В(Ф ^ г, (89) то наряду с циклом (85), (86) оно допускает цикл

хЦ, в) = 1А, вехф^), (90)

который и будем называть периодическим режимом двухкластерной синхронизации. Все циклы (90), отвечающие любым разбиениям вида (87), а также однородный цикл (85), (86), объединим в семейство и (отождествляя, естественно, циклы, отличающиеся друг от друга по в лишь на множестве нулевой меры). Аналогом Леммы 2 в данном случае является следующее утверждение.

Лемма 6. При условиях (62) каждый цикл семейства и экспоненциально орбитально устойчив. При выполнении же хотя бы одного неравенства V > 1 или сосг < 1 все циклы данного семейства неустойчивы.

Доказательство. Как и при обосновании Леммы 2, сначала обратим внимание на тот факт, что любые два цикла семейства и переходят друг в друга под действием некоторой замены из класса (88), (89). Те самым, свойства их устойчивости одинаковы, а значит, мы можем ограничиться рассмотрением лишь однородного цикла (85), (86). Добавим еще, что поскольку для абстрактного уравнения (68) очевидным образом справедлива теорема Андронова-Витта об устойчивости цикла по первому приближению, то проблема устойчивости упомянутого периодического режима сводится к анализу соответствующей системы в вариациях. Полагая в (64) г = ^0(1 + Н(1, в)) ехр(1ш01), где

Н(Ъ, з) = Н1{Ъ, з) + гН2{Ъ, 8), Н1{Ъ, 8),Н2(г, 8) е М,

и отбрасывая нелинейные по Н,Н слагаемые, приходим к выводу, что интересующее нас уравнение в вариациях имеет вид

1

К = —йо^о(Н + h) + vdl^%0(h — Н) — 2vdl^0оJ Н^в^. (91)

о

Заметим, далее, что фазовое пространство Е уравнения (91), состоящее из пар

(}ц(8 ),н2(8)): Ну (8) еь^(0,1), 3 = 1, 2,

раскладывается в прямую сумму замкнутых линейных подпространств Е1, Е2. Первое из них двумерно и задается равенством

Е1 = {(Н^в),^^)) : Ну (в) = Ну, Ну е М, 3 = 1, 2}, (92)

а второе имеет вид

1

Е2 = {(Н1(8),Н2(8)): I Ну (8^8 = 0, 3 = 1, 2}. (93)

Кроме того, нетрудно проверить, что подпространства (92), (93) инвариантны для траекторий уравнения (91).

Обозначим через Ь : Е ^ Е линейный ограниченный оператор, порожденный правой частью уравнения (91) и действующий по правилу:

(К^ ),Ко(5 )) ^ (Ке&(К,К), 1т&(К,К)),

где

1

&(к,Ъ) = -йо11(к + К) + - К) - / Н^йв, К(в) = К1(в) + ^(в).

Из приведенных выше фактов следует, что, во-первых, ЬЕ^ С Е^, ] = 1, 2; во-вторых, спектры сужений Ь|^., ] = 1, 2 совпадают с собственными значениями матриц В1 и Во , соответственно, где

+ ¿0 + у^Л - У^1 ¿0 + У^Л

В1 = -50 - , Во = -50 - - - - . (94)

у а0 + V а0 + V а\) \а0 + V ¿1 а0 - V а\)

Остается добавить, что первая из матриц (94) имеет собственные значения А4 =0, А, 2 = -2, а вторая при условиях (62) является гурвицевой. В случае же V > 1 или С0С1 < 1 у Во существует собственное значение с положительной действительной частью. Лемма 6 доказана. □

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

Перейдем теперь к несимметричной непрерывной сети

1

* = (1 + М.)>*-^М2. - ^¡Щ .)*., та

0

получающейся из (49) при условиях (61) и при т ^ Здесь, как и в дискретном случае, ее М, |е| ^ 1. Далее, по аналогии с требованиями (50), (51) предполагаем существование таких «1, в2 : 0 < < в2 < 1, что непрерывная вещественная функция в) постоянна на отрезках [0,«1] и [«2,1], а на отрезке [ вв2] нетривиально зависит от в. Последнее означает, что при Усе М множество е [ 81, 8 2] : ) = с} имеет нулевую меру Лебега.

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

, 8) = г 1^) при § е [0,81], , 8) = хо(Ь) при § е [8 2,1], (96)

где г 1^), г2(Ь) — некоторые комплекснозначные функции. В случае же отрезка 8 \ = 8 = 82 предполагаем, что для Усе С, У £ е М множество значений (« е [ в 1, в 2] : , в ) = с} имеет меру ноль. Если рассматривать 2 = , 8) как континуальный массив нелинейных осцилляторов, то в силу (96) промежутки 0 1и 8 2 = з = 1 естественно назвать зонами когерентности. Что же

касается отрезка 81 = 5 = 82, то его назовем зоной некогерентности.

Аналогом Теоремы 2 для непрерывной модели (95) является следующая

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

1

г = р(з, е) ехр[%ф( з, е) + %ю(е)£], J ф( 8, e)ds = 0,

(97)

являющийся при £ = 0 канонической химерой. Здесь непрерывные по своим переменным и аналитические по е вещественные функции р(з, е), ф( 8, е), ю(е) при е ^ 0 допускают асимптотику:

р(в, е) = ро + ) + 0( е2), ф, е) = еу^) + 0( е2), ю(е) = ю0 + е«1 + 0( е2),

(98)

где

. . И,е ) , . 1т к(8) [ , . ,

Ро = §0, Р1(э) = —£-, щ(8) =-Го-, «1 = р(s)ds

?о §о У

1т х( в )

(99)

Ф) = — 1)(«1 — Ф»,

а §,о, юо — величины (86).

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

Завершая рассмотрение непрерывной модели (95), заметим, что при е = 0 она по-прежнему инвариантна по отношению к заменам переменных (88), (89). Поэтому каноническая химера (97) порождает континуальное семейство и экспоненциально орбитально устойчивых циклов, получающихся из (97) под действием указанных замен. Как и в дискретном случае, эти периодические режимы уместно назвать химероподобными структурами.

1

Заключение

В данном разделе приводятся результаты численного анализа системы (49). Цель этого анализа — выявление возможных типов химероподобных режимов, существующих в указанной системе при фиксированных параметрах V, С1, С2, удовлетворяющих условиям (62), и при различных значениях управляющего параметра > 0.

Прежде чем перейти непосредственно к описанию численных экспериментов, остановимся на характерных особенностях изучаемой системы. Назовем два ее цикла эквивалентными, если они переходят друг в друга под действием замены переменных вида (45). Ясно, что свойства их устойчивости и происходящие с ними возможные бифуркации одинаковы. В связи с этим все эквивалентные между собой циклы уместно объединить в соответствующие семейства. Одним из таких семейств является множество химероподобных структур и, количество которых равно 2т — 1 (что при т = 100 имеет порядок 103°). Напомним, что эти химеры получаются из канонической химеры (55) под действием замен (45). Существуют и другие семейства эквивалентных периодических режимов, отличные от и и также состоящие не менее чем из 2т — 1 элементов. Все это приводит к реализации так называемого феномена флуктуационного хаоса. Суть упомянутого феномена состоит в том, что ввиду узости бассейнов притяжений отдельно

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

Перед численным интегрированием системы (49) в ней был осуществлен переход к полярным координатам

Zj = р^- ехр( г ф^-), р^- > 0, ф^ е М, 3 = 1, 2,...,т. Далее, при фиксированных параметрах с0 = 1, с1 = 3, V = 0.5, т = 100, при

№ = К 5 ^з^/п

Ц(8 ) = <( 3(2 8 - 1)

при 0 < 8 < 1/3, при 1/3 < 8 < 2/3, при 2/3 ^ й ^ 1

и при различных значениях е > 0 вычислялся отрезок траектории получившейся системы для рj, фj, отвечающий промежутку времени 0 ^ Ь ^ 1000000 и начальным условиям

Pj (0) = 0.2(4 + 0.01 ), Фj (0) = 0, ] = 1,... ,т.

На этом пути удалось выявить следующие типы химероподобных структур.

При 0 < е < 0.077 наблюдаются устойчивые циклы химерного типа, являющиеся продолжениями по параметру е циклов семейства и. На рис. 1, а и Ь при е = 0.07 для одного из химерных циклов представлены зависимости переменных рj(¿) и ф.,+1(£) — фj (£) от индекса 1 ^ з ^ т. Подчеркнем, что поскольку данный цикл имеет автомодельный вид

г j = г® ехр(гЮ), = const е С,

const е М, 3 = 1,2,... ,т,

(100)

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

0.88 0.86 0.84 0.82 0.80 0.78 0.76 0.74

а

20 40

60

80

100

0.006 0.005 0.004 0.003 0.002 0.001 0

-0.001 Ь

20

40

60

80

100

Рис. 1. а — р^', Ь — ф^+1 — ф^'; £ = 0.07 Fig. 1. а — р^', Ь — ф^+1 — ф^'; £ = 0.07

0

При последующем увеличении параметра е удается обнаружит так называемые двумерные автомодельные торы, имеющие вид

Zj = z® (t )exp( i wt ), z® (t) G C, m = const G R, j = 1, 2,...,m,

(101)

где г® (1) — некоторые периодические с периодом Т > 0 комплекснозначные функции. Добавим еще, что торы (101) заведомо не приводятся к виду (100). Такого типа химероподобные структуры будем называть квазипериодическими. На рис. 2, a, Ь изображены зависимости от 1 ^ .] ^ т отвечающих одному из этих торов функций р^ (Ь) и ф^^) — ф^- (Ь) при е = 0.08 и при фиксированном £ = 1000000. Подчеркнем, что в отличие от предыдущего случая этот химерный режим уже не является стационарным. В частности, в силу (101) отвечающие ему полярные радиусы р^ (Ь) периодичны по £ с периодом Т. График зависимости от £ на промежутке 999500 ^ ^ £ ^ 1000000 одного из этих радиусов, а именно р5(Ь), представлен на рис. 3 (точке ноль на горизонтальной оси соответствует значение = 999500).

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

1.1

1.0

0.9

0.8

20

40

60

80

100

106

-106

b

20

40

60

80

100

Рис. 2. a — pj, b — ф^+1 - ф^; е = 0.08, t = 1000000 Fig. 2. a — pj, b — ф^+i - ф^; e = 0.08, t = 1000000

0

0

0

a

0.8875 0.8850 0.8825 0.8800 0.8775 0.8750 0.8725 0.8700

Рис. 3. График рв(i) при e = 0.08 Fig. 3. Graph of рв(£) for e = 0.08

1.2

1.1

1.0

0.9

0.8

0.7

0 20 40 60

80 100

400000

300000

200000

100000

20 40 60 80 100

b

Рис. 4. a — р}, b — ф^+1 - ф^; e = 0.18, t = 1000000 Fig. 4. a — pj, b — ф^+i - ф^; e = 0.18, t = 1000000

0

a

0.86 0.84 0.82 0.80 0.78 0.76 0.74 0.72

0 200 400 600 800 1000

Рис. 5. График p5(t) при e = 0.18

Fig. 5. Graph of pg (i) for e = 0.18

Pj(t) от t. Такого рода химеры назовем турбулентными. На рис. 4, a, b для одной из турбулентных химер в случае е = 0.18 приведены зависимости pj(t) и ф^+i (i) — ф^ (t) от j при t = 1000000, а на рис. 5 представлен график p5 (t) на промежутке 999000 ^ t ^ 1000000 (точке ноль на горизонтальной оси соответствует значение t = 999000).

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

Список литературы

1. Kuramoto Y, Battogtokh D. Coexistence of coherence and incoherence in nonlocally coupled phase oscillators // Nonlinear Phenomena in Complex Systems. 2002. Vol. 5, no. 4. P. 380-385.

2. Abrams D. M., Strogatz S. H.Chimera states for coupled oscillators // Phys. Rev. Lett. 2004. Vol. 93, no. 17. P. 174102. DOI: 10.1103/PhysRevLett.93.174102.

3. Panaggio M.J., Abrams D.M. Chimera states: coexistence of coherence and incoherence in networks of coupled oscillators // Nonlinearity. 2015. Vol. 28, no. 3. P. R67.

DOI: 10.1088/0951-7715/28/3/R67.

4. Sethia G. C., Sen A. Chimera states: The existence criteria revisited // Phys. Rev. Lett. 2014. Vol. 112, no. 14. P. 144101. DOI: 10.1103/PhysRevLett.112.144101.

5. Schmidt L., Krischer K. Clustering as a prerequisite for chimera states in globally coupled systems // Phys. Rev. Lett. 2015. Vol. 114, no. 3. P. 034101.

DOI: 10.1103/PhysRevLett.114.034101.

6. Laing C. R. Chimeras in networks with purely local coupling // Phys. Rev. E. 2015. Vol. 92, no. 5. P. 050904. DOI: 10.1103/PhysRevE.92.050904.

7. Laing C.R. Chimeras in networks of planar oscillators // Phys. Rev. E. 2010. Vol. 81, no. 6. P. 066221. DOI: 10.1103/PhysRevE.81.066221.

8. Zakharova A., Kapeller M., Scholl E. Chimera death: Symmetry breaking in dynamical networks // Phys. Rev. Lett. 2014. Vol. 112, no. 15. P. 154101.

DOI: 10.1103/PhysRevLett.112.154101.

9. Omelchenko I., Zakharova A., Hovel P., Siebert J., Scholl E. Nonlinearity of local dynamics promotes multi-chimeras // Chaos. 2015. Vol. 25, no. 8. P. 083104.

DOI: 10.1063/1.4927829.

10. Omelchenko I., Omelchenko O. E., Hovel P., Scholl E. When nonlocal coupling between oscillators becomes stronger: Patched synchrony or multichimera states // Phys. Rev. Lett. 2013. Vol. 110, no. 22. P. 224101. DOI: 10.1103/PhysRevLett.110.224101.

11. Sakaguchi H. Instability of synchronized motion in nonlocally coupled neural oscillators // Phys. Rev. E. 2006. Vol. 73, no. 3. P. 031907. DOI: 10.1103/PhysRevE.73.031907.

12. Hizanidis J., Kanas V., Bezerianos A., Bountis T. Chimera states in networks of nonlocally coupled Hindmarsh-Rose neuron models // International Journal of Bifurcation and Chaos. 2014. Vol. 24, no. 3. P. 1450030. DOI: 10.1142/S0218127414500308.

13. Zakharova A. Chimera Patterns in Networks: Interplay between Dynamics, Structure, Noise, and Delay. Berlin: Springer, 2020. 233 p. DOI: 10.1007/978-3-030-21714-3.

14. Глызин С. Д., Колесов А. Ю., Розов Н. Х.Релаксационные автоколебания в сетях импульсных нейронов // УМН. 2015. Т. 70, № 3(423). С. 3-76. DOI: 10.4213/rm9659.

15. Глызин С. Д., Колесов А. Ю., Розов Н.Х. Периодические режимы двухкластерной синхронизации в полносвязных генных сетях // Дифференциальные уравнения. 2016. Т. 52, № 2. С. 157-176. DOI: 10.1134/S0374064116020035.

16. Колесов А. Ю., Розов Н.Х. Инвариантные торы нелинейных волновых уравнений. М.: Физ-матлит, 2004. 408 с.

17. Мищенко Е. Ф., Садовничий В. А., Колесов А.Ю., Розов Н.Х. Автоволновые процессы в нелинейных средах с диффузией. М.: Физматлит, 2010. 400 с.

18. Далецкий Ю.Л., Крейн М.Г. Устойчивость решений дифференциальных уравнений в банаховом пространстве. М.: Наука, 1970. 535 с.

References

1. Kuramoto Y, Battogtokh D. Coexistence of coherence and incoherence in nonlocally coupled phase oscillators. Nonlinear Phenomena in Complex Systems. 2002;5(4):380-385.

2. Abrams DM, Strogatz SH. Chimera states for coupled oscillators. Phys. Rev. Lett. 2004;93(17): 174102. DOI: 10.1103/PhysRevLett.93.174102.

3. Panaggio MJ, Abrams DM. Chimera states: coexistence of coherence and incoherence in networks of coupled oscillators. Nonlinearity. 2015;28(3):R67. DOI: 10.1088/0951-7715/28/3/R67.

4. Sethia GC, Sen A. Chimera states: The existence criteria revisited. Phys. Rev. Lett. 2014;112(14): 144101. DOI: 10.1103/PhysRevLett.112.144101.

5. Schmidt L, Krischer K. Clustering as a prerequisite for chimera states in globally coupled systems. Phys. Rev. Lett. 2015;114(3):034101. DOI: 10.1103/PhysRevLett.114.034101.

6. Laing CR. Chimeras in networks with purely local coupling. Phys. Rev. E. 2015;92(5):050904. DOI: 10.1103/PhysRevE.92.050904.

7. Laing CR. Chimeras in networks of planar oscillators. Phys. Rev. E. 2010;81(6):066221. DOI: 10.1103/PhysRevE.81.066221.

8. Zakharova A, Kapeller M, Scholl E. Chimera death: Symmetry breaking in dynamical networks. Phys. Rev. Lett. 2014;112(15):154101. DOI: 10.1103/PhysRevLett.112.154101.

9. Omelchenko I, Zakharova A, Hovel P, Siebert J, Scholl E. Nonlinearity of local dynamics promotes multi-chimeras. Chaos. 2015;25(8):083104. DOI: 10.1063/1.4927829.

10. Omelchenko I, Omelchenko OE, Hovel P, Scholl E. When nonlocal coupling between oscillators becomes stronger: Patched synchrony or multichimera states. Phys. Rev. Lett. 2013;110(22):224101. DOI: 10.1103/PhysRevLett.110.224101.

11. Sakaguchi H. Instability of synchronized motion in nonlocally coupled neural oscillators. Phys. Rev. E. 2006;73(3):031907. DOI: 10.1103/PhysRevE.73.031907.

12. Hizanidis J, Kanas V, Bezerianos A, Bountis T. Chimera states in networks of nonlocally coupled Hindmarsh-Rose neuron models. International Journal of Bifurcation and Chaos. 2014;24(3): 1450030. DOI: 10.1142/S0218127414500308.

13. Zakharova A. Chimera Patterns in Networks: Interplay between Dynamics, Structure, Noise, and Delay. Berlin: Springer; 2020. 233 p. DOI: 10.1007/978-3-030-21714-3.

14. Glyzin SD, Kolesov AY, Rozov NK. Self-excited relaxation oscillations in networks of impulse neurons. Russian Mathematical Surveys. 2015;70(3):383-452.

DOI: 10.1070/RM2015v070n03ABEH004951.

15. Glyzin SD, Kolesov AY, Rozov NK. Periodic two-cluster synchronization modes in completely connected genetic networks. Differential Equations. 2016;52(2):157-176.

DOI: 10.1134/S0012266116020038.

16. Kolesov AY, Rozov NK. Invariant Tori of Nonlinear Wave Equations. Moscow: Fizmatlit; 2004. 408 p. (in Russian).

17. Mishchenko EF, Sadovnichii VA, Kolesov AY, Rozov NK. Autowave Processes in Nonlinear Media with Diffusion. Moscow: Fizmatlit; 2010. 400 p. (in Russian).

18. Daleckii JL, Krein MG. Stability of Solutions of Differential Equations in Banach Space. Providence: American Mathematical Society; 2002. 386 p.

Глызин Дмитрий Сергеевич — родился в Ярославле (1983), окончил Ярославский государственный университет (2004). Защитил диссертацию на соискание ученой степени кандидата физико-математических наук (2006). Область научных интересов: нелинейная динамика, математическое моделирование.

Россия, 150003 Ярославль, ул. Советская, 14

Ярославский государственный университет имени П. Г. Демидова E-mail: glyzin@gmail.com ORCID: 0000-0002-0701-622X AuthorlD: 749418

Глызин Сергей Дмитриевич — родился в Новгородской области (1957), окончил Ярославский государственный университет (1979). Защитил диссертацию на соискание ученой степени кандидата физико-математических наук в ВГУ (1991) и доктора физико-математических наук в НИЯУ МИФИ (2009) в области теории нелинейных колебаний. Профессор, заведующий кафедрой компьютерных сетей. Опубликовал более 100 научных работ. Приказом Министерства образования и науки награжден нагрудным знаком «Почетный работник высшего профессионального образования». Зам. главного редактора научного журнала «Моделирование и анализ информационных систем».

Россия, 150003 Ярославль, ул. Советская, 14

Ярославский государственный университет имени П. Г. Демидова E-mail: glyzin@uniyar.ac.ru, glyzin.s@gmail.com ORCID: 0000-0002-6403-4061 AuthorlD: 35140

Колесов Андрей Юрьевич — родился в городе Старый Оскол Белгородской области (1961), окончил Ярославский государственный университет (1984). Защитил диссертацию на соискание ученой степени кандидата (1988) и доктора (1991) физико-математических наук в Математическом институте имени В. А. Стеклова РАН в области дифференциальных уравнений и динамических систем. Профессор кафедры дифференциальных уравнений. Опубликовал более 200 научных работ, среди которых 6 монографий.

Россия, 150003 Ярославль, ул. Советская, 14

Ярославский государственный университет имени П. Г. Демидова E-mail: kolesov@uniyar.ac.ru ORCID: 0000-0001-5066-0881 AuthorID: 5544