На этом графике можно видеть, как «паутина» точки x0 = 0,8 сходится к фиксированной точке, в которой пересекаются парабола — график функции f(х) — и прямая — график функции у = х. Как и следовало ожидать, этой фиксированной точкой будет точка 0,5.
Повторим описанные выше действия для другого значения параметра k. Примем его равным не 2, а 3,1. Орбита начальной точки x0 = 0,8 будет выглядеть так.
При значениях k, больших 3, происходит нечто удивительное: хотя движение по-прежнему будет оставаться правильным, орбита точки 0,8 уже не будет стремиться к какой-то одной точке. Вместо этого она будет колебаться между значениями 0,56 и 0,76. Точечный аттрактор 0,5 словно бы разделился на две точки с координатами 0,56 и 0,76. По сути, это пример орбиты с периодом, равным 2, так называемого 2-цикла, так как мы видим два точечных аттрактора. Новая паутина, которая будет порождать уже не точку, а квадрат, выглядит так.
Продолжим увеличивать значения k и рассмотрим k = 3,5. Орбита x0 = 0,8 будет выглядеть так.
Теперь орбита будет колебаться между четырьмя точками. Их координаты приблизительно равны 0,39, 0,51, 0,82 и 0,86. Это уже 4-цикл, так как одни и те же значения будут повторяться каждые четыре шага. Кажется, что с увеличением k периоды будут удваиваться: 1, 2, 4. Сначала мы наблюдали единственный точечный аттрактор, затем — два, теперь — четыре. Логично предположить, что далее их число будет равняться восьми, шестнадцати, тридцати двум и так далее. Наблюдаемая динамика уже не столь проста, однако ее по-прежнему можно назвать более или менее регулярной.
Позднее мы рассмотрим это необычное удвоение периода еще раз, а пока ограничимся тем, что изобразим новую паутину, образованную двумя основными квадратами.
И наконец, осмелимся превысить критическое значение 3,569945. Рассмотрим k = 3,9. Ситуация радикально изменится. Орбита x0 = 0,8 будет выглядеть так.
Орбита стала хаотической! В ней больше не наблюдается никаких закономерностей. Она даже не является квазипериодической, а «прыгает» с одного места на другое и кажется случайной. А что, если мы рассмотрим k = 4?
То же самое хаотическое поведение! Диаграмма-паутина будет хаотической, а представленные на ней значения будут беспорядочно колебаться между 0 и 1.
Однако орбита и диаграмма-паутина точки х0 = 0,8 — не исключение: все остальные возможные орбиты и диаграммы будут выглядеть точно так же. И вновь мы наблюдаем эффект карточной колоды.
На этом сюрпризы не заканчиваются: два различных начальных условия, близких друг к другу, определяют орбиты, которые по прошествии определенного времени будут выглядеть совершенно по-разному. Примем k = 4. Если мы хотим изучить орбиту точки а = 0,900 и по ошибке введем значение Ь = 0,901 (например, при измерении мы допустили ошибку, равную одной тысячной), то увидим, что орбиты а и b вскоре будут значительно отличаться, хотя изначально они были близки друг к другу. Орбита точки а будет образована значениями {0,900; 0,360; 0,9216; 0,2890; 0,8219; 0,5854; 0,9708…}, орбита точки b — значениями {0,901; 0,3568; 0,9180; 0,3012; 0,8419; 0,5324; 0,9958…}. Иными словами, исходная разница в одну тысячную через несколько итераций будет иметь порядок нескольких сотых. Всего за семь итераций разница увеличится в 20 раз! По прошествии определенного времени реальная и прогнозная траектории уже не будут иметь ничего общего.
И вновь мы наблюдаем эффект бабочки.
Подведем итог: изменяя значения параметра k в логистическом отображении от k = 2 до k = 4, мы показали, как система постепенно приближается к хаотическому состоянию. А где же операции растяжения и складывания, которые порождают хаос? Прямо у нас перед глазами. Логистическая функция f(х) = kx(1 — х) «растягивает» числовой интервал между 0 и 1 вследствие умножения х на k. Затем этот интервал «складывается пополам» в результате умножения kx на (1 — х) — число, меньшее единицы. Таким образом, числовой интервал растягивается и складывается, подобно подкове.
В поисках хаоса
Хотя сегодня в математике не существует четкого определения детерминированного хаоса, он рассматривается как совокупность эффекта бабочки и эффекта карточной колоды, которые мы наблюдали и в сдвиге Бернулли, и в логистическом отображении Мэя.
От какого класса динамических систем стоит ожидать хаотического поведения?
Как вы уже знаете, хаос нужно искать среди нелинейных систем — только в них действие совокупности причин может не равняться совокупному действию этих причин по отдельности и приводить к совершенно неожиданным последствиям. Также (об этом мы не упоминали) нужно искать среди неинтегрируемых систем. Система называется интегрируемой, если ее траектории или решения можно явно выразить при помощи известных функций. Интегрируемые системы (линейные и нелинейные) предсказуемы, так как известна формула, позволяющая вычислить орбиту любой точки в любой момент времени. В неинтегрируемых системах, напротив, решение нельзя представить в виде формулы, поэтому для них нельзя составить прогноз на бесконечно большой период времени. Кроме того, если мы рассмотрим такие си¬стемы с точки зрения топологии, то увидим, что траектории будут тесно сплетаться между собой.
Если мы сведем две рассмотренные выше категории воедино, то увидим, что нелинейные и неинтегрируемые системы обладают беспорядочным, непредсказуемым поведением, указывающим на присутствие хаоса. Следует заметить: даже тогда, когда хаос требует нелинейности (чтобы небольшие изменения начальных условий могли вызывать значительные изменения) и неинтегрируемости (чтобы мы не могли делать прогнозы в долгосрочном периоде), нелинейная и неинтегрируемая динамика необязательно будет хаотической. Существуют нелинейные и неинтегрируемые системы, демонстрирующие равномерное и предсказуемое поведение. Математики говорят, что эти две характеристики — нелинейность и неинтегрируемость — являются необходимыми, но не достаточными.
С другой стороны, среди нелинейных и неинтегрируемых систем выделяют два подвида: гамильтоновы системы, сохраняющие энергию, и диссипативные, которые не сохраняют энергию. Этим двум видам систем соответствуют две разновидности детерминированного хаоса, известные сегодня.
Гамильтонов хаос наблюдается в системах, сохраняющих энергию, например в системе из трех тел, изученной Пуанкаре, в звездной системе, рассмотренной Эно и Хайлсом, в моделях бильярда, описанных Адамаром и Синаем. Как мы рассказали, это хаотическое поведение возникает в силу бесконечного числа пересечений сепаратрис седловой точки, в результате которого образуется запутанная сеть траекторий. Хотя такие системы обладают очень сложной динамикой, в них отсутствуют странные аттракторы. Существует знаменитая теорема Лиувилля, согласно которой сохранение энергии препятствует возникновению аттракторов. В самом деле аттракторы — это диссипативные структуры, в которых энергия рассеивается по мере приближения системы к аттрактору.
