Об одном случае исследования принудительной синхронизации методом приближенных точечных отображений

Размещено на http://www.allbest.ru/

ОБ ОДНОМ СЛУЧАЕ ИССЛЕДОВАНИЯ ПРИНУДИТЕЛЬНОЙ СИНХРОНИЗАЦИИ МЕТОДОМ ПРИБЛИЖЕННЫХ ТОЧЕЧНЫХ ОТОБРАЖЕНИЙ

Антоновская О.Г.,

Зайцева М.Н.

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

Совершенно естественно, что наиболее доступными для исследования являются колебательные системы с малой нелинейностью, поскольку к ним в той или иной форме можно применять методы теории возмущений. Поэтому для слабо нелинейных систем имеется ряд достаточно общих асимптотических методов, применимых ко многим типичным классам колебательных систем. Основными из них являются метод малого параметра и метод усреднения, а также методы Н. Н. Боголюбова и Н. М. Крылова, представляющие собой развитие метода усреднения и получение оценок более высокого порядка [1].

Следует отметить, что до сих пор особый интерес представляет изучение систем, близких к гармоническому осциллятору (квазигармонический осциллятор) [1, С. 19-21], [2, С. 650-663]. Одним из важных достоинств такой системы является возможность использовать хорошо известные математические свойства процессов колебаний гармонического осциллятора с медленно меняющейся частотой в различного вида задачах: от задач обработки сигналов [3] до исследования неравновесных экономических систем [4].

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

Известно [5], [6, С. 5-6], что при исследовании динамики синхронизуемого осциллятора

, (1)

в котором 0 < м << 1, а 2pр - период внешней силы, методом точечных отображений [7], исследование поведения траекторий (1) может быть сведено к изучению точечного отображения Tсекущей поверхности фазового пространства в себя [5] (или секущей поверхности t = 0 в секущую поверхность ), порожденного траекториями системы. При этом с точностью до величин порядка м² точечное отображение T может быть приближено точечным отображением с функциями последования

(2)

(3)

где

, (4)

. (5)

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

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

(6)

где 0<м<<1, A>0, или, если ввести y = x, систему двух уравнений первого порядка

(7)

Задача состоит в установлении условий существования у (7) периодического решения с периодом 2р.

В этом случае точечное отображение , приближающее отображение T секущей поверхности t = 0фазового пространства x, y = x, t в секущую поверхность , порождаемое траекториями системы (7), с точностью до членов порядка имеет вид

(8)

(9)

Условия существования неподвижной точки точечного отображения (8)-(9) дают соотношения

(10)

где находится из уравнения

(11)

То есть факт существования (10) неподвижных точек определяется фактом существования корней у (11).

Детальный анализ уравнения (11) позволяет получить картину резонансных кривых при различных значениях (рис. 1-2). В самом деле, при уравнение (11) имеет единственный корень , причем при . Т.е. при резонансная кривая является разомкнутой для любого A (рис. 1). В случае для значений при существует три корня уравнения (11) , два из которых сливаются между собой при и исчезают через значение . При для малых о существует единственный корень уравнения (11), с ростом о корней становится три, а затем вновь один. При существует единственный корень уравнения (11). Кроме того, следует отметить, что резонансная кривая имеет горизонтальную касательную при , а вертикальную - в точках кривой

Рис. 1 - Вид резонансных кривых при

Рис. 2 - Вид резонансных кривых при

, (12)

которая существует только при и представляет собой эллипс с центром в точке , и с главными диаметрами соответственно. Т.е. при резонансная кривая имеет две ветви - замкнутую и разомкнутую - при . При замкнутая и разомкнутая ветви сливаются, образуя единственную разомкнутую ветвь при (рис. 2).

Устойчивость неподвижных точек точечного отображения в случае их существования определяется корнями характеристического полинома

. (13)

Корни полинома P(z) (13) будут действительными при и комплексно-сопряженными при , т.е. граница , соответствующая уходу пары корней характеристического полинома с действительной оси в этом случае будет иметь вид двух полупрямых

(14)

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

Для получаем уравнение (14).

Уравнение есть

(15)

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

, (16)

т.е. в случае ее существования представляет собой куски эллипса с центром в точке и главными диаметрами по о и по с, принадлежащие области . Заметим, что при всегда имеет точку пересечения с осью , которой соответствует .

Рис. 3 - Картина D-разбиения при

Рис. 4 - Картина D-разбиения при

Анализ существования и взаимного расположения границ (14)-(16) позволил получить картину D-разбиения при различных значениях з и малых м (рис. 3-6). Границы D-разбиения приведены с соответствующей штриховкой в сторону выхода корней характеристического уравнения из единичного круга. Однократная штриховка соответствует бифуркации корней на действительной оси, двойная - бифуркации комплексно-сопряженных корней. Область D=0 есть область устойчивости.

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

Рис. 5 - Картина D-разбиения при

Рис. 6 - Картина D-разбиения при

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

(17)

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

Факт существования неподвижных точек точечного отображения с некоторым характером устойчивости при A=const и различных о может быть установлен наложением картины поведения границ областей существования различных типов неподвижных точек на плоскости при заданном м на плоскость с резонансной кривой при заданном A.

Таким образом, во-первых, область устойчивости является ограниченной при любом конечном, хотя и малом м, все более расширяясь при . Во-вторых, при переходе на плоскости через границу N- имеет место серия бифуркаций удвоения периода и переход к хаосу (рис.7), т.е. при существует множество примыкающих друг к другу и все более сжимающихся областей, соответствующих режимам кратности 2, 4, 8, 16 и т.д., переходящих с ростом с в область хаоса, таких, что попадание резонансной кривой в соответствующую область означает существование периодического режима определенной кратности. Границы этих областей сходятся в общие точки, являющиеся точками стыковки границ .

Необходимо также отметить, что устойчивый режим, в случае его существования, имеет в фазовом пространстве отображения ограниченную область притяжения, расширяющуюся при . При дальнейшем движении по резонансным кривым при о=const в сторону увеличения с, после перехода через границу , которая, по-видимому, является границей вырождения математической модели, ограниченного притягивающего множества в фазовом пространстве не обнаружено. Устойчивость бесконечности в этом случае может быть доказана следующим образом.

Рис. 7 - Границы областей, соответствующих кратным режимам

Рассмотрим простейшую функцию Ляпунова

. (18)

Ее первая разность в силу формул (8) -(9) будет удовлетворять условию

(19)

Где

(20)

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

. (21)

а значит бесконечность всегда устойчива.

Область хаоса не является сплошной. В ней имеются “окна” [8 С. 272-275], в которых существуют режимы различной кратности. Об этом свидетельствуют данные, приведенные на рис. 8, где указаны граничные значения A, соответствующие областям кратных режимов (“окнам”). Последовательность найденных кратностей очень напоминает порядок смены кратностей неподвижных точек для одномерного отображения, установленный А.Н.Шарковским [9]. Причем переход к хаосу из этих окон может сопровождаться, а может не сопровождаться удвоением периода.

Рис. 8 - Пример структуры области хаоса

Заметим, что задача нахождения условий существования у квазигармонического осциллятора периодического решения с периодом 2р решается как задача нахождения условий существования неподвижных точек с определенным характером устойчивости у приближенно построенного точечного отображения , и выводы относительно существования и характера устойчивости неподвижных точек в приведенных выше исследованиях (а значит и соответствующих им периодических решений исходной системы) были сделаны на основе изучения свойств аппроксимирующего точечного отображения. При этом естественно встает вопрос об адекватности поведения траекторий точного и приближенного точечных отображений, т.е. о применимости полученных результатов [6, С. 81-88]. И, подобно [10], для решения этого вопроса можно воспользоваться результатами работы [5], [6, С. 123-135], и рассматривать задачу о применимости результатов приближенного исследования, оценивая степень близости приближенного точечного отображения к точному отображению.

синхронизация квазигармонический осциллятор точечный

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

1. Боголюбов Н. Н. Асимптотические методы в теории нелинейных колебаний / Н. Н. Боголюбов, Ю. А. Митропольский. - М.: Наука, 1974. - 504 с.

2. Андронов А. А. Теория колебаний / А. А. Андронов, А. А. Витт, С. Э. Хайкин. - М.: Физматгиз, 1959. - 916 с.

3. Журавлев В. М. Построение огибающей и локальной частоты стохастического процесса на основе модели осциллятора с флуктуирующей частотой / В. М. Журавлев, П. П. Миронов, С. В. Летуновский // Изв. вузов. Поволжский регион. Физико-математические науки. - 2013. -№3(27). - С. 159-169.

4. Ивинская Е. Ю. Теоретические аспекты исследования неравновесных экономических систем на основе модели гармонического осциллятора / Е. Ю. Ивинская // Теория и практика общественного развития. Экономические науки. - 2015. - № 21. - С. 57-59.

5. Антоновская О. Г. О влиянии насыщения нелинейности на результаты исследования принудительной синхронизации методом приближенных точечных отображений / О. Г. Антоновская // Математическое моделирование и оптимальное управление: Вестник ННГУ, Нижний Новгород. - 1999. - № 2(21). - С. 198-208.

6. Антоновская О. Г. Метод точечных отображений в задачах нелинейной дингамики / О. Г. Антоновская, В.И. Горюнов. - Гамбург: LAP Lambert Academic Publishing, 2013. - 140 с.

7. Неймарк Ю. И. Метод точечных отображений в теории нелинейных колебаний / Ю. И. Неймарк. - М.: Наука, 1972. - 472 с.

8. Неймарк Ю. И. Стохастические и хаотические колебания / Ю. И. Неймарк, П. С. Ланда. - М.: Наука, 1987. - 424 с.

9. Шарковский А. Н. Сосуществование циклов непрерывного преобразования прямой в себя / А. Н. Шарковский // Укр. мат. журн. - 1964. - Т.26. - № 1. - С. 61-71.

10. Антоновская О. Г. О приближенном исследовании близкого к тождественному точечного отображения плоскости в плоскость / О. Г. Антоновская // Математическое моделирование и оптимальное управление: Вестник ННГУ, Нижний Новгород. - 2004. - № 1(27). - С. 63-69.

Размещено на Allbest.ru