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

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

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

§ 4]

СИСТЕМА СО СКАЧКАМИ НА ЛИНИИ СШИВАНИЯ

 

421

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

(5) неявно определяет функцию 0 = 0(р)— вре­

мя перехода из фиксированного и на полупрямую U\. Ее пре­

дельное

значение при р, -►0

(доказательство

аналогично доказа­

тельству в приложении I)

есть

 

 

 

 

 

 

 

0

=

г In [rp/(rp —и)].

 

 

(8 )

Во втором из равенств (5) положим

0 = 0(р)

и

найдем

limwi [0(р), |х] при р-»- 0 . Это предельное значение есть

 

 

 

 

 

 

1* 1 = 0.

 

 

 

(9)

Фиксируем

теперь

и = гр.

Рассмотрим

изменение

отрезка

мо = мо(р) с уменьшением р. Для этого вычислим

 

 

 

 

*4(0) =

2P(4fc*b+l)/[r(P-v)].

 

 

Если

Ь> — 1/(4ft2) ,

то

при

малых р

величина

ио

убывает

вместе с

р и точка и = гр участвует в преобразовании S. Первое

из равенств (5)

определяет

функцию 0 ^©(р) — время

перехода

из и — гр на полупрямую

Ui. Ее предельное значение при р -►0

есть 0 =

оо (это следует из того, что 0 в равенстве (8 )

можно сде­

лать как угодно большим, если

взять и достаточно близким к

гр, и из того, что в S производная dufdQ >

0).

 

 

 

Случай 6 < — 1/ (4ft2)

рассматривается

аналогично. Величина

по возрастает с уменьшением р,

и точка и = гр участвует в пре­

образовании Т. Предельное время есть т = «>.

 

 

При

Ь= — 1/ (4ft2) имеем гр

щ. В этом случае точка и =* гр

при любом р лежит на сепаратрисе. Время доопределенного дви­ жения есть «>.

Для ц < г ( р — *j) нужно еще получить уравнения движения изображающей точки по прямой ф = фо. Для любого и это будет

у - У ( * ) = 0 .

(1 0 )

Полученная схема доопределяемых движений по верхней ча­ сти прямой ф = фо приведена на рис. 219, а . Здесь для наглядно­ сти налегающие одна на другую траектории раздвинуты в гори­ зонтальном направлении. Изображающая точка системы (2), попавшая в точку (ф = фо, у = и), и 3* r(p — 1 ), совершает мгно­ венный скачок в соответствии с формулой (7) и продолжает дви­ жение при ф > ф о в соответствии с системой (2 ).

Если гр < и < г(р — у), то изображающая точка перескакива­ ет в точку v = 0 , находится здесь в течение времени (6 ) и затем двигается при ф > фоПри и = г$ изображающая точка переска­ кивает в точку (ф = фо, у = 0 ) и остается там неограниченно долго. Это движение является вырождением при р. -►0 движения по сепаратрисе седла системы (3).

Если и < гр, то изображающая точка мгновенно перескакива­ ет в точку и\ = 0 и находится здесь в течение времени, опреде­

422

МЕТОДЫ ТЕОРИИ БИФУРКАЦИЙ В СШИТЫХ СИСТЕМАХ 1ГЛ. 19

ляемого формулой (8 ), после чего продолжает движение при

Ф < фо-

Из соображений симметрии следует, что если в описанной схеме поменять и на щ, v на ui, на — 4 и наоборот, то получим схему доопределенных движений по нижней части прямой ф = (po­

ll. $ < 0 < ч , r ^ 2 b h < 0. Состояние равновесия системы (3) при малых р — устойчивый узел, лежащий в полосе —р < х < р.

Соответствующая картина фазовых траекторий приведена

на

рис. 218,

б.

 

U

Пусть

по = wo (р)"— отрезок, отсекаемый на полупрямой

траекторией, проходящей через точку

(я = р, у = 0). При р

0

имеем ио(р)-*- г(£ —к),Ио(р)-> + оо.

Поэтому для w>r(fJ —к)

вычисления совпадают с соответствующими вычислениями п. 1 .

Если же полутраектория

системы (3) начинается

на отрезке

и

г[3([3 —к), то при малых р она целиком лежит

внутри поло­

сы

—р < х < р, и время

соответствующих доопределенных дви­

жений есть оо. Уравнением движения по прямой ф = фо, как и в случае I, будет уравнение (1 0 ) .

Полученная схема доопределяемых движений по верхней ча­ сти прямой ф = фо приведена на рис. 219, б. Изображающая точ­ ка системы (2), попавшая в точку (ф = фо, У = и), и > г(р —к)> совершает мгновенный скачок в соответствии с формулой (7) и продолжает движение при ф > фо в соответствии с системой (2 ).

При п < г ( р —ч) изображающая точка мгновенно переска­ кивает в точку (ф = фо, у = 0 ) и остается там неограниченно дол­

g 4] СИСТЕМА СО СКАЧКАМИ НА ЛИНИИ СШИВАНИЯ 423

го. Будем говорить, что эти движения ведут в состояние равно­ весия (ф = Фо, у = 0 ), являющееся вырождением при р, -*• 0 узла системы (3).

Из соображений симметрии следует, что если в описанной схе­ ме поменять и на щ, v на щ, то получим схему доопределенных движений по нижней части прямой ф = фо.

III. 0 < "г < Р, г = 2bh > 0. Состояние равновесия

системы

(3 ) есть седло, расположенное правее полосы —р < а ; <

р.

Не приводя вычислений, опишем получающуюся в этом слу­ чае схему доопределенных движений (рис. 219, в). Изображаю­

щ а я точка

системы ( 2 ) , попавшая в точку (ф = фо, у = и), и >

>

7 ),

совершает мгновенный скачок в соответствии с фор­

мулой

(7)

и продолжает движение при ф > фо в соответствии с

системой (2 ).

то

изображающая точка мгновенно пере­

Если м sg r(fl —у),

скакивает

в точку и1

= 0

и находится здесь в течение времени,

определяемого формулой (8 ), после чего продолжает движение

при ф < фо. Изображающая точка

системы (2), попавшая в точ­

ку (ф = фо, У = щ), v\ S* г(р —ч),

совершает мгновенный скачок

вверх в соответствии с равенством u\ =

V\ — г([1 — к) и продолжа­

ет движение при ф < фо в соответствии

с системой (2). Если щ <

< г(Р —Ч),

то

изображающая

точка

мгновенно

перескакивает

в точку щ = 0

и

остается

здесь в

течение времени

т = г X

X In [г[3/ (г'у + v\) ],

после чего

продолжает движение

при

ф < фо.

IV. О <

(J <

Y»

г 2bh < 0.

Состояние

равновесия

системы

(3) при малых

р — устойчивый узел,

лежащий

левее

полосы

—р < х < р. Схема доопределенных движений та же, что и в и. III.

Случай V: р < ^ < 0 , г < 0 ,

и случай VI: ^ < § < 0 , г > 0 ,

по­

лучаются соответственно из случаев IV и III заменами: р на —

и на V\, к на MI и наоборот.

 

 

 

еще шесть,

соответ­

Из имеющихся шести случаев получим

ствующих изменению знака г, заменами: р на —•{,

г

на —г, и

на

v, и\ на v\ и наоборот.

 

 

 

прямых

ф = фо

са­

Рассматривая точечные преобразования

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

П р и м е р 1. Пусть Е(ф)— функция периода 2я такая, что F (ф) = ф/я при —я < ф < я, и пусть 0 < О < 1 . Этот пример ил­ люстрирует первый случай доопределения.

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

424

МЕТОДЫ ТЕОРИИ БИФУРКАЦИЙ В СШИТЫХ СИСТЕМАХ

(ГЛ. 19

у =

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

полупрямой F{<p = —л,

v > 0 } в

полупрямую

£/{<р = л, у =

и >

0 },

осуществляемое

траекториями

системы (2).

Если обозначить

 

 

 

 

 

h± = h(n — b )v /V n ,

a 1 = V l — hl,

 

^

= /t1 /cu1,

то при 0 < hi < 1 параметрические уравнения

функции

соответ­

ствия этого преобразования запишутся в виде

 

 

 

 

 

 

(1 Q) ek i4 + (1 -f fi) ^cost) +

■—

 

sin r)

 

v = 0^ У I

s i n r ]

 

 

 

 

( 1

ft) e

^ + ( 1

— ft) ( cos Д —

sin д

 

U =

сох / я -

 

 

s m t ]

 

 

 

 

 

 

 

 

 

 

 

Ее производные —

 

 

 

 

 

 

 

 

dv

_

2 ft n j f _

 

 

 

 

 

 

 

du

 

 

v

 

 

 

 

 

 

Ц

= / л

sin ne«V>

 

-

+

.

 

 

du

 

 

 

 

 

(OjP

 

 

Асимптота —

 

 

v = и + Al/nhi.

 

 

 

 

 

 

 

 

 

Введем

параметр

a

In [(1 + £2)/(l— £2)].

Для

a > /схл кривая

и ■= u(v)

изображена на рис. 2 2 0 вместе с ломаной

 

v

v(u) =

и 2 г при

и ^

2 г,

 

 

 

0

при г ( 1

— £2 ) О <

2 г,

 

 

 

 

 

являющейся графиком функции соответствия преобразования

полупрямой ф = л,

у > 0 в себя, осуществляемого

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

ными движениями,

и построенной по формуле (7).

На рис. 220

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

рому

изображающая точка имеет остановку в

точке (ф = я,

у — 0 )

на время, определяемое первой строкой

формулы (6 ).

Легко видеть, что цикл устойчив и что более одного цикла быть не может. При увеличении г точка пересечения кривой и лома­ ной попадает в начальную точку ломаной, что соответствует влипанию цикла в сепаратрису, образующую петлю (сепаратрису, идущую из состояния равновесия (я, 0 ) в то же самое состояние равновесия). Соответствующая бифуркационная поверхность в пространстве параметров г, hu а имеет уравнения и(г\) =*=г ( 1 — £2 ),

§ 4] СИСТЕМА СО СКАЧКАМИ НА ЛИНИИ СШИВАНИЯ 425

l>(ri)=0. Исключив Т], получим

а = к1 я — arctg

)

+ Т

1+ 2 "l/nfcjT

In

V-я + V

/

2

При уменьшении г точка пересечения кривой и ломаной перехо­ дит на наклонный участок ломаной, чему соответствует превра­ щение «цикла с остановкой» в точке (я, 0 ) в «цикл без останов­

ки». Поверхность, разделяющая в пространстве параметров обла­

сти существования цикла

с

остановкой и цикла без остановки,

имеет уравнение u-(-q) = 2r,

K(TJ) = 0. Их можно переписать так:

 

"lAlsinT)

I

ft Ч

0)х cos т) — ht

ь>1е 1

 

а= — In (— cos TJ kt sin t|)

ирассматривать как параметрические уравнения сечения поверх­ ности плоскостью h\ = const.

На рис. 221 приведено сечение пространства параметров плоскостью hi = const. В области I нет циклов. В области II си­ стема имеет цикл с остановкой. В области III система имеет цикл без остановки. Качественная картина фазовых траекторий систе­

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

2 2 2 .

П р и м е р

2. Этот пример отличается от примера

1 знаками

^(ф) и г и иллюстрирует второй случай доопределения.

Уравнения

функции соответствия u = u(v), ее производных

и асимптоты получаются из соответствующих уравнений приме­

ра 1 , если считать hi = h (я + Ъ)/1/я, ©i = 1 + h\, заменить Q на —Q, тригонометрические функции на гиперболические со сходным названием и поставить минус перед второй производной.

Для получения всех значений

п и и нужно

менять TJ от 0

до °°, при этом v меняется от —

 

до

н* = К я ( 1

— Q) ((Oj + Ax),

а и меняется от °° до и* =

( 1

+

Q) (юх — /гх),

где к* и и* —

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

426

МЕТОДЫ ТЕОРИИ БИФУРКАЦИЙ В СШИТЫХ СИСТЕМАХ

[ГЛ. 19

отрезки, отсекаемые сепаратрисами седла (—я£3, 0), соответственно от полупрямых V и U. Кривая u — u(v) изображена на рис. 223 вместе с прямой v = и + 2 г, являющейся графиком функции со­ ответствия при преобразовании полупрямой ф = я, у > 0 в себя, осуществляемом доопределенными движениями.

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

попадает в точку (и*, у*). При этом предельный цикл влипает в сепаратрису седла (—я£3, 0), образующую петлю. Соответствую­

щая бифуркационная

поверхность

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

г, hi, S3 имеет уравнение

lr/(A i c o J + h j 0)х.

у* = и* +

2г или S3=

Качественная картина фазовых траекторий системы для случая, когда имеется предельный цикл, приведена на рис. 224.

П р и л о ж е н и е

I. П усть

< и <

г(Р — f). Функцию и = (т, р ),

оп­

ределяемую первым из равенств (4), не существующую при р =

0, доопре­

делим ее пределом при р -*■ 0. П олучивш аяся при этом ф ункция

 

 

 

 

 

 

(Ч4(т, р)

 

 

при р > 0 ,

 

 

 

 

 

Т’ ^

У ехр {— т/г}] при р =

0

 

 

удовлетворяет

следующим

условиям:

1)

непрерывна,

2) если

То =

г X

X In [и /(гр — в)],

то /( То, 0)

=

и, 3) производная дЦдх <

0.

 

 

Следовательно, при зафиксированном на рассмотренном интервале и

равенство и =

/( т, р)

определяет в некоторой

окрестности точки

(то, 0) од­

нозначную

непрерывную функцию т =

т„ (р)

такую, что

т* (0) =

т0. А так

как т (р) =

т* (р)

при

р Ф 0,

то lim т (р) =

т0 при р -*■ 0.

 

 

 

Пусть

теперь

в ^ г ( Р — f).

Так как

ди/дх < 0, то

функция

т =

т (р ),

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

фиксированному на

рассматриваемом

интервале и,

огра-

§4]

СИСТЕМА СО СКАЧКАМИ НА ЛИНИИ СШИВАНИЯ

427

ничена

сверху функцией т = т (ц ), соответствующей любому

и, фиксиро­

ванному

на интервале rp < и < г((1 — 'у), а предел этой ограничивающей

функции при р,-»-0 может быть сделан как угодно близким к нулю, если

взять и достаточно близким к

г(Р — -у).

 

 

П р и л о ж е н и е II. При

гР < и < г(Р — if). Это

предельное

значение

находится непосредственной подстановкой предельного значения

т (ц ). Ес­

ли и ^ г(Р — f ) , то

положим

в равенствах (4) к = ' сот. П ервое из

них при

фиксированном и и

малых р. определяет функцию к =

х (р ) такую , что

 

 

lim х (р) =

*

 

и

 

 

 

 

 

 

— In

 

 

 

 

 

 

 

ц->о

 

2

и — г (Р — У)’

 

 

 

(Доказательство

аналогично

доказательству в

приложении I.)

Положим

теперь во втором

из

равенств

(4) к — х (р )

и найдем, что

при

р

0 будет

lim г[к(р), р] =

и — r(P — if).

 

 

 

 

 

 

 

П р и м е р

3.

Пусть Е(<р)— функция

периода

такая, что

 

 

 

( —

1

при — я <

<р < 0 ,

 

 

 

 

 

^ (^)

1 +

1

при

0 < < р < я .

 

 

 

Схема доопределенных движений по прямой <р = я

приведена на

рис. 225, а (здесь масштаб по оси у уменьшен в 2h раз). Изобра­

жающая

точка, попавшая в

ь'ь

точку

(я, у), где

у> 2 Ь , со-

вершает

скачок

величиной

 

2 Ъ вниз по прямой ф = я, по­

 

сле чего продолжает движе­

 

ние при ф >

я в соответствии

 

с системой

(2). Если 6(1 —

 

—R) < у < 26, то изображаю­

 

щая точка

перескакивает в

 

точку (я, 0 ) и находится в

 

ней в течение времени

 

 

 

 

6(1 + 0)

 

2 b h In у - 6 ( 1 - О ) ’

 

после

чего

продолжает дви­

 

жение

при

ф > я

в соответ­

 

ствии

с

системой

(2). Если

 

у = b (1

— £2 ), то изображаю­

 

щая точка перескакивает

и

рис 225

точку (я, 0 ) и остается там не­

6(1 — О ), то изображающая точка

ограниченно долго. Если 0

< у <

перескакивает в точку (я, 0 ) и находится в ней в течение времени

6 (1 — О)

2 b h In 6 (1 — О) — у'

после чего продолжает движение при ф < я. Поведение изобра­ жающей точки, попавшей на прямую ф = я в нижнем полуци­ линдре, аналогично (нужно поменять у на —у и R на —й).

28*

428

МЕТОДЫ ТЕОРИИ БИФУРКАЦИЙ В СШИТЫХ СИСТЕМАХ

|ГЛ. 1»

лу,

Точка (я, 0) будет состоянием равновесия, аналогичным сед­

а траектории, проходящие через

точки

[я, 6 ( 1

—£2 )] и

[я,

—6(1 + £2)]— его сепаратрисами. Роль двух

других

сепарат­

рис выполняют траектории, выходящие из точки

(я, 0 ).

приведе­

 

Схема доопределенных движений по

прямой ср = 0

на на рис. 225, б. Изображающая точка, попавшая в точку (0, у), где у > 0 , совершает скачок величиной 26 вверх по прямой ср = 0 , после чего продолжает движение при ср > 0 в соответствии с си­ стемой (2). Если у < 0, то изображающая точка совершает ска­ чок величиной 26 вниз по прямой ср = 0 , после чего продолжает движение при ср < 0 в соответствии с системой (2) . Точка (0, 0)

и

будет состоянием равновесия, аналогичным неустойчивому узлу. Из него выходят траектории, проходящие через точки (0, у), где 0 < у < 26 или —26 < У < о . (Вертикальный участок такой тра­ ектории от точки (0 , 0 ) до точки (0 , у) изображающая точка про­ ходит скачком.)

Когда система доопределена на прямых ср = 0 и ср = я, можно проследить любое ее частное решение на любом отрезке времени

§ 4]

СИСТЕМА СО СКАЧКАМИ НА ЛИНИИ СШИВАНИЯ

429

и провести

ее качественное

исследование. Исследование

сходно

с исследованием в примере 1

, поэтому ограничимся лишь изложе­

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

Возьмем на прямой <р = 0 две полуоси: полуось и — вниз от начала координат и полуось v — вверх от начала координат. На

рис.

226

сверху изображены

траектории,

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

ото­

бражение

v (и) части полуоси

и в полуось

v (слева),

и траекто­

рии,

осуществляющие отображение и (v) части

полуоси v в

по­

луось и (справа). Графики функций v(u)

и u(v) приведены на

том

же рис. 226 внизу. Для и < щ функция

v(u)

находится

обычным способом, только с учетом того, что все v нужно сде­ лать больше на величину скачка. Участок щ ^ и ^ 112 весь пере­ ходит в одну точку vо. Для этих и имеем v(u)= const = i>o- Ана­ логичные соображения можно высказать о функции u(v).

Цикл соответствует точке пересечения линий v(u) и и (v). На рис. 227 показано изменение взаимного положения линий при изменении параметров. На рис. 228 — соответствующие превра­ щения предельного цикла (соответствующие одна другой картин­ ки на рис. 227 и 228 обозначены одинаковыми буквами). Оче­ видно, так же, как и в примере 1 , расположение траекторий рис. 228 может быть изображено на цилиндре.

Для некоторой области в пространстве параметров линии име­ ют точку пересечения на кривых участках (рис. 228, а). Цикл не имеет кусков на прямых <р = —я и <р = я. При соответствующем изменении параметров точка пересечения передвинется в точку

излома линии u(v). При этом цикл коснется прямой <р = я. При дальнейшем изменении параметров точка пересечения переходит на вертикальный участок линии (рис. 228, б) . При этом на цикле появляется вертикальный участок на прямой ф = я, который изображающая точка проходит скачком, после чего некоторое вре­ мя стоит в точке (я, 0 ), прежде чем продолжит движение по циклу. В зависимости от дальнейшего изменения параметров воз­ можны превращения двух типов: 1 ) концевая точка линии u(v) попадает на кривой участок линии v(u). При этом цикл влипает

430

МЕТОДЫ ТЕОРИИ БИФУРКАЦИЙ В СШИТЫХ СИСТЕМАХ

[ГЛ. 19

в сепаратрису состояния равновесия (я, 0 ), образующую

петлю,

и исчезает. 2) Точка пересечения переходит на горизонтальный участок линии v(u) (рис. 228, в). На цикле при этом появляется

Р и с . 228

вертикальный участок и на прямой ф = —я. При дальнейшем из­ менении параметров концевая точка линии u(v) попадает на го­ ризонтальный участок линии v(u). Цикл влипает в сепаратрису состояния равновесия (я, 0 ), образующую петлю, и исчезает.

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