Periodic solutions of a differential equation with relay nonlinearity with delay
- Authors: Bain D.D.1
-
Affiliations:
- Ярославский государственный университет им. П. Г. Демидова
- Issue: Vol 231 (2024)
- Pages: 3-12
- Section: Статьи
- URL: https://bakhtiniada.ru/2782-4438/article/view/262307
- DOI: https://doi.org/10.36535/2782-4438-2024-231-3-12
- ID: 262307
Cite item
Full Text
Abstract
For one class of second-order differential equations with relay nonlinearity and delay, orbitally stable periodic solutions are found by means of the recurrence operator, which is a suspension over some one-dimensional mapping. The analysis of this one-dimensional mapping shows that there exist domains of parameters for which exponentially orbitally stable periodic solutions exist.
Full Text
1. Постановка задачи. Рассмотрим уравнение
(1)
в котором T, b, c, d — постоянные параметры. При этом , .
Похожие системы встречаются в работах [1–3], и само уравнение (1) при различных условиях на параметры встречается в [4] в качестве предельного объекта для нелинейного осциллятора с большим параметром и в [5] в качестве модели движения судна управляемого двухпозиционным авторулевым. Мы же будем рассматривать уравнение (1) при условиях , .
Определение 1.1. Функцию мы будем называть осциллирующей, если обращается в нуль для сколь угодно больших t. Если на всяком интервале длины запаздывания осциллирующая функция обращается в нуль не более одного раза, то такая функция называется медленно осциллирующей.
Далее будем отыскивать именно медленно осциллирующие периодические решения уравнения (1), потому что, с одной стороны, при численном счёте быстро осциллирующих решений не обнаруживается, а с другой стороны, рассмотрение только медленно осциллирующих решений позволяет ввести оператор последования, при анализе которого получаются условия существования и орбитальной устойчивости периодических решений.
Обратимая замена
приводит уравнение (1) к виду
(2)
где , , и . Таким образом, уравнение (1) сводится к случаю , . В этой форме мы и будем исследовать данную систему.
2. Определение начальных условий.
Как и в [4, 5], с целью построения решений уравнения (2) определим множество начальных условий следующим образом:
(3)
Для произвольной пары обозначим через решение уравнения (2) с начальными условиями
(4)
Заметим, что при в силу (3) и (4) выполняется неравенство . Поэтому для решение совпадает с решением задачи Коши
(5)
которое имеет вид
(6)
где
(7)
Так как , то числа , вещественны и различны, при этом значения параметров b, c однозначно соответствуют значениям параметров , поэтому далее мы будем рассматривать именно , и s в качестве параметров уравнения (2), имея в виду , .
Заметим, что если в качестве начальных условий брать фиксированный v и различные , т.ч. , то соответствующие им решения будут совпадать. Поэтому, решение определяется параметром v однозначно.
Лемма 2.1. Решение не обращается в нуль при тогда и только тогда, когда
(8)
При этом знак совпадает со знаком .
Доказательство. При решение совпадает с (при они совпадают по непрерывности). Заметим, что уравнение может иметь не более одного решения, и, следовательно, функция имеет не более одного экстремума, а значит, по теореме Ролля она имеет не более двух нулей. Ввиду начальных условий (5) функция имеет нуль в точке и положительна при малых , поэтому неравенство эквивалентно тому, что решение не обращается в нуль при . Проделав преобразования неравенства получаем неравенство (8). Тот факт, что знак совпадает со знаком , также может быть доказан прямыми преобразованиями.
Пусть — наименьший положительный нуль решения , если он существует, и в противном случае. Так как при (в том числе если ), то при решение совпадает с решением задачи Коши
(9)
которое имеет вид
(10)
где
(11)
Определим множество V, состоящее из таких v, что имеет наименьший положительный нуль на интервале , формулой
(12)
Заметим, что условие является необходимым для медленной осцилляции решения , поэтому мы далее будем искать решения именно в предположении . Следующая лемма позволит нам прояснить вид множества V.
Лемма 2.2. Если , то .
Доказательство. В выражениях (6) и (10) можно, вынеся за скобку получить выражения
Заметим, что выражение в скобках положительно при положительных t. Следовательно, при любом фиксированном значения и строго возрастают по v.
Пусть . Если у решения нет положительных нулей, то при , но тогда, так как , верны неравенства
т.е. и у решения нет положительных нулей, что неверно, ведь . Значит, имеет положительный корень . При этом , так как из неравенства согласно лемме 2.1 следует, что при . Таким образом, .
Следствие. Множество V является связным множеством в , а значит, если оно не пусто, то является некоторым интервалом.
Лемма 2.3. Если , то является кратным нулём в том и только том случае, когда .
Доказательство. Пусть . Тогда и для любых верно неравенство при . Следовательно, — это минимум решения , т.е. нуль чётной кратности.
Обратно, пусть и является кратным нулём . Уравнение имеет не более одного решения , причём из условия (т.е. ) следует, что . Поэтому точка является единственным локальным минимумом , а значит, она также является и глобальным минимумом функции . Следовательно, ввиду монотонности по v, доказанной в лемме 2.2, при любых верно неравенство при , а значит, , т.е. v является верхней гранью V. Поскольку , то v является наименьшей верхней гранью V.
Для того чтобы иметь дело только с простыми нулями, будем рассматривать множество
(13)
которое, если не пусто, имеет вид
для некоторого . Ввиду леммы 2.3 для любых нуль является простым.
Вместе с будем рассматривать соответствующее множество начальных условий , заданное соотношением
(14)
3. Построение оператора последования.
Если , то известно решение при . Для построения при воспользуемся симметрией уравнения (2), которое не меняется при замене . С этой целью введём в рассмотрение одномерное отображение
(15)
Отметим, что , так как при условии нуль решения является простым.
Далее, выполним в уравнении (2) замену и сдвиг по времени . В результате придём к тому же уравнению (2), но с начальными условиями
(16)
Заметим, что и при , так что из проделанных построений вытекает включение
где S — множество (3). Значит, на последующих отрезках решение будет строиться аналогично решению с начальным условием , т.е. на множестве определён оператор последования
(17)
Заметим, что различные функции при фиксированном v переходят в одну и ту же функцию , не зависящую от . Таким образом, оператор (17) представляет собой надстройку над одномерным отображением (15), поэтому вопрос о неподвижных точках отображения П сводится к аналогичному вопросу для отображения p, а именно, всякой экспоненциально устойчивой неподвижной или периодической точке отображения p отвечает экспоненциально орбитально устойчивый цикл .
Отметим, что не всегда верно, что , или, эквивалентно, что . Условия, при которых это вложение верно, будут получены позже, но мы не будем пользоваться принципом Шаудера, так как и не замкнуты.
4. Анализ одномерного отображения.
Рассмотрим как функцию двух переменных . Ввиду формулы (10) эта функция непрерывно дифференцируема по совокупности переменных на открытом множестве . Для любого функция непрерывно дифференцируема в некоторой окрестности точки , при этом , и, так как является простым нулём, то . Тогда, по теореме о неявной функции, функция непрерывно дифференцируема и её производная находится следующим образом:
Далее, отображение можно записать в виде , а значит, оно непрерывно дифференцируемо, и его производная равна
что приводится к виду
(18)
Так как и , то знак совпадает со знаком . В свою очередь, знак совпадает со знаком .
Отображение возрастает при и убывает при , так что для случая отображение строго возрастает на всём множестве .
Определение 4.1. Областью сходимости неподвижной точки отображения будем называть множество таких , что последовательность сходится к точке .
Теорема 4.1. Пусть , определены равенствами (7) и выполнено неравенство . Тогда отображение не имеет периодических точек и, в зависимости от значений , , реализуются случаи (a)–(d), перечисленные ниже и изображённые на рис. 1.
Рис. 1. Области пространства параметров в случае .
- Если , то , т.е. отображение не определено.
- Если и , то отображение не имеет неподвижных точек.
- Если и , то отображение имеет единственную неподвижную точку, которая экспоненциально устойчива и имеет область сходимости .
- Если , то отображение имеет единственную неподвижную точку, которая экспоненциально устойчива и имеет область сходимости .
Отметим, что из (7) следует, что , так что случаи в теореме выписаны именно в предположении , однако в уравнении (2) параметры и симметричны относительно замены и друг на друга, поэтому на рис. 1 изображены аналогичные случаи для .
Доказательство теоремы. Напомним, что если , то , и . Так как , то, согласно лемме (2.1), для любых верно неравенство при . Так как , то, ввиду формулы (18), отображение является возрастающим для любых . При этом, так как отображение монотонно возрастает, то оно не имеет периодических точек. Действительно, пусть ; тогда и т. д., и для знака неравенства в обратную сторону аналогично.
Из равенств и можно получить соотношение
(19)
из которого, в частности, следует, что , так как в случае можно показать, что , что невозможно, поскольку .
(a) Пусть ; тогда можно показать, что для любых верно неравенство при . Поскольку , получаем, что при , а значит, . Следовательно, .
(b1). Пусть и . Пусть — наименьшая неподвижная точка отображения , либо, что то же самое, наименьший нуль функции . Согласно равенству (18) и ввиду неравенств и имеем
и, следовательно,
(20)
С одной стороны, при , так как и — это наименьший нуль . С другой стороны, ввиду неравенства (20), в некоторой левой окрестности точки функция принимает отрицательные значения, что является противоречием. Следовательно, не имеет неподвижных точек.
(b2), (c). Пусть теперь и ; тогда и . Пусть — неподвижная точка отображения . Выразив из равенства (19) и используя неравенство , можно показать, что . Следовательно, всякая неподвижная точка отображения является экспоненциально устойчивой, а значит, если таковая существует, то она единственна.
Пусть , определены равенством (11); тогда кусочное решение имеет вид
и монотонно стремится к значению 1, а значит нуль не пересекает; следовательно, .
Пусть ; тогда решение имеет вид (10), причём коэффициент при экспоненте с положительным показателем отрицателен, а значит, при достаточно больших кусочное решение пересекает нуль. Следовательно, для любых верно, что , а значит .
Если , то по непрерывности , но тогда ; противоречие. Следовательно, и
Далее, неподвижные точки отображения однозначно соответствуют нулям отображения . Заметим, что и .
Пусть . Тогда и, по теореме о промежуточном значении, имеет неподвижную точку, причём она единственна и экспоненциально устойчива.
Пусть теперь . Тогда , и если — наибольший корень функции , то с одной стороны, при , а с другой стороны, ; значит, в некоторой правой окрестности верно неравенство , что является противоречием. Значит, при неподвижных точек нет.
Наконец, условие может быть записано в виде
при этом, как можно проверить, из него следует условие , так что утверждение (c) доказано.
С другой стороны, условие может быть записано в виде
и можно проверить, что оно влечёт условие при . Таким образом, утверждение (b) доказано.
(d) Пусть , тогда и . Заметим,что при любых экспоненты в формуле (10) стремятся к нулю при , а значит, стремится к при . Следовательно, при некотором t кусочное решение обращается в нуль, и определено. Поскольку это верно для любых , то .
Отображение возрастает, и, как видно из равенства (19), оно ограничено сверху положительным числом . Далее, неподвижные точки отображения однозначно соответствуют нулям отображения . Заметим, что , и при , так что по теореме о промежуточном значении имеет неподвижную точку.
Так как , то, аналогично случаю (c), для любой неподвижной точки верно неравенство , а значит, эта неподвижная точка единственна, экспоненциально устойчива, причём её областью сходимости является всё множество .
На рис. 2 изображены графики отображения при различных значениях параметров, удовлетворяющих неравенству . Эти графики иллюстрируют наличие экспоненциально устойчивых неподвижных точек в случаях (c) и (d), а также характер появления неподвижной точки при переходе от случая (b) к случаю (c).
Рис. 2. Графики отображения в случае . Неподвижным точкам отображения соответствуют пересечения графиков с главной диагональю.
На рис. 3, 4 изображены графики и траектории в фазовой плоскости некоторых решений уравнения (2) с начальными условиями из . Как видно, решения стремятся к некоторым предельным циклам, причём пересечения решениями оси ординат на изображениях фазового портрета соответствуют итерациям отображения . Точки излома на фазовых портретах соответствуют переключениям запаздывающего члена: в этих точках вторая производная решения имеет скачок.
Рис. 3. Графики решений в случае (c).
Рис. 4. Графики решений в случае (d).
Теорема 4.2. Если выполнено неравенство , то отображение не имеет неподвижных точек. Точнее, в зависимости от значений , , реализуются случаи (e)–(g), перечисленные ниже и изображённые на рис. 6. [(a)]
(e) Если , то у отображения нет неподвижных точек. Верно неравенство и монотонно убывает к нулю.
(f) Если и , то у отображения есть лишь неустойчивая неподвижная точка. Верно неравенство и монотонно убывает к нулю.
(g) Если и , то у отображения нет неподвижных точек. При отображение возрастает и стремится к бесконечности при .
Данный случай отличается отсутствием устойчивых неподвижных и периодических точек. Поскольку интереса неустойчивые решения не представляют, доказательство этой теоремы мы не приводим.
Графики отображения при различных значениях параметров в случаях (b) (пунктир) и (c) (сплошная линия).
Графики отображения при различных значениях параметров в случае (d).
Графики отображения при различных значениях параметров в случаях (e) (пунктир) и (f) (сплошная линия).
Графики отображения при различных значениях параметров в случаях (f) (сплошная линия) и (g) (пунктир).
На рис. 5 изображены графики отображения при различных значениях параметров, удовлетворяющих неравенству . Эти графики иллюстрируют наличие неустойчивой неподвижной точки в случае (f), а также характер появления неподвижной точки при переходе от случая (e) к случаю (f) и при переходе от случая (g) к случаю (f).
Рис. 5. Графики отображения в случае . Неподвижным точкам отображения соответствуют пересечения графиков с главной диагональю.
Рис. 6. Области пространства параметров в случае .
Численные эксперименты, в том числе с быстро осциллирующими начальными условиями, указывают на то, что описанные в этой работе случаи исчерпывающе описывают динамику уравнения (1). Точнее, наблюдается, что при любых начальных функциях решение либо стремится к бесконечности, либо стремится к неподвижной точке, либо происходит наматывание на цикл. При этом периодических решений, кроме тех, что мы рассмотрели не наблюдается.
На обосновании этих утверждений мы не останавливаемся, так как численному исследованию уравнения (1) при произвольных параметрах будет посвящена отдельная статья.
About the authors
Danila D. Bain
Ярославский государственный университет им. П. Г. Демидова
Author for correspondence.
Email: danila.bain@yandex.com
Russian Federation, Ярославль
References
- Кащенко С. А. Сравнительный асимптотический анализ динамики автогенераторов с различными нелинейными запаздывающими связями// Фундам. прикл. мат. — 1999. — 5, № 4. — С. 1027–1060
- Дмитриев А. С., Кащенко С. А. Асимптотика нерегулярных колебаний в модели автогенератора с запаздывающей обратной связью// Докл. РАН. — 1993. — 328, № 2. — С. 174–177
- Кащенко С. А. Асимптотика релаксационных колебаний в системах дифференциально-разностных уравнений с финитной нелинейностью. I// Диффер. уравн. — 1995. — 31, № 8. — С. 1330–1339
- Глызин С. Д., Колесов А. Ю., Розов Н. Х. Явление буферности в кольцевых цепочках однонаправленно связанных генераторов// Изв. РАН. Сер. мат. — 2014. — 78, № 4. — С. 73–108
- Колесов А. Ю., Мищенко Е. Ф., Розов Н. Х. Реле с запаздыванием и его C1-аппроксимация// Тр. Мат. ин-та им. В. А. Стеклова РАН. — 1997. — 216. — С. 126–153.
Supplementary files
