Нелинейная динамика: аттракторы, бифуркации и хаос

Детерминированная система с простыми нелинейными правилами может вести себя непредсказуемо — это открытие середины XX века полностью изменило наше понимание природы. Хаос — не беспорядок, а детерминированная непредсказуемость, вызванная чувствительностью к начальным условиям.

Фазовое пространство и аттракторы

Фазовое пространство: пространство всех возможных состояний системы. Для системы с переменными x₁,...,xₙ — это ℝⁿ. Траектория системы — кривая в фазовом пространстве.

Аттрактор: множество состояний, к которому стремятся все близкие траектории. Типы:

Точечный аттрактор: система стремится к стационарному состоянию. Математический маятник с трением: стремится к положению покоя x=0.

Предельный цикл: устойчивые периодические колебания. Биологические осцилляторы (сердечный ритм, циркадные ритмы). Уравнения Лотки-Вольтерра (хищник-жертва) при определённых параметрах.

Странный аттрактор: непериодическое, ограниченное, чувствительное к начальным условиям движение. Аттрактор Лоренца, аттрактор Рёсслера.

Аттрактор Лоренца и детерминированный хаос

Эдвард Лоренц (1963) изучал упрощённую модель атмосферной конвекции:

dx/dt = σ(y − x), dy/dt = x(ρ − z) − y, dz/dt = xy − βz

Расшифровка переменных и параметров:

• x — интенсивность конвективного потока
• y — разность температур между восходящим и нисходящим потоком
• z — отклонение вертикального температурного профиля от линейного
• σ (Прандтль) = 10 — отношение вязкости к теплопроводности
• ρ (Рэлей) = 28 — тепловой градиент
• β = 8/3 — геометрический параметр

При σ = 10, ρ = 28, β = 8/3: система ведёт себя хаотически — траектория «мечется» между двумя «крыльями» бабочки аттрактора, никогда не повторяясь.

Эффект бабочки: начальные условия x₀ = (1, 1, 1) и x₀' = (1.0001, 1, 1) (отличие 0.01%) → экспоненциальное расхождение через ~30 единиц времени. Показатель Ляпунова λ₁ ≈ 0.9: начальная ошибка ε возрастает как εe^{λt}.

Бифуркации: качественные изменения поведения

Бифуркация — качественное изменение поведения системы при плавном изменении параметра. Типы:

Точка поворота (saddle-node): при изменении параметра два состояния равновесия сближаются и исчезают. Пример: популяция при снижении рождаемости ниже порога вымирания.

Вилка (pitchfork): одно стационарное состояние разветвляется на три. При нагреве магнит: ниже точки Кюри — два устойчивых состояния намагниченности (↑ или ↓). Выше — одно неустойчивое.

Пери Хопфа (Hopf): при изменении параметра стационарное состояние теряет устойчивость и рождается предельный цикл. Сердечная аритмия как переход через бифуркацию Хопфа.

Сценарий Фейгенбаума (период удвоения): логистическое отображение xₙ₊₁ = r·xₙ(1−xₙ). При r от 1 до 4: r < 3 — стабильное состояние, r ∈ [3, 3.45] — период 2, r ∈ [3.45, 3.54] — период 4, r > 3.57 — хаос. Константа Фейгенбаума δ = 4.669... — универсальная константа, одинаковая для всех систем этого класса.

Размерность Хаусдорфа и фракталы

Аттракторы хаотических систем имеют нецелую «фрактальную» размерность. Аттрактор Лоренца: размерность ≈ 2.06 (не 2D поверхность и не 3D объём).

Фракталы — объекты с самоподобием на разных масштабах. Береговая линия Норвегии: при измерении с шагом 50 км получаем одну длину, с шагом 1 км — в 2.4 раза больше. Размерность ≈ 1.52.

Применения в финансах: доходности акций — «фрактальный» временной ряд, самоподобный на разных масштабах (Мандельброт). Классический «случайный блуждание» (Броуновское движение) — частный случай с размерностью 1.5.

Численный пример: логистическое отображение

r = 3.9, x₀ = 0.2:

• x₁ = 3.9 × 0.2 × 0.8 = 0.624
• x₂ = 3.9 × 0.624 × 0.376 = 0.916
• x₃ = 3.9 × 0.916 × 0.084 = 0.300

При x₀' = 0.2001: расхождение с x₀ = 0.2 удваивается каждые ≈ 2.1 шага. Через 30 шагов начальная ошибка 0.0001 превратилась в ошибку ≈ 0.5 (весь диапазон [0,1]).

Задание: Реализуйте систему Лоренца (scipy.integrate.odeint). Для σ=10, ρ=28, β=8/3: (1) постройте 3D-траекторию аттрактора; (2) вычислите расхождение двух траекторий с начальными условиями, различающимися на 10⁻⁶, и оцените показатель Ляпунова λ₁; (3) исследуйте, что происходит при ρ=20 (докритический режим) и ρ=99 (иной хаотический режим).