Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Методы и приёмы качественного исследования динамических систем на плоскости

..pdf
Скачиваний:
4
Добавлен:
12.11.2023
Размер:
16.72 Mб
Скачать

§ 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 ) и траекторий, стремящихся

Соседние файлы в папке книги