МИНИСТЕРСТВО ОБРАЗОВАНИЯ РЕСПУБЛИКИ БЕЛАРУСЬ
Учреждение образования
«Гомельский государственный университет
имени Франциска Скорины»
Математический факультет
Кафедра дифференциальных уравнений
Качественное исследование в целом двумерной квадратичной стационарной системы с двумя частными интегралами в виде кривых второго порядка
Дипломная работа
Исполнитель:
студентка группы М-51 БРАВАЯ Е.Н.
Научный руководитель:
доцент, к. ф-м. н. ФИЛИПЦОВ В.Ф.
Рецензент:
профессор, д. ф-м. н. СТАРОВОЙТОВ Э.И.
Гомель 2003
Реферат
Дипломная работа 38 страниц, 11 источников.
Ключевые слова и словосочетания: квадратичная двумерная стационарная система, частный интеграл, парабола, гипербола, окружность, точка, характеристическое уравнение, характеристическое число, узел, седло, фокус.
Данная работа содержит результаты исследований автора, относящиеся к качественному исследованию в целом двумерной квадратичной стационарной системы.
Основным инструментом исследований является понятие частного интеграла.
Работа состоит из двух глав.
В первой главе проводится построение квадратичных двумерных стационарных систем с заданными интегралами, при этом коэффициенты интегралов выражаются через коэффициенты системы, а коэффициенты системы связаны между собой тремя соотношениями.
Во второй главе проводится качественное исследование в целом выделенных в первой главе классов систем при фиксированных значениях некоторых параметров.
Содержание
Реферат
Введение
1. Построение квадратичных двумерных стационарных систем
1.1 Построение квадратичной двумерной стационарной системы с частным интегралом в виде параболы
1.2 Построение квадратичной двумерной стационарной системы с частным интегралом в виде окружности либо гиперболы
1.3 Необходимые и достаточные условия существования у системы (1.1) двух частных интегралов (1.3), (1.13)
2. Качественное исследование построенных классов систем
2.1 Исследование системы (1.1) с коэффициентами, заданными формулами (1.28) — (1.31)
2.2 Исследование системы (1.1) с коэффициентами, заданными формулами (1.41) — (1.42)
2.3 Исследование системы (1.1) с коэффициентами, заданными формулами (1.52) — (1.53)
Заключение
Список использованных источников
Приложение А
Приложение Б
Приложение В
Введение
Известно, что в элементарных функциях и даже в квадратурах интегрируются очень немногие классы дифференциальных уравнений. В связи с этим появилась необходимость в создании такой теории, с помощью которой можно было бы изучать свойства решений дифференциальных уравнений по виду самих уравнений. Такой теорией, наряду с аналитической, и является качественная теория дифференциальных уравнений.
Впервые задача качественного исследования для простейшего случая системы двух дифференциальных уравнений с полной отчетливостью была поставлена А. Пуанкаре [7] в конце прошлого столетия. Позднее исследования А. Пуанкаре были дополнены И. Бендиксоном [3, с. 191-211] и уточнены Дж.Д. Биркгофом [4, с.175-179].
/> (0.1)
Одной из задач качественной теории дифференциальных уравнений является изучение поведения траекторий динамической системы (0.1) на фазовой плоскости в целом в случае, когда P (x,y) и Q (x,y) — аналитические функции. Интерес к изучению этой системы или соответствующего ей уравнения объясняется их непосредственным практическим применением в различных областях физики и техники.
/> (0.2)
Имеется много работ, в которых динамические системы изучались в предположении, что их частными интегралами являются алгебраические кривые. Толчком к большинству из них послужила работа Н.П. Еругина [6, с.659 — 670], в которой он дал способ построения систем дифференциальных уравнений, имеющих в качестве своего частного интеграла кривую заданного вида.
Знание одного частного алгебраического интеграла системы (0.1) во многих случаях помогает построить полную качественную картину поведения интегральных кривых в целом. Отметим ряд работ этого характера для систем (0.1), в которых P (x,y) и Q (x,y) — полиномы второй степени.
Н.Н. Баутиным [1, с.181 — 196] и Н.Н. Серебряковой [8, с.160 — 166] полностью исследован характер поведения траекторий системы (0.1), имеющей два алгебраических интеграла в виде прямых. В [10, с.732 — 735] Л.А. Черкасом такое исследование проведено для уравнения (0.2) при наличии частного интеграла в виде кривой третьего порядка. Яблонский А.И. [11, с.1752 — 1760] и Филипцов В.Ф. [9, с.469-476] изучали квадратичные системы с предположением, что частным интегралом являлись алгебраические кривые четвертого порядка.
В данной работе рассматривается система
/> (0.3)
и проводится качественное исследование в целом системы (0.3) при условии, что частным интегралом является кривая четвертого порядка, которая распадается на две кривые второго порядка, одна из которых парабола, вторая окружность или гипербола.
Работа состоит из двух глав.
В первой главе проводится построение квадратичных двумерных стационарных систем с заданными интегралами, при этом коэффициенты интегралов выражаются через коэффициенты системы, а коэффициенты системы связаны между собой тремя соотношениями.
Во второй главе проводится качественное исследование в целом выделенных в первой главе классов систем при фиксированных значениях некоторых параметров.
1. Построение квадратичных двумерных стационарных систем
1.1 Построение квадратичной двумерной стационарной системы с частным интегралом в виде параболы
Рассмотрим систему дифференциальных уравнений
/> (1.1)
Пусть система (1.1) имеет частный интеграл вида:
/>, (1.2)
где Fk (x,y) — однородные полиномы от x и y степени k.
В качестве частного интеграла (1.2) возьмем параболу вида:
F (x,y) º y+a1x2 +a2x+a3 = 0 (1.3)
Будем предполагать, что a3¹ 0, то есть парабола не проходит через начало координат.
Согласно [10, с.1752-1760] для интеграла (1.3) системы (1.1) имеет место соотношение:
/>/>/>, (1.4)
где L (x,y) = px+my+n, p, m, n — постоянные.
Тогда следуя формуле (1.4) получим равенство:
(2a1x+a2) (ax+by+a1x2+2b1xy+c1y2) + (cx+dy+a2x2+2b2xy+c2y2) = = (y+a1x2+a2x+a3) (px+my+n).
Приравнивая коэффициенты при одинаковых степенях xmyn слева и справа, получим равенства:
(2a1-p) a1= 0 (1.51)
(4b1-m) a1= 0 (1.52)
2a1c1= 0 (1.53)
(2a-n) a1+ (a1-p) a2+a2= 0 (1.61)
2a1b+ (2b1-m) a2+2b2+p= 0 (1.62)
a2c1+c2-m= 0 (1.63)
(a-n) a2-pa3n+c= 0 (1.71)--PAGE_BREAK--
a2b-a3m+d-n= 0 (1.72)
a3n= 0 (1.73)
Пусть a1¹ 0, тогда из равенств (1.51), (1.52), (1.53), (1.63) и (1.73) получаем, что
P=2a1, m=4b1, c1=0, c2=4b1, n=0 (1.8)
Из соотношений (1.61), (1.62) и (1.71) найдем выражения коэффициентов кривой (1.3) через коэффициенты системы (1.1) в следующем виде:
a1/>, (1.9)
a2/>, (1.10)
a3/>. (1.11)
Равенство (1.72) с учетом полученных выражений (1.9) — (1.11), даст условие, связывающее коэффициенты a, b, c, d, a1, a2, b1, b2:
/> (1.12)
Итак, установлена следующая теорема:
Теорема 1.1 Система (1.1) имеет частный интеграл (1.3), коэффициенты которого выражаются формулами (1.9) — (1.11), при условии, что коэффициенты системы связаны соотношением (1.12) и c1= 0, c2= 4b1, a1¹0, 2b1a-a1b¹0.
1.2 Построение квадратичной двумерной стационарной системы с частным интегралом в виде окружности либо гиперболы
Пусть теперь система (1.1) наряду с интегралом (1.3) имеет интеграл в виде:
y2+sx2+bx+gy+d=0 (1.13)
Будем рассматривать теперь систему:
/> (1.14)
Согласно формуле (1.4), где L
(x,y) = m1x+n1y+p1,
m1, n1, p1 — постоянные для системы (1.1), имеем:
(2a1-m1) s2= 0 (1.151)
(4b1-n1) s+2a1= 0 (1.152)
m1= 4b2 (1.153)
n1=8b1 (1.154)
(2a-p1) s+ (a1-m1) b+a2g=0 (1.161)
2bs+ (2b1-n1) b+ (2b2-m1) g+2c= 0 (1.162)
(4b1-n1) g+2d-p1= 0 (1.163)
(a-p1) b+cg+m1d= 0 (1.171)
bb+ (d-p1) g-n1d= 0 (1.172)
p1d= 0 (1.173)
Предположим, что кривая не проходит через начало координат, то есть d¹0.
Пусть s¹0, тогда из равенств (1.151), (1.153), (1.154) и (1.173) получаем, что
m1=4b2, n1=8b1, a1=2b2, p1=0 (1.18)
А из соотношений (1.161), (1.163) и (1.171) найдем выражения коэффициентов кривой (1.13) через коэффициенты системы (1.1) в следующем виде:
/> (1.19)
/> (1.20)
/> (1.21)
/> (1.22)
Подставляя коэффициенты s, b, g и d в равенства (1.162) и (1.172), получим два условия, связывающие коэффициенты a, b, c, d, a2, b1, b2:
/> (1.23)
/> (1.24)
Итак, установлена следующая теорема:
Теорема 1.2 Система (1.14) имеет частный интеграл (1.13), коэффициенты которого выражаются формулами (1.19) — (1.22), при условии, что коэффициенты системы связаны соотношениями (1.23), (1.24) и b1¹0, b2¹0, a1=2b2.
1.3 Необходимые и достаточные условия существования у системы (1.1) двух частных интегралов (1.3), (1.13)
В разделах 1.1-1.2 мы получили, что система (1.1) будет иметь два частных интеграла в виде кривых второго порядка при условии, что коэффициенты системы связаны соотношениями:
/>
/>(1.25)
/>
Причемb1¹0, b2¹0, a1¹0, b1a-b2b¹0.
Выражая c из первого уравнения системы (1.25), получим
/> (1.26)
Подставим (1.26) во второе и третье уравнения системы (1.25). Получим два соотношения, связывающие параметры a, b, d, a2, b1, b2:
/>
Пусть />и
/> (1.27)
Из первого уравнения системы (1.27) получим
/>
Подставляя />во второе уравнение системы (1.27), найдем
/>.
Из соотношений (1.25) при условиях (1.27) получаем, что коэффициенты системы (1.1) определяются следующими формулами:
/> (1.28)
/> (1.29)
/> (1.30)
/>, />, />, />, /> (1.31) продолжение
--PAGE_BREAK--
Равенства (1.9) — (1.11), (1.19) — (1.22) при условии, что имеют место формулы (1.28) — (1.31), дадут следующие выражения для коэффициентов интегралов (1.3) и (1.13):
a1/> (1.32)
a2/> (1.33)
a3/> (1.34)
s/> (1.35)
b/> (1.36)
g/> (1.37)
d/> (1.38)
Теорема 1.3 Система (1.1) имеет частные интегралы вида (1.3) и (1.13) с коэффициентами, определенными формулами (1.32) — (1.38), при условии, что коэффициенты системы (1.1) выражаются через параметры по формулам (1.28) — (1.31).
Пусть
/> (1.39)
Из первого уравнения системы (1.39) найдем
/>, />.
Подставляя /> во второе уравнение системы (1.39), получим равенство:
/> (1.40)
Поскольку />, то рассмотрим два случая:
/>, тогда />.
Из соотношений (1.25) при условиях (1.39) и (1.40) получаем, что коэффициенты системы (1.1) определяются следующими формулами:
/>, />, /> (1.41)
/>, />, />, />, /> (1.42)
Равенства (1.9) — (1.11), (1.19) — (1.22) при условии, что имеют место формулы (1.41) — (1.42), дадут следующие выражения для коэффициентов интегралов (1.3) и (1.13):
a1/> (1.43)
a2/> (1.44)
a3/> (1.45)
s/> (1.46)
b=0 (1.47)
g/> (1.48)
d/> (1.49)
Теорема 1.4 Система (1.1) имеет частные интегралы вида (1.3) и (1.13) с коэффициентами, определенными формулами (1.43) — (1.49), при условии, что коэффициенты системы (1.1) выражаются через параметры по формулам (1.41) — (1.42).
б) /> (1.50)
/> (1.51)
Из (1.50) найдем />:
/>
Из соотношений (1.25) при условиях (1.39) и (1.50) — (1.51) получаем, что коэффициенты системы (1.1) определяются следующими формулами:
/>, /> — любое число, /> (1.52)
/>, />, />, />/>, /> (1.53)
Равенства (1.9) — (1.11) и (1.19) — (1.22) при условии, что имеют место формулы (1.52) — (1.53), дадут следующие выражения для коэффициентов интегралов (1.3) и (1.13):
a1=0 (1.54)
a2/> (1.55)
a/>/> (1.56)
s/> (1.57)
b/> (1.58)
g/> (1.59)
d/> (1.60)
Теорема 1.5 Система (1.1) имеет частные интегралы вида (1.3) и (1.13) с коэффициентами, определенными формулами (1.54) — (1.60), при условии, что коэффициенты системы (1.1) выражаются через параметры по формулам (1.52) — (1.53).
2. Качественное исследование построенных классов систем
2.1 Исследование системы (1.1) с коэффициентами, заданными формулами (1.28) — (1.31)
Будем проводить наше исследование в предположении, что />, />, />.
Пусть мы имеем систему (1.1), коэффициенты которой определяются согласно формулам (1.28) — (1.31), тогда система (1.1) запишется в виде:
/> (2.1)
Интегральные кривые в этом случае имеют вид:
/> (2.2)
/> (2.3)
Найдем состояния равновесия системы (2.1). Приравняв правые части системы нулю и исключив переменную y, получим следующее уравнение для определения абсцисс состояний равновесия:
/> (2.4)
Из (2.4) получаем, что
/>, />, />, />.
Ординаты точек покоя имеют вид:
/>, />, />, />.
Итак, имеем точки
/>, />, />, />.
Исследуем поведение траекторий в окрестностях состояний равновесия />, />, />, />.
Исследуем точку />. продолжение
--PAGE_BREAK--
Составим характеристическое уравнение в точке />.
/>
Отсюда
/>
/> (2.5)
/>
/>
Следовательно, характеристическое уравнение примет вид:
/>=/>=0.
/>,
Или
/>.
Характеристическими числами для точки/> системы (2.1) будут
/>.
Корни /> — действительные, различных знаков не зависимо от параметра d. Следовательно, точка /> — седло.
Исследуем точку
/>.
Составим характеристическое уравнение в точке
/>.
Согласно
равенствам (2.5) характеристическое уравнение примет вид:
/>
/>,
Или
/>.
Характеристическими числами для точки /> системы (2.1) будут
/>,
то есть
/>, />.
Корни /> — действительные и одного знака, зависящие от параметра d. Если d
/> —
неустойчивый узел, если d>0, то точка
/> —
устойчивый узел.
Исследуем точку />.
Применяя равенства (2.5), составим характеристическое уравнение в точке
/>:
/>
/>
Характеристическими числами для точки
/>
системы (2.1) будут
/>,
то есть
/>, />.
Корни /> — действительные и одного знака, зависящие от параметра d. Если d — устойчивый узел, если d>0, то точка /> — неустойчивый узел.
Исследуем точку
/>.
Составим характеристическое уравнение в точке
/>.
Применяя равенства (2.5), получим:
/>
/>,
Или
/>
Характеристическими числами для точки
/>
системы (2.1) будут
/>,
то есть
/>, />.
Корни /> — действительные и различных знаков не зависимо от параметра d. Значит, точка
/> —
седло.
Исследуем бесконечно — удаленную часть плоскости в конце оси oy. Преобразование
/> [7]
переводит систему (2.1) в систему:
/> (2.6)
где />.
Для исследования состояний равновесий на концах оси y, нам необходимо исследовать только точку />. Составим характеристическое уравнение в точке/>.
/>
Получим, что
/>/>
Корни /> — действительные и одного знака. Следовательно, точка /> — устойчивый узел.
Исследуем бесконечно — удаленную часть плоскости вне концов оси oy преобразованием [7] /> Это преобразование систему (2.1) переводит в систему:
/> (2.7)
где />.
Изучим бесконечно — удаленные точки на оси U, то есть при z=0. Имеем:
/>
/>
Получаем, что />. Следовательно, состояний равновесия вне концов оси oy нету.
Теперь дадим распределение состояний равновесия системы (2.1) в виде таблицы 1.
Таблица 1.
d
/>
/>
/>
/>
∞
x=0
(-∞; 0)
седло
неуст. узел
уст. узел
седло
уст. узел
(0; +∞)
седло
уст. узел продолжение
--PAGE_BREAK--
неуст. узел
седло
уст. узел
Положение кривых (2.2), (2.3) и расположение относительно их состояний равновесия при d0 дается соответственно рис.1 (а, б).
Поведение траекторий системы в целом при d0 дается рис.4 (а, б) приложения А: Поведение траекторий системы (2.1).
Исследуя вид кривых (2), (2.3) и расположение относительно их состояний равновесия, убеждаемся, что система (2.1) не имеет предельных циклов, так как Воробьев А.П. [5] доказал, что для систем, правые части которых есть полиномы второй степени, предельный цикл может окружать только точку типа фокуса. Учитывая расположение состояний равновесия относительно кривых (1.3) и (1.13), являющиеся интегралами системы (2.1), характер состояния, заключаем, что для системы (2.1) не может существовать предельных циклов, окружающих несколько состояний равновесия.
/>
а (d
/>
б (d>0)
Рис. 1
2.2 Исследование системы (1.1) с коэффициентами, заданными формулами (1.41) — (1.42)
Будем проводить наше исследование в предположении, что
/>/>/>
Пусть мы имеем систему (1.1), коэффициенты которой определяются формулами (1.41) — (1.42). Тогда система (1.1) будет иметь вид:
/> (2.8)
Интегральные кривые в этом случае имеют вид:
/> (2.9)
/> (2.10)
Частный интеграл (1.13) в этом случае преобразовывается в две прямые (2.10)
1. Найдем состояния равновесия системы (2.8). Для этого приравняем правые части системы нулю
/>
Рассмотрим два случая:
/>
Получаем:
/>
/>
/>
Из первого уравнения найдем y:
/>
и подставляя y во второе уравнение получим:
/>
Решая это уравнение, находим:
/>.
Итак, получаем
/>, />
/>, />
Итак, получаем точки
/>, />, />, />
и прямую x=0, которая является траекторией системы (2.8).
2. Исследуем поведение траекторий в окрестностях состояний равновесия />
Исследуем точку />.
Составим характеристическое уравнение в точке />.
/>
Отсюда
/>
/> (2.11)
/>
/>
Следовательно, характеристическое уравнение примет вид:
/>
Характеристическими числами для точки /> системы (2.8) будут
/>, />.
Корни /> — действительные и различных знаков не зависимо от параметра d, значит точка /> — седло.
Исследуем точку />.
Согласно (2.11) составим характеристическое уравнение в точке />:
/>
Характеристическими числами для точки /> системы (2.8) будут
/>, />.
Корни /> — действительные и одного знака, зависящие от параметра d. Если d — неустойчивый узел, а если d>0, то точка /> — устойчивый узел.
3. Исследуем поведение траекторий в окрестности точки />.
Составим характеристическое уравнение согласно (2.11)
/>.
Характеристическими числами для точки /> системы (2.8) будут
/>, />
Корни /> — действительные и одного знака, зависящие от параметра d. Если d — устойчивый узел, если d>0, то точка /> — неустойчивый узел.
4. Исследуем поведение траекторий в окрестности точки />.
Согласно (2.11) составим характеристическое уравнение:
/>
/>
Характеристическими числами для точки /> системы (2.8) будут
/>, />
Корни /> — действительные и разных знаков не зависимо от параметра d, следовательно /> — седло.
Исследуем бесконечно — удаленную часть плоскости системы (2.8) вне концов оси oy. Преобразование [7] /> переводит систему (2.8) в систему:
/> (2.12)
где />.
Изучим бесконечно — удаленные точки на оси U, то есть при z=0. Получаем:
/>
/> продолжение
--PAGE_BREAK--
Следовательно />.
Таким образом, получаем две точки N1 (0,-1) и N2 (0,1), которые являются состоянием равновесия. Исследуем характер этих точек обычным способом.
Составим характеристическое уравнение в точке N1 (0,-1).
/>
/>
/>(2.13)
/>
Имеем:
/>
/>, />.
Корни />-действительные и различные по знаку, следовательно точка N1 (0,-1) — седло.
Исследуем точку N2 (0,1).
Согласно (2.13) составим характеристическое уравнение:
/>
/>, />.
Корни />-действительные и одного знака, значит точка N2 (0,1) — устойчивый узел.
Исследуем концы оси y с помощью преобразования [7] />. Это преобразование переводит систему (2.8) в систему:
/> (2.14)
где />.
Для исследования состояний равновесия на концах оси y, нам необходимо исследовать только точку N3 (0,0). Составим характеристическое уравнение в точке N3(0,0):
/>
/>
Корни /> — действительные и одного знака, значит точка N3 (0,0) — неустойчивый узел.
Теперь дадим распределение состояний равновесия системы (2.1) в виде таблицы 2.
Таблица 2.
d
/>
/>
/>
/>
∞
N1
N2
N3
(-∞; 0)
седло
неуст. узел
уст. узел
седло
седло
уст. узел
неуст. узел
(0; +∞)
седло
уст. узел
неуст. узел
седло
седло
уст. узел
неуст. узел
Положение кривых (2.9), (2.10) и расположение относительно их состояний равновесия при d0 дается соответственно рис.2 (а, б).
Поведение траекторий системы в целом при d0 дается рис.5 (а, б) приложения Б: Поведение траекторий системы (2.8).
Вопрос о существовании предельных циклов не возникает, так как Воробьев А.П. [5] доказал, для квадратичной системы предельный цикл не может окружать узел.
/>
а (d0)
Рис. 2
2.3 Исследование системы (1.1) с коэффициентами, заданными формулами (1.52) — (1.53)
Будем проводить наше исследование в предположении, что
/>, />/>/>.
Пусть мы имеем систему (1.1), коэффициенты которой определяются формулами (1.52) — (1.53). Тогда система (1.1) будет иметь вид:
/> (2.15)
Интегральные кривые в этом случае имеют вид:
/> (2.16)
/> (2.17)
То есть частные интегралы (1.3) и (1.13) преобразовываются в прямые таким образом, что интегральная кривая (2.16) совпадает с одной из прямых интегральной кривой (2.17).
Найдем состояния равновесия системы (2.15). Приравняв правые части системы нулю, и исключив переменную y, получим следующее уравнение для определения абсцисс состояний равновесия:
/> (2.18)
Из (2.18) получаем, что
/>, />, />.
Ординаты точек покоя имеют вид:
/>, />, />.
Итак, имеем точки
/>, />, />.
Исследуем поведения траекторий в окрестностях состояний равновесия />.
Исследуем состояние равновесия в точке />.
Составим характеристическое уравнение.
/>
Отсюда
/>
/> (2.19)
/>
/>
Следовательно, характеристическое уравнение примет вид
/>
Имеем
/>,
Или
/>.
Характеристическими числами для точки /> для системы (2.15) будут
/>.
Корни /> — комплексные и зависят от параметра d. Значит, если d — устойчивый фокус, если d>0, то точка /> — неустойчивый фокус. продолжение
--PAGE_BREAK--
Исследуем точку
/>.
Согласно (2.19) составим характеристическое уравнение в точке
/>.
Имеем
/>.
/>
Характеристическими числами для точки /> системы (2.15) будут
/>, />
Корни /> — действительные и различных знаков не зависимо от параметра d. Следовательно, точка /> — седло.
3. Исследуем точку />.
По (2.19) составим характеристическое уравнение в точке />.
Получим
/>
/>.
Решая уравнение, получим
/>,
то есть
/>, />
Корни /> — действительные и одного знака, зависящие от параметра d. Если d — неустойчивый узел, если d>0, то точка /> — устойчивый узел.
Исследуем бесконечно — удаленную часть плоскости вне концов оси oy преобразованием [7] /> Это преобразование систему (2.15) переводит в систему:
/> (2.20)
где />.
Изучим бесконечно — удаленные точки на оси u, то есть при z=0. Получаем
/>
/>
/>
Следовательно />
Итак, имеем две точки N1 (0,2) и N2 (0,-2).
Исследуем характер этих точек обычным способом. Составим характеристическое уравнение в точке N1 (0,2).
/>
/>
/> (2.21)
/>.
Следовательно
/>
/>, />
Воспользуемся параллельным переносом
/>/>
и подставим z, u в систему (2.20). Получим новую систему:
/> (2.22)
Составим характеристическое уравнение в точке N2 (0,-2)
/>
Характеристическими числами для точки N2 (0,-2), будут
/>, /> —
сложное состояние равновесия.
Для определения характера состояния равновесия воспользуемся теоремой [2, с. 196-198].
Теорема 2.1. Пусть точка (0,0) — изолированное состояние равновесия системы:
/>(2.23)
где />,/> есть полиномы от x,y начиная со второй степени, /> — решение уравнения />, а разложение функции /> имеет вид:
/>
Тогда
1) при m — нечетном и />m>0 точка (0,0) — есть топологический узел;
при m — нечетном и />m
при m — четном точка (0,0) есть седло — узел, то есть такое состояние равновесия, каноническая окрестность которого состоит из параболистического и двух гиперболических секторов. При этом
если />m
если />m>0, то отрезок отрицательной полуоси OX.
Чтобы воспользоваться теоремой, необходимо систему (2.22) привести к виду:
/>
Это можно сделать, воспользовавшись одним из следующих преобразований [2, с. 199-201]:
если />, />
если />, />, />
если />, />, />
где a, b, c, d — коэффициенты системы (2.23).
Тогда для системы (2.22) возьмем следующее преобразование:
/>
Получим
/>/>
Тогда
/> (2.24)
Найдем решение уравнения:
/>
в виде ряда по степеням Z1:
/>
/>/>
/>/>
/>/>
Следовательно
/>
Тогда
/>
Подставляя U1 в систему (2.24) получим:
/>
Отсюда
/>, />>0. продолжение
--PAGE_BREAK--
Следовательно, по теореме 2.1 получаем, что точка N2 (0,-2) — седло — узел.
Исследуем концы оси y с помощью преобразования [7] />. Это преобразование переводит систему (2.15) в систему:
/> (2.25)
где />.
Для исследования состояний равновесий на концах оси y, нам необходимо исследовать только точку N3 (0,0). Составим характеристическое уравнение в точке N3 (0,0)
/>
Соответственно характеристическими числами будут
/>
Корни /> — действительные и одного знака. Следовательно, точка N3 (0,0) — устойчивый узел.
Теперь дадим распределение состояний равновесия системы (2.1) в виде таблицы 3.
Таблица 3.
d
/>
/>
/>
∞
N1
N2
N3
(-∞; 0)
уст. фокус
седло
неуст. узел
седло
седло-узел
уст. узел
(0; +∞)
неуст. фокус
седло
уст. узел
седло
седло-узел
уст. узел
Положение кривых (2.16), (2.17) и расположение относительно их состояний равновесия при d0 дается соответственно рис.3 (а, б).
Поведение траекторий системы в целом при d0 дается Рис.6 (а, б) приложения В: Поведение траекторий системы (2.15).
Вопрос существования предельных циклов остается открытым.
/>
а (d
/>
б (d>0)
Рис. 3
Заключение
В данной дипломной работе построена квадратичная двумерная стационарная система при условии, что частным интегралом является кривая четвертого порядка, которая распадается на две кривые второго порядка, одна из которых парабола, вторая окружность или гипербола. При этом коэффициенты кривых выражаются через произвольный параметр системы.
Проведено качественное исследование системы. Найдены необходимые и достаточные условия существования у системы двух частных интегралов. В зависимости от коэффициентов были рассмотрены 3 случая. Найдены состояния равновесия трех полученных систем, которые принадлежат интегральным кривым. Исследована бесконечно-удаленная часть плоскости систем, в двух из которых доказано отсутствие предельных циклов. Выяснено поведение сепаратрис седел и построена качественная картина поведения траекторий систем в круге Пуанкаре.
Список использованных источников
Баутин Н.Н. О числе предельных циклов, появляющихся при изменении коэффициентов из состояния равновесия типа фокуса или центра // Матем. сб. — 1952. — Т.30,№1. — 458 с.
Баутин Н.Н., Леонтович Е.А. Методы и приемы качественного исследования динамических систем на плоскости. — М.: Наука, 1976. — 274 с.
Бендиксон И. О кривых, определяемых дифференциальными уравнениями. — УМН, 1941. — Вып.9. — 643 с.
Биркгоф Дж.Д. Динамические системы. М. — Л.: Гостехиздат, 1941. — 340 с.
Воробьев А.П. К вопросу о циклах вокруг особой точки типа “узел" // ДАН БССР. — 1960. — Т.4,№9. — 720 с.
Еругин Н.П. Построение всего множества систем дифференциальных уравнений, имеющих заданную интегральную кривую. — ПММ. — 1952. — Т.16, Вып.6. — с.659-670.
Пуанкаре А. О кривых, определяемых дифференциальными уравнениями. — М. — Л.: ГИТТЛ, 1947. — 839 с.
Серебрякова Н.Н. Качественное исследование одной системы дифференциальных уравнений теории колебаний. — ПММ. — 1963 Т.27, Вып.1. — 230 с.
Филипцов В.Ф. К вопросу алгебраических интегралов одной системы дифференциальных уравнений // Дифференц. уравнения. — 1973. — Т.9,№3. — 256 с.
Черкас Л.А. Об алгебраических решениях уравнения />, где P и Q — многочлены второй степени // ДАН БССР. — 1963. — Т.7,№11. — 950 с.
Яблонский А.И. Алгебраические интегралы одной системы дифференциальных уравнений // Дифференц. уравнения. — 1970. — Т.6,№10. — с.1752-1760.
Приложение А
Поведение траекторий системы (2.1)
/>
а) (d
/>
б) (d>0)
Рис. 4
Приложение Б
Поведение траекторий системы (2.8)
/>
а) (d
/>
б) (d>0)
Рис. 5
Приложение В
Поведение траекторий системы (2.15)
/>
а) (d
/>
б) (d>0)
Рис. 6