Научный журнал
Научное обозрение. Педагогические науки
ISSN 2500-3402
ПИ №ФС77-57475

ЗНАКОМСТВО С КАЧЕСТВЕННОЙ ТЕОРИЕЙ ДИНАМИЧЕСКИХ СИСТЕМ. ЭЛЕМЕНТАРНОЕ ВВЕДЕНИЕ

Булекбаев Д.А. 1 Морозов А.В. 1
1 ФГБОУ ВПО «Военно-космическая академия имени А.Ф. Можайского»
В статье рассматриваются простейшие динамические системы с параметром и . Пространство состояний таких систем одномерно, так, что все траектории располагаются на вещественной оси. Обсуждаются типы траекторий, их свойства, типичные бифуркации, формулируются критерии устойчивости и неустойчивости аттракторов (положений равновесия и циклов), делаются необходимые акценты для изложения этой теории на занятиях, приводятся подробно разобранные примеры. Во второй части статьи авторы касаются динамических систем с хаотическим поведением траекторий, такие модели были открыты сравнительно недавно – в 70-х годах прошлого века. Обращается внимание на сценарий перехода к хаосу через последовательность бифуркаций удвоения периода. Отметим, что вопросы хаоса во многих вузах не находят должного отражения. По этой причине мы приводим краткую методику изложения первого круга вопросов теории бифуркаций, включая элементы хаоса, а также начальный понятийный аппарат. Сегодня, благодаря компьютерам, студенты и даже школьники на занятиях по математике и информатике имеют возможность познакомиться с задачами, которые были поставлены и решены сравнительно недавно – 30-40 лет назад. Статья ориентирована на преподавателей математики и информатики школ и втузов, а также студентов и школьников, посещающих математические кружки.
непрерывные и дискретные динамические системы на прямой
элементарные бифуркации
методика обучения
1. Морозов А.В. Задачи с параметрами – ключ к формированию исследовательских компетенций // Современные наукоемкие технологии. 2021. № 5. С. 191–197.
2. Лерман Л.М. Лекции по обыкновенным дифференциальным уравнениям. М.-Ижевск: НИЦ РХЦ; Институт компьютерных исследований, 2016. 300 с.
3. Анищенко В.С. Знакомство с нелинейной динамикой. М.: Изд-во ЛКИ, 2006. 224 с.
4. Морозов А.Д. Математические методы в теории колебаний: учебное пособие. М. – Ижевск: Институт компьютерных исследований; НИЦ «Регулярная и хаотическая динамика», 2017. 144 с.
5. Шустер Г. Детерминированный хаос. Введение. М.: Мир, 1988. 238 с.
6. Фейгенбаум М. Универсальность в поведении нелинейных систем // УФН. 1983. Т. 141. № 2. С. 343–374.
7. Шарковский А.Н., Коляда С.Ф., Сивак А.Г., Федоренко В.В. Динамика одномерных отображений. Киев: Наукова думка, 1989. 216 с.
8. Булекбаев Д.А., Морозов А.В. Формирование и развитие навыков вычислительного эксперимента у обучающихся на примере исследования динамической системы // Труды Военно-космической академии имени А.Ф. Можайского. 2017. № 659. С. 202–209.
9. Малинецкий Г.Г. Математические основы синергетики: Хаос, структуры, вычислительный эксперимент. Изд. 7-е. М.: Книжный дом «ЛИБРОКОМ», 2012. 312 с.
10. Кузнецов С.П. Динамический хаос и гиперболические аттракторы: от математики к физике. М.-Ижевск: Ижевский институт компьютерных исследований. 2013. 488 с.
11. Леванова Т.А., Комаров М.А., Крюков А.К., Костин В.А., Осипов Г.В. Качественные и численные методы исследования динамических систем на плоскости: учебно-методическое пособие. Нижний Новгород: Нижегородский университет, 2015. 61 с.
12. Плюснина Т.Ю., Фурсова П.В., Терлова Л.Д., Ризниченко Г.Ю. Математические модели в биологии. М.-Ижевск: НИЦ: «Регулярная и хаотическая динамика», 2014. 136 с.
13. Морозов А.В., Пирожков М.А. Отображение пекаря и его траектории // Актуальные проблемы гуманитарных и естественных наук. 2017. № 04(99). С. 7–14.
14. Морозов А.В., Пирожков М.А. О канторовом множестве, сдвиге Бернулли и подкове Смейла // Теория и практика современной науки: материалы VII Международной научной конференции (г. Москва, 3–4.10.2012). М.: Изд-во «Спецкнига», 2012

Математическая модель BUL05.wmf объекта, процесса или явления может представлять собой функцию или систему функций, уравнение или их систему. Уравнения могут быть конечными (например, алгебраическими или трансцендентными), дифференциальными, интегральными, рекуррентными и др. Процессы можно описывать дискретно – итерациями отображений, либо непрерывно – траекториями дифференциальных уравнений или интегральными кривыми. Пусть модель BUL06.wmf зависит от некоторого параметра λ: BUL07.wmf. Значение параметра λ* называется бифуркационным, если при переходе через него объект изменяет качественные или количественные характеристики [1–3]. Поясним сказанное простыми примерами.

Пример. Рассмотрим модель BUL08.wmf: BUL09.wmf, с двумя параметрами BUL10.wmf, т.е. семейство квадратных уравнений. Таким образом, устанавливается биекция между точками плоскости параметров R2 и множеством квадратных уравнений. На плоскости параметров рассмотрим три множества BUL11.wmf, BUL12.wmf, BUL13.wmf. Ясно, что для точек из множества D0 чисел, которые представляют решения квадратного уравнения на вещественной оси R, нет. Как только точка (μ, λ) оказывается на параболе D1, уравнение приобретает решение BUL14.wmf. При переходе точки (μ, λ) в область D2 решение BUL15.wmf расщепляется на два BUL16.wmf, т.е. происходит бифуркация (ветвление), при этом кривую D1 естественно называть бифуркационной кривой.

Пример. Рассмотрим семейство алгебраических уравнений BUL17.wmf:BUL18.wmf. Хорошо известно, что при λ = 0 уравнение определяет конус, при λ > 0 – однополостный гиперболоид, при λ < 0 – двуполостный гиперболоид, т.е. при переходе через λ = 0 геометрический объект изменяет качество. Стало быть, значение параметра λ = 0 – бифуркационное.

Целью статьи является построение методики изложения начала качественного исследования двух динамических моделей c непрерывным и дискретным временем

BUL19.wmf,

BUL20.wmf BUL21.wmf

Несмотря на внешнюю схожесть и простоту этих моделей, вторая – оказывается труднее в исследовании, существенно содержательней и несравнимо богаче эффектами. Заметим, что в стандартных курсах дифференциальных уравнений втузов эти вопросы часто не затрагиваются. Однако в настоящее время этот материал желательно включить в программы для первого знакомства. Используемый в статье математический аппарат несложен, и мы не исключаем возможность ознакомления с ним школьников старших классов, например в рамках факультатива, учитывая, что качественные методы исследования, вычислительный эксперимент, имитационное моделирование играют в настоящее время в науке огромную роль [3; 4].

Динамические системы с непрерывным временем

Определение 1. Динамической или автономной системой на прямой называется уравнение

BUL22.wmf. (1)

Здесь x = x(t), t – время, f – непрерывная функция x; BUL23.wmf производная по t, т.е. скорость.

Уравнение (1) задает закон движения точки (изображающей точки) по прямой 0x, называемой фазовой или прямой состояний. Точки BUL24.wmf в которых f(x) = 0 называются положениями равновесия – это простейшие траектории. Интервалы BUL25.wmf …, BUL26.wmf отвечают также траекториям, по ним точки движутся либо к положениям равновесия, либо от них, либо убегают на бесконечность. Направление движения точек указывается стрелкой. При этом работает правило: если f(x) > 0 при BUL27.wmf, то стрелка на интервале ставится вправо, если f(x) < 0 при BUL28.wmf, то влево. Таким образом, фазовая прямая разбивается на траектории: положения равновесия и ориентированные интервалы. Других траекторий в этом случае нет. Заметим, что если BUL29.wmf, то интервал BUL30.wmf будет представлять собой область притяжения точки xk, т.е. окружающие его точки будут стремиться к нему. Такие положения равновесия называются устойчивыми или притягивающими. Если BUL31.wmf, то – неустойчивыми или отталкивающими. Таким образом, неравенства BUL32.wmf и BUL33.wmf выступают критериями устойчивости (неустойчивости) положений равновесия. Задача разбиения оси 0x на траектории – задача качественного исследования уравнения (1). Обращаем внимание, что мы не прибегали к операции интегрирования, и вообще, качественное исследование подразумевает изучение свойств решений (траекторий) по структуре уравнения [2; 5]. Следующий пример прояснит сказанное.

Пример 1. Рассмотрим уравнение

BUL34.wmf. (2)

Ясно, что BUL35.wmf, BUL36.wmf, BUL37.wmf – его положения равновесия. Из графика функции BUL38.wmf (рис. 1) легко заключить, что

bulekb1.tif

Рис. 1. Фазовая прямая 0x и семь траекторий на ней

x1 = 0, x3 = 2 – неустойчивые положения равновесия (иначе они называются седловыми), x2 = 1 – устойчивое (узловое), причем интервал (0,2) является его областью притяжения. Обратим внимание на неравенства в положениях равновесия BUL39.wmf, BUL40.wmf, BUL41.wmf.

Бифуркации в системах с непрерывным временем

Рассмотрим теперь уравнение с параме- тром

BUL42.wmf (3)

и поставим вопрос: как структура фазовой прямой будет меняться при изменении λ. Понятие бифуркации, обсуждаемое далее, продемонстрируем на следующем примере.

Пример 2. Рассмотрим некоторое обобщение уравнения (2)

BUL43.wmf. (4)

Так как параметр λ определяет сдвиг графика функции y = f(x, λ) вдоль оси ординат, то уравнение f(x, λ) = 0 может иметь разное количество корней. Следовательно, количество положений равновесия уравнения (4) зависит от λ. Нетрудно проверить, что при BUL44.wmf их будет три, при BUL45.wmf – два, а при BUL46.wmf – одно.

В данном примере изменения на фазовой прямой связаны с появлением (исчезновением) положения равновесия, при этом говорят о бифуркации рождения (гибели) положения равновесия. При этом бифуркационными значениями будут BUL47.wmf. Обратим внимание, что при BUL48.wmf устойчивое и неустойчивое положения равновесия сливаются, образуя, как говорят, седло-узел (или полуустойчивое положение равновесия) (рис. 2). Аналогичная картина наблюдается и при BUL49.wmf (почему?).

bulekb2.tif

Рис. 2. Неустойчивый узел и седло-узел при BUL50.wmf

Дискретные динамические системы и их траектории

Представим теперь для определенности некоторый объект (возможно процесс), состояние которого в каждый дискретный момент времени t = k определяется величиной xk = x(k). Предположим, что состояние xk в момент времени t = k полностью определяет состояние BUL51.wmf в последующий момент времени t = k + 1 по формуле

BUL52.wmf. (5)

Здесь f – заданная функция (иначе – отображение), определенная на множестве BUL53.wmf.

Определение 2. Уравнение (5), представляющее собой рекуррентную формулу первого порядка для получения итерационной последовательности, называется дискретной динамической системой с одномерным фазовым пространством.

Заметим, что численные методы решения уравнений BUL54.wmf или BUL55.wmf, основанные на идее дискретизации, также приводят к уравнениям вида (5) [5; 6] (почему?).

Определение 3. Пусть x0 – заданное число. По формуле (5) построим последовательность

BUL56.wmf (6)

Такую последовательность называют решением, а соответствующее множество точек на оси 0x – траекторией системы (5), выпущенной из точки x0. Впрочем, в дальнейшем решение (6) мы будем интерпретировать с геометрической точки зрения, т.е. как совокупность точек на прямой, и отождествлять его, таким образом, с траекторией.

Центральным вопросом в исследовании эволюции динамических систем является вопрос о финитном поведении. Иными словами: как ведут себя траектории при BUL57.wmf и что представляют собой предельные множества. Оказывается, что в отличие от непрерывных одномерных систем, в которых эволюция конкретной начальной точки заканчивается либо в положении равновесия, либо на бесконечности, структура предельных множеств в дискретных системах вида (5) может быть существенно сложнее. Это мы увидим позднее, а сейчас обратимся к следующему примеру.

Пример 3. Рассмотрим систему (5), где BUL58.wmf. Ясно, что для любой точки x0 соответствующая траектория BUL59.wmf будет монотонно возрастающей: BUL60.wmf и стремиться при k → +∞ к бесконечности. Это наглядно демонстрирует так называемая диаграмма Кёнигса-Ламерея, принцип построения которой описан на рис. 3.

bulekb3.tif

Рис. 3. Убегающая в бесконечность траектория (положений равновесия нет)

bulekb4.tif

Рис. 4. Положения равновесия

Ясно, что если BUL61.wmf, то для любой начальной точки BUL62.wmf. Нарисуйте самостоятельно диаграмму в этом случае.

Определение 4. Предположим теперь, что BUL63.wmf является решением (5). Такое решение называется стационарным, а соответствующая траектория, т.е. точка x0 на оси 0x – положением равновесия.

Положения равновесия системы (5) находятся из уравнения BUL64.wmf, корни которого суть абсциссы точек пересечения двух графиков BUL65.wmf (рис. 4).

Замечание. Легко видеть, что если система задана монотонно убывающим отображением f, то в ней возможно лишь одно положение равновесия (рис. 4а). В случае монотонно возрастающего отображения f возможно несколько (рис. 4в).

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

Пример 4. Рассмотрим систему BUL66.wmf, где BUL67.wmf. Видно, что при λ < 1 – положений равновесия нет, при λ > 1 – два, а при λ = 1 – одно (рис. 5), т.е. параметр λ = 1 – бифуркационный.

bulekb5.tif

Рис. 5. Бифуркация положений равновесия

Определение 5. Предположим теперь, что система (5) имеет решение

BUL68.wmf,…,

причем никакие два числа в наборе BUL69.wmf не совпадают. Такое решение называется периодическим, число m – периодом, а соответствующая траектория – циклом периода m. Числа BUL70.wmf называют периодическими точками периода m отображения f.

Цикл периода m будем обозначать символом BUL71.wmf.

Договоримся также о следующем обозначении: BUL72.wmf.

Тогда нетрудно проверить, что для периодических точек периода m цикла Xm имеет место m равенств: BUL73.wmf. Заметим, что цикл периода m = 1 представляет собой положение равновесия. Для нахождения цикла BUL74.wmf достаточно найти одну периодическую точку x1 периода m отображения f:

BUL75.wmf (7)

Другие компоненты периодической траектории найдутся по формулам

BUL76.wmf, BUL77.wmf, …, BUL78.wmf.

Замечание. Пусть f – непрерывное отображение и для некоторых чисел BUL79.wmf выполняется BUL80.wmf, BUL81.wmf. Это означает, что BUL82.wmf – цикл периода 2. Ясно также, что точки BUL83.wmf и BUL84.wmf расположены симметрично относительно биссектрисы y = x. В силу непрерывности f вытекает, что существует точка BUL85.wmf, для которой BUL86.wmf, т.е. x0 – положение равновесия. Таким образом, из существования цикла периода 2 следует существование положения равновесия.

Пример 5. На множестве (2, +∞) рассмотрим BUL87.wmf, где BUL88.wmf. Видно, что BUL89.wmf – положение равновесия, а f – убывающая. Кроме того, для любого x > 2 выполняется BUL90.wmf. Последнее соотношение означает, что все траектории рассматриваемой системы – циклы периода два. Такую систему можно назвать консервативной (почему?).

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

Устойчивость положений равновесия и циклов

Пусть x0 – положение равновесия (5): BUL91.wmf.

Определение 6. Положение равновесия x0 называется устойчивым, если существует окрестность (α, β) точки x0 такая, что для любой начальной точки BUL92.wmf выполняется BUL93.wmf при k → +∞. Если, кроме того, BUL94.wmf при BUL95.wmf и BUL96.wmf при BUL97.wmf, то x0 называется устойчивым узлом. Если же BUL98.wmf и BUL99.wmf, то называется устойчивым фокусом. При этом интервал (α, β) называют областью притяжения устойчивого положения равновесия x0. Если интервал BUL100.wmf, то систему (5) называют устойчивой в целом.

Определение 7. Положение равновесия x0 называется неустойчивым, если любая траектория, выпущенная из сколь угодно малой окрестности BUL101.wmf, за конечное время (конечное число шагов) покидает BUL102.wmf.

Критерии устойчивости и неустойчивости положений равновесия

Если для отображения f системы (5) выполняется неравенство BUL103.wmf, то положение равновесия x0 – устойчиво. На рис. 6 a изображен устойчивый узел, на рис. 6 b – устойчивый фокус. Если для отображения f системы (5) выполняется неравенство BUL104.wmf, то положение равновесия x0 – неустойчиво. На рис. 7 a, b изображены соответственно неустойчивый узел и неустойчивый фокус.

На рис. 8, 9 изображены графики монотонно возрастающих отображений f, порождающих дискретные системы с тремя положениями равновесия. Обратим внимание на то, что граничными точками области притяжения устойчивого положения равновесия являются неустойчивые положения равновесия.

bulekb6.tif

Рис. 6. Устойчивые узел и фокус

bulekb7.tif

Рис. 7. Неустойчивые положения равновесия

bulekb8.tif

Рис. 8. x0 – устойчивое положение равновесия

bulekb9.tif

Рис. 9. x0 – неустойчивое положение равновесия

Рассмотрим теперь цикл BUL105.wmf периода m системы (5). Напомним, что все числа, образующие цикл, различны (в противном случае мы имели бы цикл меньшего периода).

Определение 8. Цикл X называется устойчивым (притягивающим), если существует окрестность S этого цикла такая, что для любой точки BUL106.wmf траектория BUL107.wmf расщепляется на m последовательностей, сходящихся к точкам BUL108.wmf. При этом максимальная окрестность S, точки которой обладают указанным свойством, называется областью притяжения цикла X [6; 7].

Так, например, для устойчивого цикла BUL109.wmf периода 2 с областью притяжения S в зависимости от выбора x0 может быть (рис. 10)

BUL110.wmf либо BUL111.wmf

bulekb10.tif

Рис. 10. Диаграмма Кенигса-Ламерея для устойчивого цикла

Напомним, что числа x1, x2 связаны формулами f(x1) = x2, f(x2) = x1.

Критерий устойчивости цикла периода два

Рассмотрим уравнение BUL112.wmf, корнями которого являются x1, x2, и найдем производную отображения BUL113.wmf: BUL114.wmf. Учитывая, что BUL115.wmf, получим BUL116.wmf. А это означает, что отображение BUL117.wmf – монотонно возрастающее. Далее возможна альтернатива, отмеченная на рис. 11.

bulekb11.tif

Рис. 11. Пунктирной кривой отмечен неустойчивый цикл, сплошной – устойчивый

Для формулировки критерия введем определение.

Определение 9. Мультипликатором цикла BUL118.wmf называется BUL119.wmf.

Из сказанного и определения 10 вытекает критерий: цикл X является устойчивым, если BUL120.wmf, и неустойчивым, если BUL121.wmf.

Пример 6. Рассмотрим систему BUL122.wmf, где λ – параметр. Нетрудно видеть, что правая часть BUL123.wmf является монотонно убывающей на всей вещественной оси и, кроме того, BUL124.wmf является неустойчивым стационарным решением. Стало быть, x0 = λ положение равновесия. Найдем периодическую траекторию системы. Для этого рассмотрим уравнение BUL125.wmf. Откуда легко находятся корни x0 = λ и BUL126.wmf – периодические точки периода 2 и цикл BUL127.wmf. Исследуем цикл на устойчивость. Для этого вычислим производную BUL128.wmf, её значения BUL129.wmf, BUL130.wmf и мультипликатор BUL131.wmf. Учитывая неравенство BUL132.wmf и критерий устойчивости, заключаем, что цикл BUL133.wmf является устойчивым.

Введенное понятие мультипликатора распространяется и на произвольную периодическую траекторию BUL134.wmf. А именно: мультипликатором цикла X называется число BUL135.wmf. При этом остается в силе и указанный выше критерий: если BUL136.wmf, то цикл X устойчив, если BUL137.wmf, то – неустойчив.

Итак, в дискретных системах с одномерным фазовым пространством, заданных монотонными отображениями, могут наблюдаться траектории трех типов. Это положения равновесия, периодические траектории, т.е. циклы и траектории, асимптотически стремящиеся при n → +∞ либо к положениям равновесия, либо к циклам, либо уходящие на бесконечность. По существующей терминологии устойчивые положения равновесия и циклы называют аттракторами (регулярными аттракторами) [3; 4; 6].

Дискретная система, заданная треугольным отображением

В системах, заданных немонотонными отображениями f, может наблюдаться еще один тип траекторий. Этот тип связывают с фактом открытия в некоторых простых системах сложного установившегося динамического поведения. К настоящему времени уже известно множество примеров детерминированных динамических систем, в которых траектории ведут себя хаотическим, т.е. непредсказуемым образом [4; 7; 8]. Это явление называют детерминированным хаосом [6; 7]. Значимость этого открытия, сделанного в конце XX века, невозможно переоценить. Оно позволяет по-новому взглянуть на многие явления окружающего нас мира и переосмыслить наши знания о нем.

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

Начнем с примера BUL138.wmf BUL138.wmf (рис. 12).

bulekb12.tif

Рис. 12. График треугольного отображения

Видно, что функция f (треугольное отображение) задана на замкнутом промежутке [0,1] и отображает отрезок [0,1] в себя: BUL139.wmf. Таким образом, траектории, начавшиеся в отрезке [0,1], никогда в дальнейшем не покидают его. Кроме того, легко находятся две простейшие траектории. Это положения равновесия BUL140.wmf и BUL141.wmf. Учитывая неравенства BUL142.wmf, BUL143.wmf, заключаем, что оба положения равновесия неустойчивы. Из очевидного тождества BUL144.wmf для всех BUL145.wmf следует, что траектории, выпущенные из двух сколь угодно близких начальных точек BUL146.wmf, экспоненциально разбегаются друг от друга, оставаясь в пределах интервала [0,1]. Это свойство траекторий называют чувствительностью к начальным данным (экспоненциальной неустойчивостью траекторий). На рис. 13 представлена одна траектория, выпущенная из точки x0. Для того чтобы ощутить явление неустойчивости, возьмите в качестве начальной точку, близкую к x0, и постройте соответствующую ей траекторию, используя диаграмму Кенигса-Ламерея.

bulekb13.tif

Рис. 13. Траектория, выпущенная из точки x0

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

Легко видеть, например, что имеется один цикл периода 2: BUL148.wmf. Числа BUL149.wmf (периодические точки периода 2) – корни уравнения BUL150.wmf, для нахождения которых мы применили графический способ (рис. 14). Два цикла периода 3: BUL151.wmf и BUL152.wmf. Числа BUL153.wmf BUL154.wmf (периодические точки периода 3) – корни уравнения BUL155.wmf.

bulekb14.tif

Рис. 14. График отображения f 2(x)

Пусть BUL156.wmf – десятичная конечная дробь. Ясно, что каждое BUL157.wmf в силу специфики отображения f будет содержать не более n цифр после десятичной запятой. Всего же различных чисел вида 0. a1a2...an может быть 10n. Отсюда вытекает, что для выбранного x0 найдется номер BUL158.wmf такой, что BUL159.wmf, где BUL160.wmf. А это будет означать, что траектория за конечное число шагов выходит на цикл BUL161.wmf. Пусть теперь BUL162.wmf периодическая десятичная дробь. Тогда, как хорошо известно, найдутся натуральные числа s и BUL163.wmf, такие, что BUL164.wmf. Ясно также, что все последующие числа xk траектории будут представляться дробями вида BUL165.wmf, что в свою очередь также приведет к периодической траектории. Если же в качестве x0 взять иррациональное число, десятичная запись которого представляет собой бесконечную непериодическую дробь, то эффекта повторяемости уже не будет, т.е. мы будем наблюдать непериодическую (хаотическую) траекторию.

Таким образом, рациональная начальная точка x0 приводит к периодической траектории, иррациональная – к непериодической (хаотической) траектории. Однако проводить вычисления траекторий, выпущенных из иррациональных точек, в аналитическом виде мы не можем. В любом случае при вычислениях нам приходится аппроксимировать иррациональные числа рациональными числами и наблюдать, тем самым, периодические, возможно очень больших периодов, траектории. Поясним последнее. Аппроксимация иррациональных чисел рациональными числами тем точнее, чем больше десятичных знаков в его разложении мы берем. С другой стороны – с увеличением количества знаков увеличивается степень 10n – и в общем случае растет число N, т.е. возрастает период соответствующего цикла.

За последние полвека получен ряд теорем, используя которые можно определить, при каких условиях в дискретных системах возникает сложная нерегулярная динамика (детерминированный хаос). Один из известных результатов носит название теоремы А.Н. Шарковского [6; 7]. Приведем ее в формулировке математиков Т. Ли и Дж. Йорка.

Теорема. Если в дискретной системе имеется цикл периода 3, то имеются циклы любых сколь угодно больших периодов. На рис. 15 изображен один из двух циклов периода 3.

bulekb15.tif

Рис. 15. Цикл периода 3

Каскад бифуркаций и детерминированный хаос

Далее мы кратко рассмотрим пример дискретной системы BUL166.wmf, который уже стал классическим в нелинейной динамике [6]:

BUL167.wmf.

Здесь BUL168.wmf – параметр, BUL169.wmf – начальная точка траектории BUL170.wmf.

Можно сказать, что мы будем рассматривать не одну систему, а семейство и следить за динамикой в поведении ее траекторий в зависимости от λ. Оказывается, что траектории в такой системе по мере увеличения параметра λ от 0 до 4 демонстрируют удивительно разнообразное поведение.

Замечание. Квадратичное отображение BUL171.wmf впервые встречается в работах П. Ферхюльста (1845 г.) при описании динамики популяций.

Рассмотрим для начала BUL172.wmf. Тогда BUL173.wmf и система имеет единственное устойчивое положение равновесия BUL174.wmf, притягивающее все точки промежутка (рис. 16).

bulekb16.tif

Рис. 16. Траектория стремится к x0 = 0

bulekb17.tif

Рис. 17. Траектория стремится к положению

равновесия BUL179.wmf

Пусть теперь BUL175.wmf. Тогда x0 = 0 теряет устойчивость, т.к. BUL176.wmf и происходит первая бифуркация: рождается новое положение равновесия BUL177.wmf (рис. 17), и оно будет устойчивым до λ = 3, так как при этом BUL178.wmf.

Пусть BUL180.wmf. При λ > 3 происходит следующая бифуркация: положение равновесия x1 теряет устойчивость, и из него рождается цикл периода 2: BUL181.wmf. Здесь BUL182.wmf корни уравнения BUL183.wmf. Числа BUL184.wmf – называются, как уже отмечалось, периодическими точками периода 2 отображения f. Цикл X2 устойчив, так как для его мультипликатора выполняется неравенство BUL185.wmf (проверьте это). Поясним геометрически бифуркацию рождения цикла из устойчивого положения равновесия BUL186.wmf. Для этого изобразим графики отображений f(x) (сплошная линия) и f 2(x) (пунктирная линия) для следующих значений параметра: λ = 2.8 (рис. 18а), λ = 3 (рис. 18b) и λ = 3.2 (рис. 18c). Видно, что при BUL187.wmf BUL188.wmf, при λ = 3 BUL189.wmf, при BUL190.wmf BUL191.wmf. Таким образом, при λ > 3 уравнение BUL192.wmf помимо корня BUL193.wmf, приобретает два новых корня, которые и образуют цикл BUL194.wmf.

bulekb18.tif

Рис. 18. Сценарий бифуркации рождения цикла периода 2

Пусть BUL195.wmf. При BUL196.wmf происходит следующая бифуркация: цикл X2 теряет устойчивость, и из него рождается цикл периода 4: BUL197.wmf. Здесь BUL198.wmf – корни уравнения BUL199.wmf (периодические точки периода 4 отображения f). Цикл устойчив, так как для его мультипликатора BUL200.wmf.

При λ > 3.54… происходит следующая бифуркация удвоения периода с рождением устойчивого цикла периода 8.

Описанная ситуация происходит до значения BUL201.wmf…, при котором притягивающий цикл достигает бесконечно большого периода, а все циклы периодов BUL202.wmf оказываются неустойчивыми. Начиная со значения BUL203.wmf, в системе устанавливается непериодический (случайный) процесс, именуемый детерминированным хаосом. Структурировать предельное множество в этом случае не удается. В таблице приведены бифуркационные значения, при которых рождаются устойчивые циклы периода BUL204.wmf.

Период цикла

Значение бифуркационного

параметра λm

1

 
 

λ1 = 3

2

 
 

BUL205.wmf

4

 
 

λ3 = 3.544

8

 
 

λ4 = 3.56872

16

 
 

λ5 = 3.56981

32

 
 

λ6 = 3.569934

 

λ∞ = 3.569978

 

Таблица была открыта американским математиком Митчелом Фейгенбаумом – одним из основателей теории детерминированного хаоса. Для полноты приведем ещё один принадлежащий ему результат: для числовой последовательности BUL206.wmf существует предел BUL207.wmf Это число называется сегодня постоянной Фейгенбаума. Причем оказывается, что величина δ для любой дискретной системы, заданной отображением f с единственным максимумом, неизменна, т.е. является инвариантом.

Удивительным в описанной последовательности бифуркаций является не только наличие области в пространстве параметров, с характерным сложным поведением, но и переход от хаоса вновь к порядку, который наступает при значении λ = 3.83. При этом возникают два цикла периода 3: устойчивый и неустойчивый (рис. 19). Видно, что график отображения f3(x) имеет с биссектрисой x семь точек пересечения, абсциссы которых соответствуют устойчивому и неустойчивому циклам и положению равновесия.

bulekb19.tif

Рис. 19. Абсциссы точек пересечения графика с биссектрисой определяют два цикла (λ = 3.83)

На этом мы здесь остановимся, но обратим внимание, что после значения λ = 3.83 бифуркационная история не заканчивается. О её продолжении можно прочесть в книгах [4; 5; 7].

Замечание. Рассмотренная динамическая система, определяемая квадратичным отображением, несмотря на её внешнюю простоту, оказалась чрезвычайно емким математическим объектом. Ее исследованию в свое время были посвящены десятки научных публикаций, начиная с пионерских работ [6; 7]. О многих теоретических результатах и понятиях мы здесь, естественно, умолчали (эргодичность, перемешивание, корреляционная функция) [3], так как они требуют существенного погружения вглубь теории и выхода за рамки нашей ознакомительной статьи. Однако, несмотря на ограниченность материала, дальнейшее знакомство может быть продолжено в русле программирования и вычислительной математики. В этом направлении можно поставить и провести ряд интересных экспериментальных работ, например в курсе информатики и даже в средней школе.

Критерий детерминированного хаоса

В заключение сформулируем важное определение и один из критериев детерминированного хаоса.

Определение 10. Пусть BUL208.wmf – некоторая траектория динамической системы

BUL209.wmf, BUL210.wmf, λ – задано.

Число

BUL211.wmf

называется её показателем Ляпунова [3; 5].

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

Критерий: если при заданном значении λ показатель Ляпунова σ(λ) > 0, то система обладает хаотическим поведением траекторий (все траектории экспоненциально разбегаются друг от друга). Если σ(λ) < 0, то – регулярным, т.е. каждая траектория системы при n → +∞ экспоненциально стремится к своему аттрактору (положению равновесия или циклу).

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

BUL213.wmf

А так как ln2 > 0, то система обладает хаотической динамикой. Спрашивается, есть ли в данной системе циклы? Ответ – есть, но все они неустойчивы, а стало быть, не наблюдаемы в численном эксперименте.

Заключение

В курсе дифференциальных уравнений втузов изучение простых и содержательных эталонных математических моделей должно занять приоритетное место. Такие базовые модели в каждой естественно-научной дисциплине хорошо известны. Этим будет устанавливаться междисциплинарность в обучении, с одной стороны, с другой – подводить к пониманию проблем современной науки, развивать творческую активность и в целом повышать мотивацию студентов к обучению. На этих моделях необходимо демонстрировать методы, приемы и основные результаты. Там, где полный анализ модели осуществить в силу временных рамок невозможно, необходимо теоретический анализ заменить компьютерными экспериментами и тем самым, по возможности, визуализировать результаты [9–11]. Исследование нелинейных моделей при этом должно быть поставлено в центр обучения, ибо нелинейные модели играют и будут играть в дальнейшем в науке центральную роль. Студент по окончании обучения должен четко усвоить, хотя бы на качественном уровне, весь спектр возможного поведения нелинейных систем. Естественно, что обучение надо начинать с простейших динамических систем первого порядка и выработки некоторого понятийного аппарата. Это не потребует большого времени, зато сделает хороший задел для дальнейшего изучения этих вопросов, в том числе и в специальных дисциплинах. В настоящей статье намечен методический подход к изложению основных вопросов первого знакомства с динамическими системами, с опорой на конкретные примеры, набор которых можно с лёгкостью пополнить, воспользовавшись, например, источниками [12–14].


Библиографическая ссылка

Булекбаев Д.А., Морозов А.В. ЗНАКОМСТВО С КАЧЕСТВЕННОЙ ТЕОРИЕЙ ДИНАМИЧЕСКИХ СИСТЕМ. ЭЛЕМЕНТАРНОЕ ВВЕДЕНИЕ // Научное обозрение. Педагогические науки. – 2021. – № 5. – С. 10-23;
URL: https://science-pedagogy.ru/ru/article/view?id=2393 (дата обращения: 28.03.2024).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1,674