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

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

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

S 3)

БИФУРКАЦИИ НЕКОТОРЫХ ТИПОВ СЛОЖНЫХ ОСОБЫХ ТОЧЕК

171

Точнее:

 

 

 

 

 

 

 

а)

существуют е<> > 0, бо >

О такие, что в

ео-окрестности сепаратрисы L a

седло-узла

0 (х0, уо) системы

(А), идущ ей из

седло-узла в него же,

у всех

бо-близких

к системе

(А) систем

(А) сущ ествует не более

одного предель­

ного цикла;

 

 

 

 

 

 

б)

при

любом е <

е0 можно

указать б <

бо такое, что

у всякой

систе­

мы (А), б-близкой к системе (А) и такой, у которой в ео-окрестности седлоузла О нет состояний равновесия, сущ ествует предельный цикл, целиком леж ащ ий в е-окрестности сепаратрисы L 0, и притом этот предельный цикл

устойчив, если Рх (*0, yQ) + Qy (*0, У0) < 0, и неустойчив, если Рх (х0, у0) +

+<?д(жо’г'о)> 0 *

§3. Бифуркации некоторых типов сложных особых точек. Пусть для рассматриваемой системы

dx/dt = Р(х, у), dyldt = (?(х,

у)

точка 0 (хо, уо) является

(изолированной) сложной особой точ­

кой, т. е. особой точкой, для которой

 

A(*o.J/o)

К ( Х0'Уо)

К ( хо'Уо)

= 0.

 

Qm(xo^o)

<?i(Wo)

 

В этом случае всегда можно так изменить правые части этой системы, чтобы «сложная точка О» распалась на несколько осо­ бых точек, т. е. чтобы у соответствующей измененной системы

dx/dt = Р{х, у) = Р(х, у) + р(х, у), dy/dt = Q(x, y) = Q(x, y)+q(x, у)

в достаточно малой окрестности точки О было несколько осо­

бых точек. Поэтому

естественно ввести

понятие к р а т н о с т и

особой точки.

Особая

точка 0( х о,

уо) системы (А) на­

О п р е д е л е н и е .

зывается т-кратной,

если; а)

существуют

такие Ео > 0 и бо > 0,

что при всевозможных бо-добавках ранга тп у соответствующей системы (А) в Ео-окрестности О может быть не более чем тп особых точек; б) при любых б < б0 и е < е0 всегда существуют такие б-добавки ранга тп, при которых у соответствующей си­

стемы (А) в е-окрестности точки

О существует тп грубых осо­

бых точек3) .

 

 

 

 

3) Если бы среди m особых точек, на которые при сколь угодно малых

добавках разделяется данная

слож ная

особая

точка О, была бы

слож ная

особая точка (т. е. такая, для

которой

Д = 0),

то всегда можно

было бы

указать такие сколь угодно малые добавки, при которых особая точка О разделялась бы более чем на m особых точек, что противоречит условию а). Если бы среди особых точек, на которые разделяется О, была бы простая, но не грубая особая точка (т. е. такая, для которой Д > 0, а = 0), то всег­ да можно было бы указать и такие сколь угодно малые добавки, при кото­ рых все особые точки были бы грубыми. Кратность сложного состояния равновесия, очевидно, совпадает с кратностью общей точки любых двух различных изоклин, проходящих через эту особую точку.

172

БИФУРКАЦИИ ДИНАМИЧЕСКОЙ СИСТЕМЫ

[ГЛ. 10

Очевидно, особая точка кратности т может при надлежа­ щем выборе добавок разделяться и на меньшее, чем пг, число особых точек, и, в частности, могут быть такие сложные особые точки, которые при надлежащим образом выбранных, но сколь угодно малых добавках исчезают (например, седло-узел).

Т е о р е м а 1. Индекс сложной особой точки О кратности тп равен сумме индексов тех особых точек, на которые сложная особая точка О может разделяться при сколь угодно малых добавках ранга пг.

I.Разделение при малых добавках к правым частям системы

(А)сложной особой точки О с одним нулевым характеристиче­ ским корнем (т. е. особой точки, для которой А = 0, оч^О) на грубые особые точки. Такая сложная особая точка была рас­ смотрена в гл. 4. В окрестности такой особой точки, как мы ви­ дели (см. § 2 гл. 4), система может быть приведена линейным неособым преобразованием к виду

dx/dt = P(x, у) = Р2(х, у),

 

 

 

(А)

 

 

dyldt = Q (х, у) = Q2(х , у ) + Ь у ,

где Рг(х,

у),

Qz{x, у) — функции,

разложения которых по

степеням

х, у

начинаются с членов

не ниже второй степени,

и6 ^ 0 .

Вгл. 4 были введены в рассмотрение следующие функции:

а) у — ф(я), являющаяся решением уравнения

by + Q2(x, у) = 0;

б) z = i|)(ж) = Р2(х, tp(x))=Amxm+ . . . (не выписанные члены содержат х в степени выше пг), где Д„,=И=0 и тЗ * 2.

Как мы видели (см. гл. 4), четность и нечетность пг и знак Ат и Ь определяли характер рассматриваемой особой точки.

Наряду с системой (А) будем рассматривать измененную систему

 

dx/dt = Р(х,

у) + р(х,

у) *=*Р(х,

у),

(A)

 

dy/dt = Q(x,

y)+q(x,

y) = Q(x,

у).

 

 

Т е о р е м а

2.

Если

Am Ф 0, то

кратность особой точки О си­

стемы (А) есть пг4).

 

 

 

 

 

 

Т е о р е м а

3.

Если

кратность

особой точки О системы

(А)

есть пг (пг 3= 2),

то число

грубых

состояний

равновесия

(0 \,

О2, .. ., Oh) системы (А), на которые точка О может разделить­ ся при сколь угодно малых добавках р(х, у) и q(x, у) ранга пг, при тп нечетном может быть любым нечетным числом, меньшим

4)

Очевидно, m является такж е кратностью общей точки О изоклин си­

стемы

(А).

 

БИФУРКАЦИИ НЕКОТОРЫХ ТИПОВ СЛОЖНЫХ ОСОБЫХ ТОЧЕК

173

или

равным т, при т четном любым четным

числом,

мень­

шим или равным тп\ при этом-.

узла (т. е. m

нечетно и Дт < 0 ) ,

 

1) если О имеет характер

то число грубых узлов среди

особых точек Оi, . . Ok системы

(А)

на единицу больше числа грубых седел-,

нечетно и Дт > 0 ) ,

 

2) если О имеет характер

седла (т. е. m

то число грубых седел среди

особых точек

0 \, .

Ок системы

(А)на единицу больше числа грубых узлов-,

3)если О имеет характер седло-узла (т. е. m четно), то

среди особых точек Oi, ..., Ok системы (А) число узлов равно числу седел.

З а м е ч а н и е 1. В силу того, что в рассматриваемой особой точке а Ф О, ни при каких (достаточно малых) добавках среди грубых особых точек 0\, ..., Ок системы (А) не может быть фокусов и не может быть также предельных циклов, рождаю­ щихся из состояния равновесия О (это элементарно устанавли­ вается использованием критерия Дюлака).

З а м е ч а н и е 2. Из настоящей теоремы, очевидно, следует, что индекс состояния равновесия О(0, 0) системы (А), имею­ щего характер узла, равен +1, имеющего характер седла, ра­ вен — 1, и имеющего характер седло-узла, равен 0.

В случае системы (А) изоклины, имеющие наклон, отличный от нуля, не имеют особенности в точке О 5).

Нетрудно видеть, что в рассматриваемом простейшем случае сложной особой точки ее качественный характер полностью оп­ ределяется числом и характером грубых состояний равновесия, на которые она разделяется.

II.Разделение при малых добавках к правым частям систе­

мы (А) сложной особой точки О (0, 0) с двумя нулевыми кор­

нями (Д = 0,

о¥=0), для

которой | Рж(0, 0) | + | Ру (0, 0) | +

+ | Qx (0) 0) [ +

| Qy (0, 0) | ^ 0,

на грубые особые точки. Как уже

было указано в гл. 4, в этом случае система может быть приве­ дена линейным неособым преобразованием к виду

dx/dt — by + Р* (х, у) = Р (х, у), dy/dt =

Q* (х, у) = Q (х, у); (В)

как и в гл. 4, рассмотрим функции:

 

1)

функцию у=ц>(х), являющуюся решением уравнения

 

by + Р* (х, у) = 0;

 

2)

функцию

 

 

z = ф (х) = <2* (х, ф (ж)) = Дт тт

+ . ..

5) Здесь, очевидно, речь идет об особенности кривой, а не дифферен­ циального уравнения.

174 БИФУРКАЦИИ ДИНАМИЧЕСКОЙ СИСТЕМЫ [ГЛ. 10

Здесь /га 3* 2 и Ат ^ 0.

Наряду с системой (В) будем рассматривать всевозможные

измененные системы

 

 

 

 

dxfdt = Р(х,

у) = Р(х,

у) + р(х,

у),

 

dy/dt = Q(х,

y)=Q{x,

y) + q(x,

у),

( В )

 

достаточно близкие до ранга тп к системе (В).

Имеют место предложения, полностью аналогичные теоре­

мам 1 и 2.

 

4.

Если

Ат^ 0,

то кратность особой

точки О си­

Т е о р е м а

стемы (В)

есть тп.

Если

тпкратность состояния равновесия О,

Т е о р е м а

5.

то число

к грубых состояний

равновесия системы

(В) Оi, ...

..,, Ок, на которые О может разделиться при сколь угодно ма­ лых добавках р(х, у) и q(x, у) ранга тп, при тп нечетном может быть любым нечетным числом, меньшим или равным тп, при тп четном любым четным числом, меньшим или равным тп,

ипри этом:

1)если тп нечетно и особая точка О либо имеет характер узла или фокуса, либо является особой точкой с эллиптической

областью, то среди особых точек 0\, 02, . . О* системы (В) число грубых узлов или фокусов на единицу больше числа гру­

бых седел;

седла, то

среди

2) если тп нечетно и О имеет характер

особых точек 0 \, ..., Ок число грубых седел на

единицу

больше

числа грубых узлов или фокусов;

 

 

3) если тп четно, т. е. особая точка О либо является вы­ рожденной особой точкой, либо имеет характер седло-узла, то число грубых узлов или фокусов равно числу седел.

З а м е ч а н и е 1. В рассматриваемом случае о = 0 всегда можно указать такие, сколь угодно малые до ранга тп добавки,

при которых среди грубых особых точек 0\, ..., Ок системы (В) были бы фокусы и такие добавки, при которых существовали бы предельные циклы (целиком лежащие в сколь угодно малой окрестности О).

З а м е ч а н и е 2. Индекс особой точки 0 системы (В), имею­ щей характер узла, фокуса или являющейся особой точкой с эл­ липтической областью, равен +1, индекс особой точки 0, имею­ щей характер седла, равен —1 и индекс особой точки 0, имею­ щей характер седло-узла или вырожденной, равен 0.§

§ 4. Бифуркации двукратной точки, для которой А = 0 и о = 0. В предыдущих параграфах было рассмотрено расщепление при малых изменениях правых частей некоторых типов сложных состояний равновесия на грубые состояния равновесия (можно также установить расщепление на состояния равновесия, среди

§ 4] БИФУРКАЦИИ ДВУКРАТНОЙ ТОЧКИ 175

которых существуют сложные состояния равновесия меньшей

кратности). Для

сложного состояния равновесия, для которого

Д = 0, о Ф 0, все

бифуркации исчерпывающим образом описы­

ваются числом и характером состояний равновесия, на которые это сложное состояние равновесия может разделиться. Для сложных состояний равновесия, для которых А == 0, о = 0, это не так. Действительно, как было указано, из этих состояний равновесия возможно рождение предельных циклов. В настоя­ щем параграфе мы рассмотрим двукратное состояние равнове­ сия, для которого А = 0, 10 = 0 и выполнены еще некоторые дополнительные условия. Такое состояние равновесия неодно­ кратно встречается в рассматриваемых далее примерах, и зна­ ние его бифуркаций дает весьма полезную информацию при ка­ чественном исследовании конкретных систем.

Как мы уже

указывали

в

гл. 4,

состояние равновесия

0 (хо, уо), для которого выполняются условия

 

Д(яо, Уо) = 0 {хо, уо) =

О,

I Рх (х0, у0) | +

| Ру (х0, у0) | +

| Qx (х0, у0) | + | Qy (х0. у0) | ф О,

может быть приведено линейной заменой к виду

dxjdt = by + Pi (х, у) = Р (х, у),

dy'dt = Qt {х, y) = Q {х, у). (1)

Полагая

by + Р* (х. у) = v,

 

 

 

выражая из этого соотношения у через х ж v ж переходя к пе­ ременным х жv, мы придем к системе вида

dxjdt = v, dv/dt = Q* (х, v).

Возвращаясь к прежним обозначениям для переменных и за­ писывая полученную систему подробнее, мы получим

dxjdt = у, dyldt = a20x2 + аиху + а 02у2+ Q3(x,

у),

(2)

где Qo{x, у) начинается с членов степени, большей

или

равной

трем.

Нетрудно убедиться в том, что в силу предположения о дву-

кратности

состояния

равновесия а2й Ф 0. В

силу

теоремы 3 § 2

гл. 4 это

состояние

равновесия

О имеет

вид,

представленный

на рис. 105, а. Предположим, кроме того, что

ац Ф 0. Смысл

этого условия будет ясен из дальнейшего.

 

всевозможные,

Будем

наряду с

системой (2)

рассматривать

достаточно близкие к ней до ранга 3 системы. Так как каче­ ственный характер тех близких систем, у которых в окрестности

0(0, 0)

нет

ни

одного состояния

равновесия, очевиден

(рис. 105, б) , то

мы

обратимся к рассмотрению

близких

систем,

имеющих

состояния

равновесия. Такие

близкие

системы

можно,

176

БИФУРКАЦИИ ДИНАМИЧЕСКОЙ СИСТЕМЫ

[ГЛ. 10

ч

Рис. 105

tb

45

8 4] БИФУРКАЦИИ ДВУКРАТНОЙ ТОЧКИ 177

как нетрудно показать, всегда представить в виде

dx/dt = у, dy/dt =-— E2x + ny + Q2(x, у),

Q 2 ( X , у ) = а 2 о х 2+ а и х у + а о 2 у 2+ ( ? з ( х , у ) , ' '

где Q3(x, у) начинается со степеней по х и у, не меньших трех. Нетрудно видеть, что в некоторой достаточно малой окрестности состояния равновесия 0(0, 0) системы (2) у всевозможных си­ стем (3) может быть: 1) одно двукратное состояние равновесия

того же типа,

что н

у исходной

системы; 2) одно

двукратное

состояние равновесия типа седло-узел

(при е = О, р, ^

0) с устой­

чивым

узловым

сектором

при р < 0,

и неустойчивым при р > О

(рис.

105, в, г);

3) два грубых состояния равновесия: узел или

фокус

и седло

(узел

при

р Ф 0,

е Ф 0 и р2 > 4е2 и

фокус при

р ¥=О,

е=И=0 и

р2< 4

е 2) ;

4) два

состояния равновесия — слож­

ный фокус и седло при р = 0, е Ф 0 (рис. 105, д) . Непосредственным вычислением устанавливается, что первая

ляпуновская величина сложного фокуса имеет вид

•^1 = "^6 [ а 20а 11 "Ь е/ ( a ifc> Е) ] ’

где f(aik, е )— функция, зависящая от коэффициентов aik в раз­ ложении Ог(х, у), в которых г + А > 3. Так как в силу пред­ положения «гоЯц^О, то для всех достаточно близких систем

(3) «гояц ^ 0 и, следовательно,

L\

0.

Таким образом, сложный фокус, возможный у близких си­ стем (при р = 0), является устойчивым или неустойчивым в зависимости от знака выражения агояц, и из него может ро­ ждаться единственный цикл — соответственно устойчивый или неустойчивый (рис. 105, е). Однако того факта, что из сложного фокуса системы (3) может рождаться единственный предельный цикл еще недостаточно для того, чтобы сделать заключение о том, что в достаточно малой окрестности состояния равновесия О системы (3) может быть только один предельный цикл, так как предельные циклы могут еще появляться из двукратных предельных циклов, появляющихся из «уплотнения траекторий».

Если предполагать, что цикл единственный, то для различ­ ных достаточно близких к (2) систем (3) возможна, кроме ука­ занных выше, одна из качественных структур, представленных на рис. 105, ж, з.

Отметим, что знание возможных бифуркаций двукратной точ­ ки с двумя нулевыми корнями может оказаться весьма полез­ ным при качественном исследовании конкретных динамических систем.

12 н. н. Баутин, Е. А. Леонтович

178

БИФУРКАЦИИ ДИНАМИЧЕСКОЙ СИСТЕМЫ

[ГЛ. 10

§ 5. Рождение предельных циклов из особых траекторий степени негрубости выше первой. I. В § 2 было рассмотрено рождение предельного цикла из сложного фокуса первого по­ рядка, т. е. из состояния равновесия с чисто мнимыми характе­ ристическими корнями, для которого первая ляпуновская ве­ личина

L\ = ос3 Ф 0.

Здесь рассматривается случай, когда система (А) имеет слож­ ный фокус кратности т > 1, т. е. состояние равновесия с чисто мнимыми характеристическими корнями такое, что первый, не равный нулю коэффициент а,- ( i > l ) функции последования есть

cc2 m+i = ЬтФ О, т > 1.

 

Как и всюду, наряду с данной системой

 

dx/dt=P(x, у), dy/dt = Q(x, у),

(А)

имеющей сложный фокус порядка т > 1, будем рассматривать измененную систему

dxldt = Р(х, у) = Р(х, у) + р(х, у),

dyldt = Q(x, y) = Q(x, у ) + q(x, у),

где р(х, у) и

д(х,

у ) — достаточно малые добавки ранга

+ 1 .

Имеет место следующая теорема.

 

всегда можно подобрать та­

Т е о р е м а

6. При любом к ^ т

кие

(сколь угодно

малые

добавки)

р(х,

у)

и

q(x, у)

ранга

2т + \, чтобы при

переходе

от системы

(А)

к

системе

(А) из

сложного фокуса рождалось к грубых предельных циклов.

 

Точнее: пусть 0(0, 0 ) — сложный фокус порядка

m системы

(А)

(т. е.

А» =

a 2m+ 1 Ф 0, а< = 0 , i < 2m + 1 ),

тогда

существуют

е0 > 0,

о0 > 0 та­

кие,

что при любых бо-добавках ранга

+

1 у системы (А) в

ео-окрест-

ности О может существовать не

более

т предельных циклов, и,

какое бы

k

мы ни взяли, при любых е < ео, б < б0 существуют б-добавки ранга

+

1 такие, при которых у системы

(А)

в е-окрестности О существует

к предельных циклов, и все эти предельные циклы целиком лежат в е-ок­ рестности О.

Всегда можно также подобрать такие сколь угодно малые добавки ран­ га + 1, чтобы из сложного фокуса О рождался негрубый цикл любой кратности к ^ т, или s предельных циклов соответственно кратностей пи п2, . . ., п, таких, чтобы сумма их кратностей не превышала т.

Пусть теперь О есть центр системы (А). Тогда, какое бы целое к мы ни взяли, всегда можно указать систему (А), сколь угодно близкую до ранга + 1 к системе (А) и такую, что в данной сколь угодно малой окрестности О у этой системы су­ ществует к предельных циклов.

§ 5]

РОЖДЕНИЕ ПРЕДЕЛЬНЫХ ЦИКЛОВ

179

II. Предположим теперь, что система (А) имеет сложный предельный цикл Lo кратности тп выше 2 (тп> 2); т. е. если рассмотреть функцию последования на дуге без контакта I, про­ веденной через какую-нибудь точку предельного цикла L0

 

 

 

 

 

 

s

= a i s +

a.2s 2 + . . .

 

 

 

(s — параметр на дуге I),

то мы будем иметь

 

 

 

 

а х =

ехр

J [р*(ф.Ф) + <?у(ф.Ф)] dt =

1

 

 

 

 

 

 

 

о

 

 

 

 

 

 

(а; = ф(2),

г/ = ф(£)— решение,

соответствующее

предельному

циклу Lo, т — период

на L0) , eta = 0 и первый, не равный нулю

коэффициент функции последования —

 

 

 

 

 

 

 

 

 

 

а т

0 ,

т п > 3 .

 

 

 

Имеет место следующая теорема

(аналогичная теореме 6).

 

Т е о р е м а

7.

Всегда можно указать сколь угодно малые до­

бавки р(х,

у)

и

 

q(x,

у)

ранга тп такие,

чтобы при переходе от

системы (А)

к системе (А)

от замкнутой траектории Lo системы

(А) рождалось любое число к

(к^тп)

грубых

циклов

или ki

[k\ < тп) предельных циклов кратностей п{ таких, что

 

 

 

 

 

 

пх + п2 + . . . + nhi = m.

 

 

то

Точнее: если £ 0 — сложный предельный цикл кратности m системы (А),

существуют

ео > 0, б0 > 0

такие,

что при любых бо-добавках

ранга тп

в

ео-окрестности £ 0

существует

не более m предельных

циклов,

и, какое

бы ft ^ m мы

ни

взяли,

при любом

е > О можно указать б > 0

и такие

б-добавки ранга т,

р(х, у)

и q(x,

у), чтобы у соответствующей системы (А)

в е-окрестности La существовало ft грубых предельных циклов, и все эти предельные циклы лежали в е-окрестности L0 или /ц предельных циклов кратностей п^, п^, . , причем щ + л2 + ... + л* ^ т.

III. Рассмотрим случай, когда сепаратриса седла 0 (х о, уо) образует петлю, и при этом в седле

Рх (х0, Уо) +

(^о’ Уо) = О-

В этом случае возможен как

случай, когда петля устойчива,

так и случай, когда петля неустойчива, а также случай, когда все траектории, проходящие через близкие к петле точки, 8аикнуты.

Можно показать, что в этом

случае

заведомо существуют

такие сколь угодно малые до ранга 3 (или до

ранга т >

3)

добавки, что у соответствующей

системы

(А)

существует

не

менее двух предельных циклов в сколь угодно малой окрестно­ сти петли (см. [84]).

12*

ГЛАВА 11

ДИНАМИЧЕСКИЕ СИСТЕМЫ, ПРАВЫЕ ЧАСТИ КОТОРЫХ СОДЕРЖАТ ПАРАМЕТРЫ

§ 1. Возможный характер зависимости правых частей дина­ мической системы от параметров. Динамические системы, правые части которых зависят от того или другого числа параметров, всегда естественно возникают при рассмотрении различных за­ дач из приложений.

Пусть правые части динамической системы зависят от п па­ раметров, т. е. имеют вид

dx/dt Р(X)

у , Яч,

».»,

%п),

 

dy!dt = Q(x,

у, Хи

...,

К) .

(Av)

Мы будем предполагать, что правые части являются анали­ тическими функциями не только переменных х, у, но и пара­ метров Xi.

Если обратиться к рассмотрению одного из общих про­ странств динамических систем Яд (см. гл. 8), то, очевидно, при каждом наборе значений параметров Х( мы получаем точку в этом пространстве, а при всевозможных значениях параметров

Xi — в пространстве

R™

выделяется

га-мерное подпространство

динамических систем

(Ах).

 

Очевидно, что, в частности, система (Ах) может быть негру­

бой при всех значениях параметров X. Так, например, мы можем

рассматривать систему вида

 

dx

дН (x i УI

, Яп)

dt

 

ду

dy

 

дН (х, у, Я-р .. . ,

Яп)

dt

 

дх

которая является гамильтоновой системой при всех значениях параметров (т. е. с точки зрения введенной классификации — бесконечной степени негрубости)1).

■) В случае, когда мы рассматриваем всевозможные динамические си­ стемы с аналитическими правыми частями. Как уже было указано в § 8 гл. 8, если мы ограничиваемся рассмотрением только множества гамильто­ новых систем, то можно естественным образом ставить вопрос о грубости внутри этого множества, как это фактически и было сделано А. Пуанкаре.

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