CC BY
4
Известия высших учебных заведений. Прикладная нелинейная динамика. 2023. Т. 31, № 4 Izvestiya Vysshikh Uchebnykh Zavedeniy. Applied Nonlinear Dynamics. 2023;31(4)
Научная статья DOI: 10.18500/0869-6632-003054
УДК 517.9 EDN: YSXPTE
Динамика полносвязных цепочек из большого количества осцилляторов с большим запаздыванием в связях
С. А. Кащенко
Региональный научно-образовательный математический центр при Ярославском государственном университете им. П. Г. Демидова, Россия E-mail: [email protected] Поступила в редакцию 8.04.2023, принята к публикации 4.05.2023, опубликована онлайн 18.07.2023, опубликована 31.07.2023
Аннотация. Цель настоящего исследования — изучить локальную динамику полносвязных цепочек из большого количества осцилляторов с большим запаздыванием в связях. От дискретной модели, описывающей динамику большого количества связанных осцилляторов осуществлен переход к нелинейному интегродифференциальному уравнению, непрерывно зависящему от времени и пространственной переменной. Рассматривается класс полносвязных систем. Основное предположение заключается в том, что величина запаздывания в связях является достаточно большой. Это предположение открывает путь к использованию специальных асимптотических методов исследования. Выделены параметры, при которых реализуется критический случай в задаче об устойчивости состояния равновесия. Показано, что они имеют бесконечную размерность. Построены аналоги нормальных форм — нелинейные краевые задачи типа Гинзбурга-Ландау. В ряде случаев эти краевые задачи содержат и интегральные составляющие. Их нелокальная динамика описывает поведение всех решений исходных уравнений в окрестности состояния равновесия. Методы. Применительно к рассматриваемым задачам развивается методика построения квазинормальных форм на центральных многообразиях. Разработан алгоритм построения асимптотики решений, основанный на использовании квазинормальных форм для определения медленно меняющихся амплитуд. Результаты. Построены квазинормальные формы, определяющие динамику исходной краевой задачи. Получены главные члены асимптотических приближений для решений рассматриваемых цепочек. На основе приведенных утверждений выявлен ряд интересных динамических эффектов. Например, бесконечное чередование прямых и обратных бифуркаций при увеличении коэффициента запаздывания. Их отличительная особенность заключается в том, что они имеют две или три пространственные переменные.
Ключевые слова: краевая задача, динамика, запаздывание, осцилляторы, нормальная форма, устойчивость.
Благодарности. Исследование выполнено за счет гранта Российского научного фонда № 21-71-30011, https://rscf.ru/ project/21-71-30011/.
Для цитирования: Кащенко С. А. Динамика полносвязных цепочек из большого количества осцилляторов с большим запаздыванием в связях//Известия вузов. ПНД. 2023. T. 31, № 4. С. 523-542. DOI: 10.18500/0869-6632-003054. EDN: YSXPTE
Статья опубликована на условиях Creative Commons Attribution License (CC-BY 4.0).
Article
DOI: 10.18500/0869-6632-003054
Dynamics of full-coupled chains of a great number of oscillators with a large delay in couplings
S.A. Kashchenko
Regional Scientific and Educational Mathematical Center of the Yaroslavl State University, Russia
E-mail: [email protected] Received 8.04.2023, accepted 4.05.2023, available online 18.07.2023, published 31.07.2023
Abstract. The subject of this work is the study of local dynamics of full-coupled chains of a great number of oscillators with a large delay in couplings. From a discrete model describing the dynamics of a great number of coupled oscillators, a transition has been made to a nonlinear integro-differential equation, continuously depending on time and space variable. A class of full-coupled systems has been considered. The main assumption is that the amount of delay in the couplings is large enough. This assumption opens the way to the use of special asymptotic methods of study. The parameters under which the critical case is realized in the problem of the equilibrium state stability have been distinguished. It is shown that they have infinite dimension. The analogues of normal forms — nonlinear boundary value problems of Ginzburg-Landau type have been constructed. In some cases, these boundary value problems contain integral components too. Their nonlocal dynamics describes the behavior of all solutions of the original equations in the balance state neighbourhood. Methods. As applied to the considered problems, methods of constructing quasinormal forms on central manifolds are developed. An algorithm for constructing the asymptotics of solutions based on the use of quasinormal forms for determining slowly varying amplitudes has been created. Results. Quasinormal forms that determine the dynamics of the original boundary value problem have been constructed. The dominant terms of asymptotic approximations for solutions of the considered chains have been obtained. On the basis of the given statements, a number of interesting dynamical effects have been revealed. For example, an infinite alternation of direct and reverse bifurcations when the delay coefficient increases. Their distinguishing feature is that they have two or three spatial variables.
Keywords: boundary value problem, dynamics, delay, oscillators, normal form, stability.
Acknowledgements. This work was supported by the Russian Science Foundation (project No. 21-71-30011), https://rscf.ru/ project/21-71-30011/.
For citation : Kashchenko SA. Dynamics of full-coupled chains of a great number of oscillators with a large delay in couplings. Izvestiya VUZ. Applied Nonlinear Dynamics. 2023;31(4):523-542. DOI: 10.18500/0869-6632-003054
This is an open access article distributed under the terms of Creative Commons Attribution License (CC-BY 4.0).
Введение
Сначала рассмотрим цепочку из N связанных осцилляторов второго порядка с запаздывающими связями
N
из + аи2 + и2 + ^(uj, и) = ак(£ - Т), (1)
к=1
^(и, у) — достаточно гладкая нелинейная функция, имеющая в нуле порядок малости выше первого, Т > 0 — запаздывание, индекс ] меняется от 1 до N и для любого целого к значения Uk+N (¿) отождествляются с ик (^. Такого типа модели возникают во многих прикладных задачах радиофизики [1-8], лазерной физики [9-13], математической экологии [14,15], теории нейронных сетей [16-21] и см., например, [22].
Удобно значения из (¿) ассоциировать со значениями функции двух переменных и(Ь,хз), где точки хз = 1,..., N) равномерно распределены на некоторой окружности.
Предполагаем, что количество элементов цепочки N является достаточно большим:
N » 1. (2)
Это условие дает основание перейти от дискретной относительно и(Ь, х) системы (1) к уравнению с распределенными на отрезке [0, 2п] параметрами для функции и(Ь, х)
2п
д 2и ди / ди\ 1^/4/^
+ аЖ + и + М ) = 2П / Ф(5)и(Ь- Т'х + (3)
о
с периодическими краевыми условиями
и(Ь,х + 2п)= и(Ь, х). (4)
Функция Ф(з) определяет связи между элементами и, вообще говоря, зависит от параметра е. Например, в случае односторонней [17] и диффузионной [2] связях значения Ф(з) сосредоточены в окрестности одной или нескольких точек отрезка [0, 2п]. Здесь предполагаем, что функция Ф(з) от параметра е не зависит. В этом случае цепочки принято называть полносвязными.
Основное предположение, которое открывает путь к применению асимптотических методов, заключается в том, что параметр запаздывания Т в (3) является достаточно большим:
е = Т-1 < 1. (5)
В краевой задаче (3), (4) проведем замену времени Ь ^ ТЬ. В результате приходим к сингулярно возмущенной краевой задаче
2п
2д2и ди / ди\ 1 [ . . .
е 52 + еа"дЬ + " + = 2п Ф(5)и(Ь - 1,х + (6)
о
и(Ь,х + 2п) = и(Ь, х). (7)
Поставим задачу об исследовании поведения всех решений из некоторой достаточно малой и не зависящей от е окрестности нулевого состояния равновесия в (6), (7) при достаточно малых значениях е. При изучении локальной динамики (6), (7) центральное место занимает анализ вопросов устойчивости решений линеаризованного в нуле уравнения
2п
2 д2и ди 1
е —тг + еа——+ и = — / Ф(5)и(Ь — 1, х + 5)аз дЬ2 дЬ 2п ,] у ' у ' 7
(8)
с краевыми условиями (7).
Рассмотрим характеристическое уравнение для (8), (7)
е2Х2 + еаХ + 1 = /к ехр(—X), (9)
где /к — коэффициенты Фурье функции Ф(з):
2п
1
/к = 2л / Ф(5) ехр(—'¿М^, & = 0, ±1, ±2,
В том случае, когда все корни уравнения (9) при всех к = 0, ±1, ±2,... имеют отрицательные и отделенные от нуля при е ^ 0 вещественные части, решения краевой задачи (8), (7) асимптотически устойчивы, и решения (6), (7) с достаточно малыми и не зависящими от е (по норме С [о,2П] х С (В2)) начальными условиями стремятся к нулю при Ь ^ <. Если же уравнение (9) имеет корень с положительной и отделенной от нуля при е ^ 0 вещественной частью, то решения (8), (7) неустойчивы и задача о динамике (6), (7) становится нелокальной.
Здесь будет рассматриваться критический случай, когда в (9) нет корней с положительной и отделенной от нуля вещественной частью, но есть корни, которые стремятся к мнимой оси при е ^ 0. Отметим, что в случае конечной размерности критического случая методика исследования локальной динамики хорошо известна. Она опирается на метод интегральных многообразий и метод нормальных форм (см., например, [23-25]). Характерной особенностью всех рассматриваемых ниже задач является то обстоятельство, что в них реализуются бесконечномерные критические случаи, когда бесконечно много корней характеристического уравнения стремится к мнимой оси при е ^ 0. Поэтому методы интегральных многообразий и нормальных форм непосредственно неприменимы. Здесь существенно используется разработанный в [17,26-30] подход, связанный с построением бесконечномерных квазинормальных форм.
В разделах 1-3 представлены исследования наиболее важных критических случаев. В качестве основных результатов построены так называемые квазинормальные формы, нелокальная динамика которых определяет поведение всех решений исходной краевой задачи (6), (7) в малой окрестности состояния равновесия.
В разделе 4 рассмотрен существенно более сложный по сравнению с разделами 1-3 критический случай, который реализуется для несколько иной, но близкой к (6), (7) краевой задаче
,д 2 и д2
ди д
ди д
2п
2 д и ди ди 1
е2 + еа— + и + и,е— ) = — / Ф(в)
и(Ь — 1,х + в) —
—и(Ь — 1, х)
¿в, и(Ь,х + 2п) = и(Ь,х).
(10)
Содержание раздела 5 посвящено обобщению результатов предыдущих разделов на более общую краевую задачу
2 д2и ди ( ди\ 1 [
е ~о2 + еа~о1 + и + I = ^ ФЦв)и(Ь — 1,х + в)ав+
ди(Ь — 1, х + в)
о
2п
+ 2П ф2<8>-
дЬ
¿8,
(11)
и(Ь,х + 2п) = и(Ь, х). (12)
Сделаем одно упрощающее предположение относительно нелинейной функции /(и, и). Выше говорилось, что эта функция достаточно гладкая и имеет в нуле порядок малости выше первого. Чисто технически проще использовать нелинейность /(и, и) третьего порядка малости, то есть для простоты ниже считаем, что в (1)
/(и, и) = ^1и3 + Ь2и2й + Ь$ии2 + Ь^и3. (13)
Опишем здесь структуру всех следующих разделов. Она для всех них одинакова. Сначала выделяются параметры задачи, при которых наблюдается критический случай в задаче об устойчивости нулевого состояния равновесия. Затем рассматривается линеаризованная краевая задача и приводится ее характеристическое уравнение. После этого исследуется асимптотика всех тех корней характеристического уравнения, вещественная часть которых стремится к нулю при стремлении к нулю малого параметра е. Таких корней оказывается бесконечно много. На их основе для линеаризованной задачи строится совокупность специальных решений. Такие решения удается записать в форме, которая допускает использование их для анализа решений (с неизвестными амплитудами) исходной нелинейной краевой задачи. Удается определить явную форму для главного приближения (по параметру е) соответствующего решения. Условно обозначим ее здесь через eUi. Решения нелинейной краевой задачи ищем тогда в виде
x, е) = eUi + e3U3 + ....
Заметим, что отсутствие здесь квадратичных по е коэффициентов связано с тем, что в исходном уравнении нет квадратичной нелинейности. Относительно U3 заранее известно, что она периодическая по нескольким своим аргументам. Подставляя вместо u(t, x, е) приведенное выражение, для определения U3 приходим к специальной линейной неоднородной краевой задаче. Условие разрешимости этой краевой задачи в указанном классе функций позволяет выписать уравнение для неизвестных амплитуд, входящих в Ui. Получение таких уравнений и является конечной целью. Нелокальная динамика этих уравнений, их еще называют квазинормальными формами, позволяет описать локальное поведение решений исходной краевой задачи. Отметим, что выполнение условий разрешимости уравнений для U3 позволяет в явном виде определить эту функцию. Ниже будем использовать функцию U3, но формулы для этой функции для краткости изложения иногда приводить не будем.
1. Описание критических случаев
Рассмотрим характеристическое уравнение (9). Через / обозначим наибольшее из чисел |/
f = max |/1. (14)
Воспользуемся методикой работ [26-28,31]. Введем в рассмотрение величину
1, если а2 ^ 2,
Отметим, что
Y0 * а2 , 2\ 1/2 2 .
— (4 - а2} , если а2 < 2.
2
у0 ^ 1 и 0 < у0 < 1 при а < 2
у0 = 0 при а = 0. В [28] показано, что имеют место следующие утверждения. Лемма 1. Пусть выполнено неравенство
/ < Уо.
Тогда при всех достаточно малых е все корни (9) имеют отрицательные и отделенные от нуля при е ^ 0 вещественные части.
и
Лемма 2. Пусть выполнено неравенство
/ > То-
Тогда при всех достаточно малых е уравнение (9) имеет корень с положительной и отделенной от нуля при е ^ 0 вещественной частью.
Отметим, что в условии леммы 1 все решения краевой задачи (8), (7) и (6), (7) из малой и не зависящей от е окрестности нуля стремятся к нулю при Ь ^ <. Если же выполнены условия леммы 2, то задача о динамике (6), (7) является нелокальной: в окрестности нулевого решения не может быть устойчивых решений.
В дальнейшем предполагается, что имеет место критический случай: при произвольном фиксированном у1 для параметра / выполнено либо равенство
либо
/ = Уо + е2у1,
/ = — Уо — е2у1.
(15)
(16)
Тем самым исследуются краевые задачи (8), (7) и (6), (7) в критическом случае.
Рассмотрим вопрос об асимптотике всех тех корней уравнения (9), вещественные части которых стремятся к нулю при е ^ 0. Соответствующие асимптотические разложения принципиально различаются для случая, когда а0 > 2 и для случая, когда а^ < 2.
Лемма 3. Пусть выполнено условие
а2 > 2. (17)
Тогда уо = 1 и для тех корней (е) (к = 0, ±1, ±2,...) уравнения (9), вещественные части которых стремятся к нулю при е ^ 0, имеют место асимптотические равенства:
1) в случае
/ = 1 + е2у1 (18)
имеем
где
\к(е) = 2лгк + е"кк1 + е2"кк2 + о(е3), к = 0, ±1, ±2,.. Хк1 = —2лгка, Ак2 = —2п2 к2( а2 — 2) + 2лгка2 + у1
2) в случае
имеем
где
/ = —1 — е2 У1
\к(е) = т{2к + 1) + еА,к1 + е2А,к2 + о(е3), к = 0, ±1, ±2,
(19)
\к1 = — л( 2к + 1) а, \к2 = — -п2( 2к + 1) 2( а2 — 2) + гп( 2к + 1) а2 + уь Рассмотрим затем случай, когда
0 < а2 < 2.
(20)
Введем несколько обозначений. Пусть 9 = 0(е) € [0,1) такое значение, которое дополняет
— 1
до целого выражение юое 1, где
0, если а2 > 2,
®о = <; ( а2\1/2 2
1 —— 1 , если а2 < 2.
Положим
Д = (га — 2ш0)у-1 ехр ( — гО0),
Д = 2Д + Уо1 ехр (— гОо),
Д3 = (а — 2гш0) Д1у°1 ехр (— гО0).
Лемма 4. Пусть выполнены неравенства (20). Тогда для тех корней Хк(е) (к = 0, ±1, ±2,...) уравнения (9), вещественные части которых стремятся к нулю при е ^ 0, имеют место асимптотические равенства
Хк (е) = г( ш0е-1 + 6 — О0 + 2лк) + еХк1 + е2Хк2 + о(е3),
где
Хк1 = 6 — О0 + 2лк),
Хк2 = Д2(6 — О0 + 2лк)2 + Д3(6 — О0 + 2лк) + у1 ехр (— гО0), причем выполнены соотношения
Ие Д1 = 0, Ие Д < 0. (21)
2. Критический случай на нулевой моде
Здесь предполагаем, что
У0 = / и |/0| > |/к | (к = ±1, ±2,...). (22)
Исследование проведем отдельно для случаев (17) и (20). Сразу отметим, что соответствующие результаты для них существенно различны. При условии (17) решения формируются на частотах порядка 1, их назовем медленно осциллирующими. При условии (20) решения содержат частоты порядка е-1, поэтому назовем их быстро осциллирующими.
2.1. Медленно осциллирующие решения. Пусть сначала выполнено неравенство (17), то есть а2 > 2 и у0 = 1. Тогда тем корням Хк (е), о которых речь шла в лемме 3, отвечают решения Ук(¿, е) линейной краевой задачи (8), (7) и Ук(¿,е) = ехр(Хк(е)£), а значит, решением (8), (7) являются функции
то
у(г,е)= ^ 1к ехр (Хк(е)г),
к=-то
где ^к _ произвольные комплексные постоянные. Для случаев (15) и (16) это выражение можно представить соответственно в виде
то
е) = ^ 1к(х) ехр (2лгку) = £(х, у)
к=-то
е) = ^ 1к(х) ехр (лг(2к + 1)у) = £(х, у),
к=-то
где х = е2£, у = (1 — еа)£, ^к(х) = ехр ((Хк2 + о(е))х).
и
то
Решения и(Ь, у, е) нелинейной краевой задачи (6) и (7), «близкие к критическим» решениям у(Ь, е) линейной краевой задачи (8), (7), ищем в виде
и(Ь, у, е) = е1/2^(т, у) + е3/2из(т, у) + о(е2), (23)
где ^(т, у) - неизвестная вещественная функция, для которой выполнены условия:
1) в случае (15) — условие периодичности по у
£(т,у + 1) = £(т,у); (24)
2) в случае (16) — условие антипериодичности по у
£(т,у + 1) = —£(т,у). (25)
Подставим формальный ряд (23) в краевую задачу (8), (7) и будем приравнивать коэффициенты при одинаковых степенях е. Тогда, собирая коэффициенты при е3/2, получим уравнение для «3(т, у). Из условия его разрешимости приходим к выводу, что ^(т, у) является решением краевой задачи
3 = (а- — 1) I2! — а» ? + — «3 (26)
дт V 2 / ду2 ду
с краевыми условиями соответственно (24) или (25). Сформулируем итоговый результат.
Теорема 1. Пусть выполнены условия (17), (18) ((17), (19)) и Ь1 = 0. Пусть функция ^(т,у) является ограниченным при т ^ то, у € [0,1] решением краевой задачи (26), (24) ((26) и (25)). Тогда функция
и(Ь, у, е) = е£(т, у) + е3^3(т, у) (27)
при т = е2Ь, у = (1 — еа)Ь удовлетворяет краевой задаче (6), (7) с точностью до о(е3).
Из этой теоремы следует, что при сформулированных условиях построенные краевые задачи (26), (24) и (26), (25) играют роль нормальных форм для краевой задачи (6), (7). Отметим, что если 61 = 0, то есть условие теоремы 1 не выполняется, и 62 = 0, то изменения невелики. Последнее слагаемое в (26) заменяется на —Ь2?2д^/ду и формула (27) принимает вид
и(Ь, у, е) = е1/2^(т, у) + е3/2и3(т, у).
2.2. Быстро осциллирующие решения. Пусть выполнено неравенство (20), то есть а2 < 2. Корням \к(е) (к = 0, ±1, ±2,...), о которых говорилось в лемме 4, отвечают решения Ук(Ь,е) линейной краевой задачи (8), (7) Ук(Ь,е) = ехр(Хк(е)Ь). Значит и функция
то
у(Ь,е)= ^ 1к ехр(Хк(е)Ь), (28)
к=—то
где ^к — произвольные комплексные постоянные тоже удовлетворяют краевой задаче (8), (7). Учитывая асимптотические формулы для Хк (е), представленные в лемме 4, выражение (28) можно записать в виде
у(Ь,е) = Е (Ь,еЖт,у).
Здесь положено е) = ехр [(г(ю0е 1 + 6 — О0) + еД1(6 — О0))Ц ,
то
?(х,у) = ^ 1к(х)ехр(2пгку), х = е2Ь, у = (1 + геЯ^Ь,
к=-то
(х) = ехр ((Хк2 + о(е))х). Напомним, что согласно (21) значение гД вещественно.
Решения и(Ь, у, е) нелинейной краевой задачи (6), (7) в рассматриваемом случае ищем в виде формального ряда
и(Ь,у,е) = е(?(х,у)£(Ь,е) + |(х,у)£(Ь,е)) + е3^, х,у) + .... (29)
В этом выражении ?(х, у) - неизвестная, подлежащая определению комплексная функция, которая является 1-периодической по пространственной переменной у:
?(х,у + 1) = ?(х,у), (30)
а зависимость от аргумента Ь функции из - периодическая.
Подставим (29) в (6), (7) и будем приравнивать коэффициенты при одинаковых степенях е. Тогда, собирая коэффициенты при е3, приходим к уравнению относительно из. Условие разрешимости этого уравнения в классе периодических по Ь функций состоит в выполнении равенства
д? д 2^ д? 2 д?
ТГ = — + Я4т;--+ — ? о--
дх ду2 ду ду
— (3&1 + г^ + ^0^3 + зг^) ехр(—го0)?|?|2, (31)
в котором Я4 = —г(2Я2(6 — О0) + Д3), Я = Д2(6 — О0)2 + ^3(6 — О0) + У1 ехр(—гО0).
Введем обозначение. Здесь и ниже через еп = еп(60) будем обозначать такую последовательность еп ^ 0, для которой 6(еп) = 60. Сформулируем основной результат.
Теорема 2. Пусть выполнены условия (20) и (15) и пусть для произвольно фиксированного значения 60 € [0, 2п) краевая задача (31), (30) имеет ограниченное при х ^ то, у € [0,1] решение ?(х,у). Тогда на последовательности еп = еп(60) функция
и(Ь,у,еп(60)) = еп(60)[?(х,у)Е (Ь,еп(60)) + |(х,у)Ё (Ь,еп (60))] + еП(60 х,у)
при х = еп(60)^, у = (1 — еп(60)а)Ь удовлетворяет краевой задаче (6), (7) с точностью до
о(еп(60)).
Таким образом, краевая задача (31), (30) является квазинормальной формой для (6), (7) в рассматриваемом критическом случае. В силу условий (21) эта краевая задача является параболической. Структура ее решений может быть сложной. Эта краевая задача определяет главные члены асимптотических приближений решений исходной задачи краевой задачи (6), (7). Отсюда следует, что и локальная динамика (6), (7) также может быть сложной. Отметим еще, что по форме кубические нелинейности уравнения (31) сложнее, чем в классическом уравнении Гинзбурга-Ландау.
2.3. Система с малыми коэффициентами связей. Выше было показано, что критический случай для краевой задачи (6), (7) при a = 0 реализуется при f = 0. Здесь предполагаем, что коэффициенты a и f являются достаточно малыми, то есть для некоторых ai и yi имеем соотношения
a = ea1, f = еу1 (0 <s < 1). (32)
Ниже рассмотрим задачу именно такого типа, то есть с краевыми условиями (7) исследуем уравнение
2п
d2 u du s f
"dt2 + sai"dt + u + f (u, u) = — Ф(й)и(£ - s-1, x + (33)
0
Для линеаризованного в нуле уравнения характеристический квазиполином имеет вид
X2 + sa1X + 1 = sfk exp (— Xs-1).
Бесконечно много корней Xk(s) (k = 0, ±1, ±2,...) этого характеристического уравнения стремится к мнимой оси при s ^ 0 и нет корня с положительной и отделенной от нуля при s ^ 0 вещественной частью. Тем самым в задаче об устойчивости нулевого решения (33) реализуется критический случай бесконечной размерности. Для Xk(s) имеют место асимптотические представления
Xk (s) = i + sXfc! + ...,
а Xk1 является корнем квазиполинома
1 i
Xk1 + 2a1 = - 2Y1 exp (i0 - ТД^). (34)
Фигурирующая в (34) величина 0 = 0(s) € [0, 2п) дополняет значение Ts-1 до целого кратного 2п. Отметим, что квазиполином (34) имеет бесконечно много корней. Решения нелинейного уравнения ищем в виде формального ряда
u
= s1/2 (£(т, x) exp(it) + £(т, x) exp(-it)) + s3/2u3(t, т, x) + ..., (35)
в котором зависимость от аргумента Ь - 2л-периодическая. Подставим (35) в (33). Производя стандартные действия, приходим к уравнению для из. Из условия его разрешимости в классе 2п-периодических по Ь функций получаем уравнение с фиксированным запаздыванием для определения неизвестной комплексной амплитуды ^(х,х):
2п
д^ 1 г 1 [
дТ = -2а1 ^ - 2 ехр(г0) • 2П ф(5)^т - 1, х +
+2 (3ib1 - b2 + ib3 - 3b4)^|2. (36)
Сформулируем основной результат.
Теорема 3. Пусть выполнено условие (32) и для произвольно фиксированного значения 0о€[О, 2п) уравнение (36) имеет ограниченное при т ^ то решение ^(х, х). Тогда на последовательности £п = £п(0о) функция
и(Ь, х, еп) = вП/2 (?(х, х) ехр(гЬ) + |(х, х) ехр(-¿¿)) + еП/2из(Ь, х, х) при х = удовлетворяет уравнению (33) с точностью до о(еП/2).
Из этой теоремы следует, что распределенное уравнение (36) является квазинормальной формой в рассматриваемом случае. Для уравнения (36) просто исследовать вопросы о существовании и устойчивости простейших циклов вида р ехр(гох). Подробнее на этом не останавливаемся.
3. Критический случай на ненулевой моде
Здесь предполагаем, что наибольший по модулю коэффициент Фурье f функции Ф(з) имеет отличный от нуля номер ko, то есть
f = max |/fc | = |/fco ko = 0
и
Iffco I > l/fc I при k = ±ko.
Рассмотрим критический в задаче об устойчивости нулевого состояния равновесия краевой задачи (6), (7) случай, когда для некоторого 8
ffco = (Уо + £2Yi) exp(i8). (37)
Приведем асимптотические формулы для всех тех корней характеристического уравнения (9), вещественные части которых стремятся к нулю при е ^ 0.
Напомним, что характеристическое уравнение (9) получается в результате подстановки в (8) решений Эйлера
u±(t, x, е) = exp [± ikx + (е)^]. (38)
Для нахождения корней (е) приходим к уравнению
е2Х2 + еаХ + 1 = (уо + e2yi) exp (± ¿8 - X). (39)
Отсюда получаем асимптотику (е): 1) при а2 > 2 имеем Х-(е) = Х+(е) и Х+(е) = ¿8 + 2nki + еХк1 + е2Хк2 + ..где
Xfci = —¿а( 2nk + 8,
Xfc2 = 1 (2 — а2) (8 + 2nk)2 + yi + ¿а2 (8 + 2nk);
2) при а2 < 2 имеем Х±(е) = ¿(ю0/е + 6 ± 8 — Q0 + 2nk) + еХ±, + е2Х±2 + ..где
X±i = —¿(Yo exp^Qo)) ^2irno — а)К±, K± = 6 ± 8 — Qo + 2nk 2 V "fci; ню v4
Lfc = 6 ± 8 — Qo
X±2 = 1 (X±i)2 — YiY-1 — (R±)2 — aX±i — 2i®oX±i = (K±) 2 + B ±K ± + YiY-1 K± = 6 — Qo + 2nk ± 8,
= (yoexp (i(Qo ^ 8)^ — ^(а — 2irn)(Yoexp[i(Qo ^ 8)])-ij В± = —¿(а2 + 4ю2) (yo exp[i(Qo ^ 8)])-. Отметим, что Re X±i = Re = 0, Re Xk2 < 0 и Re X±2 < 0.
3.1. Случай a2 > 2. В этом случае построения в критической ситуации на ненулевой моде мало отличаются от построений раздела 2. Решения нелинейной краевой задачи (6), (7) ищем в виде
u(t, x, s) = s^^(т, y) exp [iS(1 - sa)t + ikox] + |(т, y) exp [ - iS(1 - sa)t - ikox] ) +
+s3u3(t, T,x,y) + ..., (40)
где по t,x и y зависимость периодическая, т = s21, y = (1 - sa)t. Подставляя (40) в (6), (7) и совершая стандартные действия, получаем следующий результат.
Теорема 4. Пусть ko =0 и a2 > 2 и пусть а(т, y) является ограниченным при т ^ то, y € [0,1] решением краевой задачи
д! = 1(a2 - 2)^ + (a2 + ig(a2 - 2))g + QS2(2 - a2) + Y1 + ia2s)ü - З^^2, (41)
£(т, x + 2n, y) = £(т, x, y) = £(т, x, y + 1). (42)
Тогда функция
u(t, x, s) = s^^(т, y) exp [iS(1 - sa)t + ikox] + 1(т, y) exp [ - iS(1 - sa)t - ikox] ^ +
+s3u3(t, т, x, y) (43)
удовлетворяет краевой задаче (6), (7) с точностью до o(s3).
Отметим, что при 61 = 0 и 62 = 0 последнее слагаемое в уравнении (41) заменяется на
- 62
«аа» + 3?211
1/2.
а асимптотика решений в (43) идет не по целым степеням е, а по степеням е
3.2. Случай а2 < 2. Здесь, в отличие от результатов раздела 2.2, задействованы не одна, а две цепочки корней (е) и Х- (е). Сначала введем обозначения. Положим
Е± = ехр [гкох + г(ю0е-1 + 0 - П ± 8)Ь] .
Решения нелинейной краевой задачи (6), (7) ищем в виде
и(Ь,х,е) = е(^+(х,у)£+ + |+(т,у)£+ + Г(х,у)£- + |"(х,у)£") +
+е3из (Ь, х,х,у) + ..., (44)
где зависимость от Ь, х и у — периодическая. После подстановки (44) в (6) и после стандартных действий получаем для определения неизвестных амплитуд ^±(х, у) краевую задачу параболического типа
f=- w- - '(2A± (0 - * ± »>+) f+
+ А±(0 - По ± 8)2 - В±(0 - По ± 8) + У1У-1+ + 3^±( ||2 + 2 |^|2) (61 - «2^3 - гЬ4) + гюоЬ£± |Ч± 12 (45)
с краевыми условиями
(х,у + 1) = ^±(х,у). (46)
Сформулируем полученное утверждение в виде теоремы.
Теорема 5. Пусть 0 < а2 < 2 и пусть для произвольно фиксированного значения 60 € [0, 2п) краевая задача (45), (46) имеет ограниченное при х ^ то, у € [0,1] решение ?(х, у). Тогда на последовательности еп = еп(60) функция
и(Ь, ж, еп) = еп (?+(х, у)Е + + !+(х, у)Е + + (х, у)Е- + I- (х, у)Е-) + епх, ж, у), где х = еп^, у = (1 — епа)Ь, удовлетворяет краевой задаче (6), (7) с точностью до о(еп).
4. Критические случаи в краевой задаче (10)
Принципиально новые эффекты могут возникать в ситуации, когда критические случаи реализуются одновременно на бесконечном множестве мод. Рассмотрим такую ситуацию на одном наиболее распространенном примере, когда в краевой задаче (10) имеем
Ф(з) = const = fo. (47)
Линеаризуем в нуле эту краевую задачу. В результате получим линейное уравнение
2п
2 д2u du
е + еаЖ + u = fo
1
— u(t — 1, x + s)ds — u(t — 1, x) 2n J
(48)
с периодическими краевыми условиями
и(Ь,ж + 2п) = и(Ь,ж). (49)
Характеристическое уравнение для (48), (49) имеет вид
е2Х2 + еаХ + 1 = /к ехр(—Х), к = 0, ±1, ±2,..., (50)
где
0, если к = 0,
/к Ч (51)
если к = ±1, ±2,....
Фиксируем произвольно 71 и полагаем
/0 = У0 + е2у1. (52)
При этом условии для всех, кроме одного, уравнений в (50) имеет место бесконечномерный критический случай. Следуя разработанной выше методике, решения нелинейной краевой задачи (10) ищем в виде формального ряда:
1) при условии а2 > 2 имеем у0 = 1 и
и(Ь, ж, е) = е?(х, ж, у) + е3и3(х, ж, у) + ... (53)
2) при условии а2 < 2 имеем
Yo = 0г(4 — а2)1/2
(¿,х,е) = е^^(x,x,y)exp [¿(^е i + 6 + 8 — Qo)t] +
+ !(х, x, y) exp [ — ¿(^е-1 + 6 + 8 — Qo)t] ) + + е3из (t, T,x,y) + ... (54)
и
В (53) и (54) т = е2Ь, у = (1 - еа)Ь, зависимость от х — 2п-периодическая, от у — 1-антипериодическая, от Ь в (54) - 2 л/шо-периодическая, где юо = (1 - а2/2)1/2. В силу условия (51), на функцию ^(т, х, у) необходимо наложить условие
2п
1
= 0, где = 2П / ?(т,х,у)^х. (55)
о
Подставим выражение (54) в (10). Совершая стандартные действия, получим уравнение относительно из. Из условия его разрешимости в указанном классе функций получим равенства:
1) при а2 > 2
i = (f - 0 0 - а21 + Y1? - Ч?3 - АУЛ). (56)
+ 1) = -?(х,ж,у), + 2n,y) = ?(т,ж,у); (57)
2) при а2 < 2 (в обозначениях раздела 2.2)
dT = -R2Й + R£ + ЗД - ?2д? - 2|?|2£ - M Г?2* + 2|?|2 *) (58)
ох oy2 oy oy oy V dy oy /
с условиями (57).
Приведем итоговые утверждения.
Теорема 6. 1) Пусть а2 > 2 и пусть ?(х, x, y) — ограниченное при х ^ то, x € [0, 2п], y € [0,1] решение краевой задачи (56), (57). Тогда функция
u(t,x,e) = е?(х, x,y) + e3U3(t, x,x,y)
удовлетворяет краевой задаче (48), (49) при х = e2t, y = (1 - еа)£ с точностью до о(е3).
2) Пусть 0 < а2 < 2 и пусть ?(х, x, y) — ограниченное при х ^ то, x € [0, 2п], y € [0,1] решение краевой задачи (58), (57). Тогда функция
u(t, x, е) = е^?(х, x, y) exp [¿(шое-1 + 6 + 8 - Qo)t] +
+ 1(х, x, y) exp [ - i(rnoe-1 + 6 + 8 - Qo)i]) + e3U3(t, х, x, y)
удовлетворяет краевой задаче (48), (49) при х = e2t, y = (1 - еа)£ с точностью до о(е3).
Таким образом, краевые задачи (56), (57) и (58), (57) являются квазинормальными формами для (48), (49). Отметим, что наличие в уравнениях (56) и (58) интегральных слагаемых по пространственной переменной позволяет в явном виде находить гладкие по х и y и ступенчатые по x решения. Детально эти вопросы исследовались в работах [32-34], поэтому здесь на них не останавливаемся.
5. Об одном обобщении результатов
Коротко остановимся на рассмотрении более общей краевой задачи (11), (12). Выпишем характеристическое уравнение для его линеаризации в нуле
е2А2 + еаА + 1 = / + еА/^) exp(-A), k = 0, ±1, ±2,..., (59)
где /к1 и /к2 — коэффициенты разложения функций Ф^з) и Ф2(^), соответственно, в ряд Фурье. Положим в (59) Х = гюе-1. В результате приходим к уравнению
(1 — ю2 + гаю) (/к1 + гю/к^-1 = ехр (— гюе-1). (60)
Изучим вопрос о существовании вещественного корня Ю0 в уравнении (60). Обозначим через Рк (ю) функцию
Рк(ю) = |1 — ю2 + гаю| ■ |/к1 + гю/к21-1. Эта функция неограниченно растет при ю ^ ±то. Поэтому существует такое ю0, что
Рк (ю0) = шт Рк (ю).
-то<к,ю<то>
Сформулируем один простой результат.
Лемма 5. Пусть р(ю0) < 1. Тогда при всех достаточно малых значениях е все корни уравнения (59) имеют отрицательные и отделенные от нуля при е ^ 0 вещественные части. Если же р(ю0) > 1, то при всех достаточно малых значениях е уравнение (59) имеет корень с положительной и отделенной от нуля при е ^ 0 вещественной частью.
Таким образом, критический случай в задаче об устойчивости реализуется при условии р(ю0) = 1. После этого приведенная выше методика переносится на краевую задачу (11), (12).
Заключение
Рассмотрен вопрос о локальной динамике полносвязной системы осцилляторов. Выделены критические случаи в задаче об устойчивости состояния равновесия. Показано, что они имеют бесконечную размерность. Основные результаты состоят в том, что разработан алгоритм построения специальных нелинейных краевых задач — квазинормальных форм. Их нелокальная динамика определяет асимптотику всех решений исходного уравнения в окрестности состояния равновесия.
Квазинормальными формами являются пространственно распределенные нелинейные краевые задачи, например, классические уравнения Гинзбурга-Ландау, поэтому можно сделать вывод о том, что для класса рассмотренных здесь задач характерны сложные и нерегулярные колебания.
В зависимости от величины параметра а — коэффициента при и в (1) — решения являются либо медленно осциллирующими (при а2 > 2), либо быстро осциллирующими (при а2 < 2) с асимптотически большой частотой. Квазинормальные формы из раздела 4 содержат еще одну пространственную переменную. Это, конечно, приводит к усложнению динамических свойств решений.
Показано, что количество связанных уравнений в квазинормальной форме определяется числом равных по модулю коэффициентов Фурье функции Ф(з). Важная роль принадлежит и значению аргумента 8 для соответствующих коэффициентов Ф(з).
Особо отметим, что в ряде квазинормальных форм присутствуют интегральные по пространственной переменной слагаемые от нелинейной функции. Это приводит к тому, что решения квазинормальных форм могут структурно усложниться. Например, в явном виде можно определить гладкие и периодические по времени, и ступенчатые по пространственной переменной решения и в некоторых случаях исследовать их устойчивость.
Обратим внимание, что в случае медленно осциллирующих решений для построения квазинормальных форм задействована только нелинейность ^1и3 функции /(и, и). Если 61 = 0 и 62 = 0, то происходит структурное усложнение решений и смена асимптотики решений исходного уравнения. Для быстро осциллирующих решений задействованы все коэффициенты функции /(и, и).
Динамические свойства исходной краевой задачи (3), (4) при а2 < 2 чувствительны к изменению параметров. Этот вывод следует из того факта, что, во-первых, в квазинормальной форме присутствует величина 6 = 6(е) € [0, 2п), которая при е ^ 0 бесконечно много раз пробегает все значения от 0 до 2п. Во-вторых, при различных значениях 6 динамика квазинормальной формы может различаться [35], а значит, при е ^ 0 может происходить неограниченный процесс прямых и обратных бифуркаций.
Важный вывод касается роли большого запаздывания в рассматриваемых цепочках. С одной стороны, все аналитические построения при больших T позволяют явно выделить критические случаи, найти асимптотику корней характеристических уравнений и получить асимптотические формулы для решений. С другой стороны, при T ^ 1, все квазинормальные формы содержат еще одну пространственную переменную, поэтому можно сделать вывод об усложнении динамики свойств.
Список литературы
1. Kuznetsov A. P., Kuznetsov S. P., Sataev I. R., Turukina L. V. About Landau-Hopf scenario in a system of coupled self-oscillators // Physics Letters A. 2013. Vol. 377, no. 45-48. P. 3291-3295. DOI: 10.1016/j.physleta.2013.10.013.
2. Osipov G. V., Pikovsky A. S., Rosenblum M. G., Kurths J. Phase synchronization effects in a lattice of nonidentical Rossler oscillators // Phys. Rev. E. 1997. Vol. 55, no. 3. P. 2353-2361. DOI: 10.1103/PhysRevE.55.2353.
3. Pikovsky A., Rosenblum M., Kurths J. Synchronization: A Universal Concept in Nonlinear Sciences. Cambridge: Cambridge University Press, 2001. 411 p. DOI: 10.1017/CB09780511755743.
4. Dodla R., Sen A., Johnston G. L. Phase-locked patterns and amplitude death in a ring of delay-coupled limit cycle oscillators // Phys. Rev. E. 2004. Vol. 69, no. 5. P. 056217. DOI: 10.1103/ PhysRevE.69.056217.
5. Williams C. R. S., Sorrentino F., Murphy T. E., Roy R. Synchronization states and multistability in a ring of periodic oscillators: Experimentally variable coupling delays // Chaos: An Interdisciplinary Journal of Nonlinear Science. 2013. Vol. 23, no. 4. P. 043117. DOI: 10.1063/1.4829626.
6. Rao R., Lin Z., Ai X., Wu / Synchronization of epidemic systems with Neumann boundary value under delayed impulse // Mathematics. 2022. Vol. 10, no. 12. P. 2064. DOI: 10.3390/math10122064.
7. Van der Sande G., Soriano M.C., Fischer I., Mirasso C. R. Dynamics, correlation scaling, and synchronization behavior in rings of delay-coupled oscillators // Phys. Rev. E. 2008. Vol. 77, no. 5. P. 055202. DOI: 10.1103/PhysRevE.77.055202.
8. Клинъшов В. В., Некоркин В. И. Синхронизация автоколебательных сетей с запаздывающими связями//Успехи физических наук. 2013. Т. 183, № 12. С. 1323-1336. DOI: 10.3367/UFNr.0183. 201312c.1323.
9. Heinrich G., Ludwig M., Qian J., Kubala B, Marquardt F. Collective dynamics in optomechanical arrays // Phys. Rev. Lett. 2011. Vol. 107, no. 4. P. 043603. DOI: 10.1103/PhysRevLett.107.043603.
10. Zhang M., Wiederhecker G. S., Manipatruni S., Barnard A., McEuen P., Lipson M. Synchronization of micromechanical oscillators using light // Phys. Rev. Lett. 2012. Vol. 109, no. 23. P. 233906. DOI: 10.1103/PhysRevLett.109.233906.
11. Lee T.E., Sadeghpour H. R. Quantum synchronization of quantum van der Pol oscillators with trapped ions // Phys. Rev. Lett. 2013. Vol. 111, no. 23. P. 234101. DOI: 10.1103/PhysRevLett.111. 234101.
12. Yanchuk S., Wolfrum M. Instabilities of stationary states in lasers with long-delay optical feedback // SIAM Journal on Applied Dynamical Systems. 2010. Vol. 9, no. 2. P. 519-535. DOI: 10.20347/WIAS.PREPRINT.962.
13. Grigorieva E. V., Haken H., Kashchenko S.A. Complexity near equilibrium in model of lasers with delayed optoelectronic feedback // In: 1998 International Symposium on Nonlinear Theory and its Applications (NOLTA'98). 14-17 September 1998, Crans-Montana, Switzerland. NOLTA Society, 1998. P. 495-498.
14. Kashchenko S. A. Quasinormal forms for chains of coupled logistic equations with delay // Mathematics. 2022. Vol. 10, no. 15. P. 2648. DOI: 10.3390/math10152648.
15. Кащенко С. А. Динамика цепочки логистических уравнений c запаздыванием и с антидиффузионной связью // Доклады Российской академии наук. Математика, информатика, процессы управления. 2022. Т. 502, № 1. С. 23-27. DOI: 10.31857/S2686954322010064.
16. Thompson J.M. T., Stewart H.B. Nonlinear Dynamics and Chaos. 2nd edition. New York: Wiley, 2002. 460 p.
17. Kashchenko S.A. Dynamics of advectively coupled Van der Pol equations chain // Chaos: An Interdisciplinary Journal of Nonlinear Science. 2021. Vol. 31, no. 3. P. 033147. DOI: 10.1063/ 5.0040689.
18. Kanter I., Zigzag M., Englert A., Geissler F., Kinzel W. Synchronization of unidirectional time delay chaotic networks and the greatest common divisor // Europhysics Letters. 2011. Vol. 93, no. 6. P. 60003. DOI: 10.1209/0295-5075/93/60003.
19. Rosin D. P., Rontani D., Gauthier D. J., Scholl E. Control of synchronization patterns in neural-like Boolean networks // Phys. Rev. Lett. 2013. Vol. 110, no. 10. P. 104102. DOI: 10.1103/PhysRevLett. 110.104102.
20. Yanchuk S., Perlikowski P., Popovych O. V., Tass P. A. Variability of spatio-temporal patterns in non-homogeneous rings of spiking neurons // Chaos: An Interdisciplinary Journal of Nonlinear Science. 2011. Vol. 21, no. 4. P. 047511. DOI: 10.1063/1.3665200.
21. Klinshov V., Nekorkin V. Synchronization in networks of pulse oscillators with time-delay coupling // Cybernetics and Physics. 2012. Vol. 1, no. 2. P. 106-112.
22. Клинъшов В. В. Коллективная динамика сетей активных элементов с импульсными связями: Обзор // Известия вузов. ПНД. 2020. Т. 28, № 5. С. 465-490. DOI: 10.18500/0869-6632-202028-5-465-490.
23. Hale J.K. Theory of Functional Differential Equations. 2nd edition. New York: Springer, 1977. 366 p. DOI: 10.1007/978-1-4612-9892-2.
24. Hartman P. Ordinary Differential Equations. New York: Wiley, 1965. 632 p.
25. Marsden J.E., McCracken M.F. The Hopf Bifurcation and Its Applications. New York: Springer, 1976. 408 p. DOI: 10.1007/978-1-4612-6374-6.
26. Кащенко С. А. О квазинормальных формах для параболических уравнений с малой диффузией // Докл. АН СССР. 1988. Т. 299, № 5. С. 1049-1052.
27. Kaschenko S. A. Normalization in the systems with small diffusion // International Journal of Bifurcation and Chaos. 1996. Vol. 6, no. 6. P. 1093-1109. DOI: 10.1142/S021812749600059X.
28. Кащенко С. А. Уравнение Гинзбурга-Ландау - нормальная форма для дифференциально-разностного уравнения второго порядка с большим запаздыванием // Журнал вычислительной математики и математической физики. 1998. Т. 38, № 3. С. 457-465.
29. Кащенко И. С., Кащенко С. А. Локальная динамика систем разностных и дифференциально-разностных уравнений // Известия вузов. ПНД. 2014. Т. 22, № 1. С. 71-92. DOI: 10.18500/08696632-2014-22-1-71-92.
30. Кащенко С. А. Бифуркации в окрестности цикла при малых возмущениях с большим запаздыванием // Журнал вычислительной математики и математической физики. 2000. Т. 40, № 5. С. 693-702.
31. Kashchenko S. A. Van der Pol equation with a large feedback delay // Mathematics. 2023. Vol. 11, no. 6. P. 1301. DOI: 10.3390/math11061301.
32. Grigorieva E. V., Kashchenko S. A. Rectangular structures in the model of an optoelectronic oscillator with delay // Physica D: Nonlinear Phenomena. 2021. Vol. 417. P. 132818. DOI: 10.1016/ j.physd.2020.132818.
33. Григоръева Е. В., Кащенко С. А. Локальная динамика модели цепочки лазеров с оптоэлек-тронной запаздывающей однонаправленной связью // Известия вузов. ПНД. 2022. Т. 30, № 2. С. 189-207. DOI: 10.18500/0869-6632-2022-30-2-189-207.
34. Кащенко С. А. Квазинормальные формы в задаче о колебаниях пешеходных мостов // Доклады Российской академии наук. Математика, информатика, процессы управления. 2022. Т. 506, № 1. С. 49-53. DOI: 10.31857/S2686954322050113.
35. Kashchenko I., Kaschenko S. Infinite process of forward and backward bifurcations in the logistic equation with two delays // Nonlinear Phenomena in Complex Systems. 2019. Vol. 22, no. 4. P. 407-412. DOI: 10.33581/1561-4085-2019-22-4-407-412.
References
1. Kuznetsov AP, Kuznetsov SP, Sataev IR, Turukina LV. About Landau-Hopf scenario in a system of coupled self-oscillators. Physics Letters A. 2013;377(45-48):3291-3295. DOI: 10.1016/j.physleta. 2013.10.013.
2. Osipov GV, Pikovsky AS, Rosenblum MG, Kurths J. Phase synchronization effects in a lattice of nonidentical Rossler oscillators. Phys. Rev. E. 1997;55(3):2353-2361. DOI: 10.1103/PhysRevE. 55.2353.
3. Pikovsky A, Rosenblum M, Kurths J. Synchronization: A Universal Concept in Nonlinear Sciences. Cambridge: Cambridge University Press; 2001. 411 p. DOI: 10.1017/CBO9780511755743.
4. Dodla R, Sen A, Johnston GL. Phase-locked patterns and amplitude death in a ring of delay-coupled limit cycle oscillators. Phys. Rev. E. 2004;69(5):056217. DOI: 10.1103/PhysRevE. 69.056217.
5. Williams CRS, Sorrentino F, Murphy TE, Roy R. Synchronization states and multistability in a ring of periodic oscillators: Experimentally variable coupling delays. Chaos: An Interdisciplinary Journal of Nonlinear Science. 2013;23(4):043117. DOI: 10.1063/1.4829626.
6. Rao R, Lin Z, Ai X, Wu J. Synchronization of epidemic systems with Neumann boundary value under delayed impulse. Mathematics. 2022;10(12):2064. DOI: 10.3390/math10122064.
7. Van der Sande G, Soriano MC, Fischer I, Mirasso CR. Dynamics, correlation scaling, and synchronization behavior in rings of delay-coupled oscillators. Phys. Rev. E. 2008;77(5):055202. DOI: 10.1103/PhysRevE.77.055202.
8. Klinshov VV, Nekorkin VI. Synchronization of delay-coupled oscillator networks. Phys. Usp. 2013;56(12):1217-1229. DOI: 10.3367/UFNe.0183.201312c.1323.
9. Heinrich G, Ludwig M, Qian J, Kubala B, Marquardt F. Collective dynamics in optomechanical arrays. Phys. Rev. Lett. 2011;107(4):043603. DOI: 10.1103/PhysRevLett.107.043603.
10. Zhang M, Wiederhecker GS, Manipatruni S, Barnard A, McEuen P, Lipson M. Synchronization of micromechanical oscillators using light. Phys. Rev. Lett. 2012;109(23):233906. DOI: 10.1103/PhysRevLett.109.233906.
11. Lee TE, Sadeghpour HR. Quantum synchronization of quantum van der Pol oscillators with trapped ions. Phys. Rev. Lett. 2013;111(23):234101. DOI: 10.1103/PhysRevLett.111.234101.
12. Yanchuk S, Wolfrum M. Instabilities of stationary states in lasers with long-delay optical feedback. SIAM Journal on Applied Dynamical Systems. 2010;9(2):519-535. DOI: 10.20347/WIAS. PREPRINT.962.
13. Grigorieva EV, Haken H, Kashchenko SA. Complexity near equilibrium in model of lasers with delayed optoelectronic feedback. In: 1998 International Symposium on Nonlinear Theory and its Applications (NOLTA'98). 14-17 September 1998, Crans-Montana, Switzerland. NOLTA Society; 1998. P. 495-498.
14. Kashchenko SA. Quasinormal forms for chains of coupled logistic equations with delay. Mathematics. 2022;10(15):2648. DOI: 10.3390/math10152648.
15. Kashchenko SA. Dynamics of a chain of logistic equations with delay and antidiffusive coupling. Doklady Mathematics. 2022;105(1):18-22. DOI: 10.1134/S1064562422010069.
16. Thompson JMT, Stewart HB. Nonlinear Dynamics and Chaos. 2nd edition. New York: Wiley; 2002. 460 p.
17. Kashchenko SA. Dynamics of advectively coupled Van der Pol equations chain. Chaos: An Interdisciplinary Journal ofNonlinear Science. 2021;31(3):033147. DOI: 10.1063/5.0040689.
18. Kanter I, Zigzag M, Englert A, Geissler F, Kinzel W. Synchronization of unidirectional time delay chaotic networks and the greatest common divisor. Europhysics Letters. 2011;93(6):60003. DOI: 10.1209/0295-5075/93/60003.
19. Rosin DP, Rontani D, Gauthier DJ, Scholl E. Control of synchronization patterns in neural-like Boolean networks. Phys. Rev. Lett. 2013;110(10):104102. DOI: 10.1103/PhysRevLett.110.104102.
20. Yanchuk S, Perlikowski P, Popovych OV, Tass PA. Variability of spatio-temporal patterns in non-homogeneous rings of spiking neurons. Chaos: An Interdisciplinary Journal of Nonlinear Science. 2011;21(4):047511. DOI: 10.1063/1.3665200.
21. Klinshov V, Nekorkin V. Synchronization in networks of pulse oscillators with time-delay coupling. Cybernetics and Physics. 2012;1(2):106-112.
22. Klinshov VV. Collective dynamics of networks of active units with pulse coupling: Review. Izvestiya VUZ. Applied Nonlinear Dynamics. 2020;28(5):465-490 (in Russian). DOI: 10.18500/ 0869-6632-2020-28-5-465-490.
23. Hale JK. Theory of Functional Differential Equations. 2nd edition. New York: Springer; 1977. 366 p. DOI: 10.1007/978-1-4612-9892-2.
24. Hartman P. Ordinary Differential Equations. New York: Wiley; 1965. 632 p.
25. Marsden JE, McCracken MF. The Hopf Bifurcation and Its Applications. New York: Springer; 1976. 408 p. DOI: 10.1007/978-1-4612-6374-6.
26. Kaschenko SA. Quasinormal forms for parabolic equations with small diffusion. Soviet Mathematics. Doklady. 1988;37(2):510-513.
27. Kaschenko SA. Normalization in the systems with small diffusion. International Journal of Bifurcation and Chaos. 1996;6(6):1093-1109. DOI: 10.1142/S021812749600059X.
28. Kashchenko SA. The Ginzburg-Landau equation as a normal form for a second-order difference-differential equation with a large delay. Computational Mathematics and Mathematical Physics. 1998;38(3):443-451.
29. Kashchenko IS, Kashchenko SA. Local dynamics of difference and difference-differential equations. Izvestiya VUZ. Applied Nonlinear Dynamics. 2014;22(1):71-92 (in Russian). 10.18500/0869-6632-2014-22-1-71-92.
30. Kashchenko SA. Bifurcations in the neighborhood of a cycle under small perturbations with a large delay. Computational Mathematics and Mathematical Physics. 2000;40(5):659-668.
31. Kashchenko SA. Van der Pol equation with a large feedback delay. Mathematics. 2023;11(6):1301. DOI: 10.3390/math11061301.
32. Grigorieva EV, Kashchenko SA. Rectangular structures in the model of an optoelectronic oscillator with delay. Physica D: Nonlinear Phenomena. 2021;417:132818. DOI: 10.1016/j.physd.2020. 132818.
33. Grigorieva EV, Kashchenko SA. Local dynamics of laser chain model with optoelectronic delayed unidirectional coupling. Izvestiya VUZ. Applied Nonlinear Dynamics. 2022;30(2):189-207. DOI: 10.18500/0869-6632-2022-30-2-189-207.
34. Kashchenko SA. Quasi-normal forms in the problem of vibrations of pedestrian bridges. Doklady Mathematics. 2022;106(2):343-347. DOI: 10.1134/S1064562422050131.
35. Kashchenko I, Kaschenko S. Infinite process of forward and backward bifurcations in the logistic equation with two delays. Nonlinear Phenomena in Complex Systems. 2019;22(4):407-412. DOI: 10.33581/1561-4085-2019-22-4-407-412.
Кащенко Сергей Александрович — родился в Ярославле (1953), окончил Ярославский государственный университет (1975). Защитил диссертацию на соискание ученой степени кандидата физико-математических наук в ННГУ (1976) и доктора физико-математических наук в МГУ (1989) в области теории нелинейных колебаний. Профессор, директор Объединенного института математики и компьютерных наук им. А. Н. Колмогорова. Автор монографий «Модели волновой памяти» (совместно с В. В. Майоровым) и «Релаксационные колебания в лазерах» (совместно с Е. В. Григорьевой). Опубликовал более 500 научных работ и 10 монографий. Награжден почетным званием «Заслуженный деятель науки Российской Федерации» (2020) и медалью «За вклад в реализацию государственной политики в области образования и научно-технологического развития» Министерства науки и высшего образования РФ (2023).
Россия, 150003 Ярославль, ул. Советская, 14
Ярославский государственный университет имени П. Г. Демидова E-mail: [email protected] ORCID: 0000-0002-8777-4302 AuthorID (eLibrary.Ru): 8238
