Признаки существования предельных циклов дифференциальных уравнений второго порядка Текст научной статьи по специальности «Математика»

ISSN 2074-1871 Уфимский математический журнал. Том 9. № 4 (2017). С. 3-11.

УДК 517.91

ПРИЗНАКИ СУЩЕСТВОВАНИЯ ПРЕДЕЛЬНЫХ ЦИКЛОВ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ ВТОРОГО ПОРЯДКА

М.К. АРАБОВ, Э.М. МУХАМАДИЕВ, И.Д. НУРОВ, Х.И. СОБИРОВ

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

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

Mathematics Subject Classification: 34Е

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

х" + ах' + Ьх + с\х! — \\ = 0, (1)

при определенных значениях коэффициентов а, Ь, с и параметра А, возникают предельные циклы. А именно, в работе [7] с помощью компьютерного моделирования было установлено, что уравнение (1) при а = 1,Ь = 1, с = 3/2 и А> 0 имеет предельный цикл. Более общее условие существования предельного цикла уравнения (1) получено в работе [8].

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

В настоящей работе изучается нелинейное дифференциальное уравнение второго порядка

х" + ах' + bx + с\х' — ip(x,x' )\ = 0, (2)

где а, Ь, с - вещественные числа, а функция <р(х, у) - непрерывна и удовлетворяет некоторому условию роста при |ж| + \у\ ^ ж, Найдены условия на коэфф ициенты а, Ь, си функцию <р(х,у), которые обеспечивают существование предельного цикла у уравнения (2).

M.K. Arabov, Е. Mukhamadiev, I.D. Nurov, Кн.I. Sobirov, Existence tests for limiting cycles

of second order differential equations.

©Арабов M.K., Мухамадиев Э.М., Нуров И.Д., Собиров Х.И. 2017.

Поступила 19 апреля 2016 г.

Анализ фазового портрета однородного уравнения

Сначала приведем анализ фазового портрета траекторий однородного уравнения

х" + ах' + Ьх + с|ж'| = 0, (3)

в зависимости от расположения коэффициентов (а,Ь,с) как точки трехмерного пространства Д3 (см. [8]). Еели с = 0, то (3) есть линейное уравнение второго порядка, фазовый портрет которого хорошо известен (см., папр,, |1,2, 4|), Ниже мы будем предполагать, что с > 0, Заметим, что если с < 0, то замен ой х на —х уравнение (3) сводится к такому же уравнению с коэффициентами (а,Ь, —с).

Уравнение (3) "склеивается" из линейных уравнений

х" + (а + с)х' + Ьх = 0, если х' > 0 (4)

и

х" + (а — с)х' + Ьх = 0, если х' ^ 0.

Обозначим через и корни характеристического уравнения

+ (а ± с)у + Ь = 0, соответствующего уравнениям (4) и (5):

(5)

(6)

= — 1 {(а ± с) — у/(а ± с)2 — 4&} , ^ = — 2 {(а ± с) + у/(а ± с)2 — 4Ь}

Из этих форму:: дня корней характеристических уравнений (6) следует, что в полупространстве {(а,Ь,с) : с> 0} коэффициентов уравнения (3) корни обращаются в нуль при Ь = 0 и меняют знак при возрастании Ь; корни соответственно, вещественны

при 4Ь < (а ± с)2, становятся кратными при 4Ь = (а ± с)2 и комплексно сопряженными при 4Ь > (а ± с)2.

Рис. 1. Секторы

Исходя из этих свойств корней характеристических уравнений полупространство {(а, Ь,с) : с > 0} разложим на следующие подмножества:

1. {(а, Ь,с) : Ь < 0};

2. {(а, Ь,с): 0 < 4Ъ ^ (|а| — с)2, |а| > с};

3. {(а, Ь,с): 0 < 4Ь ^ (|а| - с)2, |а| < с};

4. {(а, Ь, с) : (Н - с)2 < 4Ъ ^ (|а| + с)2};

5. {(а,Ь,с) : (|а| + с)2 < 4Ь}.

Проекция пересечений этих подмножеств с плоскостью с = сопЬ > 0 на координатную плоскость (а,Ь) приведена на рис, 1.

Отметим, что для коэффициентов (а,Ь,с), принадлежащих множествам 1-4, ненулевое решение х({) уравнения (3) с начальным значением х(10) = х0,х'(¿0) = 0 является решением линейного уравнения (4) при (£ — 10)х0 < 0 и уравнения (5) при (£ — 10)х0 > 0,

Уравнение (3) эквивалентно системе

Г ж' = У,

\у' = —аУ — Ьх — с

Если Ь = 0, то система (7) имеет единственную особую точку (0, 0), Как и в случае линейного уравнения расположение корней характеристических уравнений (6) па комплексной плоскости однозначно определяет поведение траектории кусочно-линейной системы (7) в фазовой плоскости (х,у). Как видно из ниже приводимой таблицы, у кусочно-линейных систем возникают изолированные сложные особые точки, где появляется эллиптический сектор, состоящий из траекторий, которые приближаются к особой точке как при £ ^ то так и при £ ^ —то.

7)

№

Область в пространстве параметров

Корни характери-сти ческхих) уравнения

Числа веще-

ственны и разжнх) знака

Качественная картина фазовых траекторий

Ь < 0

Седло

Числа ве-

щсствснны и 0дн01х)

знака

0 < 4 Ь < {И — с}2, |а| > с.

<

а > 0 а < 0

Узел

1

2

0 < 4 Ь

< {И — с}2,

|а| < с.

<

Если > 0, то

> 0 ^Г < 0

и < 0, (что соответствует условию с < 0). Если < 0, то < 0 ^Г > 0 и ^Г > 0, (что соответствует условию с > 0)

' \ \ n

\хч \ \ \

\ \

\ \ \

Ж

с > 0 с < 0

Узловой сектор и эллиптический сектор

{|а|— с}2 < 46 ^

< {И + с}2

Либо веще-

ственны и одного знака, и = а ± г^— комплексно со-

пряженные, либо ^+2 — вещественны и однох'о знака, и = а ± комплексно сопряженные

Ш

—

а < 0, с < 0

*2

Ш V/

а < 0, с > 0 а > 0, с < 0

Узел-фокус

46 > {|а| + с}2

Числа комплексные женные

сонря-

/ I Гг х2 X \

33 у Х1

а > 0 а < 0

Фокус (неустойчивый фокус, если а < 0, устойчивый фокус, если а > 0, центр, если а = 0)

3

4

5

В дальнейшем будем предполагать, что Ь = 0, В вышеприведенной таблице приведена классификация основных случаев поведения траекторий системы (7) па фазовой плоеко-

условия ограниченности решений на полуоси Теперь перейдем к изучению уравнения (2), Уравнение (2) эквивалентно системе

!

% У, (п\

у' = -ау - Ъхг — с\у — <р(х,у)\.

Ниже будем предполагать, что функция <р(х,у) удовлетворяет условию

lim M^i = 0. (9)

М+Ы^ \ж| + \у\

Условие (9) обеспечивает продолжимость решения системы (8) на всю числовую ось (-то, то).

Теорема 1. Пусть имеет место (9) и коэффициенты системы, (8) удовлетворяют условиям: b > 0, с > 0, а £ [min{0, 2Vb — с}, max{0, с — 2у/Ь}]. Тогда все решения системы (8) ограничены при t > 0, есл,и а > 0 и при t < 0, есл,и а < 0, т. е. а = 0 и для, любого решения (x(t),y(t)) справедливо неравенство

sup{\x(i)\ + \y(t)\ : а • t> 0} < то.

Доказательство. Пусть а > 0, Если некоторое решение (x(t),y(t)) системы (8) неогра-ничено при t > 0, то существует последовательность tk, 0 <tk < tk+i, к = 1, 2,... такая, что tk ^ то и \x(tk)\ + \y(tk)\ ^ то при к ^ то. Выберем чиела Tk i [0,tk] такие, что

\х(тк )\ + \у(тк )\ = nmax (\x(f)\ + \y(f)\).

Из неравенства \х(тк)\ + \у(тк)\ > \x(tk)\ + \y(tk)\ следует dk = \х(тк)\ + \у(тк)\ ^ то и, следовательно, Tk ^ то при к ^ то, Вектор-функции

(ик(t),vk(t)) = (х(тк + t)/dk,у(тк + t)/dk) — тк ^ t ^ 0

удовлетворяют условиям

К(t)\ + \vk(t)\ < \ик(0)\ + \vk(0)\ = 1, —тк ^ t ^ 0

и являются решением системы

= V, (10)

| v' = —av — bu — c\v — hk(i)\.

Здесь функции hk (t) определяются равенством

кк (Ь) = <^(ик (Ь),ьк (^)/(1к

и, в силу выбора чисел ¿к и условия (9), равномерно по Ь стремятся к нулю при к ^ то.

Пусть (и*, V*) - предельная точка последовательности (ик(0),ук(0)), Тогда решение системы (7), удовлетворяющее начальному условию (ж(0),у(0)) = (и*,у*), является ненулевым и ограниченным при Ь ^ 0, Но с другой стороны, в условиях теоремы, согласно п, 2 и 4, все решения системы (7) неограничены при Ь ^ 0, Полученное противоречие доказывает справедливость утверждения теоремы для случая а > 0,

Аналогично рассматривается случай а < 0, Теорема доказана.

Устойчивость изолированной особой точки

Рассмотрим общее уравнение второго порядка

х" = / (х,х'), (11)

где f (х, у) - определенная на всей плоскости (х,у) непрерывная функция. Уравнение (11) эквивалентно системе уравнений

!

ж = y,

У' = f (х,у). ^

Ниже мы будем предполагать, что любое решение системы (12) продолжимо на всю числовую ось. Как отметили выше, это предположение выполнено в частном случае, когда рассматривается уравнение (8), где функция <р(х,у) удовлетворяет условию (9),

Особые точки системы (12) лежат на оси абсцисс, т.е. имеют вид (х0, 0), где х0 является решением скалярного уравнения

f (х, 0) = 0. (13)

Пусть х0 - решение скалярного уравнения (13), Следует отметить, что исследование поведения траекторий системы (12) в окрестности точек равновесия является одной из основных задач качественной теории дифференциальных уравнений,

В дальнейшем нам понадобятся некоторые свойства, связанные с качественным поведением траектории системы (12) в окрестности особой точки (х0, 0), Для описания этих свойств используем методы функции Ляпунова [2, 3],

Лемма 1. Предположим, что функция f (х,у) в некоторой окрестности, \х — х0\ + \у\ < а, а > 0 особой точки (х0, 0) удовлетворяет условиям:

а) f (х, 0)(х — хо) ^ 0 и f (х, 0) тождественно не равна нулю в любом, интервале (-S, 0)u(0,S), 6> 0;

б) (f (х,у) — f (х, 0))у > 0.

Тогда, для любого решения (x(t),y(t)) системы, (12), отличного от стационарного решения (х0, 0), имеет .место

in0[\x(i) — Хо\ + \y(t)\] > 0. (14)

Доказательство. Предположим, что утверждение леммы не имеет места. Тогда суше-ствует решение (x(t),y(t)) системы (12), отличное от особой точки (х0, 0), и последовательность t\ < t2 < ... <tk < ... такие, что

lim \x(tk) — жо\ + \y(tk)\ = 0. (15)

Рассмотрим множество Е = {t : t > 0, \x(t) — х0\ + \y(t)\ < о\}, где 0 < о\ < min[a, \ж(0) — х0 \ + \у(0)\^. Отметим, что Е - открытое множество и в силу (15) tk G Е начиная с ^^таторого номера ^.Обозначим ч ерез (а к ,ßk) составляющие интервалы множества Е, ^^^^^^^щие числа tk- Не исключено, что интервалы (ak,ßk), соответствующие различным номерам к совпадают и ßk = ж начиная с некоторого номера к. По определению составляющего интервала имеем:

\x(t) — Х0\ + \y(t)\ < \х(ак) — Х0\ + \у(ак)\ = ох, ак <t ^ tk. (16)

Для решения (x(t),y(t)) справедливо тождество

И

- {у2/2 + G(x)} = h(x,y), (17)

где

¡•X

G(x) = — f 0)ds, h(x,y) = (f (х,у) — f (х, 0))у.

Интегрируя тождество (17) на отрезке [0,1к], имеем

уЦк)2/2 + С[хЦк)] - у(ак)2/2 — в[ак)] = [ * к(х(8),у(8)№. (18)

В силу условия б) к(х,у) > 0 при |ж — жо| + |у| < а, поэтому подынтегральная функция в правой части (18) неотрицательна. Следовательно имеет место неравенство

уЦк)2/2 + С[хЦк)] — у(ак)2/2 — С[ак)] > 0,к > ко. (19)

Так как в силу условия а) непрерывная функция С(х) положительна при 0 < |ж — ж0| < а, то у2/2 + С[х] > т > 0 при |ж — ж0| + |у| = о\. Поэтому неравенство (19) противоречит равенству (15), Полученное противоречие доказывает утверждение леммы. Лемма доказана.

Отметим, что в случае, когда функция f (х,у) имеет непрерывные частные производные

(х,у), щ(х, у) в некоторой окрестности особой точки (х0, 0), для выполнения условий а) и б) леммы 1 достаточно выполнения неравенства (х0, 0) < 0, щ(х0, 0) > 0,

Следующая лемма является обобщением леммы 1,

Лемма 2. Предположим, что функция /(х,у) в некоторой окрестности, 1х — ж0| + |у| < а, а > 0 особой точки (х0, 0) и при заданном к ^ 0 удовлетворяет условиям:

а) функция /(х, п(х — х0))(х — х0) — к2(х — х0)2 > 0, причем, тождественно не равна нулю в любом, интервале (—8, 0) м (0,^), 8 > 0;

б) а(х, к(х — Х0) + у) — ку — f (х, к(х — Х0)))у > 0.

Тогда, для любого решения (х(Ь),у(1)) системы (12), отличного от стационарного решения (х0, 0), имеет место

т0М) — + Ш1] > 0. (20)

Аналогично лемме 1 доказывается следующая

Лемма 3. Предположим,, что функция /(х,у) в некоторой окрестности, 1х — х01 + 1у1 < а, а > 0 особой, точки (х0, 0) удовлетворяет условию а) леммы, 1 и условию:

в) (I (х,у) — I (х, 0))у ^ 0.

Тогда, для любого решения (х(Ь),у(1)) системы (12), отличного от стационарного решения (х0, 0), имеет место

ш0[И*) — + 1уШ > 0. (21)

Свойство (21) решений системы (12) тесно связано со свойством устойчивости по Ляпунову стационарного решения (х0, 0), А именно имеет место

Теорема 2. Стационарное решение (х0, 0) системы (12) устойчиво по Ляпунову тогда, и только тогда, когда, для, любого решения, отличного от него, имеет место неравенство (21).

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

Продолжим изучение системы (8) предполагая, что Ь > 0,с > 0, функция <р(х,у) удовлетворяет условию (9), р(0, 0) > 0 и система (8) имеет единственную особую точку, т.е. скалярное уравнение

Ъх + с^(х, 0)| =0 (22)

имеет единственное решение х0, причем х0 = 0.

Теорема 3. Предположили что: а) коэффициенты уравнения (2) удовлетворяют условиям а > шах|0,с — 2^/Ъ}; б) функция, <р(х,у) в некоторой окрестности \х — х0\ + \у\ < а, а > 0 особой точки (х0,0) удовлетворяет условиям: ((с — а)у + с<р(х,у) — с^(х, 0))у > 0. Тогда, для любого решения (х(Ь),у(Ь)) системы (8), отличного от стационарного, и любой последовательности, Ък ^ существует подпоследовательность Ъ^ такая, что решения (х(Ь + Ъ^),у(Ь + Ъ^)) равномерно на, каждом отрезке прибллгжлются, к некоторому периодическому решению системы (8) при

Доказательство. Пусть (х(Ь),у(1)) - произвольное решение системы (8), отличное от периодического и стационарного, В силу теоремы 1 это решение ограниченно при Ь > 0, Следовательно, для любой последовательности Ък ^ последовательность решений (х(Ь + Нк),у^ + Ък)) определена, ограничена и равностепенно непрерывна на каждом отрезке, В силу теоремы Арцела существует подпоследовательность hkj такая, что решения (х\(Ь + hkj),х2(Ь + hkj)) приближаются к некоторому решению (х*(Ь),у*(Ь)) системы (8) равномерно на каждом отрезке при ] ^

В силу леммы 1, множество ш— предельных точек решения (х(1),у(Ь)) не содержит единственную особую точку системы (8), Поэтому в силу теоремы Пуапкаре-Бепдиксопа [9] решение (х*(Ь), х*2(Ь)) - периодическое. Теорема доказана. Рис, 2 иллюстрирует утверждение теоремы.

Теорема 4. Предположим, что: а) коэффициенты уравнения, (2) удовлетворяют условиям а < шт{0,2\[Ъ — с}; б) функция (р(х,у) в некоторой окрестности \х — х0\ + \у\ < а, а > 0 особой, точки (х0,0) удовлетворяет условиям: ((с — а)у + ср(х,у) — с^(х, 0))у ^ 0. Тогда, для любого решения (х(Ь),у(Ь)) системы (8), отличного от стационарного, и, любой последовательности, Ък ^ существует подпоследовательность Ък5 такая, что решения (х(Ь — Ъ^),у(Ь — Ъ^)) равномерно на, каждом отрезке прибллгжлются, к некоторому периодическому решению системы, (8) при

Рис. 2. Предельный цикл

В качестве примера рассмотрим функцию <р(х,у) = со&(х + у). Она удовлетворяет условию (9) и р(0, 0) = 1 > 0. Следующее утверждение определяет условию единственности решения уравнения (22).

Лемма 4. Уравнение ¿х + \ еов(ж)\ = 0, где d - заданное число, имеет единственное решение тогда и только тогда, когда, \d\yf 1 + х\ > 1, где х\ Е (ж/2, ж) и является решением, уравнения ооб(х) + х вт(ж) = 0 (х\ ~ 2, 798386, \/1+х{ ~ 0, 336508).

Отметим, что утверждение леммы 4 в случае Щ > 1, является следствием принципа сжимающих отображений, а в случае, когда 1 < 1 + х\ ^ \J 1 + х\ следует из анализа графиков функции z = — I cos(x) | и прямой z = dx.

Пусть коэффициенты b, с такие, что 0 < с < Ъл/\ + X i. Тогда уравнение bx+cI cos(x)| = 0, согласно лемме 4, имеет единственное решение х0 = х0(Ь/с) < 0, зависящее от дроби Ь/с. Поэтому, если коэффициенты а,Ь,с системы (8) удовлетворяют неравенствам max{0, с — 2л/Ь} < а < с(1 — sin(xo)), то выполнены все условия теоремы 3 при <р(х,у) = cos(x + у). Следовательно уравнение (2) имеет предельный цикл,

СПИСОК ЛИТЕРАТУРЫ

1. Арнольд В.И. Геометрические методы в теории обыкновенных дифференциальных уравнений. Москва-Ижевск: Редакция журнала «Регулярная и хаотическая динамика». 2000. 400 с.

2. Филиппов А.Ф. Введение в теорию дифференциальных уравнений. М.: УРСС. 2004. 239 с.

3. Баутин H.H., Леонтович Е.А. Методы и приемы качественного исследования динамических систем на плоскости. М.: Наука. 1976. 496 с.

4. Юмагулов М.Г. Введение в теорию динамических систем: учебное пособие. СПб.: Из-во «Лань». 2015. 272 с.

5. Иванов А.П. Исследование разрывных бифуркаций в негладких динамических системах. Нелинейная динамика // 2012. Т. 8. № 2. С. 231-247.

6. Каток А.Б., Хасселблат Б. Введение в теорию динамических систем. М.: МЦНМО. 2005. 454 с.

7. R.I. Leine, D.H. Van Campen Bifurcation phenomena in non-smooth dynamical system,s // European Journal of Mechanigs A/Solids. 5(2006). Pp. 595-616.

8. Мухамадиев Э.М., Нуров И.Д., Халилова М.Ш. Предельные циклы, кусочно-линейных дифференциальных уравнений второго порядка, // Уфимский математический журнал. 2014. Т. 6. № 1. С. 84-93.

9. Хартман Ф. Обыкновенные дифференциальные уравнения. М.: Мир. 1970. 720 с.

Муллошараф Курбонович Арабов,

Институт математики им. А. Джураеви АН РТ,

ул. Айни 299/1,

734063 г. Душанбе, Таджикистан E-mail: cool.araby@mail.ru

Эргашбой Мухамадиев,

Вологодский государственный университет,

ул. Ленина, 15,

160000, г. Вологда, Россия

E-mail: emuhamadiev@rambler.ru

Иехокбой Джумаевич Нуров, Таджикский национальный университет, проспект Рудаки, 17, 734000, г. Душанбе, Таджикистан E-mail: nidl@mail.ru

Хуршед Илхомиддинович Собиров, Таджикский национальный университет, проспект Рудаки, 17, 734000, г. Душанбе, Таджикистан E-mail: hurshed. sobirov@mail. ru