книги / Методы и приёмы качественного исследования динамических систем на плоскости
..pdfs 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 и 1ъ = 0 |
< 0, область |
[7] на рис. 158, 161). Область существования |
|
двух предельных |
циклов примыкает |
на интервале А,<0) < X < Хз |
к куску дискриминантной кривой и к прямой L\ = 0. На интер |
||
вале Хт < X < Я,1 |
при смещении с дискриминантной кривой в об |
|
ласть одного состояния равновесия (см. рис. 159,5—9) исчезает состояние равновесия седло-узел и остаются два предельных цикла вокруг устойчивого фокуса (на рис. 159,6—8 неустойчивый пре дельный цикл при исчезновении седло-узла возникает из замкну той траектории, образованной со-сепаратрисой седло-узла). На ин
тервале |
А,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).
