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

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

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

§ 41 «БЕЗОПАСНЫЕ» И «ОПАСНЫЕ» ГРАНИЦЫ 231

Рассмотренная граница области устойчивости, очевидно, яв­ ляется опасной.

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

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

Рассмотрим еще дополнительно поведение динамических си­ стем вблизи тех точек границы, в которых А > 0 и R = 0, где безопасная граница переходит в опасную, т. е. где первая ляпуновская величина L\ обращается в нуль. В этом случае поведе­ ние системы может быть определено знаком второй ляпуновской величины Z/2 = as ^ 0 (см. гл. 11, § 5).

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

которой

точке

Д/(Яо, ро) плоскости параметров R = 0, L\ = 0,

но Z/2 ^

0. На

рис.

120 большая точка

соответствует

точке пло­

скости

параметров,

в которой R = 0,

Ь\ =

0, Z/2 ^

0;

в точках

части линии R = 0,

обозначенной белыми

точками,

L\ > 0 (для

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

В заштрихованной части плоскости параметров R < 0, в незаштрихованной (но в которой может быть область, обозначен­ ная мелкими штрихами) R > 0. Напомним, что при переходе че­ рез часть границы R = 0, в которой L\ < 0, из заштрихованной

области

в

незаштрихованную

из

сложного

фокуса

рождается

устойчивый

предельный цикл,

а

при

переходе через

R = 0, где

L\ > 0,

из

незаштрихованной

в

заштрихованную область — не­

устойчивый цикл.

которой R = Ь\ = 0, мы

имеем

Li > 0, то­

Пусть в

точке, в

гда соответствующий

сложный

фокус

(второй степени

негрубо-

_сти) неустойчив (рис. 120,а, 121,а).

мы перейдем

по

линии

Если

в

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

R = 0 в

точки, где L\ < 0, то,

как нетрудно

показать, рассмат­

232

НЕЛИНЕЙНЫЕ ЯВЛЕНИЯ И ТЕОРИЯ БИФУРКАЦИЙ

[ГЛ. 13

ривая функцию последования (см. гл. 8, § 3), на фазовой пло­ скости из сложного фокуса второго порядка родится неустойчи­ вый предельный цикл (грубый), а фокус делается негрубым устойчивым (рис. 121,6). Если затем в пространстве параметров

Случай 1*2>С

Сличай Lz<C

Рис. 120

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

Случай 1*2 >0

Рис. 121

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

8 4] «БЕЗОПАСНЫЕ» И «ОПАСНЫЕ» ГРАНИЦЫ 233

С другой стороны, нетрудно показать, что при значениях пара­

метров в

области I I I

у

системы

вокруг

неустойчивого

фокуса

нет

предельных

циклов.

(При

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

где R = Ьх =

= О,

Ь2 >

0, на

часть

линии

R =

0, где

Ьх > 0,

циклы

не ро­

ждаются, и в силу сделанных выше замечаний не рождаются при переходе в область III.) Но тоща при движении в про­ странстве параметров из области III в область II непременно должны встретиться бифуркационные значения параметров, при которых у системы существует двукратный цикл. На рис. 120 линия в пространстве параметров, соответствующая двукратным предельным циклам, изображена штрихами. Аналогичное рас­ смотрение может быть проведено и в случае Ь2 < 0 (рис. 120,6 и рис. 122, а в).

В рассмотренном случае знак второй ляпуновской величины Ь2 играет роль, подобную знаку L\, увеличивая или уменьшая

Случай Lz<0

Рис. 122

опасность для изображающей точки быть выброшенной из ок­ рестности состояния равновесия.

Пусть Ь2> 0 (при R = L\ — 0) и пусть значения парамет­ ров достаточно близки к значениям, определяемым этими усло­ виями; тогда в достаточно малой окрестности начала координат в фазовом пространстве может быть одна из структур, изобра­ женных на рис. 121.

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

234 НЕЛИНЕЙНЫЕ ЯВЛЕНИЯ И ТЕОРИЯ БИФУРКАЦИЙ [ГЛ. 13

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

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

Пусть Z/2 < 0 (при R = L i = 0 ) , и пусть параметры опять мало изменены; тогда в достаточно малой окрестности начала координат может быть одна из структур, изображенных на рис. 122. Здесь даже нарушение опасной границы может оста­ вить изображающую точку в малой окрестности состояния рав­ новесия, если параметры достаточно близки к значениям, опре­ деляемым условиями R — L\={).

§ 5. Замечания по поводу других границ области устойчиво­ сти. Мы рассмотрели части границы области устойчивости в случае, когда эти части соответствуют негрубым состояниям рав­ новесия. Очевидно, аналогично могут быть рассмотрены гра­ ницы области устойчивости, соответствующие еще и другим си­

стемам

первой

степени негрубости п ) ; именно в нумерации

§ 6

гл. 9:

3)— системе, имеющей двукратный

предельный цикл;

4)— системе,

имеющей сепаратрису, идущую

из седла

в

седло.

 

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

место,

когда при

из­

Случай 3)

менении параметра к устойчивому предельному циклу прибли­ жается неустойчивый предельный цикл (как в рассмотренном в § 3 случае жесткого возбуждения колебаний), который сли­ вается с устойчивым циклом. Образуется двукратный цикл, ко­ торый при дальнейшем изменении параметра исчезает. Изобра­ жающая точка «срывается» и стремится либо к устойчивому со­ стоянию равновесия, либо к другому устойчивому предельному циклу.

Граница опасная. Нетрудно убедиться, что система при любой

еекачественной структуре в этом случае ведет себя необратимо.

Вслучае 4) предположим, что изображающая точка при зна­

чениях А, < А.0

двигается по

устойчивому предельному

циклу,

который при

К = А,о влипает

в сепаратрису. Очевидно,

по мере

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

13)Напоминаем, что только в случае систем первой степени негрубо­

сти в окрестности точки границы о б л а с т и в е т с т в у ю щ е м ф у н к ц и о н а л ь н о м б ы е с и с т е м ы о б р а з у ю т п л е н к у .

у с т о й ч и в о с т и в с о о т ­ п р о с т р а н с т в е не г ру -

§ 5] ЗАМЕЧАНИЯ ПО ПОВОДУ ДРУГИХ ГРАНИЦ 235

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

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

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

на кото­

рый

навивается при

t -*■ 00 несколько сепаратрис,

входящих

в границы ячеек с различными центрами притяжения,

устойчи­

выми

состояниями

равновесия или

предельными

циклами

(рис.

123,а), при возрастании параметра

стягивается к

устойчи­

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

а

5

Рис.

123

в состояние равновесия, которое делается сложным неустойчи­ вым фокусом, а затем грубым неустойчивым фокусом. При А > Ко сепаратрисы стремятся к этому фокусу, и, очевидно, нет

возможности

однозначно

указать, к центру

притяжения какой

из ячеек

(для всех

этих

ячеек фокус

является

теперь гранич­

ным) будет

стремиться изображающая

точка

(см.

[141] ) 14).

на­

Полностью аналогичная

ситуация

имеет

место также,

пример,

в

случае,

когда

изображающая точка

двигается

по

14) Таким образом, можно сказать, что мы имеем здесь «динамическую неопределенность».

236

НЕЛИНЕЙНЫЕ ЯВЛЕНИЯ И ТЕОРИЯ БИФУРКАЦИЙ

[ГЛ. 13

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

К циклу LQ приближается неустойчивый предельный цикл L, являющийся предельным для ряда сепаратрис, входящих в гра­ ницы ячеек с различными центрами притяжения. При некото­ ром значении параметра Я = Яо цикл L\ сливается с циклом LQ

иисчезает.

Врассматриваемом случае, так же как и в предыдущем, нет возможности однозначно указать поведение изображающей точ­

ки при %> Яо.

Ч А С Т Ь III

КАЧЕСТВЕННОЕ ИССЛЕДОВАНИЕ КОНКРЕТНЫХ ДИНАМИЧЕСКИХ СИСТЕМ

С АНАЛИТИЧЕСКИМИ НРАВЫМИ ЧАСТЯМИ

Г Л А В А 14

ОБЩИЕ ЗАМЕЧАНИЯ О ПРИЕМАХ КАЧЕСТВЕННОГО ИССЛЕДОВАНИЯ

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

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

Сделаем прежде всего некоторые общие замечания.

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

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

с аналитическими

(или неаналитическими) правыми

частями.

Постараемся пояснить бессмысленность таких попыток на

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

в с е в о з ­

м о ж н ы е

аналитические на некотором промежутке значений х

функции

y = f(x)

и ставится вопрос об общем универсальном

методе отыскания

(разделения) корней любой из функций

fo(x)=0.

Предположим, что рассматриваются функции f(x), аналити» ческие при всех х (—° ° < а ; < + <»), и ставится вопрос об общем универсальном методе определения числа корней любой из этих функций /о (х) — 0 на некотором конечном интервале значений х

238

 

ОБЩИЕ ЗАМЕЧАНИЯ О ПРИЕМАХ ИССЛЕДОВАНИЯ

[ГЛ. 14

(а < х < Ь ) 1).

По

сделанному

предположению

относительно

функций f(x) число

их

на конечном интервале а < х < Ь

обяза­

тельно

конечно

(но

всегда, очевидно, можно указать функцию

f (x),

у

которой

на

этом интервале любое данное число кор­

ней).

Однако функции

рассматриваемого класса

столь разнооб­

разны

 

(они зависят

от

счетного

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

напри­

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

Если бы мы обратились к вычислительным методам, то мог­ ли бы «выловить» корни, находящиеся друг от друга на расстоя­ нии, не меньшем некоторого расстояния, допускаемого точностью вычислений. Между тем в силу широты рассматриваемого класса функций среди них всегда найдется функция, у которой корни находятся друг от друга на меньшем расстоянии и количество их может быть равно любому данному числу на рассматривае­ мом промежутке.

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

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

-Р* ( х) = О,

то, как известно, для решения этой задачи существует регуляр­ ный алгоритмический метод — классический метод Штурма.

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

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

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

ВВЕДЕНИЕ

239

например для случая, когда правые части — многочлены данной фиксированной степени п, представляется имеющей смысл, од­

нако, конечно, очень далекой от решения

(такой метод

опреде­

ления числа предельных циклов был бы

в некотором

смысле

аналогичен методу Штурма).

динамических

систем,

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

правые части которых не обязательно многочлены, но зависят от конечного числа параметров.

Из сказанного выше очевидно, что не только для задачи определения числа и расположения предельных циклов, но даже для значительно более простой задачи — задачи определения числа состояний равновесия, которая сводится к определению числа общих корней пары функций

Р ( х , у ) = 0, Q(x, у) =0 ,

также можно сделать полностью аналогичные высказывания.

Вслучае, когда Р(х, у) и Q(x, у ) — многочлены данной фик­ сированной степени п, эта задача при использовании результан­ та этих многочленов, очевидно, сводится к методу Штурма.

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

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

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

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

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

240 ОБЩИЕ ЗАМЕЧАНИЯ О ПРИЕМАХ ИССЛЕДОВАНИЯ [ГЛ. 14

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

двух

параметров — бифуркационными

кривыми)

на

области

с одинаковым качественным поведением

фазовых

траекторий и

при

этом, конечно, в установлении этого

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

поведе­

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

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

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

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

метров

одновременно и

помогает

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

исследованию,

и дает

разделение на

области с

различными

качественными

структурами.

 

 

 

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

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

Один из основных вопросов качественного исследования — вопрос отыскания предельных циклов — в некоторых приклад­

ных задачах иногда удается решать

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

классическим методом исследования

нелинейных систем — ме­

т о д о м м а л о г о п а р а м е т р а .

 

Очевидно, этот метод тоже в каком-то смысле можно считать методом теории бифуркаций, так как в этом методе фактически рассматривается бифуркация от линейной (нелинейной) консер­

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