книги / Методы и приёмы качественного исследования динамических систем на плоскости
..pdf§ 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, однако для определенности мы останавливаемся на сде ланном в тексте.