Хаос, который исчезает: как отличить бурю от миража ✎

Оригинал: "Numerical analysis of dynamical systems: unstable periodic orbits, hidden transient chaotic sets, hidden attractors, and finite-time Lyapunov dimension"
· N. V. Kuznetsov, T. N. Mokaev
arXiv:1812.02201 · 2018-12-05 · arXiv non-exclusive · собственный разбор · 2 мин
Численные расчёты могут годами показывать хаос, которого на самом деле нет — траектория лишь притворяется, что никогда не успокоится.

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

Классическая система Лоренца — три уравнения, описывающие конвекцию в нагретой жидкости. При определённых параметрах она ведёт себя как маятник, который никогда не останавливается. Но при r=24 компьютер показывает иное: хаос длится, длится, а потом вдруг исчезает. Для одной начальной точки возня продолжается больше $10^7$ единиц времени, для другой — всего 18 тысяч. Обе траектории в конце концов замирают в одной и той же точке, но по пути ведут себя так, будто хаос вечен.

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

Чтобы не обмануться, учёные придумали ловушку: зажать неустойчивую периодическую орбиту, как жука пинцетом, и следить, совпадают ли показатели хаоса у стабилизированной и свободной траекторий. В системе Рёсслера они совпадают до пятого знака лишь на интервале [0, 30.4], а потом расходятся. Значит, численный хаос был иллюзией, порождённой слишком коротким временем наблюдения.

Есть и другой путь — обойтись без компьютера вовсе. Для системы Валлиса, описывающей Эль-Ниньо, удалось вывести аналитическую оценку размерности аттрактора, не запуская ни одного численного эксперимента. Формула с квадратными корнями и параметрами B и C даёт верхнюю границу: при B=102, C=3, p=0.83 размерность не превышает 2.65903. Это как измерить глубину озера, не ныряя в него, а лишь глядя на карту.

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

В системе Лоренца при r=24 хаотическое поведение может длиться более 10 миллионов единиц времени, прежде чем траектория сойдётся к устойчивому равновесию. Это сравнимо с тем, как если бы вы наблюдали бурю в стакане воды в течение года, а затем она внезапно успокоилась.

\dim_L A_V \le 3 - \frac{2(C+2)}{C+1+\sqrt{(C-1)^2+\frac{25}{4}B}}
Верхняя граница фрактальной размерности аттрактора, выраженная через параметры системы
Понятия
показатель ляпунова хаотический аттрактор число Рейнольдса гравитационная орбитальная динамика периодическая орбита число вращения уравнение Лотки-Вольтерры нелинейная динамика динамические системы стохастическая динамическая система