Работа посвящена анализу явления бифуркации рождения предельного цикла систем дифференциальных уравнений, при которых состояния равновесия меняют тип устойчивости. В частности рассмотрена бифуркация состояния равновесия типа фокус, при которой происходит появление или исчезновение предельных циклов. Исторически эта тема восходит к трудам А.М. Ляпунова, а так же А. Пуанкаре, который в 1892 г. дал метод исследования устойчивости состояния равновесия с парой чисто мнимых характеристических корней и ввёл величины, получившие впоследствии название ляпуновских. Заслуга открытия бифуркации рождения предельного цикла из состояния равновесия с чисто мнимыми корнями при изменении параметров системы и обнаружение связи этой бифуркации с ляпуновскими величинами принадлежит А. А. Андронову (1931 и 1937 гг.). Перенесение теории на системы с числом уравнений больших двух было выполнено Н. Н. Баутиным в 1941 году. Большой вклад в развитие теории бифуркации внёс математик Х. Хопф, который в 1942 г. сделал обобщения для размерности больше 2.
Задача о нахождении бифуркации рождения предельного цикла из состояния равновесия сводится к нахождению так называемых "ляпуновских величин". Для нахождения первой ляпуновской величины (L(0)) применяются формулы, приведённые в [2]. Для выполнения настоящей работы в среде Maple 17 разработана процедура, позволяющая вычислять эти значения для каждой конкретной системы уравнений.
Фазовые портреты простейших систем. Дана линейная система:
И предположим, что матрица
имеет собственные числа 
При этом будем считать, что 
Если 
одного знака, то неподвижная точка в начале координат плоскости 
(где 
- собственные вектора) называется узлом. Если 
то все траектории направлены от начало координат (к нему, если 
), то узел называется неустойчивым (устойчивым). Устойчивый узел представлен на рисунке 1. неустойчивый- рисунке 2.
Рис. 1. Устойчивый узел
Рис. 2. Неустойчивый узел
Если имеют противоположные знаки, то возникает фазовый портрет, изображенный на рисунке 3. Траектории направлены к началу координат (от него), если соответствующее собственное значение отрицательно (положительно). В таком случае начало координат называется седлом или седловой точкой.
Рис. 3. Седло
Рис. 4. Неустойчивый дискретный узел
Когда собственные числа имеют равные значения, случай соответствует узлу специального вида, который называется дикретическим узлом. Этот узел устойчив при 
неустойчив при 
Все его нетривиальные траектории являются радиальными лучами. Фазовый портрет, соответствующий точке типа дикретический узел, изображен на рисунке 4.
Если матрица не диагональная, т.е. имеет вид 
, то фазовый портрет имеет вид, представленный на картинках 5 и 6.
Рис. 5. Вырожденный неустойчивый узел
Рис. 6. Вырожденный устойчивый узел
В этом случае начало координат называется вырожденным узлом; он устойчив при 
и неустойчив при
Когда собственный значения комплексные числа, матрица имеет вид 
. Если α ≠ 0, то начало координат называется фокусом (неустойчивым при α > 0, устойчивым при α < 0). Фазовый портрет подобных систем представлен на рисунке 7.
Если α = 0, то начало координат называется центром и фазовый портрет состоит из континуума концентрических кругов. Пример фазового портрета подобной системы представлен на рисунке 8. Это единственный случай, когда в линейных системах возникает периодическое решение.
Другой тривиальный случай периодических решений для линейных систем- это неподвижные точки.
Рис. 7. Неустойчивый фокус
Рис. 8. Центр
Исследование первой ляпуновской величины. Рассмотрим систему дифференциальных уравнений вида:
Тогда значение первой ляпуновской величины рассчитывается по формуле:
В качестве примера рассмотрим систему из [1]:
Где µ - вещественный параметр, причём µ> 0 и достаточно мало.
При µ = 0 получим систему вида:
Рассмотрим состояние равновесия O (0, 0) линейного приближения данной системы. Для подсчёта в среде Maple17 была написана процедура, позволяющая вычислить значение первой ляпуновской по формуле для L(0). Тогда L(0)= π/72 то есть O (0, 0) при µ = 0 является неустойчивым фокусом. Этот вывод подтверждается фазовым портретом данной системы (рис. 9).
Рис. 9: Фазовый портрет системы (1).
Как и для двумерных систем, для трёхмерных в среде Maple 17была разработана процедура, позволяющая рассчитать значение первой ляпуновской величины для каждой конкретной системы уравнений.
Аналогично двумерному случаю, определим вид системы дифференциальных уравнений, которая и будет объектом исследования.
(1)
где
Будем говорить, что состояние равновесия системы (1) находится в точке О (0, 0). Это означает, что в данной точке Q1(x,y, z) = Q2(x, y, z) = Q3(x, y, z)= 0. Так же следует отметить, что все три функции - гладкие и имеют производные до порядка N включительно.
Для дальнейшего исследования автору необходимо было бифуркационная теорема[2, 4]. На основании анализа имеющихся источников автором была предложена собственная трактовка доказательства бифуркационной теоремы, формулировка которой приведена ниже также в авторской редакции.
Рассмотрим систему:
Наше первое предложение заключается в том, что
(а)
Второе условие состоит в том, что первая ляпуновская величина не равна нулю.
(б)
Предложения (а) и (б) означают, что О (0, 0)- нелинейный фокус при µ=0. Устойчивый при 
и неустойчивый при 
Третье предположение состоит в условии, что скорость не равна нулю.
(в)
Пусть выполнены предположения (а), (б), (в). Тогда найдётся такое значение параметра σ > 0 и такая окрестность U точки (x,y) = (0, 0), что
(I) Если 
и 
то система (1) имеет в точности один предельный цикл внутри окрестности U;
(II) Если и 
то система (1) не имеет периодических орбит внутри окрестности U.
Кроме того, предельный цикл, возникающий в случае (I) устойчив (неустойчив) при 
, и стягивается к неподвижной точке, когда 
В качестве примера рассмотрим модель Лоренца турбулентного движения жидкости. Данная модель описывается системой:
Где b, r, σ- положительные параметры. Для r ≤ 1существует единственное состояние равновесия в начале координат. При возрастании r от значения r = 1 из состояния равновесия в начале координат появляются устойчивые состояния равновесия:
и
,
которые могут поменять устойчивость, когда корни их характеристических уравнений становятся комплексными и их действительные части проходят через ноль.
Пользуясь рассуждениями, приведенными в[3], приходим к выводу, что значение первой ляпуновской величины для данной системы неотрицательно, однако точного значения не приводится.
Для определения точного значения первой ляпуновской величины в среде Maple 17 была разработана процедура, которая, позволила рассчитать эту величину для системы Лоренца при σ = 10 и β = 8/3. Получилось, что L(0) = 10.3764756164464π, что подтверждает результат, описанный в[6].
В случае, когда первая ляпуновская величина положительна, говорят, что граница области опасная. Таким образом, в модели Лоренца граница области устойчивости является опасной.
Для наглядности и подтверждения теоретических выкладок построены фазовые портреты системы Лоренца для значений параметров, для которых было рассчитано значение первой ляпуновской величины. Эти фазовые портреты приведены на рис.10,11.
Рис. 10 Фазовый портрет системы Лоренца
11 Фазовый портрет системы Лоренца
Вывод: Формула для вычисления первой ляпуновской величины хорошо известна. С целью облегчения процесса выявления знака данной величины, была составлена процедура, позволяющая точно определить значение первой ляпуновской величины. Удобные сервисы позволяют представлять полученные данные графически как для двумерных систем, так и для трёхмерных.
Полученный результат согласован с результатами других работ и подтверждается средствами качественной теории, а именно, видом фазовых портретов.