Карлос Мадрид - Мир математики. т.32. Бабочка и ураган. Теория хаоса и глобальное потепление

Можно задаться вопросом: где в этом примере выполняются операции растяжения и складывания, которые порождают хаос? Чтобы обнаружить их, нужно посмотреть, какие математические действия мы совершаем при выполнении сдвига Бернулли. Мы уже говорили, что сдвиг Бернулли представляет собой сдвиг запятой в записи десятичной дроби на одну позицию вправо с последующим удалением первой цифры полученного числа. Когда мы сдвигаем запятую, в действительности мы умножаем число на 10, то есть «растягиваем» его, а когда мы стираем первую цифру, то уменьшаем, или «складываем, сгибаем» число. И вновь мы видим магический рецепт хаоса.

* * *

СДВИГ БЕРНУЛЛИ

Символическая динамика имеет и другие интересные свойства.

1) Она не поддается компьютерным вычислениям. Так как компьютеры работают с ограниченным числом десятичных знаков в записи дробей, для них все числа представляют собой точные десятичные дроби. Следовательно, если мы запрограммируем сдвиг Бернулли, то увидим на экране компьютера, что аттрактором всех орбит (подобно орбитам всех точных дробей) будет точка 0. Ни малейшего намека на хаос.

2) Существуют периодические орбиты с произвольным периодом. Так как периодические дроби могут иметь произвольный период (например, состоящий из шести цифр: то будут наблюдаться орбиты с произвольными длинами периодов: 1, 2,3,4, 5. Математики Ли Тянь-Янь и Джеймс Йорк на основе теоремы Шарковского сформулировали знаменитую теорему, согласно которой если для непрерывной функции существует орбита с периодом 3, то для нее существуют орбиты с любым периодом. Точная формулировка теоремы звучит так: существование 3-цикла подразумевает существование n-цикла (для n — 1,2,3,4, 5…). Ли и Йорк удачно подытожили смысл теоремы в названии свой статьи: «Период, равный трем, означает хаос».

3) Адамар и Смэйл обнаружили, что символическая динамика — один из самых заметных признаков хаоса. И соленоид, и подкова Смэйла, и аттрактор Лоренца обладают символической динамикой. Если мы рассмотрим десятичные дроби в двоичной системе счисления, то сможем описать каждую траекторию аттрактора Лоренца последовательностью нулей и единиц.

К примеру, траектория 0,11000101… сначала совершит два витка вокруг правой части аттрактора (так как после запятой записаны две единицы), затем — три витка вокруг его левой части (так как за двумя единицами следуют три нуля подряд) и так далее. Применив эту символическую динамику, можно доказать существование хаоса в системе Лоренца: каждая траектория будет беспорядочно вращаться вокруг правой или левой части аттрактора.

* * *

Рассмотрим теперь логистическое отображение Мэя, которое задается следующим уравнением в конечных разностях:

х n+1= kх n(1 — х n).

Иными словами, для данного начального условия хна интервале между 0 и 1 орбита х рассчитывается путем последовательного вычисления значений функции f( х) = kx(1 — х), где k— параметр, больший 1, но меньший 4. Поведение логистической системы, названной так потому, что она используется для моделирования динамики численности определенных популяций, удивительным образом зависит от значения k. Если k меньше некоторого критического значения, которое, по оценкам, составляет 3,569945…, то траектории будут иметь правильную форму. При превышении этого критического значения траектории будут стремиться к хаосу. Эта дискретная динамическая система четко показывает, что простые математические действия могут обладать неожиданно сложными свойствами.

Функция f( х) является функцией второй степени:

f( х) = kx(1 — х) = kx— kx 2.

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

Изучим динамику логистического отображения для значений k, меньших критического, к примеру для k= 2. Примем в качестве начального условия x 0 = 0,8 и определим его орбиту с помощью калькулятора:

x 1= f(х 0) = 2 х 0(1 — х 0) = 2∙0,8∙(1 — 0,8) = 2∙0,8∙0,2 = 0,32

х 2= f(х 1) = 2х 1(1 — х 1) = 2∙0,32∙(1 — 0,32) = 2∙0,32∙0,68 = 0,4352

х 3= f(х 2) = 2х 2(1 — х 2) = 2∙0,4352∙(1 — 0,4352) = 2∙0,4352∙0,5648 = 0,49160192.

Теперь, когда мы знаем, как рассчитываются первые члены орбиты, вычислим

следующие члены напрямую: