Как отличить хаос от миража: скрытые аттракторы и конечномерные ляпуновские размерности ✎

Оригинал: "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 · собственный разбор · 3 мин
Численный анализ динамических систем сталкивается с проблемой различения настоящего хаоса и переходных хаотических множеств, которые могут существовать сколь угодно долго.

Контекст

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

Методы

Для анализа использовались классические системы Лоренца, Рёсслера и Валлиса. В системе Лоренца при параметрах r=24, σ=10, b=8/3 визуализировалось скрытое переходное хаотическое множество, траектории которого могут вести себя хаотически в течение длительного времени, прежде чем сойтись к устойчивому равновесию. Для системы Рёсслера применялась стабилизация неустойчивой периодической орбиты с помощью управления Пирагаса, что позволило сравнить поведение стабилизированной и нестабилизированной траекторий при вычислении конечновременных показателей Ляпунова. Для системы Валлиса использовался аналитический подход, основанный на методе функций Ляпунова и подходе Дуади–Остерле, что дало оценку ляпуновской размерности глобального аттрактора без численного моделирования.

Результаты

В системе Лоренца при r=24 обнаружено скрытое переходное хаотическое множество: время хаотического поведения зависит от начальных данных и может превышать $10^7$ единиц времени для точки (0,-0.5,0.5), тогда как для точки (20,20,20) хаос длится лишь около 1.8·$10^4$. Это показывает, что численное обнаружение аттрактора может быть иллюзией. В системе Рёсслера показано, что при вычислении наибольшего конечновременного показателя Ляпунова вдоль неустойчивой периодической орбиты результаты стабилизированной и нестабилизированной траекторий совпадают до 5-го знака лишь на интервале [0,30.4], после чего расходятся. Для системы Валлиса получена аналитическая оценка ляпуновской размерности глобального аттрактора: dim_L $A_V$ ≤ 3 - 2(C+2)/(C+1+√((C-1)^2+25B/4)), которая согласуется с локальными размерностями в равновесиях.

Значение

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

Развитие

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

Влияние

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

Следующие шаги

Следующим шагом может стать систематическое исследование скрытых переходных хаотических множеств в других известных системах, таких как система Рабиновича или Глуховского–Должанского. Также важно разработать методы оценки времени жизни переходного хаоса на основе аналитических свойств системы.

Ключевые проблемы

Статья напрямую связана с открытой проблемой существования скрытых аттракторов в системе Лоренца, которая является обобщением второй части 16-й проблемы Гильберта на хаотические системы. Эта проблема остаётся нерешённой, и данная работа вносит вклад в понимание того, почему её так трудно решить численно.

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

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

Ключевые числа

  • Время переходного хаоса в системе Лоренца (r=24): до 107 единиц времени
  • Интервал совпадения показателей Ляпунова с точностью 10^-5: [0,30.4]
  • Ляпуновская размерность аттрактора Валлиса при B=102, C=3, p=0.83: ≤ 2.65903
  • Период неустойчивой периодической орбиты в системе Рёсслера: ≈ 5.8811
Понятия
показатель ляпунова хаотический аттрактор число Рейнольдса гравитационная орбитальная динамика периодическая орбита число вращения уравнение Лотки-Вольтерры нелинейная динамика динамические системы стохастическая динамическая система