
книги / Методы и приёмы качественного исследования динамических систем на плоскости
..pdf392 |
МЕТОД ПОНТРЯГИНА В КУСОЧНО-СШИТЫХ СИСТЕМАХ |
[ГЛ. 18 |
|
Кривая 1р2 = гМЛДо) для любых Яо, отличных от нуля, |
ухо |
дит в бесконечность при h -*■°°, но будет лежать внутри полосы сколь угодно малой ширины е на любом заданном интервале
изменения h, если выбрать Хо достаточно малым. |
\pi(/i) + |
||||||
Проследим |
за изменением |
числа |
нулей |
функции |
|||
+ тЫ^До) = ф(ЛДо) |
на интервале 0 |
< Л < ° ° |
при возрастании |
||||
параметра Хо от нуля. |
е < |
I 'ф1 т т (Л*) I. Если Ао |
выбрано |
так, что |
|||
Пусть е < |
I ipi (0) I, |
||||||
на интервале 0 < h < h* |
кривая 'фг(ЛДо) лежит внутри полоски |
||||||
шириной е (—е < ф 2 |
< 0 ), то |
число |
нулей функции ipi + тр2 по |
сравнению с ipi не может измениться ни за счет изменения знака ipmin в точке минимума, ни за счет изменения знака \|3i(0 ). С дру гой стороны, как бы ни было мало Хо, функции г|ц + фг и г|ц име
ют при достаточно большом h разные знаки Ш т |
= 2 лу0» |
\ h ~ * оо |
|
lim (TJJJ + ф2) = |
— оо). Очевидно поэтому, что при возрастании Хо |
||
h - ю о |
) |
|
нуль функции ф] + фгТак |
от нуля из бесконечности появляется |
|||
|
/ |
/ |
0 , то при этом у системы |
как в нуле, очевидно, будет % + |
ф2< |
(8 ) из бесконечности появляется устойчивый предельный цикл.
Если ао |
и |
[} выбраны так, |
что \|3i,mi n < 0 , i|Ji(0)>0 и Яо>0 |
|||||
и мало, то система (8 ) |
будет |
иметь три предельных цикла (два |
||||||
устойчивых п один неустойчивый). |
|
|
||||||
§ 3. Автоколебания синхронного мотора. Рассмотрим систему |
||||||||
d(p/dt = |
у, |
dy/dt = Т —simp — р (а + р cos 2<р — if sin Ц>)у, |
||||||
|
|
|
|
а, р, |
к, Т > |
О, |
|
|
где р — малый параметр |
[59, 60, 67]. |
|
|
|||||
В качестве фазового пространства будем рассматривать по |
||||||||
лосу, заключенную между прямыми <р = —я |
и ср = я; точки этих |
|||||||
прямых, имеющие одинаковые координаты у, отождествляем. |
||||||||
Полагая |
(рис. 203) |
— 2 я -1<р — 2 , |
— я ^ ф < 1 — я/2 , |
|||||
|
|
|
||||||
sin ф « |
s (ф) |
|
2 я _1 ф, |
— я/2 |
< ф < я / 2 , |
|||
|
|
|
— 2 я _1ф -f- 2 , |
я / 2 |
ф я; |
|||
|
|
|
— 4я_1ф — 3, |
— я ^ ф <: — я /2 , |
||||
cos ф « |
с (2 ф) |
|
4я_1ф + 1, |
- я/ 2 ^ ф ^ О, |
||||
— 4я_1ф + 1, |
0 |
^ ф ^ 2 /л, |
||||||
|
|
|
||||||
|
|
|
|
4я-1ф — 3, |
|
«Ss Y • |
получим систему, близкую к кусочно-линейной консервативной системе
dq>/dt = у = Р, dy/dt = Т - s(<p)- IL[<X + $с(2<р)~ 4s((p)]y = Q (1)
(р — положительный параметр).
§ 31 |
АВТОКОЛЕБАНИЯ СИНХРОННОГО МОТОРА |
393 |
|
|
1. Консервативная система. |
|
|
|
d<p/dt = у, |
dy/dt = Т — s(q>). |
(2) |
Система (2) имеет замкнутые фазовые траектории при Т < 1, Замкнутые фазовые траектории системы (2) охватывают состоя ние равновесия. Другие траектории представляют собой спирали, накручивающиеся на фазовый цилиндр (рис. 204).
Сшивание траекторий консервативной системы внутри полосы
(—я, я) |
(кусков эллипсов |
и гипербол) |
осуществляется |
на пря |
||||||||||
мых ф = ±я/2. |
Сшивание |
траекторий |
системы |
(1) |
происходит |
|||||||||
также и на прямой ф = 0 . |
|
|
|
|
|
|
|
|
||||||
Консервативная система (2) имеет интеграл |
|
|
|
|
|
|||||||||
н ( ф, »)— £ |
+ |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
F |
(ф + я + ip) |
+ 2лТ |
|
|
|
. |
. |
|
Я |
|||
|
|
|
|
Я < |
ф |
< |
— - у , |
|||||||
+ |
|
|
|
|
|
|
|
|
|
Я _ |
_я |
, |
||
|
|
|
|
|
|
|
= |
й’ |
у < |
ф |
< |
т |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
--7 Г |
Ф~ л + "г) |
|
|
|
■f < Ф < Я . |
|
|||||||
|
|
|
|
|
|
|
|
|
||||||
Замкнутым |
кривым |
соответствует интервал |
изменения |
k |
от |
|||||||||
|
I |
|
|
|
(состояние равновесия типа |
центр) |
до |
h — 0 |
||||||
h — ----5 - я (1 — Г)2 |
||||||||||||||
(петля сепаратрисы седла). Точка пе |
|
|
|
|
|
|
||||||||
ресечения |
петли |
сепаратрисы с |
осью |
|
|
|
|
|
|
|||||
ф лежит |
в |
полосе |
(0 , я /2 ) , если |
2 — |
|
|
|
|
|
|
||||
—У2 < Г < 1 |
; в полосе |
(—я/2 , 0), |
если |
|
|
|
|
|
|
|||||
3 —2 У2 |
< Г < |
2 |
—V2 ; |
в |
полосе |
(—я, |
|
|
|
|
|
|
—п/2), если 0 < Г < 3 —2У2. В зависи мости от Т изменяется число сшива ний для периодического движения. В дальнейшем ограничимся случаем, когда периодическое движение сшива
ется из траекторий не более чем трех полос ( Т > 3 —2У2). При этом будут исчерпаны все возможные случаи разбиения фазового пространства на траектории.
394 |
МЕТОД ПОНТРЯГИНА В КУСОЧНО-СШИТЫХ СИСТЕМАХ |
[ГЛ. 18 |
2. Состояния равновесия. Седловая величина. Система (1)
имеет при Т < 1 два состояния равновесия: Oi(nTI2,0) и 0г(я — яГ/2, 0). Состояние равновесия 0\ будет неустойчивым фокусом, если для точки 0 \ будет
(^ ф + = — р (а + Р — 2РТ — уТ) = — рог > 0, а < 0,
устойчивым фокусом, если а > 0, и центром, если а = 0. В послед
нем легко убедиться: |
при условии а = а + р —2$Т —уТ = 0 си |
||
стема (1 ) на полосе 0 |
< ф < я / 2 |
может быть представлена в виде |
|
а, |
( г — И |
“ - » 1 1 2 1 1 + *)■'] |
|
dq> |
|
у |
|
и проинтегрирована. |
|
|
величина |
Состояние равновесия Ог — всегда седло. Седловая |
|||
будет иметь значение |
|
|
|
|
(^*ф + Qy)i = — F 7- |
(3) |
Существенно, что выражения £*Ф+ Qy для фокуса и для сед ла имеют одинаковое значение. Как будет показано, это опре деляет существенные особенности структуры разбиения фазового пространства и пространства параметров.
3. Предельные циклы. Если систему (1) записать в виде
dyldt = Ну, dy/dt = — Яф + р9 (ф, у), |
(4) |
где ф(ф, !/) = [—а — рс(2 ф)+фз(ф)]1/, то значения константы |
ho, |
выделяющие кривые Снаконсервативной системы, вблизи которых
при |
малом р будут |
предельные циклы системы (4), определя |
ются как корни уравнения ф(й) = 0 , где |
||
Ф ( h ) |
= f J Я у d<f> d y = |
—j j [a + Рс(2ф) — у«(ф)] cfy d y = |
= — J fa + pc (2ф) — ys (ф)] у dq>. (5) ch
Чтобы выполнить интегрирование, нужно сделать определен ные предположения о расположении петли сепаратрисы консер вативной системы на плоскости и тем самым о числе кусков, из которых сшивается кривая
а) Пусть |
2 —У 2 < Г < 1 . |
При изменении h |
в |
интервале |
||
—я(1 — Г)2^ |
< /г< —я(1 — Г)2/4 |
кривые |
Ch будут |
эллипсами, |
||
целиком расположенными |
в |
полосе |
(0 , я/2 ) , |
а |
интервалу |
(—л/4) (1 —Т)2< h < 0 будут соответствовать замкнутые кривые, сшитые из кусков эллипса и гиперболы, расположенных соответ ственно на полосах (0 , я /2 ) и (я/2 , я ) .
396 МЕТОД ПОНТРЯГИНА В КУСОЧНО-СШИТЫХ СИСТЕМАХ [ГЛ. 18
для Г > 0 ,5 будет
_ JL (i _ Tf < - ■ £ ( 1 - T f < - - f ( Г 2 - AT + 2 ) < 0 .
Выполняя |
интегрирование |
формуле (5) для |
случая |
0,5 < |
|||||||||||
.< Т < 2 —V2, получим |
|
|
— a2 < f t < |
— а2 /2 , |
|
||||||||||
|
|
|
— oF0(h), |
|
|
||||||||||
|
ijf(ft) = |
— oFi(h), |
|
— а2 / 2 < |
ft < |
— ft2 /2 |
, |
/7\ |
|||||||
|
|
|
F . W - o F ^ h ) , |
|
- f t 2 /2 < f c < 0 . |
|
|
|
|
||||||
Здесь ft2 = |
~ - л (Т2— АТ + 2), |
F0(h) |
и Fi(h) |
имеют |
значение, |
||||||||||
указанное в |
(6 ), и Fi(h) |
будет |
|
|
|
|
|
|
|
|
|
||||
Ft {К) = |
2 (2 а2 - |
ft2) / |
2 Л + |
ft2 + -i-(2A + |
ft2 / |
2 — |
|
|
|
|
|
||||
|
|
|
|
|
— 2(к + а2)У 2а2— Ьг(лЦ — 2 arcsin |
|
—— |
||||||||
Имеем также |
|
|
|
|
|
|
|
|
|
|
|
|
V 2(h+a*) |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
F, (h) = 4 (2h + |
ft2/ |
2 - 2 / 2 а 2 |
- |
b2 л — 2 |
arcsin |
/ |
2 а2 - ь 2 \ |
||||||||
|
|
|
|
|
|
|
|
|
|
|
/ |
2 |
(fe+ |
а2)/ |
|
Fl (h) |
2 / |
2 fe + |
ь2 |
|
|
|
|
|
|
|
|
|
|
|
|
М -а2 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
F2( - |
ft2 /2) = |
F't ( - |
ft2 /2) = F; ( - |
Ь212) = 0, |
F%( - |
ft2 /2) > 0, |
|||||||||
Из последнего следует, что / ’г(Л) на |
интервале |
—b2/ 2 < h ^ 0 |
есть функция, принимающая положительные значения, монотон но возрастающая и обращающаяся в нуль в единственной точке h = - ft2/2 .
Выражения (7) определяют ф(Л) как непрерывную функцию, сшитую из трех кусков с различными аналитическими представ лениями на интервалах (—а2, —а2/2 ), (—а2/2 , —ft2 /2 ) и (—ft2/2 , 0 ). Функция г|з(Л) на интервале (—а2, —ft2 /2 ) имеет знак, противо
положный знаку о, обращаясь при а = 0 в нуль |
тождественно. |
На интервале (—ft2 /2,0) функция ty(h) состоит из |
двух слагае |
мых, одно из которых положительно, а другое имеет знак, про тивоположный знаку о. Выбором о значение ф(й) на всем интер вале (—ft2 /2 , 0 ) может быть сделано как положительным, так и отрицательным.
Проследим за корнями функции ф(й) и изменением качест венной структуры разбиения фазового пространства на траекто рии при возрастании о от отрицательных значений. При о < 0 функция ф(й) на всем интервале (—а2 , 0 ) будет положительной.
6 3] АВТОКОЛЕБАНИЯ СИНХРОННОГО МОТОРА 397
Предельных циклов нет. Фокус неустойчивый. Структура раз биения фазового пространства изображена на рис. 205,1.
При о = 0 |
функция O|J(A) тождественно обращается в нуль |
на интервале |
(—а2, —Ь2/2 ) и сохраняет положительное значение |
на интервале (~Ь2/2,0). Существует континуум замкнутых кри вых, окружающих состояние равновесия типа центр. Петля се паратрисы не может существовать, так как о|з(0 )¥=0 ; (о-сепарат- риса седла скручивается с границы области, заполненной замк
нутыми кривыми; «-сепаратриса седла |
уходит в |
бесконечность |
||||||||||
по верхнему полуцилиндру (рис. 205,7—2). |
|
о |
достаточно |
|||||||||
При |
о > 0 |
фокус становится |
устойчивым. Если |
|||||||||
мало, |
то |
на |
интервале |
(—а2, —Ь2/2) будет |
а|з(й)<0 , |
но сохра |
||||||
нится |
тр (0) > 0. |
На |
интервале |
(—Ь2/2,0) |
будет |
существовать |
||||||
h = h\ — корень |
функции т|з(h), |
соответствующий |
н е у с т о й ч и |
|||||||||
во м у |
(тр' (hi) > 0 ) |
предельному циклу, возникшему из границы |
||||||||||
области, заполненной замкнутыми кривыми (рис. 205,2). |
||||||||||||
При возрастании |
о кривая ар (А) |
будет опускаться — при этом |
||||||||||
корень |
о|:(h) |
будет |
возрастать — и |
для |
достаточно |
больших h |
||||||
функция а|з(й) становится отрицательной на всем |
интервале |
|||||||||||
(—а2, |
0). Предельных |
циклов |
нет. Фокус |
устойчив. |
Структура |
398 |
МЕТОД ПОНТРЯГИНА В КУСОЧНО-СШИТЫХ СИСТЕМАХ |
[ГЛ. 18 |
||
разбиения фазового пространства |
однозначно |
определяется |
||
(рис. 205, 4). |
переходе |
от структуры |
||
|
Исчезновение корня h = hi при |
рис. 205,2 к структуре рис. 205, 4 может происходить либо при возрастании ho до нуля (этому соответствовало бы влипание не у с т о й ч и в о г о предельного цикла в петлю сепаратрисы седла), либо при слиянии его с другим корнем функции i|i(h). Не про водя подробного исследования поведения функции г|>(/г) при раз ных о, можно по знаку седловой величины заключить, что при
возрастании о реализуется именно последний случай. |
(3), |
имеет |
||||
Седловая |
величина, определяемая |
выражением |
||||
при о > 0 отрицательное значение и, |
следовательно, |
к |
петле |
|||
сепаратрисы |
может стянуться или |
из |
нее появиться |
лишь |
||
у с т о й ч и в ы й предельный цикл. Отсюда |
следует, |
что |
обраще |
ние в нуль величины г|>(0 ), соответствующее возникновению петли сепаратрисы, должно предшествовать исчезновению корня
h = h\. При |
перемене знака ф(0) |
(когда ф(0) становится |
отри |
|||||
цательным) |
появляется |
второй |
корень |
h — Аг функции |
ф(й), |
|||
убывающий |
с возрастанием о и |
соответствующий |
у с т о й ч и |
|||||
в о м у |
(г|/(Й2 ) < 0 ) |
предельному циклу. При дальнейшем |
возра |
|||||
стании |
о корни hi |
и йг |
сближаются, |
сливаются |
и исчезают. |
Описанному процессу соответствует на фазовой плоскости смена
структур, представленных на рис. 205: рис. 205, 2—3 |
(ф(0) = |
|
= 0)— возникла петля сепаратрисы; рис. 205,5 (AI </»2 |
< 0 ) — |
|
из |
петли сепаратрисы возник устойчивый предельный |
цикл; |
рис. |
205,5—4 (hi = h2) — возник двойной полуустойчивый пре |
дельный цикл; рис. 205,4 (ф (/г)< 0 ) — циклов нет.
Приведенные простые рассуждения дают строгое доказатель ство существования области пространства параметров, для точек которой система (1 ) имеет по крайней мере два предельных цик ла. Непосредственное исследование функции г|>(/г) позволяет убедиться, что г|>(Л) не может иметь более двух корней. Это исследование, однако, довольно громоздко [67].
Г Л А В А 19
КАЧЕСТВЕННОЕ ИССЛЕДОВАНИЕ СШИТЫХ СИСТЕМ МЕТОДАМИ ТЕОРИИ БИФУРКАЦИИ
§ 1. Кусочно-линейная система с тремя параметрами. Рас смотрим дифференциальное уравнение [71]
х + а [1 — $F' (х) ] i + F(x) = 4 .
Полагая
F W = |
"ЗГ х |
при |
я ____ л |
|
||
2 |
^ ^ ^ 2 |
’ |
||||
|
||||||
F = |
— ~k х + 2 |
при |
я |
______ 3 |
: |
|
2 |
2 |
|||||
|
и вводя новые переменные и параметры, приходим к системе вида
dx/dt = у, |
|
при |
я |
_ |
_я |
|
||
dy/dt — — х — 2hxy + |
а |
2 |
|
^ |
' 2 |
’ |
||
dx/dt = у, |
|
при |
я . |
. |
3 |
( 1) |
||
dy/dt = х — 2h2y + а |
— я |
2 |
^ |
х ^ |
2 |
: |
||
|
Будем рассматривать фазовую цилиндрическую поверхность склеенной системы, развернутую на часть плоскости, соответст
вующей |
неравенствам —я / 2 ^ х ^ |
|
< Зя/2 (рис. 206). Прямая |
х = я/2 |
|
разбивает |
рассматриваемую |
часть |
плоскости на области 1 и 2, в каж дой из которых фазовые траектории определяются соответственно линей ными системами. Прямые х = —я / 2 И 1 = Зя/2 отождествляются. Склеен ная система имеет два состояния равновесия: Oi(a, 0 ) и Ог(я — а, 0 ).
Точка 0{ — устойчивый фокус |
при |
0 < /г < 1, устойчивый узел при |
1 . Точка 0 2 — всегда седло, |
сепаратрисы которого определяются уравнениями
У ~ ( — ^ j ± V + 1 ) [# — (я —- о )].
400 МЕТОДЫ ТЕОРИИ БИФУРКАЦИЙ В СШИТЫХ СИСТЕМАХ [ГЛ. 19
Обозначим через Si, 5г и S3 полупрямые х = —я /2 , х = я / 2 И 1 = Зя/2, соответствующие значениям у > 0.
Фазовые траектории системы осуществляют точечные преоб
разования полупрямой |
S i в S 2 (преобразование |
L i ) |
и полупря |
||
мой S 2 |
в S 3 (преобразование L 2 ) . Пусть S i , |
S 2 |
и S 3 |
— ординаты |
|
точек |
соответствующих |
полупрямых; т/ам, |
О/002 — времена про |
бега изображающей точки (см. рис. 206) через области 1, 2, со ответствующие преобразованиям L \ и L 2 .
Величины S i , S 2 , т, 0 принимают положительные или нулевые значения. Интегрируя линейные уравнения (1) в областях 1 и 2, обычным образом получаем параметрические уравнения для функций соответствия.
Для области 1 получаем
Si (т, К |
а) = |
|
( |
|
|
А |
т |
+ |
ctg( - тJ ++ hj),а ) К |
||
S%(т, h |
а) = |
Si (т, |
hu |
d), |
|
|
|
|
(Li) |
||
Здесь cOj = V |
|
i — h\, |
ki = hi/d)i |
и 0 < h < i . |
Выражения |
для |
|||||
производных могут быть представлены в виде |
|
|
|
|
|||||||
dS1 |
. P2hiT ll. |
|
|
|
|
|
|
|
|||
dSl |
-e |
5 i, |
|
|
|
|
|
|
|||
£ f i _ |
|
|
|
|
|
|
|
|
|
|
|
dS% |
|
|
|
|
|
|
|
|
|
|
|
В случае |
h > |
1 |
выражения |
для |
Si (т, hi, |
а) |
и |
Sz(T, hi, |
а) и |
производных получим, если в правых частях полученных выра
жений заменим sin т и ctg т соответственно |
через sh т и |
cth т. |
|
Параметр т меняется в пределах от т = 0 |
(5 2 (0 ) = 5 I (0) = 00) |
||
до значения т = т*, при котором Si |
или S2 обращаются в |
нуль. |
|
Тот или иной случай при 0 < h < 1 |
реализуется в зависимости |
||
от знака выражения |
|
|
|
(я/ 2 — а) Л л — (я/ 2 + а). |
|
(2 ) |
|
При h > 1 всегда реализуется случай 5г(т*) = 0, S i (т*)> 0 . |
Кри |
||
вая 5 1 = 5 1 (5 2 ) имеет асимптоту 5 1 = |
5 2 + 2 я/&1 . Некоторые |
воз |
можные виды кривой, соответствующие преобразованию Li, изо бражены на рис. 207.
Рис. 207, а соответствует случаю (я/2 — а) ек*п — (я/2 + а) <С 0 й, в частности, всегда реализуется при малых hi. Промежуточ ный между изображенными случай, когда кривая проходит че рез начало координат, соответствует т* = я и обращению в нуль выражения (2 ).