CC BY
10
Бифуркации в динамических системах.
Детерминированный хаос. Квантовый хаос
Известия высших учебных заведений. Прикладная нелинейная динамика. 2021. Т. 29, № 6 Izvestiya Vysshikh Uchebnykh Zavedeniy. Applied Nonlinear Dynamics. 2021;29(6)
Научная статья УДК 517.929
DOI: 10.18500/0869-6632-2021-29-6-869-891
Влияние связи на динамику трех осцилляторов с запаздыванием
А. А. Кащенко
Центр интегрируемых систем, Ярославский государственный университет имени П. Г. Демидова, Россия E-mail: a.kashchenko@uniyar.ac.ru Поступила в редакцию 15.06.2021, принята к публикации 5.07.2021, опубликована 30.11.2021
Аннотация. Цель настоящего исследования - построить асимптотику релаксационных режимов системы дифференциальных уравнений с запаздыванием, описывающей три диффузионно связанных генератора с нелинейной финитной запаздывающей обратной связью в предположении, что множитель перед функцией обратной связи является достаточно большим. Также целью является изучение влияния связи между осцилляторами на нелокальную динамику рассматриваемой модели. Методы. Мы строим асимптотику решений рассматриваемой модели с начальными условиями из специального множества. По асимптотике решений получаем оператор сдвига по траекториям, переводящий множество начальных функций в множество того же типа. Главная часть этого оператора описывается конечномерным отображением. Изучение его динамики позволяет уточнить асимптотику решений исходной модели и сделать выводы о ее динамике. Результаты. Из вида построенного отображения следует, что при положительных параметрах связи у исходной модели, начиная с некоторого момента времени, все три генератора имеют одинаковую главную часть асимптотики - генераторы «синхронизируются». При отрицательных значениях параметра связи возможны как неоднородные релаксационные циклы, так и нерегулярные режимы. Описана связь этих режимов с режимами построенного конечномерного отображения. Заключение. Из результатов работы следует, что на динамику рассматриваемой модели принципиальным образом влияет значение параметра связи между генераторами.
Ключевые слова: запаздывание, нелокальная динамика, асимптотика, релаксационные колебания.
Благодарности. Исследование выполнено при финансовой поддержке гранта Президента РФ № МК-1028.2020.1.
Для цитирования: Кащенко А. А. Влияние связи на динамику трех осцилляторов с запаздыванием // Известия вузов. ПНД. 2021. T. 29, № 6. С. 869-891. DOI: 10.18500/0869-6632-2021-29-6-869-891
Статья опубликована на условиях Creative Commons Attribution License (CC-BY 4.0).
© Кащенко А. А., 2021
869
Article
DOI: 10.18500/0869-6632-2021-29-6-869-891
Influence of coupling on the dynamics of three delayed oscillators
A. A. Kashchenko
Center of Integrable Systems, Demidov Yaroslavl State University, Russia E-mail: a.kashchenko@uniyar.ac.ru Received 15.06.2021, accepted 5.07.2021, published 30.11.2021
Abstract. The purpose of this study is to construct the asymptotics of the relaxation regimes of a system of differential equations with delay, which simulates three diffusion-coupled oscillators with nonlinear compactly supported delayed feedback under the assumption that the factor in front of the feedback function is large enough. Also, the purpose is to study the influence of the coupling between the oscillators on the nonlocal dynamics of the model. Methods. We construct the asymptotics of solutions of the considered model with initial conditions from a special set. From the asymptotics of the solutions, we obtain an operator of the translation along the trajectories that transforms the set of initial functions into a set of the same type. The main part of this operator is described by a finite-dimensional mapping. The study of its dynamics makes it possible to refine the asymptotics of the solutions of the original model and draw conclusions about its dynamics. Results. It follows from the form of the constructed mapping that for positive coupling parameters of the original model, starting from a certain moment of time, all three generators have the same main part of the asymptotics — the generators are "synchronized". At negative values of the coupling parameter, both inhomogeneous relaxation cycles and irregular regimes are possible. The connection of these modes with the modes of the constructed finite-dimensional mapping is described. Conclusion. From the results of the work it follows that the dynamics of the model under consideration is fundamentally influenced by the value of the coupling parameter between the generators.
Keywords: delay, nonlocal dynamics, asymptotics, relaxation oscillations.
Acknowledgements. Research funded by the Council on grants of the President of the Russian Federation (MK-1028.2020.1).
For citation: Kashchenko AA. Influence of coupling on the dynamics of three delayed oscillators. Izvestiya VUZ. Applied Nonlinear Dynamics. 2021;29(6):869-891. DOI: 10.18500/0869-6632-2021-29-6-869-891
This is an open access article distributed under the terms of Creative Commons Attribution License (CC-BY 4.0).
Введение
Системы дифференциальных уравнений с запаздыванием играют важную роль во многих областях науки: в медицине и физиологии [1-4], биологии [5, 6], физике [7-10]. Например, они моделируют барорефлекс [1,3,4], генные сети [2], системы передачи информации [10]. Отдельным классом систем дифференциальных уравнений с запаздыванием являются системы дифференциальных уравнений с финитной нелинейностью.
Рассмотрим дифференциальное уравнение с запаздыванием и финитной нелинейностью
и = -ш + ХР(и(Ь - Т)). (1)
Здесь функция и является скалярной, параметры V, X, Т - положительные, нелинейная функция Р является финитной. Это уравнение является математической моделью генератора с нелинейной запаздывающей обратной связью с ЯС-фильтром нижних частот первого порядка [7,8]. Генераторы, описываемые моделью (1), применяются при изготовлении гидролокаторов и усилителей класса Э [8]. Также модель (1) описывает многие биологические процессы, в которых переменная состояния и распадается со скоростью, пропорциональной и в настоящее время, и создается со скоростью, зависящей от значения и некоторое время в прошлом (см. таблицу 1 в работе [5] и литературные ссылки в ней). Кроме того, динамика моделей с финитной нелинейностью представляет общенаучный интерес [6,11,12].
В работах [13-18] аналитически изучалась динамика модели (1). Были найдены сложные периодические режимы [16,17], глобальный аттрактор [15], показано явление мультистабильности [14], и изучены свойства решений при достаточно больших значениях параметра Т [13,18].
В работах [19-22] изучалась динамика двух связанных генераторов вида (1)
( и 1 + и\ = №(щ(Ь - Т)) + - И1), [ й2 + П2 = №(и2^ - Т)) + у(и,1 - П2)
при достаточно больших значениях параметра X и разной величине связи между осцилляторами у (параметр V был взят равным единице, поскольку это значение мы можем получить в исходной модели с помощью подходящих нормировок). Было показано, что при положительных значениях связи (у > 0) порядка единицы с некоторого момента времени оба уравнения имеют одинаковую главную часть асимптотики решения, при уменьшении положительной связи происходят бифуркации, в результате которых рождаются устойчивые неоднородные релаксационные циклы, возникает мультистабильность, а при отрицательной связи порядка единицы -1/2 < у < 0 были найдены синхронные неоднородные режимы.
В настоящей работе мы рассмотрим динамическую систему, описывающую три диффузионно связанных осциллятора вида (1)
( щ + щ = ХР(щ(г - Т)) + у(щ-1 - 2щ + щ+1), (з = 1, 2, 3), ^
^ и0 = из, и4 = щ.
Здесь и^ (] = 1, 2, 3) - скалярные переменные, параметры X и Т положительные, у отлична от нуля, а знакопеременная ограниченная нелинейная кусочно-гладкая функция Р является финитной, то есть вне отрезка [-р,р] (где р - некоторая положительная константа) данная функция тождественно равна нулю. Кроме того, от функции Р(и,) на отрезке и € [-р,р] мы требуем выполнения условий
^ (р) = Р (-р) = 0,
Р(и,) = 0 за исключением конечного числа точек, (3)
если Р(и*) = 0, то существует п € N такое, что Р(п\и*) = 0.
Ключевое предположение состоит в том, что параметр X является достаточно большим
X » 1.
В работе мы будем строить асимптотику решений модели (2) при условиях (3), при достаточно больших значениях параметра X, при различных ненулевых и независящих от X значениях параметра связи между осцилляторами у для всех начальных условий из некоторого множества 5. Кроме того, мы будем изучать, как зависит динамика модели (2) от значений параметра у.
Как и в работах [20-22], изучение нелокальной динамики модели (2) при условии X ^ 1 будет проводиться с помощью специального метода большого параметра [9]. Вначале будет выбрано множество начальных условий 5. Для всех начальных условий из этого множества будет построена асимптотика решений при X ^ (см. раздел 1). Далее будет построен оператор сдвига по траекториям, переводящий множество начальных условий в множество такого же типа, и найдены главные части этого оператора при положительных и отрицательных значениях параметра у (см. разделы 2 и 3). Эти главные части будут иметь вид конечномерных отображений. Будет исследована их динамика и установлена связь динамики построенных отображений с динамикой исходной системы при положительных значениях параметра у (см. раздел 2) и при отрицательных значениях параметра у (см. раздел 4).
1. Построение асимптотики решений
Запишем систему (2) в более удобном виде:
= X
Р(«1(4 - Т))
Р(и2(1 - Т)) Р(и3(1 - Т))
+ А
где
А =
-2у - 1 у у
у -2у - 1 у у у -2у - 1
(5)
Рассмотрим линейную часть системы (4) - систему
и = Аи, (6)
где и = («1,«2, ^з)т.
Выпишем решение системы (6) для произвольных начального момента времени т и начальных условий и(т) = (и1(т),и2(т),и3(т))
т
здесь
л*
и(г) = еА(г-т)и(т),
( е— + 2е-(1+з^ е.- - е-(1+3у)* - е-(1+3у)* \
- е-(1+3у)< е— + 2е-(1+3у)* - е-(1+3у)*
^ - е-(1+3у)< - е-(1+3у)Л + 2е-(1+3у)* у
(7)
Заметим, что если все компоненты решения и системы (4) находятся вне полосы ^ (¿)| < р (] = 1, 2, 3) на протяжении времени больше или равного величине запаздывания, то на следующем отрезке времени длины запаздывания система (4) будет иметь вид (6). Из этого факта и вида решения линейной системы (формула (7)) следует, что при у < -1/3 у системы (4) существуют такие начальные условия, находящиеся вне полосы |и^ (з)| < р (] = 1, 2, 3) при в € [-Т, 0), что ни одна из компонент решения с этими начальными условиями никогда не попадет внутрь полосы 1из (-Ю! < Р. Более того, эти решения будут либо стремиться к константе (при у = -1/3), либо стремиться к бесконечности (при у < -1/3). Нас интересуют релаксационные решения, поэтому везде далее предполагаем, что выполняется неравенство у > -1/3.
Рассмотрим произвольные начальные условия системы (4) из фазового пространства С[_г,0](^3), удовлетворяющие условиям ^^(з)| ^ р, Uj(0) = Ьур, где 8 € [-Т, 0), |bj| ^ 1 для всех ] = 1, 2, 3. При у > -1/3 существует такой момент времени ¿1 ^ 0, что все компоненты решения на полуинтервале £ € [¿1 - Т, ¿1) находятся вне полосы ^ | < р, а на некотором интервале ненулевой длины, начинающемся с Ь = ¿1, хотя бы одна из компонент решения щ(1) находится внутри полосы |ui(t)| < р (это следует из того, что при у > -1/3 решение линейной системы (6) при £ ^ стремится к нулю, поэтому не могут все три его компоненты быть вне полосы |uj(¿)| < р на бесконечном промежутке времени).
Рассмотрим отрезок [¿1 - Т, ¿1] и введем на нем множество начальных условий для системы (4). Обозначим это множество через Б(х1,х2,х3). В Б(х1,х2,х3) содержатся тройки (и1(з),и2(з),и3(в)) непрерывных на отрезке в € [¿1 - Т, ¿1] функций. Эти функции по модулю не меньше чем р на полуинтервале в € [¿1 - Т,Ь1):
^(в)| ^ р при 8 € [¿1 - Т,^), 2 = 1, 2, 3,
(8)
а в точке 8 = Ь1 модуль хотя бы одной из функций равен р, а модули остальных функций больше или равны р:
Uj(Í1 ) = Xjp, \xj| > 1 (i = 1, 2, 3), 3i e {1, 2,3} : \Xi| = 1.
(9)
При этом существует такая положительная константа 5, что при t G (ti ,t1 + 5) выполняется неравенство (£)| < р (это требование означает, что хотя бы одна из компонент решения при t = t1 попадает внутрь полосы (£)| < р).
Из вида решений линейной системы (7) нетрудно найти, что в случае |х11 = 1 параметры х1, х2, х3, фигурирующие в описании множества S(х1 ,х2,х3), при у > 0 будут принадлежать множеству
Х1 = 1, |Ж2 | > 1, |^3 | > 1, Xs < 2 + 1/y - Х2, Х1 = -1, |Х2 | > 1, |хз | > 1, Хз > -2 - 1/y - Х2,
а при -1/3 < y < 0 - множеству
(10)
Х1 = 1, |х2 | > 1, |х3 | > 1,
Хз > 2 + 1/y - х2,
х1 = -1, |х2 | > 1, |хз | > 1,
хз < -2 - 1/y - х2.
(11)
Если же минимальным из модулей Х{ будет |ж21 (или |хэ |), то мы получим соответствующее множество, поменяв местами х1 и х2 (или х1 и х3, соответственно) в формулах (10) и (11). На рис. 1 изображены множества (10) и (11).
Построим асимптотику всех решений системы (4) с начальными условиями из множества £(х1, х2, х3) при X ^
Будем интегрировать систему (4) методом шагов.
Пусть £ Е , ^ + Т]. На рассматриваемом отрезке выполняются условия Р (иj (£ — Т)) = 0 (^ = 1, 2, 3), поэтому система (4) при £ Е \Ъ1 ,Ь1 + Т] принимает вид (6), и ее решение имеет вид (7), где т = ¿1, и(г1) = (^1 (¿1), и2(¿1), и3(¿1 ))т = (х1р, х2р, х3р)т.
Рис. 1. a - Часть множества (10) при х\ = 1, у = 0.5; b - часть множества (10) при х\ = — 1, у = 0.5; c - часть множества (11) при х\ = 1, у = —0.15; d - часть множества (11) при х\ = —1, у = —0.15. Легенда: горизонтальная ось - Х2, вертикальная ось - хз
Fig. 1. a - Part of the set (10) at x\ = 1, у = 0.5; b - part of the set (10) at x\ = —1, у = 0.5; c - part of the set (11) at x\ = 1, у = —0.15; d - part of the set (11) at x\ = —1, у = —0.15. Legend: horizontal axis - X2, vertical axis - хз
Хотя бы одна из компонент и3 на промежутке £ € + 8) П (¿1, ^ 1 + Т) находит-
ся в полосе (-р,р), поэтому на следующем шаге длины запаздывания нелинейность в системе (4) не будет тождественно нулевой, более того, она будет отлична от нуля на интервале Ь € € (¿1 + Т, ¿1 + 8 + Т) П (¿1 + Т, ¿1 + 2Т).
Пусть г € [¿1 + Т,г1 + 2Т]. Функции Р(и^(Ь - Т)) а = 1, 2, 3) известны с прошлого шага, поэтому на текущем промежутке времени систему (4) можно рассматривать как линейную неоднородную систему обыкновенных дифференциальных уравнений. Отсюда получаем, что решение системы (4) на отрезке Ь € [¿1 + Т, ¿1 + 2Т] имеет вид
( и1(*) \ V )
= е
/ х1р \
Х2Р
V %3'Р )
+ ХО(Ь),
(12)
где
&(г) =
( 91(1) \ 92(1)
V 9з(1) )
92{Ь,Х1,Х2,Хз)
\ дз(1,х1,х2,хз) ) ¿1+т
0А(г-з)
( Р(и1(в - Т)) \
РЫ8 - Т)) V Р(из(8 - Т)) )
йв.
Пусть выполняются следующие два предположения.
Предположение 1. Для каждого г (г = 1,2,3) количество точек ¿* € [¿1 + Т, ¿1 + 2Т], при которых §1(1*,х1,х2,х3) = 0, конечно. Для каждой точки ¿* такой, что д^*,х1,х2,х3) = 0,
д3дг (1,Х1 ,Х2,Хз)
существует ] € N такой, что
дР
г=г*
= 0.
Предположение 2. Для всех г (г = 1,2,3) верно дг(11 + 2Т, х1,х2,хз) = 0.
Тогда в Предположении 2 при достаточно больших X > 0 в точке Ь = ^+2Т верно равенство
и(*1 +2Т)= Х(С(*1 +2Т)+ о(1)). (13)
На следующем отрезке £ € [¿1 + 2Т, ¿1 + 3Т], как и ранее, будем рассматривать систему (4) как систему неоднородных обыкновенных дифференциальных уравнений. Верен следующий результат.
Лемма 1. При выполнении Предположений 1 и 2 на отрезке £ € [¿1 + 2Т, ¿1 + 3Т] асимптотика решения системы (4) при достаточно больших X > 0 будет иметь вид
и
(г) = X (еА(^1-2Т) + о(1)Щь + 2Т) =
X
1
= _(ег1+2т- - ^ЭД^т-})^! +2Т)+ 92(11 +2Т)+ дз(Ь + 2Т)) ^ 1 j +
( 91(и + 2Т) \ + Хе(1+з^)(*1+2тI д2(Ь + 2Т) I + о(Х).
V <7з(*1 + 2Т) )
(14)
Формула (14) верна потому, что при условии выполнения Предположений 1 и 2 главную часть асимптотики решения модели (4) задает решение линейной задачи (6) с начальными условиями (13), а нелинейность дает вклад, имеющий порядок о(Х) при X ^ (этот факт доказывается полностью аналогично Лемме 1 в работе [21]).
Изучим асимптотику решения модели (4) при £ > 3Т. Пока все три компоненты находятся вне полосы | ^ р, система (4) имеет вид (6) и главная часть ее решения имеет вид (14). Если же какая-то из компонент решения попадает внутрь полосы | < р в точке Ь = ¿о > ¿1+ 2Т, тогда на следующем шаге Ь Е [¿о + Т, ¿о + 2Т] нелинейность Р станет ненулевой и может изменить главную часть асимптотики решения (а может и не изменить). Изменится ли главная часть асимптотики решения, зависит от значений выражений
Ог(*) = £ д3 (¿1 + 2Т) + е-3^1-2Т) ( 3дг(¿1 + 2Т) — £ д3 (¿1 + 2Т) ], (15)
(г = 1, 2, 3) в окрестности точки ¿о.
Заметим, что при £ Е [¿1 + 2Т, ¿о] решение системы (4) в терминах функций Ог будет иметь вид
ш(г) = - (Сг(1) + о(1)) е-(^1-2Т), I = 1, 2, 3. (16)
3
Существуют две принципиально различные ситуации, когда одна (или несколько) функций щ попадает внутрь полосы 1щ1 ^ р в точке £ = ¿о > ¿1 + 2Т.
1. Второй множитель в формуле (16) при одном или нескольких г из множества {1, 2, 3} в некоторой точке из асимптотически малой при - ^ окрестности точки £ = ¿о равен нулю.
2. При всех г из множества {1, 2, 3} второй множитель в формуле (16) отличен от нуля в некоторой не зависящей от - окрестности точки £ = ¿о, а третий множитель асимптотически мал при - ^ в окрестности точки £ = ¿о.
Заметим, что, во-первых, оба случая не могут выполняться одновременно при достаточно больших значениях параметра - (второй множитель зануляется при £ не зависящих от -, а третий множитель становится асимптотически малым при асимптотически больших по - значениях ¿), а, во-вторых, существуют такие функции Р и параметры у, х\, Х2 и хз, что случай 1 не имеет места (например, это верно в случае у > 0,
Ж1 — ^2 — ^з = 1 и функции Р(и), принимающей только положительные значения на интервале и Е (—р,р)). Пусть функция Р и параметры у, х\, Х2 и хз такие, что случай 1 реализуется. Тогда верен следующий результат.
Лемма 2. Пусть хотя бы одна из компонент решения системы (4) попадает внутрь полосы 1щ1 < р в точке £ = ¿о > ¿1 + 2Т и формула (14) верна на отрезке £ Е [¿1 + 2Т, ¿о]. Если в асимптотически малой окрестности точки £ = ¿о существует точка, в которой второй множитель в формуле (16) равен нулю, тогда главный член асимптотики решения задачи (4) на отрезке £ Е [¿о + Т, ¿о + 2Т] будет иметь вид (14).
Доказательство. Для начала рассмотрим наименьшее значение ¿о > ¿1 + 2Т, при котором одна или несколько компонент решения попадают в полосу 1щ1 < р.
Во-первых, заметим, что, если второй множитель в формуле для щ (16) равен нулю, то при достаточно больших - > 0 существует такое, что Gí(t*) =0 и |£о — | = о(1), при этом
= 0(1) при - ^
Во-вторых, отметим, что каждое из уравнений Ог(Ь) = 0 (г = 1, 2, 3) имеет не более одного корня > ¿1 + 2Т (более того, три уравнения = 0, С2 (¿) = 0 и С3(£) = 0 имеют в
совокупности не более двух корней, удовлетворяющих неравенству £ > ¿1 + 2Т). В-третьих, отметим, что, если Gí(t*) = 0, то = 0.
Зафиксируем произвольное небольшое положительное значение 8 так, чтобы при £ = ¿о - 8 все три компоненты решения были вне полосы |«г| < р (оно существует в силу Предположения 2 и условия ¿о > t1 + 2Т). Тогда на отрезке £ € [¿0 - 8, ¿0 - 8 + Т] решение задается формулой (14), и, следовательно, на этом отрезке для функций С г выполняется Предположение 1 (с заменой временного отрезка [¿1 + Т,^ + 2Т] на [¿0 - 8, ¿0 - 8 + Т]). Если в точке £ = ¿0 - 8 + Т функции Сг будут отличны от нуля (то есть для них будет выполнено Предположение 2), то главную часть асимптотики решения системы (4) на отрезке [¿0-8+Т, ¿0-8+2Т] задает решение линейной системы (6) с начальными условиями в точке Ь = ¿0 - 8 + Т (доказательство этого факта полностью повторяет доказательство Леммы 1). Если же какая-то из функций С г в этой точке равняется нулю, то отступим назад на малое, но фиксированное значение 81, тогда в точке ¿0-8-81+2Т > ¿0-8+Т выполняется Предположение 2, следовательно, на отрезке [¿0-8-81+Т, ¿0-8-81+2Т] для решения системы (4) верна формула (14). Поскольку корней уравнений Ог(Ь) = 0 не более двух, то действуя так же, мы можем показать, что данная асимптотика верна при £ € [¿0-8-81+2Т, ¿0-8-81+3Т], поэтому она верна £ € [¿0,^ + 2Т]. □
Заметим, что в силу того, что суммарно корней уравнений Ог(Ь) = 0 не более двух, то возможно существование еще одного значения ¿0, для которого выполняется первый случай. Тем не менее и для него Лемма 2 верна (доказательство полностью аналогично приведенному выше).
Для дальнейших рассуждений нам понадобится ввести обозначение для времени выхода из полосы 1щ1 < р в случае 1 (если дважды наблюдался случай 1, тогда возьмем максимальное из этих времен, а если его не было, то возьмем в качестве него значение ¿1 + 2Т). Обозначим этот момент времени как Заметим, что при X ^ величина Ьь имеет оценку 0(1).
Из Леммы 2 вытекает следующее утверждение.
Следствие 1. При всех 1>1ь все компоненты щ решения системы (4) находятся вне полосы |иг| < р до тех пор, пока не произойдет случай 2.
Следствие 2. При всех £ > ¿1 + 2Т, пока не произойдет случай 2, решение системы (4) имеет вид (14).
Пусть произошел случай 2. Перепишем формулу (14) в другом виде. Главная часть асимптотики решения системы (4) на промежутке Ь € [¿1 + 2Т, ¿0] - решение линейной системы (6), где у матрицы А собственные значения -1 кратности один (соответствующий собственный вектор (1,1,1)т) и -(1 + 3у) кратности два (соответствующие собственные векторы (-1/2, -1/2,1)т и (1, -1/2, - 1/2)т). Поэтому решение системы (4) записывается в виде
( М*) \
V ^зсо)
(
= X
21
-1/2
-1/2 1
+ ¿2
1
-1/2 -1/2
\
+ 0(1)
е(*1+2Т-)(1+з7) +
+ X(zз + о(1)) I 1 I е(*1+2Т(17)
где
= К+2Т) - 92(11 + 2Т)), 22 = 2(91(Ь +2Т) - <72(*1 + 2Т)),
д1(11 + 2Т) + ^(¿1 + 2Т) + дз(11 + 2Т)).
2з =
(18)
Рассмотрим две ситуации: у > 0 и -1/3 < у < 0.
1
1
Пусть у > 0. Предположим, что выполнено условие невырожденности
23 = 0
и произошел случай 2. Тогда существуют положительные константы Мю, Ш1 и М1 такие, что в некоторой не зависящей от X окрестности точки £ = ¿0 при всех г = 1, 2, 3 выполняются неравенства
М = I Ог(1) + о(1) | | (X/3)e(tl+2Т-tо) | < М0, 0 <Ш1 < |Gг(t)+ о(1)| < Мъ
Следовательно, для положительных у выполняется неравенство < М2, где М2 -
некоторая положительная константа. Отсюда мы получаем, что при достаточно больших X
¿0 - ¿1 = (1 + о(1))ЫX, и что в окрестности точки £ = ¿0 решение системы (4) представимо в виде
(20)
( иЛ*) \ V из{Ъ) )
1
= X(zз + о(1))
е
(41+2Т—)
(21)
1
Теперь рассмотрим случай отрицательных значений параметра связи: -1/3 < у < 0. Предположим, что выполнено следующее условие невырожденности:
I
-+ 22
-- ц) • I"1 - -2>=0
(22)
(при выполнении этого условия ни одна из компонент вектора - первого слагаемого в формуле (17) не будет нулевой) и произошел случай 2.
Тогда существуют такие константы М0, тз и Мз, что в некоторой не зависящей от X окрестности точки £ = ¿0 при всех г = 1, 2, 3 выполняются условия
^ = 1 (Сг(1)+ о(1))ез^-1-2Т) 1 1 (X/3)e(1+3Y)(íl+2Т-о) 1 ^ М0, 0 < тз < 1 (Сг(1)+ о(1))ез^-1-2Т) 1 < Мз.
Таким образом, в окрестности Ь = ¿0 мы получаем ^е^3^^—1—о) | < М4, где М4 - некоторая положительная константа. Отсюда при достаточно больших X вытекает оценка на ¿0:
к - 11 = ((1 + 3у)-1 + о(1))ЫX. Следовательно, в окрестности точки £ = ¿0 решение системы (4) имеет вид
(23)
( М*) \
V из(^ )
= X
( - 2 + г2 + °(1) \
2 22 т
-----+ 0(1)
2 2 22
V 21 - ~2 + 0(1) )
е
(41+2Т—¿)(1+3у)
(24)
Из формул (20) и (23) следует, что ¿0 - > Т. Кроме того, из формул (21) и (24) видно, что, если для какого-то индекса г из множества {1, 2, 3} верно |uг(tо)| = р, то существует положительная константа 8 > 0 такая, что |ui(t)| < р при всех £ € (¿0, ¿0 + 8).
1
Таким образом, существует ¿2 (оно равно значению ¿о из случая 2), для которого выполняются следующие соотношения
(1 + о(1))1п-, у > 0,
(25)
((1 + 3у)-1 + о(1))1п-, —1/3 <у< 0, К(в + ¿2)1 >Р, (з = 1, 2, 3), 5 Е [—Г, 0), (26)
1щ(¿2)1 = %р, | ^ 1 а = 1, 2, 3), Зг : = 1. (27)
Из приведенных выше рассуждений вытекает следующий результат. Лемма 3. При £ из промежутка [¿1 + 2Т, ¿2] решение системы (4) имеет вид (14).
Сопоставляя формулы (8), (9) и (26), (27), легко видеть, что в точке ¿2 мы вернулись к исходной ситуации - множество 5(ж1,ж2,жз) преобразовалось в множество такого же типа -Б(ж^ ж2, жз). Поэтому существует оператор последования, который переводит исходное множество начальных условий Б(х1,х2,хз) в множество Б(ж1,ж2,жз). Если, начиная с точки £ = ¿2, проделать все построения, описанные выше, и при этом выполнятся условия Предположений 1 и 2 с заменой ¿1 на ¿2 и для новых значений параметров Хг, то мы получим асимптотику решения системы (4) на отрезке [¿2,£з], потом на отрезке [¿з,£4] и так далее. Таким образом с помощью оператора последования можно получить отображение на параметры х?, х?, х?, где п = 2, 3,... По динамике этого отображения мы можем судить о динамике исходной задачи, поскольку параметры х? (г = 1,2,3) определяют решение системы (4).
Перейдем к нахождению вида операторов последования и построению отображений на параметры х? (г = 1, 2, 3, п = 2, 3,...) в случаях у > 0 и —1/3 < у < 0.
2. Построение отображения и динамика исходной системы в случае положительной связи
В этом разделе мы построим отображение на параметры (ж?, х?, ж™) в случае положительной связи ( > 0) и на основе его динамики сделаем вывод о динамике исходной системы.
Введем обозначения. Пусть на п-й (п ^ 2) итерации построений, описанных в предыдущем разделе, мы получили значения 9г(Ь? + 2Т) (г = 1, 2, 3). Тогда введем значения г?, г?, г? по правилу
= \{дз(1п + 2Т, х?, х?, х?) — 92(1п + 2Т, х?, х?, ж?)),
= 3 + 2Т, х?, X?, х?) — 92(Ьп + 2Т, х?, X?, ж?)), (28)
= 1(91^? + 2Т, х?, х?, х?) + 92(1? + 2Т, х?, х?, ж?) + 9з(1? + 2Т, х?, х?, ж?)).
Формула (28) получается из (18) заменой ¿1 на Ь? и ж1, ж2, жз на ж?, ж?, ж?. Ниже мы уточним возможные значения ж?, ж?, х? при п ^ 2 в случаях у > 0 (см. раздел 2) и —1/3 < у < 0 (см. раздел 3) и в результате получим более информативные представления для величин г?, г? и г? при п ^ 2.
Пусть
= 0, (29)
t2 — 11 =
и Предположения 1 и 2 выполняются для значений х', х', х' и tn для всех натуральных п ^ 2. Тогда, действуя, как описано в предыдущем разделе, мы получим, что при положительной связи между генераторами для значений Ui(tn+i) выполняются соотношения
^(Wi) = Ц*П + o(1))e(t"+2T-t"+l), (i = 1,2,3).
Поэтому в случае у > 0 значения tn (п = 1, 2, 3,...) удовлетворяют соотношению
Wi - tn = (1 + o(1))\n X (30)
при X ^ +те.Тогда по асимптотике решения мы получаем главную часть значений параметров
™n+1. •ьг
х'+1 = kn + о(1),
= kn + о(1), (31)
х'+1 = kn + o(1),
где
kn = sign(^n) = sign (^g1(tn + 2T, xn, xn, x')+
+g2(tn + 2T, x', X', xn) + gs(tn + 2T, x', X', X')) .
Заметим, что согласно определению множества начальных условий модуль хотя бы одной из величин х'+1 в точности равен единице, а остальные модули больше или равны единице. Кроме того, из формулы (31) следует, что при всех п ^ 2 верно х' = х' + о(1), х' = х' + о(1).
Напомним, что все эти соотношения будут выполняться, только если на всех предыдущих и на текущей итерации для полученных значений параметров х' и tn будут выполняться Предположения 1 и 2. Сформулируем их для этих значений параметров. Сначала заметим, что при п ^ 2
t
gn(t,kn-1,kn-1,kn-1) = J es-tF(kn-1peT+t"-s)ds + o(1),
tn+T
и что
i i J es-tF(kn-1peT+tn-s)ds = J es-iF(kn-rpeT-s)ds, tn+T T
где £ = t — ¿п.
Поскольку в Предположении 1 время £ изменяется на промежутке \Ъп + Т,Ьп + 2Т], то при любом п переменная I изменяется на отрезке [Т, 2Т], а переменная кп-\ может принимать всего два значения: 1 и —1. Отсюда следует, что Предположение 1 для всех п = 2, 3,... выглядит одинаково (может меняться лишь кп, но оно принимает не более двух значений).
Таким образом, из выполнения следующего утверждения следует выполнение Предположения 1 при всех п = 2, 3,...
Предположение 3. Количество точек ¿* Е [Т, 2Т] таких, что Н(к, ¿*) = 0, конечно. Если к(к, ¿*) = 0, то существует ] Е N такое, что выражение дз^ к = 1 или к = —1, а функция к(к, ¿) имеет вид
отлично от нуля. Здесь
t=t*
h(к, t) = j es-i F( kpeT-S)ds. T
При выполнении Предположения 3 и условия невырожденности
h(1, 2T)h(-1,2Т) = 0 (32)
(это условие гарантирует нам выполнение Предположения 2 при п = 2, 3,...) мы получаем, что решение с начальными условиями из множества S имеет следующую асимптотику на отрезке
[tn, ¿n+1].
Пусть t G [ín, tn + Т], тогда
Uí(t)=pkn-1 é"-t + о(1), i = 1, 2, 3. (33)
При t G [ín + Т, tn + 2Т] имеем
Ui(t) = l(h(kn-1, t - ín) + o(1)), (34)
а при t G [ín + 2Т, tn+1]
Ui(t) = \(h(kn-1,2Т) + o(1))é"+2T-t, i = 1,2,3. (35)
Отметим, что формулы (33), (34), (35) получаются соответственно из формул (7), (12), (14) при подставлении в них значений хП = kn-1 + о(1) и замены í1 на tn.
Заметим, что при выполнении Предположения 3 и условия (32), начиная с п = 2,
kn = sign( z3n ) = sign(h( kn-1,2Т)).
Из построений, приведенных выше, вытекает следующее утверждение.
Теорема 1. Пусть y > 0 и для начальных значений х\, х^, х^ выполнены Предположения 1, 2 и неравенство (19). Пусть выполняются Предположение 3 и неравенство (32). Тогда для каждого достаточно большого X существует такое положительное значение t2, зависящее от (ж1,^1,^1)'1, что при t > t2 решение уравнения (4) будет иметь асимптотику (33), (34), (35).
Заметим, что из формул (33), (34), (35) следует, что главные части асимптотики всех трех компонент решения совпадают для любого значения t > t2, поэтому Теорему 1 можно трактовать следующим образом: все три осциллятора «синхронизируются» при t > ¿2.
3. Построение отображения и асимптотика решений исходной системы в случае отрицательной связи
В этом разделе мы построим отображение на параметры (х?, х?, х?) в случае отрицательной связи (—1/3 < у < 0) и уточним асимптотику решений исходной системы. Вначале приведем основной результат. Для его формулировки нам потребуется несколько обозначений.
Обозначим через 2п вектор
( — ? + % + 0(1) ^
2п =
п
—И. — £2. + 0(1) 2 2
п 2
V Щ — -2 + 0(1) )
а через Е™ - его г-ю координату (г = 1, 2, 3). Кроме того, обозначим через wn величину
wn = шт |Еп|.
¿€{1,2,3}
(36)
Сформулируем основной результат.
Теорема 2. Пусть —1/3 < у < 0 и для значений хп,хп,хп при п = 1,2,...,ц выполнены Предположения 1 и 2. Пусть для п = 1,2,... ,д выполняется неравенство
(" ? + -)(— I — !)(— §)=*■
(37)
Тогда для каждого достаточно большого X при п = 2,3 ... ,д решение уравнения (4) будет иметь следующую асимптотику: при £ е [1п, Ьп + Т]
( М*) \
V из(1) )
(
р_ Wn
п 1
п2 1 п 1
+*п-1
п 1
п2 1 п 1
п- 1
+ 0(1) + (1) + (1)
\
^ „-¿)(1+3у)
(38)
при £ е ^п + т, Ьп + 2Т]
( М*) \ V из(г) }
= X
- «)
р -М -
р
Wn р
Wn р
Wn
п 1
р -
2
п- 1 1
+ ¿Г1
п 1
2
е( 1п+т - 8)(1+3у) е( 1п +т- 5)(1 + 3у)^
г(1 п+т - 8)(1+3у)^
йв + о(Х), (39)
при Ь е [Ьп + 2Т, £п+1 ]
( М*) \
V м*))
= X
—1/2
—1/2 1
+ гп
1
—1/2 —1/2
+ о(1)
¿1 п+2Т - й)(1+3у) +
/
+ Х( гп + о(1))
1
1 1
(I „+2Т -1)
■ (40)
п
1
1
2
2
2
2
1
2
2
п
1
Здесь значения параметров г} и х^ определяются по формуле (28), а значения г?, х? и г? при п = 2,3 ... ,д принимают вид
2Т
т
\шп\ 1 2 )
? 1
—г (- (—
р_
и?
? 1
? 1
)е(1+зУ)(т 3 + 0(1)
(41)
2Т
2? =
- е
т
(1+зу)(8-2Т)
V
? 1
+ 2? -1
—(- (—
? 1
? 1
)е(1+зУ)(т+ 0(1)
(42)
2Т
/?_1 ?_1
[" (I —V К+^-'О+? (I (—V+--'У'+^-'О +
т
+(—^ 2
? 1
2
э(1+зу)(Т
^Х,-2^ + 0(1). (43)
Доказательство. Предположим, что выполняются условия невырожденности (37) (общий вид х? приведен в формуле (28)) и Предположения 1 и 2 для параметров ж?, ж?, х? и Ь? при всех п = 1, 2,..., д, где д - некоторое натуральное число. Тогда, действуя, как описано в разделе 1, мы получим, что при отрицательной связи между генераторами для значений и(Ъ?+1) (где п < д) выполняются равенства
( «1(^+1) \ ?+1) V из(г?+1) у
=-
Отсюда вытекает, что
1 — % + + о(1) Х
? ?
— Ъ — £2. + 0(1)
2 2
?
— у +о(1)
^+1 — и = ((1 + 3у)-1 + о(1)) 1п -
Э(гп+2Т-*„+1)(1+зу)
(44)
при - ^ Поскольку в точке ¿?+1 минимальный из модулей ж?+1 равен единице, то из
формулы (44) следует, что отображение на параметры ж?+1 имеет вид
? +1
? +1
= ¿(" I+* ■
и? 1
= - — ^ — +о(1) ,
? +1
и?
— ( г? — ^ +о
ш 1 1
(45)
Заметим, что «в главном» величины ж?+1 зависят только от двух параметров: х? и х? (а зависимость от х? дает вклад о(1)). Таким образом, при условии, что Предположения 1, 2 и условие (37) будут выполняться при всех п = 1, 2, 3,..., д, мы получаем представление величин ж?+1 через величины г? и х? при п ^ д. Подставив это представление в формулы (7), (12), (14), получим
1
2
2
2
2
1
2
1
2
2
2
?
з
2
1
?
2
з
формулы (38), (39), (40), соответственно, для асимптотики решения на отрезке £ € [Ьп, Ьп+1] при п = 2, 3,..., д. А формулы (41), (42), (43) получаются из формулы (28) при подстановке в нее
¿п— i
п— 1
значений хП, выраженных через г'{ "и г'2 Теорема доказана.
Введем обозначение уП ( % = 1, 2) для нормированных значений величин г':
уП
1
Wr,
□
(46)
Тогда в новых обозначениях величины x™+i примут вид
уП
x?+i = + УП + 0(1),
п+1 = - 4 - 4+о(1),
х.
2
2
х
П+ = уп - 4 + 0(1),
(47)
Поскольку хП+1 выражаются через величины у'П и уП с точностью до о(1) при X ^ то, начиная с п = 2, значения х', х', х' будут находиться на расстоянии о(1) от плоскости У1Оу2 (плоскости, натянутой на векторы (-1/2, -1/2,1)т и (1, -1/2, —1/2)т - см. формулы (40) и (46)). Более того, из формул (36) и (46) следует, что не все значения в плоскости у-\_Оу2 будут приниматься (из-за нормирующего множителя).
Выясним, какие значения могут принимать величины у'П и у'П. Из того, что вектор 2п имеет три компоненты, следует, что возможны три случая: когда минимальным среди модулей компонент вектора 2п является модуль первой, второй или третьей компоненты. Решая три системы неравенств, в результате получаем объединение 12 лучей:
У2 = -2 - У1, где У1 ^ -2 или yi ^ 0, у2 = 2 - yi, где yi ^ 0 или yi ^ 2, У2 = 2 + 2yi, где yi ^ -2 или yi ^ 0, У2 = -2 + 2yi, где yi ^ 0 или yi ^ 2, У2 = -1+ yi/2, где yi ^ -2 или yi ^ 2, У2 = 1+ yi/2, где yi ^ -2 или yi ^ 2.
(48)
Из приведенных выше рассуждений вытекает следующее утверждение.
Теорема 3. При п = 2, 3,... ,д величины х хх' могут принимать только те значения, для которых значения у' и у' находятся вблизи множества (48).
Часть множества (48) (все лучи продолжаются до бесконечности) изображена на рис. 2.
Рис. 2. Часть множества (48) Fig. 2. Part of the set (48)
п
4. Связь динамики «главной части» отображения (41), (42), (46) с динамикой исходной системы при -1/3 < у < 0
Заметим, что в силу геометрии множества (48) (оно не является выпуклым) применение принципа Шаудера [23] для доказательства теоремы о соответствии между периодическими режимами отображения (41), (42), (46) и релаксационными циклами исходной системы (4) невозможно. Однако возможно провести численный эксперимент и сформулировать гипотезы о взаимосвязи динамики отображения (41), (42), (46) и системы (4). В данном разделе мы опишем результаты численного анализа «главной части» отображения (41), (42), (46) («главной частью» мы называем отображение, которое получится после отбрасывания о-малых в формулах (41), (42)) и сопоставим их с результатами численного счета системы (4).
Остановимся на динамике «главной части» отображения (41), (42), (46). Отметим, что, если после отбрасывания о(1) в формулах (41), (42) мы выразим ,г™-1 и через у™-1 и у22-1, то получим представление для «главной части» отображения (41), (42), (46) в виде отображения на переменные у2 и у'2. Это отображение будет иметь вид (46),
2Т ^
= 2 е(1+3у)(в-2Т)
т
2 1 2 1
-р(р( - ^ - У^у(1+ШТ-^8, (49)
2Т
О Г г „.2 — 1
= 2 / ^ш*-2^
р(Р( - ^ + У2^)е(1+^)(г-в))
т
/ / ч 2 -1 ч 2 -1 \ м
- Ь_ - (50)
Везде далее мы будем рассматривать «главную часть» отображения (41), (42), (46) как отображение на переменные у™ и у2 вида (46), (49), (50).
Перечислим основные особенности динамики этого отображения.
1. Для любой функции Р, удовлетворяющей условиям, поставленным в начале статьи, и набора исходных параметров реализуется одна из возможных ситуаций:
1) у отображения есть три цикла (у 1, у2):
(2, 0) ^ (-2, 0), (0, 2) ^ (0, -2), (-2, -2) ^ (2, 2); (51)
2) у отображения есть шесть неподвижных точек (у1,у2):
(2,0), (0,2), (-2, -2) (52)
и
(-2,0), (0, -2), (2,2); (53)
3) у отображения есть три неподвижные точки (у1,у2): тройка (52) или тройка (53).
Это утверждение несложно проверить, подставив приведенные выше значения (у1, у2) в формулы (46), (49), (50).
Других циклов и неподвижных точек ни аналитический, ни численный анализ отображения не выявили.
2. Если мы возьмем значения у™ и такие, что два из трех аргументов функции Р в формулах (41), (42) на отрезке 5 Е [Т, 2Т] будут находиться вне полосы [-р,р], то на следующем шаге переменные и у^'+1 попадут (в зависимости от функции Р и значений параметров) в одно из состояний равновесия (52) или (53) или в один из циклов (51). Это следует из того, что Р = 0, когда аргумент функции Р находится вне отрезка [-р,р].
3. Стоит заметить, что при подстановке в формулу (47) первых трех (см. формулу (52)) (или вторых трех (53)) неподвижных точек получаются циклические перестановки одного и того же набора ж™'1. Поэтому, в силу симметричности исходной системы (4), мы получаем, что эти тройки неподвижных точек соответствуют одному и тому же режиму в исходной модели (с точностью до переобозначения номеров осцилляторов). Следовательно, эти режимы будут иметь одинаковые динамические свойства. То же самое верно и для циклов (51).
4. При численном анализе были найдены три принципиально разные ситуации:
1) у «главной части» отображения (41), (42), (46) все сосуществующие состояния равновесия (52), (53) или циклы (51) «устойчивы на множестве (48)» . Мы называем режим «устойчивым на множестве (48)», если все точки из малой окрестности этого режима на лучах (48) под действием итераций отображения не выходят из этой малой окрестности;
2) у «главной части» отображения (41), (42), (46) сосуществуют «устойчивые на множестве (48)» состояния равновесия (52) (или (53)) и неустойчивые (во всех просчитанных малых окрестностях этих состояний равновесия есть точки, убегающие из окрестности) состояния равновесия (53) (или соответственно (52));
3) у «главной части» отображения (41), (42), (46) все сосуществующие состояния равновесия (52), (53) или циклы (51) неустойчивы.
Был выполнен численный анализ системы (4). Ниже представлены основные результаты численного счета.
I. Как и говорит Теорема 3, значения ж™'1 (полученные на основе численного интегрирования разностного аналога системы (4)) оказались близки к плоскости у1 Оу2, а соответствующие значения у™, у% оказались близки к множеству (48) (рис. 3).
-150 ОООР
-4 .
b
У2
4 -
У 1
Рис. 3. a - Пример решения системы (4) при t Е [11000,12000]; b - итерации соответствующих значений (у™, у2) на плоскости у\Оу2 при t Е [0,12000]. Легенда: a - и\ - синий цвет, и2 - серый цвет, и3 - оранжевый пунктир. Значения параметров: \ = 10000, у = -0.15, Т =1, р =1
Fig. 3. a - Example of solution of system (4) at t Е [11000,12000]; b - iterations of the corresponding values (у™, y2 ) on the plane yiOy2 at t Е [0,12000]. Legend: a - щ - blue line, u2 - gray line, из - orange dashed line. Values of parameters: I = 10000, у = -0.15, T = 1, p = 1
a
II. Мы будем говорить, что режим системы (4) соответствует какому-то состоянию равновесия (или циклу) «главной части» отображения (41), (42), (46), если при численном счете разностного аналога системы (4) мы получаем режим, у которого значения у™ и у^ при всех п близки к этому состоянию равновесия (или циклу) «главной части» отображения (41), (42), (46).
В случае выбора параметров системы (4), соответствующих пункту 4.1), у исходной системы (4) при численном счете наблюдались релаксационные циклы, соответствующие устойчивым на множестве (48) режимам «главной части» отображения (41), (42), (46). Кроме этого, для некоторых функций в случае, когда сосуществовали три «устойчивые на множестве (48)» цикла периода два «главной части» отображения (41), (42), (46), у исходной системы наблюдались циклы, у которых параметры yf, полученные на основании численного интегрирования разностного аналога системы (4), находятся близко к множеству (48), но соответствующего близкого цикла у «главной части» отображения (41), (42), (46) нет (рис. 4).
III. В случае выбора параметров системы, соответствующих пункту 4.2), у исходной системы (4) при численном счете наблюдались релаксационные циклы, соответствующие устойчивым на множестве (48) режимам «главной части» отображения (41), (42), (46) (рис. 5, а), и для любых начальных условий (при которых производился численный счет разностного аналога системы (4)), начиная с некоторого п £ N, значения у™ и у'£, полученные численно, были близки к какому-то одному из устойчивых на множестве (48) состояний равновесия «главной части» отображения (41), (42), (46). То есть наблюдалась «сваливаемость» решений в окрестность режимов, соответствующих устойчивым на множестве (48) состояниям равновесия (41), (42), (46) (рис. 5, b). При этом термин «сваливаемость» нельзя заменить на устойчивость по Ляпунову, поскольку, если мы возьмем начальные условия для системы (4), близкие к начальным условиям, соответствующим какому-то устойчивому на множестве (48) состоянию равновесия, но не находящиеся в точности на множестве (48), то можем «свалиться» в окрестность цикла, соответствующего другому устойчивому на множестве (48) состоянию равновесия.
IV. В случае выбора параметров системы (4), соответствующих пункту 4.3), у исходной системы (4) наблюдались нерегулярные релаксационные режимы (полученные численно значения
Рис. 4. a - Релаксационный цикл системы (4) при t £ [14000,15000]; b - значения у" и у2 при п £ [400, 500]. Легенда: a - и\ - синий цвет, и2 - серый цвет, и3 - оранжевый пунктир; b - у" - синий цвет, у% - красный цвет. Значения параметров: X = 12000, у = -0.2, Т = 0.5, р =1
Fig. 4. a - Relaxation cycle of system (4) at t £ [14000,15000] b - values y" and y2 at n £ [400, 500]. Legend: a - ui - blue line, u2 - gray line, u3 - orange dashed line; b - y" - blue, y2 - red. Values of parameters: X = 12000, Y = -0.2, T = 0.5, p = 1
5000 0
-5000 -10000 -15 000 -20000 -25 ООО -30000
щ
30 ОООt 20000 10000 о
-10000 -20000 -30000 -40000 b
40
SOO
Рис. 5. a - Релаксационный цикл системы (4), соответствующий устойчивому на множестве (48) состоянию равновесия отображения (41), (42), (46). b - «Сваливаемость» к циклу системы (4), соответствующему устойчивому на множестве (48) состоянию равновесия отображения (41), (42), (46). Значения параметров: X = 10000, y = -0.15, Т = 1, р = 1
Fig. 5. a - Relaxation cycle of system (4), corresponding to the equilibrium state of the map (41), (42), (46) that is stable on the set (48). b - Transition to the cycle of the system (4), corresponding to the equilibrium state of the map (41), (42), (46) that is stable on the set (48). Values of parameters: X = 10000, y = -0.15, T = 1, p =1
a
у 1, у2 долгое время показывали нерегулярную динамику, рис. 6) или релаксационные циклы, у которых параметры у™, полученные на основании численного интегрирования разностного аналога системы (4), находятся близко к множеству (48), но соответствующего близкого цикла у «главной части» отображения (41), (42), (46) нет.
V. Если брать начальные условия для системы (4) такими, что получившиеся значения у\ и у2 попадут в точности в одну из точек (52), (53) (но сам исходный вектор (х\,х2,х3) не обязательно находился в плоскости параметров (у\, у2)), то при численном счете получим режим, соответствующий (в зависимости от значений параметров системы и функции Р) одному из циклов (51) или одной из неподвижных точек (52) или (53). Это верно даже в случае, когда параметры системы (4) таковы, что выполняется случай пункта 4.3) (то есть даже для неустойчивого
Рис. 6. a - Нерегулярный релаксационный режим (4) при t £ [0,1000]; b - значения у™ и при t £ [0,15000]. Легенда: a - щ - синий цвет, и2 - серый цвет, из - оранжевый пунктир; b - у™ - синий цвет, у% - красный цвет. Значения параметров: X = 10000, y = -0.25, Т = 0.7, р = 1
Fig. 6. a - Irregular relaxation mode (4) at t £ [0,1000] b - values у™ and y2 at t £ [0,15000]. Legend: a - щ - blue line, U2 - gray line, из - orange dashed line; b - у™ - blue, y2 - red. Values of parameters: X = 10000, y = -0.25, T = 0.7, p = 1
Рис. 7. a - Релаксационный цикл системы (4), соответствующий неустойчивому циклу (0, 2) ^ (0, -2) «главной части» отображения (41), (42), (46), при t € [14500,15000]; b - итерации значений у"на плоскости у1Оу2 при i € [0,15000]. Легенда: a - щ - синий цвет, и2 - серый цвет, и3 - оранжевый пунктир. Значения параметров: X = 10000, Y = -0.2, Т = 0.5, р =1
Fig. 7. a - Relaxation cycle of system (4), corresponding to the unstable cycle (0, 2) ^ (0, -2) of the "main" part of map (41), (42), (46); b - values (^1,^2) on the plane yiOy2 at t € [0,15000]. Legend: a - ui - blue line, u2 - gray line, u3 - orange dashed line. Values of parameters: X = 10000, y = -0.2, T = 0.5, p =1
состояния равновесия или цикла «главной части» отображения (41), (42), (46) есть множества начальных условий системы (4) такие, что решение системы с этими начальными условиями «свалится» к режиму, соответствующему рассматриваемому режиму «главной части» отображения (41), (42), (46)) (рис. 7).
Выводы
В работе исследована нелокальная динамика модели, описывающей три связанных осциллятора с нелинейной запаздывающей обратной связью. Была построена асимптотика решений с начальными условиями из множества S, и с помощью специального метода большого параметра было получено конечномерное отображение, которое помогает уточнить значения параметров, фигурирующих в формулах для асимптотики решения. Было показано, что при положительной связи между осцилляторами (у > 0), начиная с некоторого момента времени, все три осциллятора имеют одинаковую главную часть асимптотики («синхронизируются»). В случае отрицательной связи между осцилляторами (-1/3 < у < 0) было построено множество (48), которое помогает описать главную часть асимптотики решений системы (4), найдены различные релаксационные циклы исходной системы, описана их связь с циклами построенного конечномерного отображения (41), (42), (46), а также были найдены нерегулярные релаксационные режимы системы (4), показаны условия для их существования.
Список литературы
1. Ringwood J. V., Malpas S. C. Slow oscillations in blood pressure via a nonlinear feedback model // Am. J. Physiol. Regul. Integr. Comp. Physiol. 2001. Vol. 280, no. 4. P. R1105-R1115. DOI: 10.1152/ajpregu.2001.280.4.R1105.
2. Глызин С. Д., Колесов А. Ю., Розов Н.Х. Явление буферности в кольцевых генных сетях // Теоретическая и математическая физика. 2016. Т. 187, № 3. С. 560-579.
DOI: 10.4213/tmf9052.
3. Karavaev A. S., Ishbulatov Y. M., Ponomarenko V.I., Prokhorov M.D., Gridnev V.I., Bezruch-ko B. P., Kiselev A. R. Model of human cardiovascular system with a loop of autonomic regulation of the mean arterial pressure // J. Am. Soc. Hypertens. 2016. Vol. 10, no. 3. P. 235-243.
DOI: 10.1016/j.jash.2015.12.014.
4. Karavaev A. S., Ishbulatov Y.M., Prokhorov M.D., Ponomarenko V.I., Kiselev A. R., Runnova A. E., Hramkov A. N., Semyachkina-Glushkovskaya O. V., Kurths J., Penzel T. Simulating dynamics of circulation in the awake state and different stages of sleep using non-autonomous mathematical model with time delay // Front. Physiol. 2021. Vol. 11. P. 612787.
DOI: 10.3389/fphys.2020.612787.
5. an der Heiden U., Mackey M. C. The dynamics of production and destruction: Analytic insight into complex behavior// J. Math. Biol. 1982. Vol. 16, no. 1. P. 75-101. DOI: 10.1007/BF00275162.
6. Erneux T. Applied Delay Differential Equations. New York: Springer-Verlag, 2009. 204 p. DOI: 10.1007/978-0-387-74372-1.
7. Дмитриев А. С., Кислов В. Я. Стохастические колебания в радиофизике и электронике. М.: Наука, 1989. 280 с.
8. Kilias T., Kelber K., Mogel A., Schwarz W. Electronic chaos generators - design and applications // International Journal of Electronics. 1995. Vol. 79, no. 6. P. 737-753.
DOI: 10.1080/00207219508926308.
9. Кащенко С. А., Майоров В. В. Модели волновой памяти. М.: Книжный дом «Либроком», 2009. 288 c.
10. Ponomarenko V.I., Prokhorov M. D., Karavaev A. S., Kulminskiy D. D. An experimental digital communication scheme based on chaotic time-delay system // Nonlinear Dyn. 2013. Vol. 74, no. 4. P. 1013-1020. DOI: 10.1007/s11071-013-1019-0.
11. Lakshmanan M., Senthilkumar D. V. Dynamics of Nonlinear Time-Delay Systems. Berlin, Heidelberg: Springer-Verlag, 2011. 313 p. DOI: 10.1007/978-3-642-14938-2.
12. Преображенская М. М. Дискретные бегущие волны в релейной системе уравнений типа Мэки-Гласса с двумя запаздываниями // Теоретическая и математическая физика. 2021. T. 207, № 3. С. 489-504. DOI: 10.4213/tmf10038.
13. Mallet-Paret J., Nussbaum R.D. Global continuation and asymptotic behaviour for periodic solutions of a differential-delay equation // Annali di Matematica Pura ed Applicata. 1986. Vol. 145, no. 1. P. 33-128. DOI: 10.1007/BF01790539.
14. Losson J., Mackey M. C., Longtin A. Solution multistability in first-order nonlinear differential delay equations // Chaos. 1993. Vol. 3, no. 2. P. 167-176. DOI: 10.1063/1.165982.
15. Krisztin T., Walther H.-O. Unique periodic orbits for delayed positive feedback and the global attractor// J. Dyn. Diff. Equat. 2001. Vol. 13, no. 1. P. 1-57. DOI: 10.1023/A:1009091930589.
16. Stoffer D. Delay equations with rapidly oscillating stable periodic solutions // J. Dyn. Diff. Equat. 2008. Vol. 20, no. 1. P. 201-238. DOI: 10.1007/s10884-006-9068-4.
17. Krisztin T., Vas G. Large-amplitude periodic solutions for differential equations with delayed monotone positive feedback // J. Dyn. Diff. Equat. 2011. Vol. 23, no. 4. P. 727-790.
DOI: 10.1007/s10884-011-9225-2.
18. Kashchenko I., Kaschenko S. Normal and quasinormal forms for systems of difference and differential-difference equations // Communications in Nonlinear Science and Numerical Simulation. 2016. Vol. 38. P. 243-256. DOI: 10.1016/j.cnsns.2016.02.041.
19. Kashchenko A.A. Non-rough relaxation solutions of a system with delay and sign-changing nonlinearity // Nonlinear Phenomena in Complex Systems. 2019. Vol. 22, no. 2. P. 190-195.
20. Кащенко А. А. Релаксационные циклы в модели двух слабо связанных осцилляторов со знакопеременной запаздывающей обратной связью // Теоретическая и математическая физика. 2020. T. 202, № 3. С. 437-446. DOI: 10.4213/tmf9806.
21. Kashchenko A.A. Dependence of dynamics of a system of two coupled generators with delayed feedback on the sign of coupling // Mathematics. 2020. Vol. 8, no. 10. P. 1790.
DOI: 10.3390/math8101790.
22. Kashchenko A. A. Relaxation modes of a system of diffusion coupled oscillators with delay // Communications in Nonlinear Science and Numerical Simulation. 2021. Vol. 93. P. 105488. DOI: 10.1016/j.cnsns.2020.105488.
23. Schauder J. Der Fixpunktsatz in Funktionalraumen // Studia Mathematica. 1930. Bd. 2, Nr. 1. S. 171-180. DOI: 10.4064/sm-2-1-171-180.
References
1. Ringwood JV, Malpas SC. Slow oscillations in blood pressure via a nonlinear feedback model. Am. J. Physiol. Regul. Integr. Comp. Physiol. 2001;280(4):R1105-R1115.
DOI: 10.1152/ajpregu.2001.280.4.R1105.
2. Glyzin SD, Kolesov AY, Rozov NK. Buffering in cyclic gene networks. Theor. Math. Phys. 2016;187(3):935-951. DOI: 10.1134/S0040577916060106.
3. Karavaev AS, Ishbulatov YM, Ponomarenko VI, Prokhorov MD, Gridnev VI, Bezruchko BP, Kiselev AR. Model of human cardiovascular system with a loop of autonomic regulation of the mean arterial pressure. J. Am. Soc. Hypertens. 2016;10(3):235-243.
DOI: 10.1016/j.jash.2015.12.014.
4. Karavaev AS, Ishbulatov YM, Prokhorov MD, Ponomarenko VI, Kiselev AR, Runnova AE, Hramkov AN, Semyachkina-Glushkovskaya OV, Kurths J, Penzel T. Simulating dynamics of circulation in the awake state and different stages of sleep using non-autonomous mathematical model with time delay. Front. Physiol. 2021;11:612787. DOI: 10.3389/fphys.2020.612787.
5. an der Heiden U, Mackey MC. The dynamics of production and destruction: Analytic insight into complex behavior. J. Math. Biol. 1982;16(1):75-101. DOI: 10.1007/BF00275162.
6. Erneux T. Applied Delay Differential Equations. New York: Springer-Verlag; 2009. 204 p. DOI: 10.1007/978-0-387-74372-1.
7. Dmitriev AS, Kislov VY. Stochastic Oscillations in Radiophysics and Electronics. Moscow: Nauka; 1989. 280 p. (in Russian).
8. Kilias T, Kelber K, Mogel A, Schwarz W. Electronic chaos generators - design and applications. International Journal of Electronics. 1995;79(6):737-753.
DOI: 10.1080/00207219508926308.
9. Kashchenko SA, Maiorov VV. Models of Wave Memory. Moscow: Book House «Librokom»; 2009. 288 p. (in Russian).
10. Ponomarenko VI, Prokhorov MD, Karavaev AS, Kulminskiy DD. An experimental digital communication scheme based on chaotic time-delay system. Nonlinear Dyn. 2013;74(4):1013-1020. DOI: 10.1007/s11071-013-1019-0.
11. Lakshmanan M, Senthilkumar DV. Dynamics of Nonlinear Time-Delay Systems. Berlin, Heidelberg: Springer-Verlag; 2011. 313 p. DOI: 10.1007/978-3-642-14938-2.
12. Preobrazhenskaya MM. Discrete traveling waves in a relay system of Mackey-Glass equations with two delays. Theor. Math. Phys. 2021;207(3):827-840. DOI: 10.1134/S0040577921060106.
13. Mallet-Paret J, Nussbaum RD. Global continuation and asymptotic behaviour for periodic solutions of a differential-delay equation. Annali di Matematica Pura ed Applicata. 1986;145(1): 33-128. DOI: 10.1007/BF01790539.
14. Losson J, Mackey MC, Longtin A. Solution multistability in first-order nonlinear differential delay equations. Chaos. 1993;3(2):167-176. DOI: 10.1063/1.165982.
15. Krisztin T, Walther HO. Unique periodic orbits for delayed positive feedback and the global attractor. J. Dyn. Diff. Equat. 2001;13(1):1-57. DOI: 10.1023/A:1009091930589.
16. Stoffer D. Delay equations with rapidly oscillating stable periodic solutions. J. Dyn. Diff. Equat. 2008;20(1):201-238. DOI: 10.1007/s10884-006-9068-4.
17. Krisztin T, Vas G. Large-amplitude periodic solutions for differential equations with delayed monotone positive feedback. J. Dyn. Diff. Equat. 2011;23(4):727-790.
DOI: 10.1007/s10884-011-9225-2.
18. Kashchenko I, Kaschenko S. Normal and quasinormal forms for systems of difference and differential-difference equations. Communications in Nonlinear Science and Numerical Simulation. 2016;38:243-256. DOI: 10.1016/j.cnsns.2016.02.041.
19. Kashchenko AA. Non-rough relaxation solutions of a system with delay and sign-changing nonlinearity. Nonlinear Phenomena in Complex Systems. 2019;22(2):190-195.
20. Kashchenko AA. Relaxation cycles in a model of two weakly coupled oscillators with sign-changing delayed feedback. Theor. Math. Phys. 2020;202(3):381-389.
DOI: 10.1134/S0040577920030101.
21. Kashchenko AA. Dependence of dynamics of a system of two coupled generators with delayed feedback on the sign of coupling. Mathematics. 2020;8(10):1790. DOI: 10.3390/math8101790.
22. Kashchenko AA. Relaxation modes of a system of diffusion coupled oscillators with delay. Communications in Nonlinear Science and Numerical Simulation. 2021;93:105488.
DOI: 10.1016/j.cnsns.2020.105488.
23. Schauder J. Der Fixpunktsatz in Funktionalraiimen. Studia Mathematica. 1930;2(1):171-180 (in German). DOI: 10.4064/sm-2-1-171-180.
Кащенко Александра Андреевна - родилась в Ярославле (1990). Окончила с отличием математический факультет Ярославского государственного университета имени П. Г. Демидова (2012). Защитила диссертацию на соискание учёной степени кандидата физико-математических наук по специальности «Дифференциальные уравнения, динамические системы и оптимальное управление» (2015, ЯрГУ). Работает на кафедре математического анализа Ярославского государственного университета в должности доцента. Научные интересы - динамика моделей, описываемых дифференциальными уравнениями с запаздыванием.
Россия, 150003 Ярославль, Советская, 14
Ярославский государственный университет имени П. Г. Демидова E-mail: a.kashchenko@uniyar.ac.ru ORCID: 0000-0003-3823-9351 AuthorlD: 815088
