Особые точки дифференциальных систем второго порядка
ГРОДНЕНСКИЙ ГОСУДАРСТВЕННЫЙ
УНИВЕРСИТЕТ
Особые точки дифференциальных систем второго порядка
Беленик Татьяна Васильевна
Гродно,2012
Введение
1. Особые точки
2. Геометрическая классификация особых точек 13
3.Исключительные направления.
Поведение интегральных кривых в нормальной
области
4. Пример
Заключение
Литература
Математический анализ –ветвь
математики, оформившаяся в 18 столетии
и включающая в себя две основные
части: дифференциальное и интегральное
исчисления. В моей работе будем
рассматривать часть
Дифференциальное уравнение — уравнение, связывающее значение некоторой неизвестной функции в некоторой точке и значение её производных различных порядков в той же точке
Исследование особых точек систем дифференциальных уравнений на плоскости представляет актуальную и важную задачу качественной теории дифференциальных уравнений и их приложений в физике, технике, химической кинетики, биологии, медицине, экологии и других областях. Основополагающие результаты о приведении траекторий в окрестности о собой точки были получены А.Пуанкаре, А.И.Ляпуновым, И.Бендиксоном. И.Бендиксоном для аналитических систем был разработан метод расщепления сложной особой точки, позволяющий конечным числом шагов определить топологический тип расположения траекторий (с точностью до различения центра и фокуса).
Изучение особой точки дифференциальных уравнений, т. е. по существу изучение поведения семейств интегральных кривых в окрестности особой точки, составляет один из разделов качественной теории дифференциальных уравнений и играет важную роль в приложениях, в частности в вопросах устойчивости движения .
В моей курсовой характеризуются возможные топологические типы поведения кривых в окрестности изолированной особой точки. Изложение методов расщепления особенностей основывается на использовании нормальных областей. Расскажем о исключительных направлениях, а так же рассмотрим пример решив его с помощью методом полярных координат.
1. Особые точки
Пусть дана система уравнений
тогда x=,y= -ее особая точка, т.е.
X(,)=0, Y()=0 (1)
Предположим, что эта особенность не исчезает при рассмотрении вместо системы одного уравнения
=,
Далее мы будем предполагать, что точка(,) –изолированная особая точка, т.е. вокруг (,) можно описать столь малую окрестность, что внутри нее нет других особых точек.
Поставим следующее проблемы:
1.Характеризовать возможные топологические типы поведения интегральных кривых в окрестности изолированной особой точки.
2.Дать аналитические критерии, которые позволили бы, исходя из функций X(x,y) и Y(x,y), конечным или счетным числом операций определить, какой тип расположения интегральных кривых имеет место для данного уравнения, причем при решении этой задачи, в некоторых случаях, мы будем определять тип расположения интегральных кривых с точностью до инвариантов линейного преобразования.
Для решения этих проблем нам потребуется наложить ряд ограничений на функции X(x,y) и Y(x,y). Именно, считая для простоты x=, y=, предположим, что уравнение (1) имеет вид
= ,
где (x,y) и (x,y) - однородные полиномы степени n и m, а функции (x,y) и (x,y) соответственно суть o и о (), где, , т.е., иными словами, мы требуем, чтобы x=,y= было нулем не ниже первого порядка, а также Y(x,y) и X(x,y) были достаточное число раз дифференцируемыми в окрестности изучаемой особой точки.
Основной для классификации встречающихся случаев служит классификация, предложенная Пуанкаре для случая системы однородных линейных уравнений
=Ax+By;
или одного уравнения
= ,
где
Приведем эту классификацию.
Замена z= приводит нас к уравнению
. ()
Рассмотрим квадратное уравнение
Могут быть четыре случая:
1. Корни этого уравнения действительны и различны.
2. Уравнение имеет кратный корень.
3. D=0, т.е уравнение линейное.
4. Корни комплексны.
Рассмотрим каждый случай в отдельности.
1-й случай. Корни этого уравнения действительны и различны. Тогда они дают два решения дифференциального уравнения ():
или соответствующие два решения уравнения (2):
Перепишем уравнение () в виде
.
Разделяя переменные и интегрируя, получим:
.
:
;
тогда имеем
или
Так как - то окончательно:
Здесь могут представиться два подслучая:
- имеют различные знаки;
- имеют одинаковые знаки.
Введем на плоскости x,y косоугольные координаты:
И рассмотрим семейство интегральных кривых в новой системе координат. Мы имеем:
Или, обозначив через , получим:
Если а), то показатель отрицателен, и в семействе интегральных кривых будет только две кривые, проходящие через начало: Каждая из остальных остается внутри координатных углов.
В этом случае особая точка (0,0) называется седло (рис.1)
Если б), то показатель , и каждая кривая семейства при любом проходит через начало координат. Причем ,если , то все интегральные кривые, кроме прямой , кроме , касаются оси . На плоскости (x,y) это же явление будет наблюдаться не по отношению к осям координат, а по отношению к прямым . Такое расположение интегральных кривых называется узлом (рис.2).
2-ой случай. Уравнение (3) имеет кратный корень . Тогда дифференциальное уравнение () примет вид
x
И общий интеграл получит выражение:
.
Все кривые этого семейства входят в начало координат и там касаются прямой (вырожденный узел) (рис.3).
3-й случай. D=0.
Этот случай распадается на ряд подслучаев.
А) Общий подслучай: С≠0 и В-С≠0.
Исследуемое решение имеет вид
x=P+Qz,
где
P=, Q=.
Отсюда получаем:
x=,
или
.
Из последнего вида общего интеграла имеем:
.
Если Q>0, то при x, т.е. имеем узел. Все интегральные кривые , кроме интегральной кривой х=0, касаются прямой y=-.
Пусть теперь Q<0. Этот случай разбивается на три:
) -1<Q<0; ) Q=-1; ) Q<-1.
)Все интегральные кривые входят в начало; , т.е все интегральные кривые касаются оси у в начале координат (вырожденный узел).
) Q=-1; -+Px. Интегральные прямые параллельны прямой y=+Px. Особая точка – исчезающая (вырожденное седло). (рис.4).
) Q<-1; y=. Единственными интегральными кривыми, входящими в начало, будут: x=0 и y= Все остальные кривые – типа гипербол, т.е имеем седло.
Б) Первый особый случай;
Исследуемое уравнение принимает вид:
x.
Если А≠0, то z=, где Р; отсюда y=.
Все кривые этого семейства проходят через начало и касаются в нем оси y;
Поэтому в рассматриваемом случае имеем узел (вырожденный).
Если А=0, то уравнение превращается в x=0. Отсюда х=0 и z=С, т.е y=Cx. Интегральные кривые входят в начало по всем направлениям (дикритический узел) (рис. 5).
С) Второй особый случай: С=0,В≠0. Из условий С=0, D=0 вытекает, что рассматриваемая система превращается в систему
Или
Т.е. x=const. Все интегральные кривые параллельны оси y. Однако, чтобы выяснить «динамическую » картину, т.е зависимость решения от параметра t, будем интегрировать систему (4), не сводя ее к одному уравнению.
Прежде всего заметим, что прямая Ax+Bx=0 есть особая линия, т.е она сплошь состоит из особых точек.
Приступаем к интеграции системы (4). Пусть x=x(0),y= при t=0, причем A+B≠0, т.е. точка () не является особой.
Имеем:
x= при любом t, и .
Отсюда
=.
Следовательно, при tкоордината y, монотонно изменяясь, приближается к - ,т.е. точка () приближается к точке, лежащей на «особой прямой».
На (рис.6) изображен примерно ход интегральных кривых.
4-й случай. Рассмотрим теперь случай, когда уравнение (3) имеет комплексные корни:
Дифференциальное уравнение () можно записать в этом случае так:
x
Разделяя переменные и интегрируя, получим общий интеграл:
или
Преобразуем плоскость (x,y), положив
и введем на плоскости () полярные координаты:
тогда общий интеграл (2̎ ) примет очень простой вид:
, где k=.
Интегральные кривые уравнения (2̎ ) после преобразования плоскости (x,y) превратились в семейство логарифмических спиралей. Если k>0,то при ,т.е. при неограниченном возрастании числа оборотов по направлению часовой стрелки кривая асимптотически приближается к точке (0,0). Если k<0, то мы имеем то же самое при , т.е. при вращении против часовой стрелки. Если k=0, то решения суть концентрические окружности.
При переходе обратно к плоскости (x,y) картина меняется лишь незначительно: логарифмические спирали остаются кривыми, которые закручиваются вокруг особой точки, асимптотически к ней приближаясь. Особая точка этого рода называется фокусом (рис.7).
Окружности также остаются замкнутыми кривыми, превращаясь в эллипсы. В этом случае особая точка называется центра (рис.8).
п.2.Геаметрическая классификация особых точек.
Пусть точка О - изолированная особая точка. Различаем прежде всего два класса особых точек: устойчивые и неустойчивые.
Особую точку называем устойчивой по Биркгофу, если существуют замкнутые интегральные кривые произвольно малого диаметра, окружающие особую точку, во всех остальных случаях точку называем неустойчивой.
Для того чтобы выяснить дальнейшую структуру точек и другого типа ,докажем несколько общих предложений.
Теорема 1. В достаточно малой окрестности изолированной особой точки не может содержаться замкнутых траекторий, не заключающих внутри себя эту особую точку.
Теорема 1 позволяет провести анализ возможного поведения интегральных кривых в окрестности устойчивой особой точки.
Пусть О – устойчивая особая точка и пусть S(О,r) – столь малый круг, описанный около особой точки, что внутри него и на границе нет других особых точек.
Рассмотрим какой-либо радиус этого круга. Обозначим через М верхнюю границу множества таких точек на этом радиусе, через которые проходят замкнутые интегральные кривые, содержащиеся целиком внутри или на границе круга S(O,r).Через точку М, как это непосредственно вытекает из теоремы о непрерывной зависимости от начальных условий, проходит тоже замкнутая кривая L,содержащая О.Будем рассматривать только интегральные кривые, лежащие внутри кривой L. Отметим на радиусе совокупность точек, лежащих внутри L или на ней самой, через которые проходят замкнутые кривые, очевидно, окружающие точку О. Множество точек этих кривых, если добавить к нему точку О, будет замкнутым множеством, которое мы обозначим через F. Допустим, что существует точка (), лежащая внутри L и не содержащая в F. Движение x=x (t,), y=y (t,), оставаясь все время в кольцеобразной области между двумя замкнутыми кривыми, окружающими точку О, имеем все свои предельные точки как при t,так и при t, на ограничивающих эту область замкнутых кривых и, следовательно, представляет собой спираль, навивающуюся на эти кривые назовем центром-фокусом (рис..9). Следует выделить особо тот случай, когда внутри L заключаются только замкнутые кривые; подобное расположение, как уже упоминалось, Пуанкаре назвал центром.
Переходим к изучению окрестности неустойчивой точки.
Теорема 2. Если изолированная особая точка неустойчива, то всегда существует полутраектория, имеющая эту особую точку единственной своей предельной точкой.
Теорема 3. Около неустойчивой особой точки всегда можно найти столь малую окрестность, что каждая полутраектория будет либо входить в особую точку, либо покидать окрестность через промежуток времени.
Итак, все траектории,
которые могут наблюдаться в
достаточно малой окрестности
неустойчивой особой точки,
- Параболические –одним концом входящие в особую точку, другим –выходящие за границу окрестности.
- Гиперболические или седловые –в обе стороны выходящие за границу окрестности.
- Эллиптические –обоими концами входящие в особую точку.
Соответственно, эти же названия мы присвоим областям, за полненным траекториями того или иного сорта.
Теорема 4. Точки из достаточно малой окрестности изолированной особой точки, лежащие на эллиптических и гиперболических траекториях, если таковые имеются, заполняют множества, содержащие внутренние точки, причем лежащие на гиперболических заполняют области.
Итак, пусть дана особая точка О; опишем около нее малую окрестность S(O,r), удовлетворят условиям теоремы 3.
Рассмотрим точки границы этой области. Выделим сначала множество таких точек, через которые проходят интегральные кривые, хотя бы одним концом входящие в особую точку. Множество этих точек на основании теоремы 4 замкнутое. В этом замкнутом множестве содержится замкнутое множество Ф, состоящие из точек, через которые проходят траектории, обоими своими концами входящие в особую точку. Наконец, через точки смежных интервалов проходят гиперболические кривые. Они заполняют некоторую счетную совокупность областей, расположенных внутри рассматриваемой окрестности. Среди этих областей будут области двух типов: истинные гиперболические области, т.е. такие, которые на своей границе имеют особую точку, и ложные гиперболические области; легко видеть, что истинных гиперболических областей может быть лишь конечное число и что все ложные гиперболические области находятся от особой точки на расстоянии . Топологическая структура гиперболической области весьма проста, именно интегральные кривые, заполняющие эти области, могут быть топологически отображены на семейство отрезков, расположенных на системе параллельных прямых.
Множество Ф состоит из точек, определяющих эллиптические области. Но следует отметить, что, во-первых, не все эллиптические области определяются этими точками, во-вторых, топологическая структура эллиптической области может быть весьма сложной, и эллиптические области могут принадлежать к бесчисленному множеству различных топологических типов(рис.10).
3 .Исключительные направления. Поведение интегральных кривых в нормальной области.
Для проведения дальнейшего анализа окрестности особой точки введем понятие исключительного направления. Для простоты будем предполагать, что исследуемая особая точка находится в начале координат.
Определение . Направление, определяемое полярным углом , назавем исключительным направлением, если существует последовательность точек такая, что и сверх того, если и направлением поля в точке , то .
Если некоторая интегральная кривая входит при или при в особую точку с определенными направлением касательной, то это направление является исключительным. Следовательно, если мы сумеем определить все исключительные направления, то мы сможем найти все направления, по которым интегральные кривые могут входить в начало координат.
Теорема 1.Если на числовой оси нет интервалов, целиком заполненных исключительными направлениями, то каждая полутраектория, содержащаяся в некоторой: достаточно малой окрестности особой точки, либо входит в него с определенной касательной, либо представляет собой спираль, неограниченно приближающуюся к особой точке.
Теорема 2.Пусть есть угол, не заключающий исключительного направления. Тогда можно найти столь малое , что на всякой интегральной кривой, начинающейся в секторе , , переменная меняется монотонно вместе с t и интегральная кривая выходит из этого сектора.
Введем теперь определение нормальной области.
Определение. Некоторый достаточно малый сектор с центром в особой точке и радиуса R называется нормальной областью, если:
1.он заключает одно и только одно исключительное направление;
2.внутри и на границе
направление поля не
Обозначим радиусы, ограничивающие сектор, через .
Будем рассматривать интегральные кривые, пересекающие
В точках пересечения рассмотрим касательные к этим кривым и отметим на них направления, по которым кривые входят в рассматриваемый сектор. Тогда есть три типа нормальной области:
1-й тип. На касательные к интегральным кривым направлены «к центру» (рис.11а).
2-й тип. На касательные к интегральным кривым направлены «от центру» (рис.11b).
3-й тип. На одном из радиусов касательные направлены «к центру», а на другом «от центра» (рис.11c).
Изучим возможные поведение интегральных кривых в каждом из этих случаев.
Радиусы будем называть боковыми сторонами нормальных областей, а дугу окружности ее задней стенкой.
Теорема 3.Если интегральная кривая при возрастании или убывании t входит в нормальную область первого типа, то она при дальнейшем n продолжении войдет в начало координат.
Теорема 4. Если интегральная кривая входит в нормальную область второго типа через боковую сторону, то она при дальнейшем продолжении выходит из этой области; однако на задней стенке существует точка или целая дуга, состоящая из точек, через которые проходят интегральные кривые, входящие в начало координат при или при
Итак, второй тип приводит к возможности существования двух случаев, которые мы назовем случай единственности и не единственности.
Теорема 5. Если нормальная область принадлежит третьему типу (например, на входящие кривые направлены от центра, а на к центру), то либо все интегральные кривые, входящие в область, из нее выходят, либо на отрезке можно будет найти отрезок OS, через все точки которого проходят интегральные кривые, входящие в начало координат.
Итак, третий тип области тоже приводит к двум случаям, которые называть будем случаем отсутствия входящей кривой и случаем не единственности.
Пусть дано уравнение
(4.1)
и пусть особая точка, является изолированной, т.е начало координат О(0,0).
Представим
Данное уравнение решим с помощью полярных координат :
x=rcos
y=rsin
Перейдем к полярным координатам r и .
Получаем
имеем
r (4.2)
Характеристическое уравнение
Исключительные направления находятся из tg
Отсюда имеем три корня: .
Для простого корня уравнения =0 и причем имеем . Следует
Теорема: Если в (4.2)
то существует единственная кривая системы (4.1), касающаяся в точке О луча .
Значит , луч можно заключить в нормальную область второго типа, из нашей теоремы вытекает существование единственной кривой, касающейся в начале координат луча . Аналогично образом доказывается существование единственной кривой, касающейся в точке О луча ,
Для обозначим угол через t: откуда
и с помощью тригонометрических преобразований получим:
найдем теперь
Подставляя в наше получим
(Аsint+Bcost)(Acost-Bsint)(
=(*
*((3A+2B)sint+(3B-2A)cost=
=(3A+2B)cost+AB(3A+2B)tsint+…+
=6
Посчитав получим 6
Отсюда АВ(2А+2В)≠0 (находится в первой четверти).
Аналогично для . Подставив и посчитав получим
-108-24+54≠0 следовательно АВ(2А+2В)≠0 (находится в третей четверти ).
Фазовый портрет системы (4.1) в окрестности особой точки О имеет вид седла
ЗАКЛЮЧЕНИЕ
Особая точка
где Р и Q - непрерывно дифференцируемые функции.
Предполагая, что особая точка расположенна в начале координат и используя Тейлора формулу, можно представить уравнение выше в виде:
где P(x, у) и Q(x, у)- бесконечно
малые по отношению к .Характер поведения
интегральных кривых около особой точки зависит от корней l1 и
l2 характеристического уравнения.
В окрестности особой точки система может обладать очень сложной динамикой. В зависимости от собственных значений матрицы, различают четыре типа невырожденных особых точек линейных систем: узел, седло, фокус, центр. Кроме основных типов появляются так называемые сложные: седло-узел, центро-фокус и т.д. Они образуются как результат слияния точек покоя основных типов.
Мы имеем классы устойчивых и неустойчивых точек. Ясно, что наличие исключительных направлений, заключается в одну из областей 1 или 2, исключает возможность устойчивой особой точки. Но вообще наличие исключительного направления не препятствует устойчивости точки, даже наличию центра. Отсутствие исключительных направлений приводит, вообще говоря, к двум возможным случаям: или особая точка есть фокус, либо устойчивая особая точка. Если же исключительные направления имеются в конечном числе, тогда легко видеть, когда могут возникать в окрестности начала координат эллиптические, гиперболические и параболические области.
Параболическая возникает при областях: типа1,2,3.(рис. а, б, в)
Гиперболическая – в том случае, если, рассматривая два средних исключительных направления, мы убеждаемся, что оба они могут быть заключены в область типа 3 (или 3)(рис.г).
Эллиптические получаются в том случае, если два соседних исключительных направления могут быть заключены в область типа 1(рис. д).
Наконец, если правая часть характеристического уравнения тождественно обращается в 0, то только по конечному числу направлений в особую точку могут входить интегральные кривые.
Литература