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

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

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

§ 3] СИСТЕМЫ, ЗАВИСЯЩИЕ ОТ ПАРАМЕТРА 131

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

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

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

грубыми.

В гл. 8 дается точное определение грубой динамической си­ стемы и при этом уточняется смысл слов «малые изменения динамической системы».

Грубость динамической системы именно н можно считать тем

свойством, которое мы выше назвали с у щ е с т в е н н о й

н е к о н ­

с е р в а т и в н о с т ь ю 3). До сих

пор мы все

время

говорили

лишь о динамических системах,

правые части

которых — анали­

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

§ 3. Измененные системы. Системы, правые части которых зависят от параметра. Прежде чем переходить к определению грубой системы, понятию, являющемуся основным в дальнейшем, мы приведем некоторый необходимый вспомогательный материал.

3) Следует отметить, что хотя многие понятия теории бифуркаций ди­ намических систем переносятся на случай многомерных систем, но в этом случае все значительно сложнее, и имеют место совсем иные факты (см. [111, 25, 24, 137-141]).

9*

132

ДИНАМИЧЕСКИЕ СИСТЕМЫ ТЕОРИИ КОЛЕБАНИЙ

[ГЛ. 7

 

Всюду в дальнейшем наряду с заданной системой

 

 

dx/dt = P(x, у), dy/dt = Q(x, у),

(А)

которую мы оудем предполагать определенной в некоторой огра­ ниченной области G (или замкнутой ограниченной области 5), будем рассматривать также другие системы вида

dxjdt = Р(х,

у) = Р(х,

у) + р(х,

у),

(А)

dy/dt = Q(x,

y) = Q(x,

y)+ q(x,

у).

 

Систему (А) мы будем называть исходной

системой,

отличные

от (А) системы (А) — измененными системами.

у) назы­

Функции р = Р(х, у ) - Р ( х , у) и q = Q(x,

 

y ) - Q (x ,

ваются добавками к правым частям системы

(А). В дальнейшем

всегда предполагается, что правые части являются в рассматри­ ваемой области_аналитическимп функциями 4) .

В области G каждая из систем (А) и (А) задает свое век­ торное поле. Синус угла между направлением векторного поля, заданного системой (А), п направлением векторного поля, задан­ ного системой (А) в каждой точке, дается выражением

Q(*’ V)p (*’ V ) - P

У)

(2)

У Р 2 (*, У) + Q* (*, у) У Р2 (X, у) +

Q2 (X, у)

 

Очевидно, в точках, в которых

Q(x, у)Р(х, у) — Р(х, y)Q(x, у)> О,

угол 0 положителен, в точках, в которых

Q(x, у)Р(х, у ) - Р ( х , y)Q(x, у)< О,

этот угол отрицателен, а в точках, где

Q(x, у)Р(х, у) — Р(х, y)Q(x, у) = О,

направления поля систем (А) и (А) совпадают или прямо про­ тивоположны. В том частном случае, когда во всех точках пло­ скости (плп рассматриваемой области)

 

Q(x, у)Р(х, у ) - Р ( х , y)Q(x, у)> 0 ,

 

мы

будем говорить, что система (А) дает поворот поля

систе­

мы

(А) (или просто поворот поля) на неотрицательный или не­

положительный угол.

 

 

 

 

Предположим, что рассматривается динамическая система,

правые части которой зависят

от некоторого

параметра,

 

 

dx/dt = Р(х, у, ц),

dy/dt = Q(x,

у, ц),

(Ац)

4) Понятие грубости динамической системы может быть введено при значительно более широких предположениях относительно правых частей системы (А), именно при предположении, что правые части имеют лишь непрерывные частные производные (см. § 8 гл. 8 ).

§ 4] ТЕОРЕМЫ О ЗАВИСИМОСТИ РЕШЕНИЯ 133

причем эту систему при некотором частном значении параметра, например при р = 0 , т. е. систему

dxldt = P(x, у, 0), dyldt = Q(х, у, 0), (А0)'

мы будем принимать за исходную систему (А). Тогда изменен­ ной системой будет система (А„) при |Х^ О, и она, очевидно, может быть получена из исходной системы (А) с помощью до­ бавок

Р(х,

У,

У, 0) = р(х,

у,

|х),

Q(x,

У, I

У, 0) = д(х,

у,

|х).

В дальнейшем мы часто будем рассматривать добавки, линейно зависящие от р, т. е. будем наряду с данной системой (А) рас­ сматривать измененную систему вида

dxldt = P(x, у) + цр(х, у), dyfdt = Q(х, y)+ \iq(x, у).

Отметим еще один встречающийся в дальнейшем частный слу­ чай измененной системы, именно

 

dxldt = P (x, у)±|х<?(х, у),

dy/dt — Q(х, у)Ч=цР(х, у).

Нетрудно видеть, пользуясь формулой (2)

или выражением для

tg 0 ,

которое нетрудно получить,

что эта

система

дает поворот

поля системы (А) на постоянный угол,

тангенс

которого ра­

вен

±|х.

 

 

 

В этом частном случае измененной системы ее состояния рав­ новесия совпадают с состояниями равновесия системы (А) (хотя характер их может быть отличен от характера состояний равно­ весия системы (А)).

Действительно, нетрудно видеть, что при любом |х мы можем иметь одновременно

Р(х, y)± \iQ (x, у) = 0, <?(х, у)+ цР (х, у) = 0,

лишь когда одновременно

Р{х, у) = 0 и Q(x, у) = 0.

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

§ 4. Основные теоремы о зависимости решения от изменения правых частей динамической системы5). В настоящем параграфе излагаются основные теоремы, касающиеся изменения решения системы дифференциальных уравнений, рассматриваемого на ко­ нечном промежутке значений t, при изменении правых частей системы. На эти теоремы опирается все дальнейшее изложение.

6) См. [ИЗ, 116, 130, 134].

134 ДИНАМИЧЕСКИЕ СИСТЕМЫ ТЕОРИИ КОЛЕБАНИЙ [ГЛ. 7

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

логической) структуры разбиения

на траектории

в целом6).

Пусть динамическая

система

(А)

определена в

некоторой

ограниченной замкнутой

области

G и

наряду с

ней

рассматри­

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

Т е о р е м а 1 (о н е п р е р ы в н о й з а в и с и м о с т и р е ш е ­ н и я от и з м е н е н и я п р а в о й ч а с т и и н а ч а л ь н ы х у с л о ­ вий). Пусть

x = y ( t — t'); х*,у*),

У= — t0; х*, Ро)

решение системы (А), определенное при всех значениях

т< t < T ,

и f1 , t2какие-нибудь числа между х и Т, удовлетворяющие неравенству t i < t 0< ti.

Тогда при любом е > 0 существует б > 0 такое, что при усло­

вии, что в области G

 

\ р ( х , у ) \ < б ,

\ q ( x , y ) l < 6

и, кроме того,

 

U o — ^ о | < б .

\Уо— У о \ < Ь ,

решение системы (А), соответствующее начальным значениям

^о> хо 1 Уо:

х= 4>(t — t0,x*, yo),

у= ^ { t — t0, x* , y* ),

определено при всех значениях

t, fi < К

и

при всех этих

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

 

 

| Ф(1 — t0,x*, У*) — ф(* — t0’xo,y*) I <

е,

| $ (t — i0. *о*. Уо) —

(t t0, х*0, у*) | <

е.

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

dx/dx = P(x, у, р), dyldx = Q(x, у, р),

6) Это аналогично тому, как знание структуры траекторий в малом в окрестности отдельной неособой точки не позволяет судить о качественной структуре траекторий в целом (см. § 13 гл. 1).

9 4]

ТЕОРЕМЫ О ЗАВИСИМОСТИ РЕШЕНИЯ

135

а, следовательно, решение этой системы зависит от р,:

x = q > (t-t0, х0, уо,

М-)»

у = ф (£ -£ 0, хо,

Уо,

М-),

то функции ф(£ — to, Хо, уо, р) И

(£ — to, ХО, Уо, р) являются

непрерывными функциями ц.

 

 

Теорема 1 может быть сформулирована в следующей геомет­ рической форме:

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

Наряду с теоремой 1 основную роль в

дальнейшем

играет

также следующая теорема, уточняющая по

сравнению

с

теоре­

мой 1 характер близости решений систем (А) и (А^

в

случае,

когда близки не только правые части систем (А)

и (А),

но и их

частные производные до порядка к.

 

 

 

Пусть по-прежнему решение системы (А) определено при

значениях f: fi

t ^ h.

 

 

 

 

б > 0 такое, что

Т е о р е м а 2.

Д ля всякого г > 0 существует

если в области G выполняются неравенства

 

 

 

\р{х,

у)|

< 6 ,

\q(x,

у ) 1 < б ,

 

 

 

4 - * о

| < 6 ,

| г / о —

* / о * | < 6 ,

 

 

дпр

дпд

 

 

 

 

 

 

дх*дуп~1 < 6,

dxidyn~ i

< б,

п = 1 , 2 , . . . , k;

i = 1 , 2 ,

, п,

то решение системы (А)

 

 

 

 

 

 

x = q ( t — t0, Хо ,Уо),

у = ф (t — t0, Хо, у*)

 

определено при всех значениях t, t\ < t ^ £г, и при всех этих значениях выполняются неравенства

йПф ( г -

го’ х * ’ 2'о)

5Пф

( г -

г0. V ' i ' o )

■дхЪдП~'у0

 

дх'дп~1у0

°п* ( г

- г о, х;,у*0)

а ” ф

( ^

- г 0, * 0*, у *

дх10дп~ \

Предположим теперь, что правые части рассматриваемой ди­ намической системы содержат параметр р, так что система име­ ет вид

dxldt = Р(х, у, р ) , dyldt = Q(x, у, ц).

136

ДИНАМИЧЕСКИЕ СИСТЕМЫ ТЕОРИИ КОЛЕБАНИЙ

(ГЛ. 7

Предположим, кроме того, что функции Р(х, у, ц), Q(x, у,

р) —

аналитические функции также и параметра р.

 

 

Т е о р е м а 3. Если Р(х, у, р),

Q(x, у, р) — аналитические

функции своих аргументов, то и функции

 

 

s = <p(f-f0, so, у0, р),

 

 

у = ty(t — to, so,

уо, р)

 

также являются аналитическими в окрестности всякой системы значений t — to, хо, уо, р, для которой они определены1).

§ 5. Грубость динамической системы и теоремы о непрерыв­ ной зависимости решения от изменения правых частей. На осно­ вании приведенных теорем мы можем утверждать, что на любом конечном замкнутом промежутке значений (на котором опреде­ лено решение исходной системы) при малых изменениях правых частей решение измененной системы мало отличается от реше­ ния исходной системы.

Однако на основании этих теорем нельзя сделать никаких заключений о неизменности поведения траектории на неограни­ ченном интервале значений t и тем более о неизменности харак­ тера разбиения на траектории в целом.

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

dx/dt = 2 s, dy/dt = у,

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

С другой стороны, рассмотрим систему

 

dx/dt = у, dy/dt = —s,

(3)

у которой все траектории замкнуты (начало координат является состоянием равновесия типа «центр»). Рассмотрим наряду с этой системой измененную систему

dx/dt = цх — у, dy/dt ~ \ху — х.

(4)

Все траектории этой системы, кроме состояния равновесия,— спирали (состояние равновесия 0 (0 , 0 ) есть фокус).

7) Из теоремы 3 могут быть, как следствие, получены как утверждения теоремы 1 , так и утверждения теоремы 2 .

S 5]

ГРУБОСТЬ ДИНАМИЧЕСКОЙ СИСТЕМЫ

137

Хотя в силу теоремы 1 на любом конечном промежутке зна­ чений t при достаточно малом р витки спирали системы (4) сколь угодно близки к соответствующей замкнутой траектории системы (3), но очевидно, что при сколь угодно малых р =И=О топологические структуры разбиений у систем (3) и (4) раз­ личны. Таким образом, требование неизменности всей качествен­ ной картины траекторий при малых изменениях правых частей в целом непосредственно не вытекает из приведенных теорем о непрерывной зависимости решения от изменения правых ча­ стей и требует специального рассмотрения. Мы проведем это рассмотрение в следующей главе.

Г Л А В А 8

ГРУБЫЕ ДИНАМИЧЕСКИЕ СИСТЕМЫ

§ 1. Определение грубой динамической системы. Мы будем предполагать, что у всех рассматриваемых динамических систем

dx/dt = P(x, у), <dyldt = Q(x,

у)

(А)

правые части определены в некоторой области

W

плоскости (х, у)

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

Будем наряду с данной фиксированной системой (А) рас­

сматривать всевозможные измененные системы

 

dx/dt = Р(х, у), dy/dt = Q(x, у),

(А)

правые части которых также определены и апалитичны в обла­

сти

W.

 

считать

измененную

систему

{А)

б л и з к о й к систе­

ме

Будем

(А)

в

замкнутой

области

G, целиком

(вместе

с

границей)

лежащей в W (в которой определены системы

(А)

и

(А)), если

в каждой

точке

М(х,

у)

замкнутой

области G

не

только функ­

ции Р(х, у) и Q(x, у)

близки соответственно к функциям Р(х,у)

и Q(x,

у

но

и

первые

производные

от функций Р(х, у) и

 

 

]' (х,у)

и

Ру(х,у), Qx(x,y)

и

Qy(x,y) соответственно

близки к производным от функций Р(х,

у)

и

Q(x,y): Рх (х,у) и

(Х>у) и Qv(x,y)2)*. В соответствии с этим мы будем говорить, что система (А) мало меняется, если наряду с систе­ мой (А) рассматриваются всевозможные измененные системы (А), близкие к системе (А) в указанном смысле.

Впервые определение грубости динамической системы было дано (см. [2 , 3, 1 1 ]) при некотором дополнительном предполо­ жении относительно множества рассматриваемых динамических

■) Понятие грубой динамической системы [3, 13, 26, 144] имеет смысл также и при значительно более общих предположениях относительно пра­ вых частей (см. [13] и § 8 настоящей главы).

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

§ 1) ОПРЕДЕЛЕНИЕ ГРУБОЙ ДИНАМИЧЕСКОЙ СИСТЕМЫ 1 3 9

систем^Именно, дополнительно предполагалось, что граница об­ ласти G, в которой рассматривается система (А), является цик­ лом без контакта для траекторий этой системы, т. е. простой гладкой замкнутой кривой С, не имеющей контактов (не касаю­ щейся траекторий системы (А)). Очевидно, тогда кривая С яв­ ляется циклом без контакта также и для траекторий всякой системы (А), достаточно близкой к (А). Хотя это предположе­ ние сильно ограничивает класс рассматриваемых динамических систем, но при этом смысл понятия грубости системы сохраня­ ется, а определение грубости значительно прощением при общих предположениях относительно границы области G.

О п р е д е л е н и е

I. Динамическая система (А) называется

грубой

замкнутой

области G, граница которой есть цикл без

контакта),

если для

любого е > 0 можно указать 6 > 0 такое,

что для всевозможных измененных систем (А), правые части которых Р (х, у) и Q(x, у) удовлетворяют в области G условиям

IР (х, У) Р (х, у) | <

б,

| Q (х, у) Q (х, у) | < б,

 

IК

{х, у) — Р'х {х, у) I <

б,

I Ру (х, у) Ру (х, у) I <

б,

I &

(X, y) — Qx(x,y) | <

б,

I Q'y (х, у) Q'y (х, у) I <

б,

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

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

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

Для формулировки этого общего определения грубости введем неко­ торую вспомогательную терминологию. Пусть, как и выше, система (А)

и измененная система (А) рассматриваются в замкнутой ограниченной об­

ласти G. Пусть дано некоторое б >

0.

_

 

системе (А),

1. Измененная система (А) называется б- б л и з к о й в G

к

если во всех точках области G выполняются неравенства

 

 

I Р (х , у) Р (*, у) I <

б;

| Q (х, у) Q (х, у) | <

6;

 

| р 'х (х >у ) —

р 'х (*. у) | <

б;

| Р у (*> у) — Р у (*, у ) | <

б;

| Q'x (*> у) —

Q x (*. у ) | <

б;

| Q'y (*. у) — <?«(*. у) | <

в.

Предположим теперь, что система (А) рассматривается в некоторой

замкнутой области Я , целиком вместе с границей лежащей в области G,

а данная измененная система (А) —в некоторой замкнутой области В , также целиком вместе с границей лежащей в области G (так что ни одна

граничная для областей И или И точка не является граничной для об­ ласти G).

140

ГРУБЫЕ ДИНАМИЧЕСКИЕ СИСТЕМЫ

[ГЛ. 8

2. Разбиение области Н на траектории системы (А) называется е-тож-

дестеенным разбиению области Н на траектории системы (А), если сущест­

вует топологическое отображение замкнутых областей Н и Н, при котором

траектории систем (А) и (А) отображаются друг в друга, и при этом со­ ответствующие точки находятся на расстоянии, меньшем е (при этом вся­

кая точка области Н находится в е-окрестности некоторой точки области Н). Рассматривая исходную систему (А) и измененные системы (А), оп­ ределенные в замкнутой области G, мы будем говорить о грубости систе­ мы (А) — по самому смыслу этого понятия — не во всей области G, а в не­

которой (произвольной) замкнутой области Go, целиком содержащейся в открытой области G.

Будем при этом предполагать, что граница области G0 является про­ стой замкнутой кривой3) (но теперь уже эта граница может и не быть

целиком без контакта).

 

называется грубой

О п р е д е л е н и е _Г. Динамическая система (А)

в замкнутой

области Go с: G, если существует замкнутая область Н, цели­

ком содержащаяся в G (HczG)

и целиком содержащая G0

(G0czH ), в ко­

торой выполняются следующие

условия: при любом е > 0

можно указать

6 > 0 такое,

что, какую бы систему (А), б-близкую £

области G к систе­

ме (А), мы

ни взяли, существует замкнутая область Н cz G,

разбиение_ко-

тороп на траектории системы (А) е-тождественно разбиению области Н на траектории системы (А).

В приведенном определении может вызвать недоумение рассмотрение

вспомогательных областей Н жН. Непосредственно представляется естест­ венным следующее определение: система называется грубой в замкнутой

области G0 cz G, если при любом е > 0 можно указать б > 0 такое, что, ка­ кую бы систему (А), б-близкухо_в области G к системе (А), мы ни взяли,

существует замкнутая область G0, разбиение которой на траектории систе­

мы (А) е-тождественно разбиению области G0 на траектории системы (А). Однако нетрудно видеть, пользуясь введенным ниже понятием грубой и негрубой траекторий, что это определение не запрещает наличия не­

грубых траекторий (негрубых состояний равновесия, негрубых предельных

циклов), лежащих на границе области G0. А это, очевидно, не соответст­ вует содержанию понятия грубости. Данное в тексте определение с введе­

нием вспомогательных областей Н и Н выделяет системы, полностью адек­ ватные интуитивному понятию грубой системы.

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

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

(А) в е-окрестности каждого состояния равновесия системы (А) будет лежать одно и только одно состояние равновесия и при этом

3) Приводимое ниже определение грубости динамической системы не изменится, если_сделать и более общие предположения относительно гра­

ницы области G0, однако для определенности мы останавливаемся на сде­ ланном в тексте.

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