Портал | Содержание | О нас | Авторам | Новости | Первая десятка | Дискуссионный клуб | Чат Научный форум --> Первая десятка "Русского переплета" Темы дня:

Президенту Путину о создании Института Истории Русского Народа. |Нас посетило 40 млн. человек | Чем занимались русские 4000 лет назад?

| Кому давать гранты или сколько в России молодых ученых?

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

В. С. АНИЩЕНКО

Саратовский государственный университет им. Н.Г. Чернышевского

ВВЕДЕНИЕ

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

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

Еще более важной проблемой является анализ устойчивости сложных, многокомпонентных систем. Наблюдая за эволюцией живой и неживой природы, мы можем подметить одно интересное свойство: развитие той или иной сложной системы всегда сопровождается потерей устойчивости некоторыми режимами ее функционирования и рождением новых, устойчивых. Одни структуры гибнут, рождаются новые, которые видоизменяются, совершенствуются и затем вновь уступают место новым. Изменения могут накапливаться плавно, а могут происходить скачком в виде катастроф. Формирование новых структур всегда сопровождается потерей устойчивости (даже разрушением) предшествующих. И здесь скрыта важная проблема - проблема перехода системы из одного режима функционирования в другой режим, отличающийся принципиально. Предшествующий режим потерял устойчивость. Но что при этом происходит? Система выбирает новый устойчивый режим, который может наследовать некоторые свойства предыдущего, а может быть и резко отличным. В таких случаях говорят о бифуркациях динамических систем.

Приведенные рассуждения являются качественными и приобретают вполне определенный смысл лишь в том случае, когда удается перевести их на формальный язык математики. Основы строгой математической теории устойчивости были заложены в трудах русского математика А.М. Ляпунова около ста лет назад. Развитие качественной теории и теории бифуркаций динамических систем связано с именами российских ученых А.А. Андронова, В.И. Арнольда и их учеников.

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

ЛИНЕЙНЫЙ АНАЛИЗ УСТОЙЧИВОСТИ

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

где x(t) - переменная состояния, F - некоторая функция состояния, характеризующая закон эволюции, m - параметр системы. Если задано начальное состояние x(t0), то существует единственное решение уравнения (1), которое предсказывает будущее состояние x(t) для любых t > t0 . Если число переменных состояния равно двум (или более), то моделью будет система двух (или более) уравнений:

Число параметров также может быть больше, чем один.

В связи с тем что проблема устойчивости связана с анализом реакции системы на малое возмущение ее состояния, на первом этапе она может быть исследована в рамках линейного приближения. Поясним это. Пусть x0(t) есть некоторое частное решение уравнения (1). Устойчивость этого решения (состояния) мы хотим исследовать. Введем в рассмотрение переменную y(t), которая задает малое отклонение от частного решения:

y(t) = x(t) - x0(t)

(здесь x(t) - возмущенное решение).

Наша задача состоит в исследовании эволюции во времени малого возмущения y(t), которая подчиняется уравнению (1). Разложим функцию F в степенной ряд в окрестности частного решения x0(t):

Производные функции F должны вычисляться в точках, соответствующих частному решению.

Перепишем уравнение (1) для возмущения y(t) с учетом (4):

где

Слагаемые F(y) включают все члены с yn (n $ 2), то есть учитывают нелинейные добавки. По определению, переменная y(t) есть малое отклонение от частного решения. Поэтому нелинейными членами в уравнении (5) в первом приближении можно пренебречь.

Таким образом, для эволюции малого возмущения мы получаем линейное уравнение

Рассмотрим пример. Пусть динамическая система задана уравнением

, a > 0, b > 0.

Найдем стационарные состояния этой системы и исследуем их устойчивость. В стационарном состоянии изменений во времени нет, значит, и мы получаем

Рассмотрим уравнение для возмущений (7) применительно к первому стационарному состоянию

Решением уравнения (10) будет y = exp (lt). Возмущение y экспоненциально затухает во времени (l есть отрицательное число). Это означает, что состояние устойчиво. Так как второе состояние отличается от первого только знаком, то решение уравнения (10) в этом случае будет экспоненциально нарастающим во времени. Стационарное состояние неустойчиво.

Достаточно простая идея предсказания устойчивости по линейному приближению оказалась весьма плодотворной. Используя математический формализм, можно обобщить результат (7) на случай двух и более переменных состояния. Например, в случае N = 2 уравнение (7) примет вид

где

Если одномерное уравнение (1) описывает эволюцию исключительно в окрестности стационарных состояний, то уравнение (2) может иметь в качестве решения не только стационарные, но и периодические решения. С увеличением размерности исходной системы (1) в общем случае усложняются и типы возможных решений. Это создает определенные проблемы в исследовании устойчивости: ведь для решения уравнений для возмущений типа (7) необходимо знать частное решение x0(t)! С применением современных компьютеров эти трудности легко преодолимы.

БИФУРКАЦИИ ДИНАМИЧЕСКИХ СИСТЕМ

Если внимательно присмотреться к окружающей нас природе, то можно сделать следующее интересное наблюдение. Жизнь на планете Земля возможна лишь благодаря тепловому излучению Солнца, которое служит источником энергии. Летом северное полушарие получает световой энергии больше, чем зимой. И картина летней природы при этом заметно отличается от зимней. Давайте рассмотрим в качестве примера объем воды в озере. Количественной мерой привносимой солнечной энергии является температура воды. Летом вода в озере теплая и можно купаться. С наступлением осени температура воды постепенно уменьшается. Купаться уже не хочется, однако вода и при более низкой, но плюсовой температуре остается водой! Глубокой осенью верхний слой воды в озере остывает до нулевой температуры и вдруг превращается в лед! Далее и при - 20?C лед остается льдом. Что же произошло? При прохождении температуры через нуль вода резко изменила свои свойства: она из жидкого состояния перешла в твердое. И не плавно, а скачком.

Если рассматривать температуру воды как некий параметр системы, то хорошо известно, что с изменением параметра вода резко меняет свои свойства при переходе через 0?С, через 100?C, когда вода превращается в пар. Есть и другие особые значения температуры воды. Оказывается, что большинство интересных физических задач при их математическом описании приводят к дифференциальным уравнениям, зависящим от одного или нескольких параметров.

Рассмотрим в качестве другого примера уравнения колебаний обыкновенного маятника или (что с математической точки зрения полностью идентично) параллельного RLC-контура:

Уравнение (12) содержит два параметра: a - параметр затухания, характеризующий трение, и w0 - параметр, определяющий частоту колебаний. Если потери энергии отсутствуют (параметр затухания a = 0), то решением уравнения (12) будут гармонические незатухающие колебания. При малом трении 0 < a < 1 движение системы будет колебательным с амплитудой, которая уменьшается во времени по экспоненциальному закону. Наконец, при достаточно большом трении (a > 1) движение маятника будет апериодическим, затухающим во времени. Уже в этом простом примере выделяются два особых значения параметра a = 0 и a = 1, отклонения от которых качественно меняют свойства системы.

Изменение параметра в уравнении может вызвать потерю устойчивости одного состояния (или режима функционирования) системы и переход ее в другое, отличное от первого состояние. Это явление, как сказано выше, называется бифуркацией (от слова раздвоение), а значение параметра, при котором оно происходит, - точкой бифуркации. Состояние системы ниже точки бифуркации и выше ее при изменении параметра все-таки меняется. Ясно, что вода при температуре + 3 и + 22?C - это разные состояния. Но при этом вода остается водой! В математике и физике существует понятие грубости или структурной устойчивости. Суть этого понятия в том, что при малом изменении параметра грубая система хоть и изменяет в деталях режим функционирования, но не принципиально. С этой точки зрения для грубых систем переход через точку бифуркации означает смену одного структурно устойчивого режима на другой. При этом в точке бифуркации система не является грубой: малое изменение параметра в ту или иную сторону приводит к резким изменениям состояния.

Давайте вернемся к нашему примеру с устойчивостью стационарных состояний в системе (8). Мы условились, что в уравнении (8) параметры a и b положительны. Устойчивость определяется знаком производной правой части уравнения (8) в стационарной точке, то есть знаком величины l (10). При положительных значениях параметров a и b эта производная всегда отлична от нуля. А что, если мы будем уменьшать значение параметра a? Как видно из (10), при a = 0 (независимо от величины b > 0) величина l обращается в нуль, возмущение y не нарастает и не затухает. Мало того, при a = 0 в системе (2) оба стационарных состояния как бы сливаются в одно (x = 0)! Далее, если a < 0, то стационарных состояний нет вовсе. Действительно, в этом случае (j - мнимая единица), то есть при a < 0 становятся чисто мнимыми.

Приведем теперь результаты математического анализа этой бифуркации, которая известна как бифуркация "двукратное равновесие". Вновь рассмотрим уравнение (8). Пусть x0(a) есть грубое состояние равновесия, то есть l(a) ? 0. Это означает, что при малой вариации параметра a состояние x0(a) продолжает существовать как устойчивое или неустойчивое стационарное состояние.

При некотором значении параметра a = a* собственное число l(a*) в стационарном состоянии может обратиться в нуль:

Для реализации бифуркации "двукратное равновесие" необходимо, чтобы вторая производная была отлична от нуля:

Для выполнения условий (13) и (14) в общем случае необходимо, чтобы исходное уравнение в правой части включало как минимум квадратичное нелинейное слагаемое, как в нашем примере (8).

Если условия (13) и (14) выполнены, то x0 есть двукратный корень исходного уравнения (8).

Значение параметра a*, при котором выполняется условие (13), является точкой бифуркации. До точки бифуркации a < a* мы имеем два стационарных состояния. В точке бифуркации a = a* они сливаются в одно, далее при a > a* этих состояний в системе не будет. В нашем случае (8) a* = 0.

Результаты можно представить графически (рис. 1).

"МЯГКИЕ" И "ЖЕСТКИЕ" БИФУРКАЦИИ. КАТАСТРОФЫ

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

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

Другой пример бифуркации качественно представлен на рис. 3. При m < m* (рис. 3, 1 ) шарик находится в устойчивом стационарном состоянии. При этом существует еще одно, неустойчивое состояние (помечено звездочкой на рис. 3, 1 ). В точке бифуркации m = m* устойчивое и неустойчивое состояния сливаются в одно (рис. 3, 2 ). Далее они исчезают, и система выбирает новый режим (например, как это показано на рис. 3, 3 ), который существенно отличается от предыдущего и не находится в непосредственной близости от исходного режима. Такой тип бифуркаций называют жестким, и именно жесткие бифуркации явились предметом анализа в теории катастроф.

Рассмотренный выше пример бифуркации "двукратное равновесие" в системе (8) представляет собой типичный пример жесткой бифуркации, который качественно проиллюстрирован на рис. 3.

ЗАКЛЮЧЕНИЕ

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

ЛИТЕРАТУРА

1. Арнольд В.И. Теория катастроф. М.: Изд-во МГУ, 1983.

2. Постон Т., Стюард Я. Теория катастроф и ее приложения. М.: Мир, 1977.

3. Анищенко В.С. Сложные колебания в простых системах. М.: Наука, 1990.

Рецензент статьи Б.С. Бокштейн

* * *

Вадим Семенович Анищенко, доктор физико-математических наук, профессор, зав. кафедрой радиофизики и нелинейной динамики Саратовского государственного университета, заслуженный деятель науки Российской Федерации, член-корреспондент РАЕН и Международной академии информатизации. Лауреат премии по физике Фонда им. А. Гумбольдта. Область научных интересов - нелинейная динамика и статистическая радиофизика. Автор более 250 научных работ, семи научных монографий.