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

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

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

s 2]

Э Л Е К Т Р И Ч Е С К А Я Ц Е П Ь С Т У Н Н Е Л Ь Н Ы М ДИ О ДО М

3 0 1

структура рис. 159,5 с двумя предельными циклами, между ко­ торыми нет состояний равновесия. На интервале 1 < А < А (1> при некотором А = А<0) предельные циклы сливаются в двойной полуустойчивый предельный цикл (рис. 159,4) и затем исчезают. По­ следовательность структур вдоль нижней ветви дискриминантной кривой представлена на рис. 159,4—10.

3.5. Структуры разбиения фазового пространства и бифурка­ ции внутри дискриминантной кривой в области трех состояний равновесия. Дискриминантная кривая, представленная в парамет­ рическом виде уравнениями (3), может рассматриваться как оги­ бающая семейства прямых о XXQ — cp(£o)= 0 в плоскости (А,, о). Полупрямые

о —Ххо —ф (хо) = О, А «S —cp'(zo),

(16)

касающиеся в точке А = —ф'(£о) дискриминантной кривой, одно­ кратно покрывают область, ограниченную линией симметричных структур и нижней ветвью дискриминантной кривой, при хо, из­ меняющемся от точки перегиба до минимума характеристики

ф(ж) (при

b ^

Захо ^

Ъ+ 1/Ь2— Зас).

Движение в

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

параметров

А,

о вдоль

полупрямых

(16) от точки

касания со­

ответствует для системы (1) повороту против часовой стрелки изоклины о —Хх у = 0 вокруг седла, возникшего в точке х = хо при разделении седло-узла на седло и узел.

Рассмотрим бифуркации, осуществляющиеся при движении по полупрямым (16), касающимся дискриминантной кривой на ин­

тервале Ai < А < Аг.

Здесь возникнут

как бифуркации состояний

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

бифуркации

сепаратрис и предельных

циклов.

 

 

При уменьшении А состояние равновесия седло-узел внутри

устойчивого предельного цикла разделяется на седло и неустой­ чивый узел, который при дальнейшем уменьшении А превраща­ ется в фокус (рис. 160,44).

В точках пересечения рассматриваемой полупрямой (16) с пря­ мыми L\ = 0 и 4,2 = 0 происходят бифуркации состояний равнове­ сия: при уменьшении А сначала из фокуса х\ (рис. 160,10) и за­ тем из фокуса Х2 рождаются неустойчивые предельные циклы (фо­ кусы становятся устойчивыми) и возникает структура с тремя пре­ дельными циклами (рис. 160,5). Так как при А = 0 предельных циклов нет = о — интегральная прямая, качественная структура эквивалентна структуре рис. 160,4), то рассуждениями, аналогич­ ными проведенным в п. 3.4, находим, что при убывании А до ну­ ля должны осуществиться следующие бифуркации сепаратрис: возникновение петли сепаратрисы вокруг верхнего фокуса, вокруг нижнего фокуса, возникновение большой петли, содержащей внутри два состояния равновесия. Так как седловая величина по­

ложительна {Рк + Qy= — ф' (х0) — 1 = А — 1, где А — координата

302

КАЧЕСТВЕННОЕ ИССЛЕДОВАНИЕ ДИНАМИЧЕСКИХ СИСТЕМ [ГЛ. 16

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

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

§ 2] ЭЛЕКТРИЧЕСКАЯ ЦЕПЬ С ТУННЕЛЬНЫМ ДИОДОМ 303

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

Для полупрямых (16), касающихся нижней границы дискри­ минантной кривой на интервалах (1, Я,(0>), (Х<0>, А,ш ), (А,(1), Х<2)), (Х(2\ Я[), бифуркации будут аналогичными, но число их от ин­ тервала к интервалу будет уменьшаться за счет того, что неко­ торые бифуркации уже произошли при движении вдоль дискри­ минантной кривой.

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

метричных

структур

и

на

дискриминантной

кривой.

Все

три

бифуркационные

кривые,

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

трем

типам

петель

сепара­

трис,

пересекаются

в

точке

X —Х+ на линии симметрич­

ных структур (рис. 161). Они

заканчиваются

в точках X =

= Х(2>, X ='Ха> и X = 1 на ди­

скриминантной

кривой.

В

точках X =*Х(2> и X = Х(1> осу­

ществляются

структуры

с

петлями сепаратрис.

 

В точ­

ке X — 1 сепаратрисы вырож­ денного седло-узла (см. гл. 4)

нужно рассматривать как вырождение петли сепаратрисы во­ круг верхнего фокуса, стягивающейся вместе с фокусом в одну точку. Кривая двойных циклов проходит между точкой X = X* на линии симметричных структур и точкой X = Х<0> на дис­ криминантной кривой слева от кривой, на которой осуществляет­ ся большая петля.

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

3.6. Структуры разбиения фазового пространства и бифурка­ ции вне дискриминантной кривой (в области одного состояния

304 КАЧЕСТВЕННОЕ ИССЛЕДОВАНИЕ ДИНАМИЧЕСКИХ СИСТЕМ [ГЛ. 16

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

нантной кривой

в области между

прямыми L\ = 0 и = 0

< 0, область

[7] на рис. 158, 161). Область существования

двух предельных

циклов примыкает

на интервале А,<0) < X < Хз

к куску дискриминантной кривой и к прямой L\ = 0. На интер­

вале Хт < X < Я,1

при смещении с дискриминантной кривой в об­

ласть одного состояния равновесия (см. рис. 159,59) исчезает состояние равновесия седло-узел и остаются два предельных цикла вокруг устойчивого фокуса (на рис. 159,68 неустойчивый пре­ дельный цикл при исчезновении седло-узла возникает из замкну­ той траектории, образованной со-сепаратрисой седло-узла). На ин­

тервале

А,1 < X < Хз при

переходе из области L\ < 0 в область

L\ > 0

(с убыванием а)

из фокуса появляется второй неустойчи­

вый предельный ци кл(аз>0

при Х<Хз). Бифуркационная кри­

вая (штриховая на рис. 158,

161), соответствующая слиянию ус­

тойчивого и неустойчивого предельных циклов, начинается в точ­ ке X = Хз на прямой L\ = 0 (здесь аз = 0) и пересекает дискри­ минантную кривую при X = Х<0>, выделяя некоторую окрестность дискриминантной кривой и прямой L\ =>0, для точек которой есть одно устойчивое состояние равновесия и два предельных цикла (область [2] на рис. 158, 161). При переходе из области L\ < 0 в область Ь\ > 0 при X > Хз устойчивый цикл стягивается к фокусу (аз < 0) и возникает структура без предельных циклов (область [5] на рис. 158, 161). Границами области без предельных циклов служат кусок дискриминантной кривой (для 0 ^ X < Хт ), кривая

двойных циклов (для А,10) ^ Х < Х з ) и прямая L\ = 0 (для X

Хз).

3.7.

Разбиение пространства параметров. Разбиение

простран­

ства параметров на области различной качественной структуры по обе стороны линии симметричных структур L = 0 приведено на рис. 161. Соответствующие различным областям грубые структу­ ры разбиения фазового пространства (обозначенные теми же но­ мерами) представлены на рис. 160. Жирными линиями изобра­ жены сепаратрисы и предельные циклы, штриховой — неустойчи­ вые предельные циклы. Устойчивые состояния равновесия —чер­ ные точки, неустойчивые — светлые.

Разбиение на рис. 161 соответствует предположению, сделан­ ному в п. 3.4, об отсутствии предельных циклов для структур в точке X = 1 дискриминантной кривой. Если предположение не вы­ полняется, то кривая двойных циклов (штриховая на рис. 161)‘

будет проходить не через точку

X = Х(0> >

1,

а через

точку X =

= Х(0> < 1, и тогда исчезнут области [2]

и

[5], а на

рис. 160

соответственно структуры 2 ж3.

 

 

 

 

§ 3]

ДВУМЕРНАЯ МОДЕЛЬ ДИНАМИКИ ТВЕРДОТЕЛЬНОГО ЛАЗЕРА

305

 

§ 3. Двумерная модель динамики твердотельного лазера

[119].

Система дифференциальных уравнений

 

 

Ж = С (ге- - ^ Т Т - i ) ™ = P(rn,n),

 

 

Of = noi — (т + 1) п = Q(т , п)

(1)

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

По физическому смыслу задачи рассматриваются область фа­ зового пространства т > 0 и параметры, удовлетворяющие услови­ ям р > 0 и G > 1: noi > 0, гео2 > 0.

Координаты состояний равновесия находим, приравнивая нулю правые части системы (1):

(2)

Ли — (тп + 1) п = 0.

Исключая из этих уравнений п, для абсцисс состояний равновесия получаем

Ординаты же состояний равновесия однозначно определяются абс­ циссами:

п = noi/(m + 1 ) .

(4)

Из (3) следует, что система (1) при всех рассматриваемых зна­

чениях параметров всегда имеет одно состояние

равновесия 0\

с координатами

(5)

Ш\ = 0, щ = ге01

и не может иметь более трех состояний равновесия. Кроме m = 0,

уравнение (3)

имеет еще два корня

(Гj —

тгз = ~2~ (»oi

п°2

^

4 " ( ”OI ~ - п 02— 1 — j - j 2 + - у (п01 — п02— 1) (6)

(первый индекс 2 соответствует знаку плюс перед радикалом, вто­ рой — 3 — знаку минус), являющихся в случае, когда они дей­ ствительны, абсциссами двух других состояний равновесия систе­ мы (1): Оч и 0%. В дальнейшем, изучая разбиение пространства параметров на области с различной качественной структурой, ес­ тественно рассматривать характер разбиения квадранта плоскости (noi, «0 2 ), соответствующего «oi>0, « 0 2 > 0 , при различных зна-

20 в . Н. Баутин, В. А. Леонтович

3 0 6

КАЧЕСТВЕННОЕ ИССЛЕДОВАНИЕ ДИНАМИЧЕСКИХ СИСТЕМ

[ГЛ. 16

чениях р и G. На плоскости (reoi, п02)

точкам прямой

 

 

п<п — п02 — 1 =

0,

(7)

очевидно, соответствуют системы (1), имеющие двукратное состоя­ ние равновесия с равной нулю координатой т, а точкам кривой (параболы)

A23 = 4 " (reoi---- F " 02_1 — т ) + "«г — 1) = ° (8)

— системы, имеющие двукратное состояние равновесия, получаю­ щееся от слияния состояний равновесия 0 2 и 0 3. Исследуя урав­

нение

параболы обычными методами, нетрудно видеть, что

при р <

1 она расположена вне рассматриваемого квадранта плос­

кости (тгоь «ог), а при р > 1 — в этом квадранте.

Кроме того, как нетрудно видеть, при р < 1 корень тз отрица­ телен, и, следовательно, у системы (1) при р < 1, при значениях

reoi, П02 выше прямой (7)— два состояния равновесия,

а ниже —

одно — 0 \ 2)*.

1 для области значений тгоь пог,

Дгз <

0 система

В случае р >

(1) имеет одно состояние равновесия — Оi, а для

области значе­

ний, где Дгз > 0

(ниже прямой (6))— три состояния равновесия

0 1 , О2, Оз. Прямая (7) в этом случае, очевидно, также соответству­ ет системам с кратными состояниями равновесия. Общая точка N параболы (8) с прямой соответствует динамической системе с трех­ кратным состоянием равновесия и делит прямую (7) на две части Д12 и Д1 3 1 2 соответствует слиянию 0\ и 0 2, а Д13 — слия­ нию 0\ и Оз). На рис. 162—164 цифрами /, II, III указаны соот­ ветственно области плоскости параметров, при которых система имеет одно, два и три состояния равновесия с координатой m > 0.

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

A = P'mQ n-PnQ 'm =

= G [ro” -

( (pm”+i)*

+ г а ~ Т ?М :Т ~

1) (та+ 1)] =

0’ (9)

о = Pm + <?n =

— (m. + 1) +

+

И —

l] = 0,

 

 

 

 

(10 )

2) При переходе через прямую (6), для точек которой система имеет двукратное состояние равновесия, двойное состояние равновесия не исчезает (как в общем случае), а опять разделяется на два, но у одного из них координата m делается отрицательной.

§ 3]

ДВУМЕРНАЯ МОДЕЛЬ ДИНАМИКИ ТВЕРДОТЕЛЬНОГО ЛАЗЕРА

307

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

Для того чтобы получить на плоскости (поь Ног) кривую, со­ ответствующую о = 0 (смене устойчивости), нужно исключить т и п из уравнений

а{т, п)= 0, Р(т, п) = 0, Q(m, п) = 0.

(11)

Однако такое исключение весьма сложно, и проще получить па­ раметрические уравнения этой кривой. Исключая п из

Р(тп, п) = 0, Q(m, п) = О,

получаем

 

 

 

 

 

 

 

n0i(pm +

1)— n02(m + 1) —(ртге+ 1) (m + 1) = 0

(12)

и, полагая в (10)

n = n0il(m + 1), с помощью

(10) и

(12)

выра­

жаем noi и по2 через m

(пг — параметр):

 

 

 

 

(m + l)2 (P m + l)

, _ , а

_

(m +

1) (р/ге +

I)2

(13)

----------Gpm---------

+ ™ + 1 >

« 0 2 =

----------Gpm---------

 

Кривую на плоскости (поь ног), определенную равенствами

(13),

или, что то же,

уравнениями (11),

будем обозначать через S.

Исследование характера этой кривой при различных р и G даны в приложении I. Некоторые основные случаи ее характера и рас­ положения относительно кривой Дгз даны на рис. 162—164. Из исследования в приложении I следует, что: 1) кривая S имеет

ветви, уходящие в бесконечность при тге-^0 и m

и

при

тп-+ 0 имеет асимптоту (штриховая прямая на рис. 164);

2)

при

р > 1 кривая S имеет одну общую точку М с кривой Дгз, соответ-

20*

308 КАЧЕСТВЕННОЕ ИССЛЕДОВАНИЕ ДИНАМИЧЕСКИХ СИСТЕМ [ГЛ. 16

ствующую значению т = т. Нетрудно убедиться на основании

выражений (9) и (10), что:

 

 

 

A) Состояние равновесия 0\ — седло, когда т\ — ног — 1 > 0,

и устойчивый узел, когда Hoi — ног — 1 < 0 .

 

 

 

Б) В случае р < 1 состояние равновесия

0 2 , для которого

0 =

0, всегда сложный фокус. В случае р >

1

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

сия,

в котором о = 0 , есть сложный фокус

0 2

для т < т

(т. е.

для

ветви AM кривой S) и седло 0з для т > т (т. е. для

ветви

МВ кривой S ) .

 

 

 

B)

При р > 1 для значений параметров, соответствующих точ­

ке М,

система имеет двукратное состояние равновесия, для кото­

рого о = 0. Это состояние равновесия рассмотрено в § 4 гл. 10.

В точке М кривая S касается кривой Ага (это вытекает из ска­

занного в § 3 гл.

1 1 и может быть также получено непосред­

ственно, если для параболы, как и для кривой S, найти парамет­

рические уравнения). Точкой М кривая

(Агз) разделяется на две

части. Точки одной

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

системе (1),

имеющей седло-

узел с устойчивой

узловой областью,

а точки

другой части —

с неустойчивой узловой областью.

На

рис. 162—164 штрихами

указаны

области, где состояние

равновесия 0 2 соответственно ус­

тойчиво

(о < 0) и неустойчиво

( а > 0 ) .

Очевидно,

при

переходе

в плоскости (ио1 , йог) через кривую S у

системы

(1)

возможно

рождение предельных циклов, на котором мы остановимся ниже. Предельные циклы. Покажем сначала, что существует зам­ кнутая область в фазовом пространстве, содержащая все состоя­ ния равновесия, внутрь которой входят траектории системы. Дей­ ствительно, рассматривая прямые Гх: та = 0 , Г2 :и = ио1 , можно показать, что все траектории входят внутрь полосы, ограниченной

этими прямыми.

Рассматривая поле направлений на прямой

Г3 : т = 0i, где

0 1 — некоторая положительная постоянная, для

Hoi ^ 1 получа­

ем, что траектории пересекают эту прямую сверху вниз. Таким

образом, для

Hoi ^ 1 существование замкнутой области показано,

а для Hoi ^ 1

необходимо еще рассмотреть поле

направлений на

прямой

Г4 :

тп = — G(noi — 1)н + 0 2 . Выбирая

0 2

достаточно

большим,

можно убедиться, что на интервале 1

< Hi

noi траек­

тории пересекают эту прямую справа налево, т. е. по-прежнему существует область, внутрь которой входят все траектории си­ стемы (1 ).

Принимая во внимание тот факт, что существуют значения параметров, при которых система имеет одно состояние равно­ весия, лежащее на интегральной прямой тп= 0 , а также прини­ мая во внимание, что состояние равновесия 0з — седло, можно утверждать:

а) если предельные циклы существуют, то они обязательно ох­ ватывают состояние равновесия 0 г;

• 3]

ДВУМЕРНАЯ МОДЕЛЬ ДИНАМИКИ ТВЕРДОТЕЛЬНОГО ЛАЗЕРА

309

чая

б) предельных циклов в области значений параметров /, вклю­

границы — прямую (7) и кривую (8), быть не может;

 

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

Мы покажем, что в некоторых случаях у системы может су­ ществовать два предельных цикла. Рассмотрим выражение для

первой ляпуновской

величины

(полученное

в приложении II)

L{m) = l!&+*)... ' К

v

'

4 д У д

т (pm + I) 2

где /(то) = р(р + Gp — 1)то2 + 2(р — 1 — Gp)m + р + 1.

1)р < 1. В этом случае, очевидно, существует не более двух

состояний равновесия, седло лежит на прямой то = 0, и, очевид­ но, его сепаратрисы не могут образовать петлю (рис. 165).

Возможны два случая:

а) р — 1 + Gp < 0. Уравнение /(то)= 0 не имеет положитель­

ных корней,

L ( m ) < 0

при всех то.

б) р — 1

+ Gp > 0.

Граница области устойчивости S состоит

из двух участков, на

одном из которых L < 0, на другом L > 0;

 

 

Р и с .

165

в точке Мз, граничной для

этих участков, L = 0. В точках, где

К О , при переходе из области,

где о < 0, в область, где о > 0,

у системы из фокуса

рождается

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

а в точках, где L > 0,

при

обратном переходе из области, где

о > 0 , в области, где о < 0 ,

рождается неустойчивый предельный

цикл.

Особенности в поведении системы в окрестности тех значений, при которых L = 0, можно выяснить, устанавливая знак второй ляпуновской величины Ьч. Однако мы можем рассмотреть воз­ можности, которые здесь могут иметь место и без вычисления этой величины.

310

КАЧЕСТВЕННОЕ ИССЛЕДОВАНИЕ ДИНАМИЧЕСКИХ СИСТЕМ [ГЛ. 16

 

Если, выйдя из достаточно близкой к прямой

(7)

точки, в ко­

торой, как мы уже говорили выше

(см. п. б)),

нет

предельных

циклов, мы перейдем из области о <

0 в область а > 0 через учас­

ток кривой S, где L < 0, то при этом у системы рождается устой­ чивый предельный цикл. Пока мы не пересечем вновь кривую S, этот устойчивый предельный цикл сохранится — ему некуда деть­ ся (если в области о > 0 есть точки, соответствующие двукрат­ ным циклам, то тогда во всяком случае число устойчивых циклов все время будет на единицу больше, чем неустойчивых). Если за­ тем мы пересечем границу S на части ее, где L > 0, переходя из области о > 0 в область а < 0, то при этом у системы из фокуса рождается неустойчивый предельный цикл, устойчивый же сохра­ няется, и таким образом теперь у системы есть устойчивый и не­ устойчивый предельные циклы (или одинаковое число и тех и других).

Если, не пересекая кривую S, мы вернемся в область вблизи прямой (7), где нет предельных циклов, то мы непременно долж­ ны пройти через значения параметров, при которых эти предель­ ные циклы исчезают. Следовательно, должна быть бифуркация, при которой устойчивый и неустойчивый циклы сливаются, обра­ зуя четнократный цикл, который затем исчезает. Непременно су­ ществует, следовательно, в области параметров, где а > 0, бифур­ кационная кривая >ST , упирающаяся в точку Мз, соответствующая четнократному циклу (см. рис. 162).

Возможные качественные структуры на фазовой плоскости

представлены на рис. 165.

S есть общая точка

2)

р > 1. В этом случае у кривых Дгз и

М (предположим, что она лежит на части S,

где L < 0). При

двукратное состояние равновесия, для которого о — 0. Бифуркации такой точки описаны в гл. 10. В частности, при

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

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