Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Учебники / НЕЛИНЕЙНАЯ ДИНАМИКА СЛОЖНЫХ.pdf
Скачиваний:
640
Добавлен:
30.03.2022
Размер:
28.89 Mб
Скачать

http://profbeckman.narod.ru/

Таким образом, зная след и якобиан, легко определить тип устойчивости неподвижных точек и найти условия всех локальных бифуркаций. Можно показать, что неподвижная точка ( x a ,

y=0) является неустойчивой, а неподвижная точка ( x a , y=0) - устойчивой.

Анализ бифуркаций двумерного отображения удобно выполнить, используя плоскость параметров след – якобиан матрицы монодромии. Устройство этой плоскости в общем случае показано на рис. 23. На нем представлен треугольник устойчивости неподвижной точки двумерного отображения, границы которого, в общем случае, заданы линиями 1-S+J=0 ( =+1), 1+S+J=0 ( =-1) и J=1 (| |=1) . Таким образом, для устойчивой неподвижной точки двумерного отображения имеют место три бифуркации коразмерности 1 на границах треугольника устойчивости: касательная бифуркация (μ=+1), бифуркация удвоения периода (μ=-1) и бифуркация Неймарка-Сакера (|μ|=1). А так же три бифуркации коразмерности 2 в вершинах треугольника: μ1=-1 и μ2=-1 (резонанс 1:2), μ1=+1 и μ2=+1 (резонанс 1:1) и μ1=-1 и μ2=+2 (fold-flip

бифуркация).

Рис. 23. Треугольник устойчивости неподвижной точки двумерного отображения на плоскости след S - якобиан J матрицы монодромии.

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

16.5 Бифуркация рождения пары устойчивых замкнутых траекторий.

Пусть при μ<0 в системе существует устойчивый цикл x0(t,μ) один мультипликатор которого равен +1, а все остальные мультипликаторы лежат внутри единичного круга. При переходе значений параметра μ через точку μ=0 один простой мультипликатор цикла переходит через точку +1 единичной окружности, что соответствует пересечению мнимой оси одним простым вещественным показателем Флоке линеаризованной системы или одним вещественным собственным значением матрицы D(μ). Поэтому нормальная форма такой бифуркации цикла совпадает с одномерной нормальной формой бифуркации особой

точки типа вилки:

=

Õ̇

. Циклу, очевидно, соответствует нулевое решение этого

уравнения.

 

 

 

 

В результате бифуркации цикл x0(t,μ) теряет устойчивость (но не исчезает), и

одновременно рядом с

ним на расстоянии

 

в направлении собственного

вектора, соответствующего показателю,

проходящему через мнимую ось, рождается пара

 

4 = ±

 

устойчивых предельных циклов. На рис. 24

родившиеся

циклы лежат в параллельных

исходному циклу плоскостях, так как вектор и направлен вертикально.

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

цикла, а другой вне него.

Рис. 24. Рождение из устойчивого цикла (а)

двух других

устойчивых

циклов (б)

(бифуркация вилки предельных циклов).

В случае

размерности

фазового

пространства

m>2,

цикл,

потерявший

устойчивость,

становится

седловым.

Бифуркация является мягкой.

В качестве примера рассмотрим систему трех дифференциальных уравнений

̇ = −

+ (

+ 1)(1 −

) ,

http://profbeckman.narod.ru/

,

̇ =

 

+

+ (

+ 1)(1 −

− )

 

T

̇ =

Õ −

 

 

 

 

 

Эта система имеет своим решением особую точку O=(0,0.0) , устойчивую при μ<-1. При μ>-1 из нее в результате бифуркации Андронова-Хопфа рождается устойчивый предельный цикл x0(t)=(cost, sint,0)T, лежащий в плоскости переменных (x1,x2). Линеаризованная на цикле система примет вид:

̇ = −2

( + 1) cos æ−

 

(

+ (

+ 1) sin2 æ) −

 

 

sin æ+ ,

,

 

̇ =

(

− ( + 1) sin2 æ) − 2

(

+ 1) sin

æ+

cos æ+

 

 

̇ =

Õ −

 

1

 

1

2

2 1 2

 

 

 

 

 

 

где разложения функций f (y ,y ) и

f (y ,y ) в ряды в точке (0,0) начинаются с членов второго

порядка. Заменой переменных y(t) =Q(t)z(t) с 2π/ периодической матрицей

 

 

 

( ) = ( ̇( ),

( ),(0 0 1) ) =

 

 

− sin

cos

0

̇= +

( ,

 

координатах:

1

),

 

 

,

 

 

 

 

 

0

0

перейдём к системе, записанной в связанных с циклом

 

cos

sin

0

̇ = −2(

+ 1).

+ (

,

)

 

 

 

 

 

 

 

 

 

̇ =

 

 

 

 

 

 

 

 

 

 

 

 

 

Линейная

часть

последней

системы

имеет

действительную

постоянную матрицу

G(t, )=G(μ)=Β(μ) с нулевым первым столбцом, т.е. координата z1 однозначно определяется координатами z2 и z3 системы, лежащими в плоскости S, перпендикулярной плоскости цикла. Поэтому устойчивость цикла и его возможные бифуркации определяются исключительно двумерной системой уравнений

̇ = −2(

+ 1).

+ ( ,

),

3

3

 

̇ =

2

2

с диагональной матрицей D(μ) = diag(-2(μ+1),μ). Заметим,

записанной в координатах u =z

и u =z

 

что диагональные элементы матрицы Ω(μ) являются показателями Флоке предельного цикла исходной системы. Цикл является устойчивым в интервале изменения параметра -1<μ<0. При μ > 0 нулевая особая точка и цикл системы становятся неустойчивыми. Происходит рождение двух

новых устойчивых особых точек u3=z3=x3=

 

 

, лежащих на отрезке, перендикулярном плоскости

цикла. Этому соответствует рождение

двух устойчивых циклов, лежащих в параллельных

 

 

расстоянии =

 

от неё. Бифуркация рождения из

плоскости исходного цикла плоскостях на

 

 

циклов

предваряет субгармонический и

одного устойчивого цикла двух других

 

устойчивых

 

 

 

 

гомоклинический каскады бифуркаций рождения

устойчивых циклов в системе Лоренца.

 

 

 

16.6 Транскритическая (обмена устойчивостью между циклами) бифуркация.

В этом случае при μ<0 наряду с устойчивым циклом x0(t, μ) в системе существует также и неустойчивый (седловой) цикл x1(t, ), лежащий на расстоянии и u=μ от устойчивого цикла в направлении собственного вектора, соответствующего показателю Флоке устойчивого цикла (простому вещественному собственному значению матрицы D(μ)), проходящему через мнимую ось. При переходе значений параметра μ через точку μ=0 один простой мультипликатор цикла, соответствующий указанному показателю Флоке, переходит через точку +1 единичной окружности, а циклы обмениваются устойчивостью - устойчивый цикл x0(t, μ) становится неустойчивым (седловым), а

неустойчивый (седловой) цикл x1(t,μ) ~ устойчивым (рис. 25).

Рис.

25.

Бифуркация

обмена

устойчивостью между циклами.

 

Нормальная форма такой бифуркации цикла совпадает с одномерной нормальной формой транскритической бифуркации

 

 

 

http://profbeckman.narod.ru/

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

 

(на рис. 25 направление вектора и совпадает с вертикалью).

Бифуркация

является мягкой. Если направление вектора и совпадает с горизонталью, то

 

̇= Õ −

 

при этом оба цикла должны лежать в одной плоскости (один внутри другого). В результате бифуркации внутренний устойчивый, например, цикл становится неустойчивым, а внешний неустойчивый устойчивым циклом.

16.7Бифуркация исчезновения (рождения) пары замкнутых траекторий.

Вэтом случае при μ<0 наряду с устойчивым циклом x0(t,μ) существует также и неустойчивый седловой цикл x1(t,/ ). При возрастании значений параметра μ эти циклы сближаются друг с другом и при μ = 0 сливаются в один вырожденный цикл. При μ>0 оба цикла исчезают (рис. 26). Поэтому нормальная форма такой бифуркации цикла совпадает

с одномерной нормальной формой седло-узловой бифуркации особой точки:

 

(на рис. 26 направление вектора и совпадает с вертикалью). Бифуркация

является жесткой

 

̇= +

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

̇ = −

В качестве примера рассмотрим систему трех дифференциальных уравнений

 

 

 

 

 

 

+ (

+ 1)(1 −

) ,

,

 

 

 

 

 

 

̇ =

.

 

+ ( + 1)(1 −

)

 

 

 

 

 

 

 

̇ =

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При -1<μ< 0 исходная система кроме устойчивого предельного цикла

 

 

 

 

− ) , sin(

− ) ,

 

− )

, лежащего в плоскости

= −

неустойчивый

от него в

 

 

 

( )

= (

 

)

, имеет также и ( ) = (

î (

 

 

 

î (

+

− ) , sin( + − ) , −

, лежащий на расстоянии

2

 

предельный цикл

 

 

 

. Заметим, что движение по

 

 

 

 

 

плоскости

 

 

циклам происходит с разными

частотами. При μ=0 оба

цикла сливаются, образуя лежащий√ √ в плоскости х =0 неустойчивый цикл

,

=

 

 

 

 

 

3

 

 

 

исчезающий при >0. Анализустойчивости цикла также сводится к анализу устойчивости нулевой

особой точки двумерной системы

 

 

 

 

 

 

 

 

̇ = −2( +.

1)

 

+

(

,

),

 

 

 

 

 

 

 

 

 

̇ =

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если рассмотреть процесс изменения параметра μ от больших значений к меньшим, то описываемая бифуркация означает, что при μ=0 внезапно возникает замкнутая траектория (цикл), которая при дальнейшем уменьшении значений параметра расщепляется на две замкнутые траектории, одна из которых устойчивая, а другая нет.

Другими словами, происходит внезапное одновременное рождение устойчивого и седлового предельных циклов. Каскад бифуркаций этого типа, порождающих бесконечную последовательность пар различных устойчивых и седловых циклов субгармонического и гомоклшшческого каскадов бифуркаций, можно наблюдать в системе Лоренца и практически во всех других нелинейных системах дифференциальных уравнений с хаотической динамикой

Рис. 26. Седло-узловая бифуркация предельных циклов.

16.8 Бифуркация удвоения периода цикла

Эта бифуркация может произойти только в фазовом пространстве размерности m>2. Она связана с переходом через мнимую ось слева направо одного комплексного показателя Флоке α(μ) исходного цикла x0(t, μ) при μ=0. При этом соответствующий вещественный мультипликатор цикла при μ=0 имеет вид: .

Переход одного комплексного показателя Флоке цикла х0(t.μ), имеющего период Т, через мнимую ось, что эквивалентно переходу мультипликатора цикла через точку -1 единичной окружности, означает, что цикл x0(t,μ) теряет устойчивость (но не исчезает), и

http://profbeckman.narod.ru/

одновременно возникает другой устойчивый цикл, имеющий ту же амплитуду и удвоенный период 2Т (рис. 27).

 

 

Рис. 27. Суперкритическая бифуркация удвоения периода

 

 

цикла: С0-исходный цикл; С - цикл удвоенного периода после

 

 

бифуркации.

 

 

 

 

 

 

 

Бифуркация удвоения соответствует второй границе

 

 

зоны устойчивости неподвижной точки, когда =-1.

 

 

Суперкритический вариант этой бифуркации получил

 

 

широкую известность благодаря той роли, которую она играет

 

 

при возникновении хаотического аттрактора по сценарию

 

 

Фейгенбаума.

 

 

 

= − − .

Его неподвижные точки

Рассмотрим отображение

 

 

 

 

 

4 + 4

 

 

 

 

 

 

x * и эволюцией ее собственного значения, которое можно

Будем следить за точкой, =1−2 ±

2

= −1 ±

+ 1

 

вычислить следующим образом:

 

 

 

= −(2

 

+ 1) = −1 + 2 − 2 + 1 = 1 − 2 + 1

(26)

 

 

 

 

 

 

При -1< <0 1<1 и неподвижная точка x1* устойчива. В точке =-3/4 1=0. Это означает, что при переходе через это значение характер сходимости меняется от монотонного при -1< <-3/4 к знакопеременному при >-3/4. При =0 1=-1. Этот случай соответствует потере устойчивости неподвижной точкой. Траектория начинает расходиться, причем отклонение меняет свой знак на каждой итерации. Если бы отображение было линейным, то это было бы единственным событием. Однако наличие

нелинейности обогащает картину бифуркации.

 

 

Рассмотрим дважды примененное отображение (26):

 

= − ( − − ) − ( − −

)

(27)

 

 

Определенные выше неподвижные точки отображения (27) будут таковыми и для дважды примененного отображения, при этом траектория в их окрестности уже не осциллирует. Однако расчёт собственного значения в точке =0 даёт значение +1!

В дважды примененном отображении критической ситуации 1=-1 соответствует совсем другое событие, 1,2=+1. При =0 имеет место вилообразная бифуркация. x1* теряет устойчивость, однако в её окрестности рождаются две новые устойчивые неподвижные точки.

Рис. 28. Бифуркация удвоения.

С точки зрения исходного отображения обе родившиеся точки принадлежат одной и той же траектории – циклу периода 2, т.е. отвечают условию

xn+2=xn.

Таким образом, потеря устойчивости точкой x1* при =0 приводит к рождению устойчивого цикла удвоенного периода. По этому признаку данный случай

получил название «бифуркация удвоения».

Функция последования отображения (27) качественно не меняется в момент бифуркации =0, но неподвижная точка теряет устойчивость и при >0 возникает замкнутая траектория из двух неподвижных точек.

Субкритическая бифуркация удвоения.

Отметим, что условие 1=-1 описывает локальную динамику системы (27) в точке x*1 , тогда как наличие и устойчивость цикла удвоенного периода определяются

http://profbeckman.narod.ru/

свойствами отображения в некоторой окрестности x*1. Другой вид нелинейности может привести к иному сценарию развития бифуркации. Например, отображение вида

xn+1=-xn(1+ )-2xn3

(28)

имеет неподвижную точку x*1=0, которая также претерпевает бифуркацию

удвоения при =0, так как 1=-(1+ )-6(x12)=- -1.

этом случае претерпевает не

Однако дважды примененное отображение в

суперкритическую, а субкритическую вилообразную бифуркацию. Это означает, что в отображении (28) до точки бифуркации при <0 существует неустойчивый цикл периода 2, который «влипает» в x*1 в момент бифуркации удвоения при =0. Далее траектория уходит от потерявшей устойчивость точки x*1 в бесконечность.

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

Примером может служить динамика численности биологической популяции, если наблюдение за её численностью производится, например, один раз в год.

Пример 14. Рассмотрим бифуркацию удвоения периода в системе трёх дифференциальных уравнений

̇ = − −

− 2

− (1 − − ) ,

,

 

 

 

 

̇ =

+

− 2

− (1.

− − )

 

 

 

 

 

̇ = 2

Эта

+ +

− 1

 

 

 

 

 

, 0)

, лежащий в плоскости

 

система имеет предельный цикл

система имеет вид:

переменных (х12). Линеаризованная на цикле

 

( ) = (cos ,

¾

 

 

 

̇ = 2

cos −

( − sin2 ) −

( ¾

+ 2cos ) +

 

 

̇ =

( +

sin2 ) + 2

sin +

(cos − 2sin ) +

 

 

 

1 1̇ = 2

cos + 2 sin + 2

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

где разложения функций f (y , y2) и f2(y1, у2)

в ряды в точке (0,0)

порядка. Заменой y(t)=Q(t)z(t) с 2π/ν- периодической матрицей Q(t) приведем исходную систему (к связанным с циклом координатам

̇= +

( ,

),

, ),

̇ = 2

− 2

+ (/

̇ = 2

+ 2

 

 

Анализ устойчивости цикла сводится к анализу устойчивости нулевой особой точки

двумерной системы

 

, ),

̇ = 2

− 2

+ (,

̇ = 2

+ 2

 

 

матрица D(μ) линейной части которой имеет комплексно сопряженные собственные значения 2μ± 2i. Следовательно, при μ<0 нулевая особая точка и цикл x0(i) устойчивы. При μ>0 в системе в результате бифуркации Андронова-Хопфа рождается устойчивый предельный цикл периода Τ π, что соответствует рождению устойчивого двумерного тора с частотами, приближенно равными и и 2. Четкие переходы от резонансных к нерезонансным торам можно наблюдать при малых μ > 0 и варьировании параметра .

Рассмотренная бифуркация рождения устойчивого тора T2 является супер критической (мягкой) бифуркацией. Наряду с ней, как и в случае бифуркации Андронова-Xопфа. может происходить также и субкритическая (жесткая) бифуркация. В этом случае при μ< 0 система имеет устойчивый цикл, лежащий внутри неустойчивого тора. При μ>0 тор влипает в устойчивый цикл, который после этого становится неустойчивым. Происходит исчезновение аттрактора, и. следовательно, бифуркация является кризисом.

Случай бифуркации удвоения периода цикла, которому соответствует переход при μ=0 одного комплексного собственного значения матрицы Ε(μ) через мнимую ось важен с точки зрения приложений. При этом все остальные собственные значения должны иметь отрицательные вещественные части.

Пример 15. Рассмотрим систему трёх дифференциальных уравнений

̇ = − − Ñ − ( − 1) + 1 (1 − − ),

,

̇ = + 2 (1 − /4) − ( − 1) (1 − −

)