|
|
Статьи Соросовского Образовательного журнала в текстовом формате
Рассматриваются основные вопросы качественной теории двумерных динамических систем в элементарном изложении. Приводятся основные сведения о циклах одномерных отображений. Цель статьи - вызвать интерес ушкольников и учителей к этой сложной математической теории.
ЭЛЕМЕНТАРНОЕ ВВЕДЕНИЕВ КАЧЕСТВЕННУЮ ТЕОРИЮ
И ТЕОРИЮ БИФУРКАЦИЙ ДИНАМИЧЕСКИХ СИСТЕМ
В. Н. БЕЛЫХ
Волжская государственная академия водного транспорта, Нижний Новгород
ВВЕДЕНИЕ
Наша конечная цель - познакомиться с основами математической теории хаоса в динамических системах без случайных воздействий. Явление динамического хаоса в настоящее время обнаруживают практически во всех областях знаний, где возможно математическое моделирование, ему посвящено огромное число научных статей, книг и отдельных международных журналов, таких, как "Chaos". Ясно, что представить сложную математическую теорию, не углубляясь в математику, даже в самом элементарном изложении, - задача не простая, и одной статьи здесь недостаточно. Поэтому разговор о самом явлении хаоса мы будем вести и дальше, а здесь будет представлено введение в качественную теорию и теорию бифуркаций динамических систем, составляющих основу языка, на котором говорят современные специалисты по хаотической динамике в различных областях знаний.
Любые системы, состояние которых изменяется во времени, называют динамическими. В математике под динамическими системами понимают дифференциальные уравнения и отображения.
Это не случайно, поскольку в математике большинство уравнений динамики как фундаментальных - Ньютона и Гамильтона в механике и Максвелла в электродинамике, - так и феноменологических, например Ходжкина-Хаксли в биофизике нейрона, Лотки-Вольтера в экологической задаче хищник-жертва и Леонтьева в модели развивающейся экономики, записываются в форме дифференциальных уравнений, а компьютерные представления последних - в форме разностных уравнений, то есть отображений.
Основы качественной теории и теории бифуркаций динамических систем были заложены в трудах великого французского ученого Анри Пуанкаре, который первым понял, что можно, не интегрируя дифференциальных уравнений, представить все основные качественные особенности поведения его решений.
Мы рассмотрим простейшие динамические системы: системы дифференциальных уравнений второго порядка и одномерные отображения, а затем в качестве предисловия к следующей статье обсудим возможности применения представленной информации к теории хаоса.
К сожалению, даже в популярной форме ответить на вопрос о том, что такое хаос в динамической системе и как он возникает, невозможно без основных понятий из теории дифференциальных уравнений и отображений. Поэтому наберемся терпения и постараемся освоить азы теории динамических систем. Зато потом это окупится: разобравшись в азах, мы сможем понять главное о хаосе.
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ
Начнем с самого простейшего уравнения первого порядка
Слева в этом уравнении - производная по времени (независимой переменной), справа - непрерывная гладкая функция переменной x и постоянного параметра l. Решить уравнение (1) означает найти функцию
x = S (t, x0 , l),
которая при подстановке в (1) обращает это уравнение в тождество. Величина x0 в (2) есть начальное условие, удовлетворяющее соотношению
x0 = S (0, x0 , l).
Константа x = c обращает производную в нуль, и если она одновременно служит решением уравнения
F (x, l) = 0,
то x = c есть решение уравнения (1), называемое состоянием равновесия.
Теперь если F > 0, то x(t) возрастает, а если F < 0, то x(t) убывает в силу (1). В этом вся динамика уравнения (1): корни уравнения (3) - состояния равновесия С1 , С2 , _, а на интервалах (Ci , Ci + 1) происходит либо рост, либо убывание x(t) в зависимости от знака F (скорости). Рис. 1 иллюстрирует изменение x(t), указанное стрелками.
Состояния равновесия C1 и C3 устойчивы, а C2 неустойчиво. На рис. 1 видно, что если производная , то состояние равновесия устойчиво, а если , то неустойчиво.
Следует заметить, что аналитически решения уравнения (1) можно представить в виде неопределенного интеграла
Однако для качественного анализа этот интеграл не нужен, более того, такие интегралы, как правило, не берутся.
Пример 1. Линейное уравнение
имеет решение x(t) = x0elt, определяющее экспоненциальный рост (убывание), если l > 0 (l < 0) соответственно. Вспомним, что уравнение (4) определяет динамику цепной реакции (l > 0) и распада ядра (l < 0). Единственное состояние равновесия уравнения (4) x = 0 устойчиво при l < 0 и неустойчиво при l > 0.
Пример 2. Нелинейное уравнение
имеет два состояния равновесия при l > 0 и не имеет состояний равновесия при l < 0 (dx / dt < 0 всюду и x(t) убывает всегда). Правильность этого утверждения легко проверить, нарисовав график F (x) = l - x2 и расставив стрелки, как на рис. 1. Значение l = 0 в этом примере соответствует бифуркации рождения (гибели) двух состояний равновесия с ростом l (при убывании l) соответственно. Важнейшее понятие - бифуркация - было введено в конце прошлого века Анри Пуанкаре.
Рассмотрим теперь систему дифференциальных уравнений второго порядка
с непрерывными гладкими функциями F1 и F2 , зависящую от параметра l, решения которой есть система двух функций
Плоскость (x, y) для дифференциальных уравнений (6) называют фазовой плоскостью, а график функции (7) как параметрически заданной кривой (параметр - время t) - траекторией. Траектории пересекаться не могут по теореме единственности решения (7) с заданной начальной точкой. Изменяя начальную точку (x0 , y0), получаем семейство траекторий на фазовой плоскости, которое с отмеченным стрелками направлением изменения точки (7) при увеличении времени образует фазовый портрет системы (6). Состояния равновесия системы (6) аналогично уравнению (1) определяются системой алгебраических уравнений
Пример 3. Линейная система дифференциальных уравнений
составленная из двух независимых уравнений первого порядка вида (4), очевидно, имеет решение
x = x0elt, y = y0egt.
Состояние равновесия O(0, 0) может быть трех типов:
а) при l < 0, g < 0 точка O - устойчивый узел;
б) при l > 0, g > 0 точка O - неустойчивый узел;
в) при l < 0, g > 0 (или l > 0, g < 0) состояние равновесия O называют седлом.
Фазовые портреты в этих трех случаях изображены на рис. 2.
Заметим, что состояния равновесия являются отдельными траекториями, так что пересечения, реализующиеся для системы (9) и изображенные на рис. 2, суть не пересечения траекторий, а пересечения так называемых интегральных кривых, составленных из трех траекторий, одна из которых - состояние равновесия.
Пример 4. Рассмотрим систему
заданную в полярных координатах (j, r), связанных с декартовой системой координат формулами x = = r cos j, y = r sin j, где j - угол, r - полярный радиус, l и w - параметры. При этом, очевидно, r $ 0. Снова уравнения в (11) независимы, и мы имеем решение
r = r0elt, j = wt + j0 ,
где (j0 , r0) - начальная точка. Состояние равновесия r = 0 (в декартовых координатах это O(0, 0)) есть устойчивый (неустойчивый) фокус при l < 0 (l > 0) соответственно. Фазовые портреты для этих случаев изображены на рис. 3.
Заметим, что к системе (11) (при l < 0) приводит задача о затухающих колебаниях маятника или колебательного контура (по поводу последнего см., например, мою статью в последнем издании "Физической энциклопедии", т. 2). При этом полярный радиус есть экспоненциально убывающая амплитуда, угол j - полная фаза, а w - частота колебаний.
Пример 5. Нелинейная система в полярных координатах
c разделенными переменными исследуется с помощью рис. 1 с важной особенностью: устойчивое состояние равновесия r = l > 0 первого уравнения соответствует замкнутой траектории системы (13) на плоскости (x, y), называемой предельным циклом. При l <0 dr / dt < 0 и цикл отсутствует (рис. 4).
Предельный цикл был введен А. Пуанкаре, впрочем как и все основные понятия качественной теории динамических систем, предмет которой - не интегрируя систему (6), нарисовать ее фазовый портрет и проследить за его изменениями при изменении параметров (предмет теории бифуркаций). В 30-е годы предельный цикл был интерпретирован как математический образ автоколебаний А.А. Андроновым, основателем Нижегородской школы нелинейной теории колебаний. К уравнениям (13) в простейшем случае приводит задача о колебаниях томпсоновского (лампового) генератора. Возникновению колебаний соответствует бифуркация рождения цикла из состояния равновесия (рис. 4), открытая в 30-е годы Андроновым и обобщенная в начале 40-х годов Хопфом на случай систем дифференциальных уравнений более высокого порядка.
Пример 6. Нелинейная система
при | l | < 1 имеет следующую последовательность фазовых портретов: при l < 0 производная dr / dt отрицательна и система (14) имеет устойчивый фокус, как и в примере 5; при l = 0 dr / dt обращается в нуль при r = 1, то есть в системе происходит рождение предельного цикла. Траектории приближаются к нему снаружи и удаляются от него внутри. Поэтому такой цикл называют полуустойчивым. При l > 0 этот цикл разваливается на два цикла: один - устойчивый снаружи и другой - неустойчивый внутри. Стало быть, имеет место бифуркация рождения двух циклов через появление полуустойчивого цикла. Проанализировать эту бифуркацию можно с помощью исследования первого уравнения в (14) при использовании рис. 1.
Пример 7. Система вида
при | l | < 1 помимо неустойчивого состояния равновесия O(0, 0) имеет два состояния равновесия P1(r1 = 1, j1 = arcsin l) и P2(r2 = 1, j2 = p - arcsin l). Точка P1 есть устойчивый узел, а P2 - седло. При | l | = 1 состояния равновесия P1 и P2 сливаются, образуя так называемый седло-узел, исчезающий при | l | > 1. Поскольку при l > 1 производная dj / dt > 0 (при l < -1 dj / dt < 0), то есть угол j(t) растет, система (15) имеет устойчивый предельный цикл. Последовательность изменения фазовых картин в этом случае изображена на рис. 5. На рисунке видно, что при l = 1 система (15) имеет траекторию G, выходящую из седла-узла и возвращающуюся в него же. Такие траектории Пуанкаре назвал гомоклиническими. Как оказалось впоследствии, такие траектории играют важнейшую роль в теории хаоса в многомерных динамических системах без случайных воздействий.
Таким образом, мы получили бифуркацию рождения предельного цикла в ситуации, когда исчезает седло-узел, имеющий гомоклиническую траекторию.
Пример 8. Пусть нелинейная система (6) в квадрате Q = (| x | # d, | y | # d) совпадает с линейной системой (9) и имеет в нем фазовый портрет, соответствующий седлу (рис. 2, в) при l < 0, g > 0, а вне этого квадрата имеет фрагмент фазового портрета, изображенный на рис. 6.
В этом случае (при некотором l = lG) существует гомоклиническая траектория G , выходящая из седла и возвращающаяся в него же при t ?. Тот факт, что для | l | > g при изменении l, скажем, в сторону l > lG происходит бифуркация рождения устойчивого предельного цикла (рис. 6, б ), был изучен А.А. Андроновым и Е.А. Леонтович в начале 50-х годов.
Мы привели частные примеры, которые с помощью уравнения первого порядка объясняют различные случаи изменения фазовых портретов системы двух уравнений при изменении параметра. Но, оказывается, эти примеры отражают общую ситуацию качественного изменения фазовых картин любой динамической системы (6), поскольку все главные деформации системы (6) общего вида (со сколь угодно сложными гладкими функциями F1 и F2) происходят только так, как в наших примерах. Говоря "главные деформации", мы имеем в виду, что если в системе (6) рождается сразу много состояний равновесия или предельных циклов, то малыми изменениями F1 и F2 это рождение можно разделить на последовательность элементарных бифуркаций, происходящих по схеме, указанной примерами 2, 5-8, то есть путем поочередного рождения пар состояний равновесия (пример 2) и предельных циклов либо через бифуркацию Андронова-Хопфа (пример 5), либо через бифуркацию полуустойчивого цикла (пример 6), либо через бифуркацию гомоклинической траектории седло-узла (пример 7) или седла (пример 8). Других элементарных бифуркаций в диссипативных системах не бывает. Так что наши примеры не являются частными, а системы, рассмотренные в них, как отражающие общие бифуркационные свойства называют нормальными формами.
Приведем теперь несколько важнейших понятий теории динамических систем общего характера.
Первое. Все сказанное относится к так называемым диссипативным системам (исключая, конечно, тривиальные случаи линейных систем типа (1)), траектории которых при t ? не уходят в бесконечность, а стремятся к ограниченному предельному множеству, называемому аттрактором. Легко проанализировать, что таким множеством в случае системы (6) могут быть состояния равновесия и предельные циклы.
Второе общее понятие, важнейшее в теории динамических систем, - это понятие грубости (или структурной устойчивости), введенное в конце 30-х годов А.А. Андроновым и Л.С. Понтрягиным и играющее центральную роль в современной теории динамических систем. Смысл этого понятия состоит в следующем. Если при малом изменении параметров системы (6) или целиком правой части (функций F1 , F2 и их производных) вид фазового портрета топологически остается тем же самым, то система называется грубой. Из примеров следует, что качественные изменения фазовых портретов происходят при бифуркациях состояний равновесия (примеры 2 и 5), предельных циклов (примеры 5-8) и гомоклинических траекторий (примеры 7 и 8).
Третий общий вывод состоит в том, что хаоса, то есть явления, когда решение системы дифференциальных уравнений (6) ведет себя как случайная функция, в системах второго порядка не бывает.
Заметим, что в математической теории систем второго порядка существует множество нерешенных задач. Достаточно упомянуть о 16-й проблеме Гильберта (ее второй части) о числе предельных циклов в зависимости от порядка полиномов в правой части системы (6). По решению этой проблемы в последние годы получено много интересных результатов, но проблема остается нерешенной.
ОДНОМЕРНЫЕ ОТОБРАЖЕНИЯ
Используемые в математике последовательности действительных чисел x(0), x(1), x(2), _, x(n), _, рассматриваемые как функция дискретного аргумента x(n): Z R, могут задаваться уравнением
x(n + 1) = f (x(n)),
связывающим последующее значение x(n + 1), обозначаемое обычно , с предыдущим x(n), обозначаемым просто x. Уравнение (16) называется разностным и, когда под n понимают дискретное время, представляет собой динамическую систему. Последовательность x(0), x(1), x(2), _, x(n), _ называют траекторией, удовлетворяющей начальному условию x(0). Иначе говоря, мы имеем функцию f : R1 R1, отображающую действительную ось в себя, поэтому уравнение (16) называют также одномерным отображением. Действие отображения f геометрически изображают с помощью диаграммы Ламерея (рис. 7), на которой отрезками изображают переход от аргумента x(n) к значению функции x(n + 1) и указывают переход от значения функции к новому значению аргумента (через биссектрису). Итак, мы представили простейшую динамическую систему как объект, доступный любому старшекласснику.
Пример 9. Арифметическая прогрессия задается отображением сдвига
x(n + 1) = x(n) + d,
где начальное условие x(0) есть ее первый член, а константа d - ее разность.
Пример 10. Геометрическая прогрессия со знаменателем q задается отображением сжатия (| q | < 1) или растяжения (| q | > 1) вида
Нуль является неподвижной точкой отображения (18) как решения уравнения x = qx. При | q | < 1 эта точка устойчива: к ней стремятся точки траекторий, а при | q | > 1 она неустойчива: точки траекторий от нее удаляются. Заметим, что нулевой первый член в геометрической прогрессии выбрасывается, поскольку порождает бесконечную последовательность нулей, которая для отображения (18) есть отдельная траектория - неподвижная точка.
Пример 11. Отображение вида
- "сжатие плюс сдвиг" - имеет неподвижную точку
как решение уравнения x = qx + d и имеет такое же семейство траекторий, как и отображение (18), если перенести начало координат на диаграмме Ламерея в точку (x*, x*).
Из приведенных примеров легко понять, что если функция f (x) в (16) монотонно возрастает, то отображение может иметь только неподвижные точки и траектории, стремящиеся к ним или удаляющиеся от них. Если функция f (x) монотонно убывает, то отображение может иметь еще и циклы периода 2, то есть траектории, удовлетворяющие условию x2 = f (x1), x1 = f (x2), x1 ? x2 . Например, в случае (18) при q = -1 все точки, кроме нуля, являются точками циклов периода 2. В общем случае, когда функция f может менять свою монотонность, отображение (16) может иметь циклы любого периода k, то есть траектории, удовлетворяющие условиям
x2 = f (x1), x3 = f (x2), _, xk = f (xk - 1),
xk + 1 = x1 , xi ? x1 , i = 2, 3, _, k
(внутри цикла точки не могут совпадать, иначе получим цикл меньшего периода). Легко, например, получить цикл периода 3 для так называемого логистического отображения
переводящего единичный отрезок [0, 1] в себя при l # 4. Итак, аттрактором невзаимно однозначного отображения может быть предельное множество, включающее различное число циклов, в том числе и бесконечный их набор.
Известный украинский математик А.Н. Шарковский в 1964 году доказал теорему о том, что если отображение (16) имеет цикл периода 3, то оно имеет циклы всех периодов.
Теперь давайте представим себе общую картину поведения траекторий нелинейного отображения. Пусть функция f (x) в (17) удовлетворяет условию
x > f > 0, x > D; x < f < 0, x < - D,
как показано на рис. 7, a. Тогда все траектории попадают в область | x | < D и остаются в ней навсегда. В случае когда в этой области существуют циклы всех периодов (а именно они составляют основу аттрактора, то есть множества, к которому стремятся все траектории отображения), то появляется чувствительность к выбору начальных условий: стоит слегка изменить начальную точку, как в финале (n ?) траектория может оказаться в окрестности другого цикла.
Итак, мы пришли к первому признаку хаоса: изменение начальных условий влечет изменение финального поведения. Но это еще не все и не главное, поскольку периодическое движение, каким бы сложным оно ни было, обладает порядком. Для хаоса нужна неустойчивость в финальном поведении.
Пример 12. Рассмотрим отображение
имеющее разрыв в точке x = 0 и представляющее собой растяжение плюс сдвиг по обе стороны от этой точки (см. рис. 7, б ). Это отображение имеет две неустойчивые неподвижные точки (считаем 1 < q < 2)
и интервал (- d, d ), на который попадают все траектории из области | x | < x* и остаются на нем навсегда. Заметим, что условие q > 1 означает растяжение, то есть неустойчивость отображения в обеих областях непрерывности (x < 0) и (x > 0), а условие q < 2 принято для того, чтобы выполнять неравенство d < x* и отрезок [- d, d] был притягивающим, содержащим аттрактор.
Итак, мы получили, что отображение имеет притягивающий отрезок, не имеющий внутри себя ни одной устойчивой траектории: устойчивая неустойчивость! Таким образом, отображение (21) при 1 < q < 2 демонстрирует настоящий хаос. Я.Г. Синаю, выдающемуся российскому математику, первому удалось построить вероятностную теорию для динамических систем, простейшим представителем является (21). Вероятность в такие системы вносится через случайное распределение начальных условий. Вопрос этот непрост и требует отдельного обсуждения, тем более на элементарном уровне.
Мы закончим наше обсуждение ответом на вопрос, откуда берутся отображения как динамические системы. Во-первых, как уже упоминалось, при замене производной в (1) разностью
уравнение (1) сводится к отображению (16) с функцией
f (x(n)) = x(n) + hF (x(n)).
Так поступают, когда приближенно решают дифференциальные уравнения с помощью компьютера.
Второй путь получения отображений - отображения последования Пуанкаре, смысл которого в том, что для системы (6) выбирается секущая, например y = 0, находится время возвращения траектории, начинающейся на этой секущей, на нее же из уравнения 0 = S2(t, x, 0, l), и это время t(x) подставляется в первую функцию в (7): . Полученное выражение и есть одномерное отображение вида (16).
Пример 13. Пусть в примере 8 динамическая система (6) в квадрате Q задана линейной системой (9) с решением (10) в случае седла, то есть при l < 0, g > 0. Построим отображение последования Пуанкаре отрезка p = {x = d, 0 # y # d} на полуось p ' = = {x = d, y # d}, считая, что вне квадрата Q траектории системы индуцируют отображение последования L отрезка q = {y = d, 0 # x # d} на полуось p ' вида
Получим отображение T : p q p'. Для этого запишем начальные и конечные условия для решения (10) в виде
x(0) = d, y(0) = y0 , d = y0e gt, x(t) = x0e lt,
где t - время движения по траектории от отрезка p до отрезка q. Решая совместно уравнения (22) и (23), получаем искомое отображение последования
где обозначено
Задача. Провести исследование отображения (24) при различных значениях параметров m, a, n и тем самым проверить утверждение примера 8.
Рассмотренные нами вопросы более подробно освещены в научной литературе.
РЕКОМЕНДУЕМАЯ ЛИТЕРАТУРА
1. Андронов А.А., Витт А.А., Хайкин С.Э. Теория колебаний. 2-е изд. М.: Физматгиз, 1959. 926 с.
2. Андронов А.А., Леонтович Е.А., Гордон И.И., Майер А.Г. Теория бифуркаций динамических систем на плоскости. М., 1967.
3. Динамические системы-5 / Под ред. В.И. Арнольда. Итоги науки и техники. Современные проблемы математики. М.: ВИНИТИ, 1986.
* * *
Владимир Николаевич Белых, доктор физико-математических наук, профессор, зав. кафедрой высшей математики Волжской государственной академии водного транспорта (Нижний Новгород), заслуженный деятель науки Российской Федерации. Специалист по теории нелинейных колебаний и качественной теории динамических систем. Автор более 150 статей и одной монографии по теории динамических систем.
