книги / Методы и приёмы качественного исследования динамических систем на плоскости
..pdf§ 8] |
СВОДКА СВЕДЕНИЙ |
81 |
Рис. 44 |
Рис. 45 |
\У
Рис. 48 |
Рис.' 49 |
в Н. Н. Баутин, Е. А. Леонтович
82 ИССЛЕДОВАНИЕ СОСТОЯНИЯ РАВНОВЕСИЯ [ГЛ. 3
Ниже перечислены все возможные типы состояний равновесия
с не равными нулю действительными |
частями |
характеристиче |
|||||||||
ских корней и приведены схематические рисунки |
расположения |
||||||||||
траекторий в их окрестности. |
|
|
|
|
|
|
|
|
|||
При этом для недикритического узла и |
для |
седла рисунки |
|||||||||
приводятся как в случае, когда рассматриваемая |
система имеет |
||||||||||
канонический вид (когда направления, по |
которым к состоянию |
||||||||||
равновесия |
стремятся |
траектории, совпадают |
с |
направлением |
|||||||
осей координат), так и в общем |
случае (т. е. |
в |
случае, |
когда |
|||||||
система не имеет канонического вида, так что направления |
к\ |
и |
|||||||||
& 2 могут быть любыми). |
|
|
|
|
|
|
|
|
|||
I. Узел |
(характеристические корни |
Ал и Аг |
действительны |
и |
|||||||
одинаковых знаков, т. е. А1 А2 |
> 0 ; |
Д > |
0, о2 |
— 4Д > 0). |
|
|
|||||
А. Невырожденный узел |
(Ал Ф Аг): |
|
|
|
|
|
|
||||
а) |
устойчивый: Ал < |
0, Аг<0, |
т. е. о < 0 |
(рис. 41, 44); |
|
|
|||||
б) |
неустойчивый: Ал > 0, Аг > |
0, т. е. о > |
0. |
|
|
|
|
(Рисунки даются только для устойчивого узла, в случае неус тойчивого узла надо переменить направление стрелок. При этом рис. 41 соответствует случаю, когда си стема имеет канонический вид, а
рис. 44 соответствует общему виду.)
Б. Вырожденный узел (Ал = Аг = А,, но в канонической форме р Ф 0, Д > 0,
02 _ |
4Д = 0): |
а) |
устойчивый: А < 0 (рис. 43, |
^4 5 );
хб) неустойчивый: А > 0 .
|
(Рис. 43 соответствует случаю, ког |
|||
да система имеет канонический вид, а |
||||
рис. 45 соответствует общему виду.) |
||||
|
В. Дикритический |
узел |
(Ai = Аг =А |
|
и |
р = |
0 ): |
0 (рис. 42); |
|
|
а) |
устойчивый: А < |
||
|
б) |
неустойчивый: А > 0 . |
|
|
II. Седло (характеристические корни Ai и Аг действительны и |
||||
разных знаков, т. е. AiАг < |
0, |
либо Ai > 0 , Аг < 0, |
либо Ai < 0, |
Аг>0, Д < 0 ) изображено на рис. 46 и 47 |
(рис. 46 соответствует |
|||
случаю, системы |
в |
каноническом виде |
при Ai > 0 , |
А г< 0 , |
рис. 47 — общему случаю). |
|
|
||
III. Фокус (характеристические корни комплексные сопря |
||||
женные, т. е. Д > |
0, а2— 4Д < 0): |
|
|
|
а) устойчивый: а < |
0 (а < 0) (см. рис. 48 и 49 для |
устойчи |
вого фокуса в случае канонического вида системы: рис. 48 соот ветствует случаю Р > 0, рис. 49 — случаю р < 0; рис. 50 соответ ствует случаю, когда система имеет общий вид);
б) неустойчивый: а > 0 ( а > 0 ).
§ 8] |
СВОДКА СВЕДЕНИЙ |
83 |
П р и м е р 1. |
|
|
х = у, у = х(а2 — х2)+Ъу. |
|
|
Состояние равновесия |
0(0, 0)— седло. Определим |
направления |
сепаратрис в седле. Уравнение для нахождения углового коэф фициента сепаратрис в седле имеет вид к2— Ък — а2= 0 , откуда
Ли = 6 / 2 ± Гб2/4 + а2.
П р и м е р 2. |
|
|
|
|
х = —х (2 + у), у = х + $у. |
|
|
Состояние равновесия 0(0, 0)— седло. Уравнение для |
определе |
||
ния направлений |
сепаратрис в седле: |
(Р + 2 )Л + 1 = |
0 , откуда |
к = —1/(2 + Р). |
Нетрудно видеть, что |
второе значение |
к есть <», |
а сепаратриса с наклоном к = <х>есть прямая х = 0 .
Г Л А В А 4
КАЧЕСТВЕННАЯ СТРУКТУРА ОКРЕСТНОСТЕЙ НЕКОТОРЫХ СЛОЖНЫХ СОСТОЯНИЙ РАВНОВЕСИЯ
§ 1. Направления, в которых траектории стремятся к слож ному состоянию равновесия. При исследовании сложных состо яний равновесия иногда бывает весьма существенно знание на правлений, в которых траектории могут стремиться к этому со стоянию равновесия 1).
Рассмотрим динамическую систему
dx/dt = P (х, у), dy/dt = Q(x, у),
для которой начало координат является сложным состоянием равновесия, так что
Р (0 , 0 ) = 0 , m |
о )= о , |
|
|||
0 |
0 |
) |
Ру (0,0) |
|
|
д __ К ( , |
|
|
= 0. |
(1) |
|
Q'x (0 , 0 ) |
|
||||
Qy(0, 0 ) |
|
Предположим, что разложения правых частей в ряд по степеням х, у ъ окрестности точки 0 (0 , 0 ) имеют вид
Р(х, у) = Рт(х, у)+<р(х, у),
Q(x, y)=Q m(x, y)+<b(x, у),
где т > 1, Рт(х, у) и Qm(x, у ) — однородные многочлены, состоя щие из всех членов т-то порядка соответствующих разложений, а функции ф(а;, у) и ф(я, у )— ряды, состоящие из членов более высоких порядков. При этом мы считаем, что многочлены Рт(х, у) и Qm(x, у) одновременно не равны тождественно нулю (в противном случае мы бы взяли т ' > т).
Рассмотрим выражение
xQm(x, у )— уРт(х, у), |
(2 ) |
а также выражение
k ) - k P m( 1, к ) ,
') Методы исследования сложных особых точек (метод Бендиксона [143] и метод Фроммера [132]) опираются на рассмотрение траекторий, стремящихся в определенном направлении.
s 1] |
НАПРАВЛЕНИЯ СТРЕМЛЕНИЯ ТРАЕКТОРИЙ |
85 |
которое мы получим из предыдущего, если поделим его на |
1 |
|
и введем обозначение |
|
|
|
к = у/х. |
|
Мы будем также рассматривать выражение
cos 0 Qm(cos 0, sin 0) — sin 0 Pm(cos 0, sin 0 ).
Имеет место Т е о р е м а 1. Всякая полутраектория системы
dx/dt = Pm(x, у)+ <p(z, у), dy/dt = Qm{x, у)+ it>(a;, у)
(Pm и Qm не равны тождественно нулю), стремящаяся к состоя
нию равновесия |
О (0 , 0 ), либо является спиралью, |
стремящейся |
к О при t +°° |
(или t->----°°), либо стремится к |
О в опреде |
ленном направлении 0 *.. |
|
При этом:
I. Если хоть одна из траекторий системы является спиралью, стремящейся к О при t +°° (или t —°°), то все траектории, проходящие через точки некоторой окрестности состояния равно
весия О, являются такими же спиралями |
(т. е. точка есть устой |
||||
чивый или неустойчивый «фокус высшей сложности»). |
|||||
II. Если выражение (2) |
не обращается тождественно в нуль, |
||||
то наклоны к, с которыми |
траектории |
стремятся к |
состоянию |
||
равновесия О, удовлетворяют уравнению |
|
(3) |
|||
Qm(i, |
k * ) - k * P m(i, к*)= 0, |
||||
или, иначе, направления |
0 *, с |
которыми траектории |
стремятся |
||
к О, удовлетворяют уравнению |
|
|
|
||
cos 0* Qm(cos 0*, sin 0*)— sin 0* Pm(cos 0*, sin 0*) = O. (4) |
|||||
III. Если |
|
|
|
|
|
xQm(x, y )— yPm(X, y)= 0 |
|
||||
и, следовательно, |
|
|
|
|
|
P m ( x , y) = xQm-i (x, y), |
Qm(x, y)= yQm-\(x, |
y), |
где Qm-i(x, y ) — некоторый не равный нулю тождественно одно родный многочлен степени тп— 1 , то, какое бы направление 0 , не удовлетворяющее уравнению
Qm-1 (cos 0*, sin 0 *) = 0 , |
(5) |
мы ни взяли, существует в точности одна полутраектория, стре мящаяся к О в направлении 0. Для особого же направления 0*, удовлетворяющего (5), может оказаться, что не существует ни одной полутраектории, стремящейся к О в этом направлении 0 *, либо есть конечное число таких траекторий, либо, наконец, таких траекторий может существовать бесчисленное множество.
З а м е ч а н и е . Если существует траектории, стремящаяся к состоянию равновесия О с определенным наклоном к*, то этот
СТРУКТУРА СЛОЖНЫХ СОСТОЯНИЙ РАВНОВЕСИЯ |
[ГЛ. 4 |
наклон, согласно сформулированной теореме 1 , является действи тельным корнем уравнения Qm(1, к )— кРт(1, к) = 0. Однако если это уравнение имеет действительные корни, то это еще не означает, что существуют траектории, стремящиеся к О с этим наклоном: возможны случаи, когда при этом все траектории яв ляются спиралями или замкнутыми траекториями.
§ 2. Сложное состояние равновесия (особая точка) с нуле выми характеристическими корнями2). В настоящем параграфе мы приведем результаты исследования одного простейшего типа сложных особых точек.
Рассмотрим систему
dx/dt = Р(х, у), dy/dt = Q(x, у),
где Р(х, у) и Q(x, у) — аналитические функции, не имеющие об щего множителя, отличного от постоянного.
Пусть начало координат является сложным состоянием равно весия этой системы, т. е. мы имеем
Р( 0 , 0 ) = 0 , Q(0 , 0 ) = 0 ,
д = р;(о,о) |
(о, о) |
(о, о) |
<?; (о, о) |
и, следовательно, хотя бы один из характеристических корней этого состояния равновесия равен нулю. Мы будем рассматривать здесь такие сложные состояния равновесия, когда в разложениях по степеням х жу функций Р(х, у) и Q(x, у) хотя бы один из линейных членов не равен нулю, т. е. когда
1^(0,0)1 + |р;(0,0)1 + |&(0,0)1 + |&(0,0)М 0.
Рассмотрим наряду с величиной Д(0, 0) величину
о = Р'х(0,0) + Qy(0,0).
Среди состояний равновесия, для которых выполняется условие (6 ), естественным образом выделяются два случая в зависимости от того, что имеет место: а(0, 0)Ф 0 или а(0, 0) = 0. Так как характеристическое уравнение имеет вид
Я2 — сгЯ + А = 0,
2) Качественный характер таких состояний равновесия был рассмот рен методом Бендиксона (см. 0143, 60, 70]) и методом Фроммера (см. [1, 132, 133]). Мы приводим здесь величины, определяющие характер этих со стояний Ьавновесия, однако способом, указанным в работе [70] (см. также монографии [12, 13]), так как этот способ их введения значительно более естествен при рассмотрении бифуркаций этих состояний равновесия в гл. 10.
§ 2] |
СОСТОЯНИЕ РАВНОВЕСИЯ С НУЛЕВЫМИ КОРНЯМИ |
87 |
|
то, очевидно, в случае, когда о |
0 , только один характеристиче |
||
ский корень равен нулю, второй |
же равен о. В случае, |
когда |
о= 0 , оба характеристических корня равны нулю.
I. А(0, 0) = 0, о (0, 0 )^ 0 . В этом случае существует неособое линейное преобразование (см. § 2 гл. 3), с помощью которого
система в окрестности начала 0 (0 , 0 ) может |
быть представлена |
||||||||||||
в следующем каноническом виде: |
|
|
|
|
|
|
|
||||||
|
|
dx/dt = Р*(х, у), |
dyldt = by + Q*(x, у), |
|
|
(7) |
|||||||
где Ь Ф 0, а разложения по степеням х, у |
функций |
Р* (х, у) и |
|||||||||||
Q*(x, у) |
начинаются с членов не менее чем второго порядка. |
|
|||||||||||
Введем в рассмотрение функцию |
|
|
|
|
|
|
|
||||||
|
|
|
|
у = ц>(х), |
|
|
|
|
|
|
|
||
являющуюся решением уравнения3) |
|
|
|
|
|
|
|
||||||
|
|
|
|
by + Q *{x,y) = 0 . |
|
|
|
|
|
|
|
||
Подставим функцию у = у{х) в Р* (х, |
у) |
и |
введем |
обозна |
|||||||||
чение |
|
|
|
ф(.х)= Р * {х, ф(х)). |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|||
Так как мы предположили, |
что функции |
Р(х, |
у) |
и |
Q(x, |
у) |
|||||||
не имеют |
общего |
множителя, |
отличного |
от |
постоянного, |
то |
|||||||
Р* (х, ср (х )) не |
может быть |
тождественно равна |
нулю и, следо |
||||||||||
вательно, |
в разложении функции ф(я) по степеням |
х |
заведомо |
||||||||||
будут отличные от нуля члены. |
|
|
|
|
|
|
|
|
|||||
Таким образом, мы можем написать |
|
|
|
|
|
|
|
||||||
|
|
|
ф {х) = Р* {х, ф{х)) = Атхт+ ..., |
|
|
|
|
|
|||||
где т 3* 2 |
(так как разложение Р*(х, у) по степеням х |
и у |
на |
||||||||||
чинается |
с членов не ниже второй степени и Ат ¥= 0). Число т, |
||||||||||||
очевидно, |
х а р а к т е р и з у е т |
к р а т н о с т ь |
о б щ е й |
т о ч к и 4) |
|||||||||
кривых |
|
Р * (х ,у )= 0 |
и |
by + Q*(x, у)= 0 . |
|
|
|
|
|||||
|
|
|
|
|
|
||||||||
Т е о р е м а |
2. |
Состояние |
равновесия 0(0, |
0), |
для |
которого |
|||||||
А(0, 0 ) = 0 |
и о(0, |
0 )^ 0 , может иметь следующий качественный |
|||||||||||
характер: |
|
|
|
|
|
|
0 ); |
|
|
|
|
||
а) характер седла (при тп нечетном и Дт > |
|
|
|
|
|||||||||
3) Такая |
функция всегда существует, так как |
для |
уравнения by + |
||||||||||
+ Q*(х, у) = 0 выполняются условия существования неявной функции. |
|||||||||||||
Действительно, при х = 0, у = 0 будет |
|
|
|
|
|
|
|
||||||
|
|
by + Q* (х, *)'- 0, |
-щ- (by + Q* (х, у)) = ЬФ 0. |
|
|
|
4) Кроме того, Дт является ляпуиовскоп величиной, соответствующей одному, равному нулю характеристическому корню (см. гл. 6 ).
88 |
СТРУКТУРА СЛОЖНЫХ СОСТОЯНИЙ РАВНОВЕСИЯ |
[ГЛ. 4 |
б) |
характер узла (при тп нечетном и Дт < 0 ) ; причем при |
|
Ъ< 0 |
узел устойчивый, а при Ь > 0 — неустойчивый; |
|
в) |
состояние равновесия с одним узловым сектором и двумя |
|
седловыми (при тп четном и любом знаке Дт ). При Ъ< |
0 узловой |
|
сектор устойчивый, при Ъ> 0 — неустойчивый. Кроме |
того, если |
ЬДт < 0 , то траектории узлового сектора стремятся к О (при t -*■ +оо или при t —оо в зависимости от знака Ь) слева от оси у
(рис. 51, а), а если ЬДт > 0 , то справа от оси у (рис. 51, б)5). Состояние равновесия в случае а) мы будем называть слож
ным седлом, в случае в)— седло-узлом, а в случае б)— сложным узлом6).
Нетрудно видеть (см. §1), что, когда система имеет канони ческий вид (7), существует два направления, по которым траек тории могут стремиться к рассматриваемому состоянию равнове сия, это: 0, я и я/2, Зя/2.
На рис. 51, а, б представлен седло-узел в случае, когда систе ма приведена к каноническому виду (7).
Очевидно, в случае, когда в рассматриваемых координатах х жу система не имеет канонического вида, направления, в кото рых траектории стремятся к началу координат, могут быть от личны от направления осей. Такой случай представлен для слу чая седло-узла на рис. 51, в.
В дальнейшем особый интерес для нас будет представлять случай тп = 2 , который мы будем называть случаем простейшего двукратного седло-узла. В случае тп > 2 будем называть состоя ние равновесия сложным седло-узлом.
5)Хотя на рис. 51 оси координат нарисованы, но они не относятся к качественной структуре состояния равновесия и не надо придавать значе ния деталям взаимного расположения траектории и осей.
6)Геометрически в рассматриваемом случае сложное седло и сложный узел ничем не отличаются от простого седла и узла.
§5) |
СОСТОЯНИЕ РАВНОВЕСИЯ С НУЛЕВЫМИ КОРНЯМИ |
89 |
II. Д(0, 0 ) = 0 , о(0, 0) = 0. В этом случае, очевидно, оба ха рактеристических корня состояния равновесия равны нулю.
В рассматриваемом случае система линейным неособенным преобразованием приводится к виду
dx/dt = у + P t (х , у), dy/dt = Q\ (х, у),
где Р * (х , у) и Q* (х, у) — аналитические функции, разложения ко торых по степеням х ж у начинаются с членов не менее чем второго порядка.
Рассмотрим следующие функции:
1 ) функцию у = у(х), являющуюся решением уравнения
у + Р * (х,у) = 0;
2 ) функцию у = $ (х), определяемую формулой
■ф(ж) = ^*(ж, ф(ж));
эта функция заведомо не равна нулю тождественно (в силу предположения об отсутствии общих множителей, отличных от постоянных у правых частей рассматриваемой системы), поэтому в разложении -ф(л:) по степеням х заведомо будут отличные от нуля члены, и мы можем написать
■ф(ж)= Q*(x, <р(ж))= ahxk + ...; акФ 0;
3) функцию
|
а (х) = Р*х (х, ф (х)) + Qly (х, ф (х)) = |
Ъпхп + . . . |
||
Функция |
а (я), |
в отличие от ф(ж), |
может |
тождественно обра |
щаться в |
нуль. |
Рассмотрим сначала |
случай, |
когда о(х)¥=0, так |
что при некотором п Ъп Ф 0. Числа к, п и коэффициенты ак и Ьп характеризуют качественную структуру особой точки. При этом число к характеризует кратность общей точки изоклин (см. ч. II).
Имеют место следующие теоремы. |
|
|
|||
Т е о р е м а |
3. Пусть к четное, к = 2т. |
|
|||
Тогда: |
1) |
в случае, если п < т, особая точка 0(0, 0) имеет |
|||
качественный характер седло-узла |
(рис. 52); |
полутраектория, |
|||
2 ) |
в случае, когда |
т, |
существует одна |
||
стремящаяся к О при t -*■ — |
и одна полутраектория, стремящая |
||||
ся к О при t -*■ +°°, все остальные траектории и при возрастании, |
|||||
и при убывании t выходят из |
окрестности О (г. е. |
окрестность |
|||
особой точки О состоит из двух седловых секторов). |
|
Такое состояние равновесия мы будем называть вырожденным седло-узлом (рис. 53).
Отметим, что рис. 53 выполнен при условии Ьп > 0 и а%т< 0, а рис. 52 соответствует случаю Ьп > 0, а^т < 0.
Т е о р е м а 4. Пусть k = 2m + 1 — нечетное число и агт + 1 Ф 0, и пусть у = bn + 4(m + 1 ) а2т+1.
90 |
СТРУКТУРА СЛОЖНЫХ СОСТОЯНИЙ РАВНОВЕСИЯ |
|
[ГЛ. 4 |
||||
Тогда: 1 ) если а2т+\ > |
0, |
то особая точка О (0, 0) |
имеет каче |
||||
ственный характер седла |
(рис. 54); |
|
|
|
|||
2 ) |
если а2т+\ < 0 , то особая точка имеет: |
|
|
п = |
|||
а) |
характер фокуса или центра при п > т, а также при |
||||||
= т и •у < 0 [92]; |
|
п |
четное и при |
этом |
п < т |
или |
|
б) |
характер узла, если |
||||||
п = т и у 3* 0 ; |
|
|
(эллиптическую) область, две со |
||||
в) |
одну замкнутую узловую |
||||||
провождающие ее узловые |
области и одну |
седловую область |
|
|
|
|
|
|
|
Рис. 54 |
|
|
|
(рис. 55), если п — нечетное число и при этом п < тп или |
га = |
m |
||||||||
и у > 0. |
|
|
|
|
0; в случае |
Ьп < О |
||||
Рис. |
54 и 55 выполнены при условии Ъп > |
|||||||||
|
|
|
расположение |
|
траекторий |
получается |
||||
|
|
|
симметричным |
отображением |
относитель |
|||||
|
|
|
но оси х. |
|
|
|
|
|
|
|
|
|
|
Нетрудно видеть, что в случае рас |
|||||||
|
|
|
сматриваемых |
состояний равновесия урав |
||||||
|
|
|
нение для определения направления, по |
|||||||
|
|
|
которому траектории стремятся к состоя |
|||||||
|
|
|
нию равновесия: |
|
|
|
|
|
||
|
|
|
bk2 + (a — d)k — с = 0, |
(8 ) |
||||||
|
|
|
имеет двукратный нулевой корень (так |
|||||||
|
|
|
как мы имеем |
а = d = с = О, Ъ Ф 0). Все |
||||||
|
Рис |
|
стремящиеся |
к |
состоянию |
равновесия |
с |
|||
|
5 5 |
о п р е д е л е н н ы м |
н а п р а в л е н и е м |
|||||||
оси х |
(см. |
рис. |
траектории стремятся |
к |
нему, касаясь |
|||||
52—55). Однако |
если состояние |
равновесия |
есть фокус или центр, то имеем, очевидно, случай, возможность которого была указана: когда, несмотря на наличие дейст вительных корней уравнения (8 ) и траекторий, стремящихся