🦋 Атрактор Лоренца

Система Лоренца (Едвард Лоренц, 1963) — це три пов'язан ODE, що моделюють спрощену атмосферну конвекцію. При класичних параметрах (σ=10, ρ=28, β=8/3) траєкторії ніколи не повторюються, але залишаються в обмеженій області — дивний атрактор. Зна менитий ефект метелика: дві траєкторії з різницею початкових умов 10⁻⁵ розходяться по-різному, та обидві трасують одну й ту саму форму крила метелика. Перетягуйте для обертання. 🇬🇧 English

Параметри

Траєкторія

Точок0
x—
y—
z—
ρ≈28 → класичний метелик; ρ<24.74 → стабільно.

Рівняння Лоренца

dx/dt = σ(y − x)    dy/dt = x(ρ − z) − y    dz/dt = xy − βz
Інтегрується методом РК4 (dt ≈ 0.005). Два крила атрактора відповідають двом нестійким рівновагам у (±√(β(ρ−1)), ±√(β(ρ−1)), ρ−1). Показник Ляпунова λ₁ ≈ 0.91 для класичних параметрів — траєкторії розходяться зі швидкістю e^(0.91t).

Про цю симуляцію

Ця симуляція чисельно інтегрує систему Лоренца — dx/dt = σ(y − x), dy/dt = x(ρ − z) − y, dz/dt = xy − βz — методом Рунге–Кутти четвертого порядку (RK4) з фіксованим кроком dt ≈ 0.005. На кожному кадрі стан просувається на кілька кроків RK4 (задається повзунком «Швидкість»), і нова точка (x, y, z) додається до кільцевого буфера на до 20 000 точок, що малюється як згасаючий, кольоровий 3D-слід через Three.js. За класичних параметрів σ=10, ρ=28, β=8/3 траєкторія ніколи не повторюється й не заспокоюється, обертаючись навколо двох нестійких точок рівноваги і малюючи знаменитий двокрилий силует метелика — підручниковий приклад дивного атрактора.

🔬 Що показано

Одна точка рухається у 3D фазовому просторі за рівняннями Лоренца, її шлях малюється як безперервно зростаючий слід. Два «крила» відповідають двом нетривіальним нерухомим точкам (±√(β(ρ−1)), ±√(β(ρ−1)), ρ−1); траєкторія кружляє навколо однієї з них, непередбачувано перемикається на іншу і ніколи не замикається сама на себе — візуальна ознака детермінованого хаосу.

🎮 Як користуватися

Перетягуйте полотно, щоб обертати вигляд (OrbitControls із демпфуванням); повзунки σ, ρ і β змінюють форму атрактора в реальному часі. «Швидкість» задає кількість кроків RK4 на кадр анімації. «Скинути» перезапускає з (0.1, 0, 0); «Autorotate» вмикає/вимикає повільне автоматичне обертання камери. Кнопка «Info» відкриває вікно з повними рівняннями й описом керування.

💡 А ви знали?

Едвард Лоренц знайшов цю систему випадково у 1961 році, повторно запускаючи модель погоди: надрукувавши проміжне значення як 0.506 замість повного 0.506127, він отримав під час повторного введення геть інший прогноз. Ця похибка округлення виявила чутливу залежність від початкових умов — так виник термін «ефект метелика».

Поширені запитання

Що насправді означають σ, ρ і β?

Вони походять зі спрощеної моделі атмосферної конвекції Лоренца: σ (число Прандтля) — це відношення в'язкості рідини до теплопровідності, ρ (число Релея) вимірює прикладену різницю температур, яка рухає конвекцію, відносно порогу, необхідного для її початку, а β пов'язане з фізичним співвідношенням сторін конвективного шару. Значення «класичного хаосу», що використовуються за замовчуванням, — σ=10, ρ=28, β=8/3, взяті безпосередньо з оригінальної статті Лоренца 1963 року.

Чому це називається «дивний атрактор»?

Атрактор — це множина станів, до яких система прямує з часом; він «дивний», бо має фрактальну (нецілу) розмірність, а траєкторія на ньому ніколи не повторюється й не перетинає сама себе, при цьому назавжди залишаючись в обмеженій області. У цій симуляції видно, як слід нескінченно обвиває обидва крила, ніколи точно не повторюючи свій шлях, — саме це і є визначальною поведінкою дивного атрактора.

Як ця симуляція демонструє ефект метелика?

Симуляція щокадру інтегрує ті самі три пов'язані рівняння методом RK4 з фіксованої початкової точки (0.1, 0, 0). Оскільки система Лоренца має додатний показник Ляпунова (λ₁ ≈ 0.91 за класичних параметрів), дві траєкторії, що стартують на нескінченно малій відстані одна від одної, розходяться експоненційно — приблизно як e^(0.91t), подвоюючи відстань кожні ≈0.76 одиниці часу — навіть попри те, що обидві залишаються в межах одного й того ж атрактора у формі метелика.

Чому інтегрування виконується методом RK4, а не простим методом Ейлера?

Система Лоренца хаотична й дуже чутлива до чисельної похибки, тож метод низького порядку, як-от Ейлера, швидко накопичував би похибку і спотворював довгострокову форму траєкторії. Метод Рунге–Кутти четвертого порядку скасовує члени похибки аж до O(dt⁴) на кожному кроці, тож із крихітним кроком dt ≈ 0.005 змодельований атрактор точно відповідає справжньому неперервному розв'язку, залишаючись достатньо швидким для роботи в реальному часі в браузері.

Що станеться, якщо знизити ρ нижче приблизно 24.74?

При ρ ≈ 24.74 (для σ=10, β=8/3) система проходить через субкритичну біфуркацію Хопфа. Нижче цього порогу дві симетричні точки рівноваги стають стійкими, і замість хаотичного метелика траєкторія закручується всередину й осідає в одній з них — жодного хаосу, жодного атрактора, лише проста збіжність. Вище цього порогу точки рівноваги нестійкі, і траєкторії відштовхуються на дивний атрактор.