Научная статья на тему 'Сечения Пуанкаре при описании поведения нелинейных систем'

Сечения Пуанкаре при описании поведения нелинейных систем Текст научной статьи по специальности «Математика»

CC BY
1074
254
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
СЕЧЕНИЕ ПУАНКАРЕ / POINCARЙ SECTION / НЕЛИНЕЙНАЯ ДИНАМИКА / NONLINEAR DYNAMICS / ХАОС / CHAOS / ПРЕДЕЛЬНЫЙ ЦИКЛ / LIMIT CYCLE / СТРАННЫЙ АТТРАКТОР / STRANGE ATTRACTOR / ВЫЧИСЛИТЕЛЬНЫЕ МЕТОДЫ / COMPUTATIONAL METHODS

Аннотация научной статьи по математике, автор научной работы — Кондратьев Александр Сергеевич, Ляпцев Александр Викторович

Неавтономные системы исследуются путём рассмотрения их карт Пуанкаре. Карты Пуанкаре используются для выявления периодических и субгармонических решений и для изучения систем, решения которых проявляют хаотическое поведение.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Nonautonomous systems are studied by considering a Poincare map for the system. A Poincare map is used to detect periodic and subharmonic solutions and to study systems whose solutions exhibited chaotic behavior.

Текст научной работы на тему «Сечения Пуанкаре при описании поведения нелинейных систем»

УДК 53.072

Кондратьев Александр Сергеевич, Ляпцев Александр Викторович

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

Аннотация

Неавтономные системы исследуются путём рассмотрения их карт Пуанкаре. Карты Пуанкаре используются для выявления периодических и субгармонических решений и для изучения систем, решения которых проявляют хаотическое поведение.

Ключевые слова: сечение Пуанкаре, нелинейная динамика, хаос, предельный цикл, странный аттрактор, вычислительные методы.

Метод фазовой плоскости широко используется для анализа поведения автономных систем, позволяя с помощью фазовых диаграмм выявить многочисленные важные черты поведения нелинейных систем [1, 2]. Например, замкнутые фазовые траектории соответствуют периодическим во времени решениям. Ситуация становится гораздо более сложной в случае неавтономных систем, когда расшифровка фазовых диаграмм представляет собой очень трудную задачу. Один из возможных путей в этом случае -это использование карт (или сечений) Пуанкаре [3, 4]. Такие карты позволяют заменить изучение неавтономных систем исследованием динамических систем, определенных положениями на фазовой плоскости в равномерно распределенные моменты времени, такие как г = 2тг, п = 0, 1, 2, ... Достоинства данного метода становятся особенно впечатляющими, когда рассматриваются нелинейные системы, для которых неизвестны точные решения. В этом случае фазовые траектории рассчитываются численными методами, например, методом Рунге-Кут-та. Существуют пакеты программ, позволяющих проводить построение карт Пуанкаре для подобных систем.

Принцип построения карт Пуанкаре для неавтономных систем, описываемых дифференциальным уравнением второго порядка,

© Кондратьев А.С., Ляпцев А.В., 2012

состоит в следующем. В фазовом пространстве исследуемой системы рассматриваются поперечные сечения фазовых траекторий плоскостями t = const. Обычно в качестве временного интервала этих сечений берется период внешнего воздействия на систему. Точки пересечения фазовых траекторий с плоскостями t = const указываются на одной плоскости (x, y), которая называется сечением Пуанкаре. В результате возникает удобная для наглядного представления двумерная картина, которая содержит последовательность точек пересечения фазовых траекторий с различными плоскостями t = const. Это своего рода «дискретная орбита», показывающая (x, y) координаты начальной точки спустя один, два и так далее периодов внешнего воздействия на систему.

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

x(t) + w2(x(t) - x0) = F cos t, (1)

где F- положительная амплитуда приведенной вынуждающей силы, частота которой положена равной единице. Общее решение уравнения (1) при w Ф 1 дается выражением:

F

x(t) = x0 + A sin( wt + y) + —— cos t, (2)

w -1

где A и y- произвольные постоянные. Для скорости с помощью (2) имеем:

F

u(t) = wA cos( wt + y) —-— sin t, (3) w -1

Заметим далее, что дифференциальное уравнение второго порядка для неавтономной системы (1) можно свести к системе трех дифференциальных уравнений 1 -го порядка для автономной системы:

X (t) = i(t), i) (t) = F cos( j(t)) -w2( x(t) - x0), (4) j(t) = 1.

где j- циклическая переменная (угол). Решения системы (4) даются выражениями (2) и (3), к которым добавляется уравнение:

j(t) = t, (5)

Переменные x(t), i (t), j(t) можно рассматривать как некоторые криволинейные координаты в трехмерном фазовом пространстве, параметрически задающие некоторую кривую в этом пространстве. Особенно наглядная кривая получается при использовании цилиндрических координат r, z, j. Именно, если положить z = i(t), r = x(t), j = j(t), то траектория движения в трехмерном фазовом пространстве представляет собой линию, «намотанную» на тор. На рис. 1 приведены траектории для w = 20 (рис. 1 а) и w= 1/20 (рис.1 б).

Движение точки в трехмерном фазовом пространстве можно представить как сово-

купность движений по двум окружностям с частотами £ = 1 и о (рис. 1 а). Сечение Пуанкаре получается как сечение тора полуплоскостью (рис. 1 б).

Далее мы будем выбирать секущую плоскость соответствующую углу ф = 0, что соответствует моментам времени г = 2кп, п = 0, 1, 2, ... Мы будем рассматривать частный случай х0 = 0, тогда карта Пуанкаре определяется выражениями:

хп = х(2пп) = А $т(1п(Оп +у) + ——,

о -1

ип = и(2жп) = соА со$,(2жсоп + у), (6) п = 0, 1, 2, ...

На рис. 2 и 3 показаны первые одиннадцать (п = 0, 1, 2, ... 10) положений точек карты Пуанкаре (хп, ип) на фазовой плоскости (х, и) для различных значений о. Для простоты выбраны значения А = Е = 1, ф = 0. Такие графики называют сечениями Пуанкаре.

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

К

x

w = 2

u

V

v8 ,

x

w = 3

' 4 1л v"3,2y

В

4 _ 1

v"3, _

u

x

1

(Q = -2

щим образом: пусть при w = 2 функция x (t) является решением (1) c начальными условиями x(0) = 1/3, u(0) = 2 и пусть X(t) = x (t + 2p). Поскольку сечение Пуанкаре состоит из единственной точки (1/3, 2),

то X(0) = x(2p) = 1 и X(0) = x(2p) = 2 . Таким образом функции x (t) и X(t) имеют одинаковые начальные значения при t = 0 и то же справедливо для производных этих функций. Далее, поскольку cos t является 2p периодической функцией, то справедливо:

X (t) + (О2 X (t) = = X(t + 2p) + w2 x(t + 2p) =

= cos(t + 2p). Это означает, что функции x (t) и X(t) должны совпадать на всем интервале [0, ¥), поскольку справедлива теорема единственности решения дифференциального уравнения при заданных начальных условиях [3].

Таким образом, х (?) = X (?) = х (? + 2р) для всех ? > 0: функция х (?) является 2р-перио-дичной. С помощью таких же аргументов можно убедиться в том, что при со = 1/2 имеется решение с периодом 4р, которое осциллирует между двумя точками на фазовой плоскости, показанными на сечении Пуан-

, п=0

(0 = -Д

Рис. 3

каре. При ю = 1/3 имеется решение с периодом 6р, обращающееся между тремя точками на сечении Пуанкаре, при со = 1/4 имеется 8р- периодическое решение, обращающееся между четырьмя точками. Эти решения называются субгармоническими.

А вот случай с = 42 приводит к совершенно иному результату (рис. 3). На сечении Пуанкаре видно, что точки на фазовой плоскости не повторяются (линии, соединяющие точки проведены, чтобы был понятен порядок следования этих точек при увеличении значения п). Может быть, рассмотрено слишком мало точек? Оказывается, что благодаря тому, что 42 - число иррациональное, любое новое число п приводит к новой точке на сечении Пуанкаре. Все эти точки лежат на эллипсе. Действительно при а = 42, А = Е = 1, р = 0 имеем: хп = 8ш(2л/2рп) +1, и =42 С08(2л/2рп), п = 0, 1, 2,...

Отсюда следует, что каждая точка с координатами (хп, ип) лежит на эллипсе, определяемом уравнением:

(х -1)2 + ^ = 1.

2

Мы рассматривали 2р периодические решения уравнения (1), поскольку вынуждающая сила имеет период, равный 2р. При этом обнаружились субгармоники 2 • 2р, 3 • 2р, 4 • 2р при значениях со, равных 1/2, 1/3, 1/4 соответственно. При иррациональ-

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

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

x(t)+2bx(t)+W (x(t) - x0) = F cos t, (7) где b - положительная константа. При условии b2 < 4ю2 общее решение уравнения (7) имеет вид:

x(t) = A exp(-bt) sin(Vw2 - b2 ■ t + j) +

F

1 ■ sin(t + q), (8)

д/(ю2 -1)2 + 4b2

где tgd =

Q)2-1 2b !

а A и j-произвольные по-

стоянные. Построим карту Пуанкаре для значений г = 2рп, п = 0, 1, 2, ... Выберем для расчета значения Ь =0,11, с = А = Г = 1,

j = 0. Поскольку tg в =

ю2-1 2b

= 0, то выбе-

рем значение 6 = 0. Результаты расчета соответствующей карты Пуанкаре приведены на рис. 4.

Точками показаны соответствующие значения п = 0, 1, 2, ..., а кружком на оси и обозначено значение хп = 0, ип » 4,545, к которому стремятся хп и ьп после некоторого достаточно большого числа итераций. Можно предположить, что существует 2р- периодическое решение, соответствующее выбранным значениям А и (. Это действительно так, как следует из выражения (8) при А = 0.

Имеет место существенное различие между сечениями Пуанкаре для уравнений (1) и (7). На рис. 2 и 3 положение всех точек зависит от выбранных начальных значений (А = 1, р = 0). На рис. 4 положение первых нескольких точек зависит от начальных условий, однако положение предельной точки от них не зависит. Это характерная черта уравнений, содержащих затухание, причем

Рис. 4

для более сложных уравнений, чем (7), характерны не отдельные предельные точки, а более сложные объекты. Так, для уравнения:

X(t) + 0,22X(t) + x(t) = cos t + cos(V2t) (9)

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

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

x0 2> U0 4

x0 = -1, U0 = 1

Рис. 5

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

Проиллюстрировать хаотическое поведение системы можно на примере уравнения Дуффинга [5]:

х(?)+2Ьх(?) - х(?)+х3 (?) = ^ 8ш(т0 .(10)

о. в

D.S О 4 D.2 D -0.2 -0.4 -0.6 -0.8

F = 0,2

F = 0,28

0.2 0.4 О.Б 0.8

F = 0,29

F =0,37

Рис. 6

-is-■-1-1-

-2-1012

9 Г:9

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

1

Рис. 7

Поскольку решение уравнения (8) не может быть записано в явном виде, то построение карты Пуанкаре придется вышолнить путем численного решения для заданных начальных значений и нахождения значений для х(2тп!%) и и(2тп/у) = х(2тп/у). Действительно, поскольку внешняя сила 2тп!у-периодична, то следует искать 2т/ у- периодические решения. На рис. 6 показаны сечения Пуанкаре при Ь = 0,15 и у = 1,2 для случаев Г = 0,2, Г = 0,28, Г = 0,29, Г = 0,37.

Для наглядности точки сечения Пуанкаре нанесены на фазовые траектории, что позволяет понять последовательность следования этих точек во времени.

Видно, что при увеличении значения Г характер карт Пуанкаре изменяется. При Г = 0,2 существует 2т!у - периодическое решение. При Г = 0,28 появляется субгармоника с периодом 4т!у а при Г = 0,29 и Г = 0,37 существуют субгармоники с периодами 8т!у и 10т!у соответственно. Все драматически преображается при Г = 0,5. Теперь уже нет ни 2т! у - периодического решения, ни субгармоник. Сечение Пуанка-

Рис. 8

ре для этого случая (n = 500) показано на рис. 7 б, а на рис. 7 а приведена фазовая траектория. Как видно из рисунков, фазовая траектория в этом случае не несет большой информации, при увеличении числа точек фазовая траектория практически полностью заполняет некоторую область пространства. В то же время сечение Пуанкаре имеет характерную структуру. При этом, если исключить несколько точек с малыми значениями n, которые соответствуют начальному режиму с установлением колебаний, структура точек сечения Пуанкаре не зависит от начальных условий. Сложный вид этой структуры свидетельствует о том, что это действительно странный аттрактор [3]. Отметим, что вид сечения Пуанкаре изменяется, если построить различные сечения, например, взять t = (2кп + к 12)1 у, n = 0,1,2,...

Рассмотрим странный аттрактор для вынужденных колебаний физического маятника, описываемого уравнением:

x(t) + 2bx (t) + sin( x(t)) = F cos t, (11)

в котором x(t) соответствует углу между осью симметрии маятника и вертикалью. На рис. 8 показано сечение Пуанкаре для значений F = 2,7 и b = 0,11, полученное при использовании примерно 5000 точек. Отметим определенную цикличность картины -линии, начинающиеся при x = -к, являются продолжениями линий, заканчивающихся при x = к.

Сечения Пуанкаре для уравнения Дуф-финга и физического маятника демонстрируют не только странный аттрактор, но и явление детерминированного хаоса. Малые изменения в начальных условиях приводят

-2-1012

к огромным изменениям поведения решений уравнений. Как писал еще А. Пуанкаре, «.. .Может случиться, что небольшие различия в начальных условиях вызовут огромные изменения в протекании явления. Маленькая ошибка в прошлом приводит к чудовищной ошибке в дальнейшем. Предсказания становятся невозможными» [3].

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

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

Исследование решений нелинейных уравнений показывает, что сечение Пуанкаре, соответствующее странному аттрактору, имеет специфическую структуру, описываемую понятием фрактал и имеющую дробную размерность, лежащую в интервале между значением 1 (длина линии) и 2 (площадь поверхности) [6]. Именно сечение Пуанкаре позволяет наглядно показать, что решение соответствует странному аттрактору. Для сравнения приведем фазовую диаграмму и сечение Пуанкаре для слегка измененного уравнения (9):

х(?) + 0,22 X (?) + х(?) = 0.2 008 ? + 008(^2?)

после некоторого периода времени, соответствующего режиму установления (рис. 9 а), и сечение Пуанкаре для п = 500 (рис. 9 б).

б)

. 9

Как видно из рисунка, фазовая диаграмма подобно фазовой диаграмме рис.7 а практически полностью заполняет некоторую область фазового пространства. Однако сечение Пуанкаре показывает, что в данном случае аттрактор есть линия, а именно эллипс.

В заключение отметим, что предельные циклы в виде странных аттракторов характерны не только для неавтономных систем. Разнообразные примеры систем различной природы, предельным циклом для которых является странный аттрактор, можно найти в [6]. Одна из простейших механических систем, описывающая автоколебания вследствие сил сухого трения, предельным циклом которой является странный аттрактор, рассмотрена в статье [7]. Схема автоколебательного устройства изображена на рис. 10. На стержне А закреплен груз В, который скользит по ленте транспортера, движущейся со скоростью V (устройство, похожее на самописец). Стержень может поворачиваться (угол поворота р) и удерживается в рав-

" У

с

Рис. 10

/

г----

J*

-0.8-1-1-1-1-1-

-1.5 -1 -0.5 0 0.6 1 1.5

новесном (симметричном) положении пружиной. Стержень с пружиной закреплен на «каретке» С, которая может без трения скользить по направляющим, перпендикулярным ленте. Автоколебательная система имеет две степени свободы, однако в силу трансляционной симметрии координата каретки у не входит в систему уравнений. Вследствие этого, колебания системы описываются системой из 3 дифференциальных уравнений 1-го порядка, в которой переменными, зависящими от времени, являются р, ю = р, Уу = у. Такая система уравнений является минимальной системой, для которой возможны хаотические движения. Численный эксперимент показывает [7], что в случае, когда сила упругости пружины нелинейно зависит от угла и сила сухого трения зависит от относительной скорости тела и ленты, в системе могут реализовываться как периодические колебания, так и хаотические колебания, в зависимости от скорости движения ленты.

В данном случае фазовое пространство трехмерное, и сечение Пуанкаре является множеством точек пересечения фазовой траектории любой поверхностью, в частности, плоскостью. С учетом симметрии удобно, например, рассмотреть сечение плоскостью, определяемой равенством р= 0. Такое сече-

б) -=i-■-■-■-■-■-

0.6 ■ ,

□ .4 ■

0.2 - . f- о -0.2 ■ „ -0.4 ■ -0.6 ■

-0.8-'-'-'-'-'-

-15 -1 -0.5 0 0.5 1 15 omega

. 11

ние Пуанкаре для хаотических колебаний приведено на рис. 11 а, а для периодических колебаний на рис. 11 б. Как видно из рисунка, периодическому движению соответствует система из 6 точек, симметричная относительно поворота графика на 180°, что и должно быть вследствие симметрии физической системы. Сечение Пуанкаре для хаотических колебаний построено при расчете приблизительно 100000 точек и, тем не менее, не является строго симметричным относительно поворота графика. Это - еще одно свойство странного аттрактора - несмотря на симметрию уравнений, полностью симметричная картина сечения Пуанкаре не может быть получена в каком-либо пределе, хотя визуально такая симметрия будет наблюдаться при увеличении числа точек.

Так же, как и в случае построения сечений Пуанкаре при периодическом воздействии на систему, вид сечений Пуанкаре для автономной системы зависит от выбора той или иной переменной и значений этой переменной, для совокупности которых производится построение. В вышеприведенном примере можно, например, построить карту Пуанкаре, в которой по осям отложены переменные ри Vy, а точки карты Пуанкаре вычисляются при значениях w= 0.

Литература

1. Бордовский Г.А., Кондратьев A.C., Чоудери А.Д.Р. Физические основы математического моделирования. М.: Academia, 2005.

2. Кондратьев A.C., Ляпцев A.B. Физика, задачи на компьютере. М.: Физматлит, 2008.

3. Kent Nagle R., EdwardB. Saff. Fundamentals of Differential Equations. Addison- Wesley Publishing Company, 1993.

4. BlanchardP., Devaney R.L., HallG.R. Differential equations. Brooks/Coll, 2002.

5. Кондратьев А. С., Ляпцев A.B. Компьютерное моделирование при изучении физики. Вынужденные колебания нелинейного осциллятора // Компьютерные инструменты в образовании, 2005. № 2. С. 66-71.

6. Гритченко B. Т., Маципура B. Т., Снарский А.А. Введение в нелинейную динамику: Хаос и фракталы. М.: Изд. ЛКИ, 2007.

7. Ляпцев A.B. Странный аттрактор в простейшей механической системе // Компьютерные инструменты в образовании, 2010. № 6. С. 57-66.

Abstract

Nonautonomous systems are studied by considering a Poincare map for the system. A Poincare map is used to detect periodic and subharmonic solutions and to study systems whose solutions exhibited chaotic behavior.

Keywords: Poincare section, nonlinear dynamics, chaos, limit cycle, strange attractor, computational methods.

Кондратьев Александр Сергеевич, академик РА О, доктор физико-математических наук, профессор кафедры методики обучения физике РГПУим. А.И. Герцена, Копйга1б125@тай. ги,

Ляпцев Александр Викторович, доктор физико-математических наук, профессор, заведующий кафедрой методики обучения физике РГПУ им. А.И. Герцена. ирт_ено@таИ. ги

© Наши авторы, 2012. Our authors, 2012.

i Надоели баннеры? Вы всегда можете отключить рекламу.