Девлог №79 – Хвиля 59: дивні атрактори, вовк-вівця та Курамото

Нові симуляції

🌀 Дивні атрактори — 3D хаотична топологія

Що робить атрактор "дивним"?

Дивний атрактор — це множина, до якої еволюціонує динамічна система, але де сусідні траєкторії розходяться експоненційно — ознака детермінованого хаосу. Ключові властивості: (1) обмежений рух, (2) чутливість до початкових умов (додатний показник Ляпунова) та (3) фрактальна структура (нецілочислова розмірність Гаусдорфа).

Симуляція реалізує чотири класичні системи за допомогою інтегратора Рунге-Кутти 4-го порядку (RK4). Точки зберігаються в кільцевому буфері сліду та рендеряться відрізками з альфа-змішуванням. Проста перспективна проєкція з двома кутами обертання (rotX, rotY) дозволяє повне 3D-дослідження через перетягування мишею чи дотик.

Чотири атрактори

• Лоренц (1963) — ẋ=σ(y−x), ẏ=x(ρ−z)−y, ż=xy−βz. При σ=10, ρ=28, β=8/3: метелик із показником Ляпунова λ≈0.905 та фрактальною розмірністю ≈2.06.
• Реслер (1976) — ẋ=−y−z, ẏ=x+ay, ż=b+z(x−c). Простіший за Лоренца; атрактор нагадує складену стрічку. λ≈0.071 при a=b=0.2, c=5.7.
• Томас (1999) — ẋ=sin(y)−bx, і циклічно. Один параметр b контролює дисипацію; при b≈0.19 виникають три частки з 3-кратною симетрією.
• Гальворсен — циклічний квадратичний зв'язок: ẋ=−ax−4y−4z−y². Створює елегантну потрійну спіральну структуру з 3-кратною обертальною симетрією.

Кольорові режими

Три режими розкривають різну геометрію атрактора:

• Швидкість — відтінок HSL відображається на миттєву швидкість (повільно=синій, швидко=червоний). Підсвічує зони "швидкого згину" атрактора.
• Час — безперервно циклічний відтінок. Ефективно для виявлення розривів періодичності.
• Глибина — відтінок відображається на z-координату 3D сліду. Розкриває топологічне шарування.

🐺 Хижацтво вовк-вівця — просторова агентна екологія

Від диференціальних рівнянь до окремих агентів

Класичні рівняння Лотки-Вольтерри описують динаміку хижак-жертва гладкими диференціальними рівняннями. Ця симуляція замінює їх окремими агентами на сітці 50×50, де кожна вівця та вовк приймають власні локальні рішення на кожному тику. Цикли глобального масштабу виникають із суто локальних взаємодій — жодні рівняння явно не кодують коливання.

Правила агентів

• Трава відростає з фіксованою швидкістю кожен тик, обмежена 100 одиницями енергії на клітину.
• Вівці відчувають близьких вовків (радіус 5) і тікають; інакше рухаються випадковим блуканням. Вони їдять траву на своїй поточній клітині, втрачають 1 енергію/тик і розмножуються, коли енергія перевищує поріг (імовірність 4%/тик).
• Вовки відчувають близьких овець (радіус 8) і рухаються до найближчої; інакше рухаються випадковим блуканням. Вони з'їдають одну вівцю за тик (якщо перебувають в одній клітині), втрачають 3 енергії/тик і розмножуються при 70 енергії (імовірність 2.5%/тик).
• Агенти гинуть, коли енергія ≤ 0.

Графік часового ряду популяції записує історію за 350 тиків. Типова поведінка: вівці перевищують норму, коли трави багато → популяція вовків зростає з затримкою → вівці зазнають краху → вовки голодують → вівці відновлюються. Просторова структура створює прихистки — ділянки, де вівці виживають хвилі вовків — які підтримують цикли, а не призводять до вимирання.

🔄 Синхронізація Курамото — фазовий перехід у зв'язаних осциляторах

Модель

Йосікі Курамото (1975) запропонував одну з найелегантніших моделей у нелінійній науці. Кожен осцилятор i має фазу φᵢ ∈ [0, 2π), що еволюціонує як:

dφᵢ/dt = ωᵢ + (K/N) Σⱼ sin(φⱼ − φᵢ)

де ωᵢ береться з гауссового розподілу з розкидом σ, а K — глобальна сила зв'язку. Природні частоти заморожені (закалений безлад), тому без зв'язку кожен осцилятор просто рухається з ωᵢ. Член зв'язку sin притягує кожен осцилятор до фаз його сусідів.

Параметр порядку та фазовий перехід

Ступінь синхронності вимірюється комплексним параметром порядку:

r e^(iψ) = (1/N) Σⱼ e^(iφⱼ)     →     r ∈ [0, 1]

r = 0 означає, що фази рівномірно розподілені (некогерентні); r = 1 означає, що всі фази однакові (повністю синхронізовані). Точний розв'язок Курамото для лоренцівського розподілу частот дає критичний зв'язок:

K_c = 2γ   (Lorentzian)    K_c ≈ 2σ√(2/π)  (Gaussian)

Нижче K_c: r → 0 при N → ∞ (некогерентна фаза). Вище K_c: r безперервно зростає з нуля — фазовий перехід другого роду. Симуляція використовує гауссовий розподіл та оцінює K_c ≈ 2σ як швидкий орієнтир.

Спробуйте це

• Почніть з K=0 і спостерігайте, як фази рівномірно розповсюджуються. Гістограма плоска, r ≈ 1/√N (шум скінченного розміру).
• Повільно збільшуйте K вище K_c ≈ 2σ. Спостерігайте, як r зростає, а гістограма формує пік — кластер зафіксованих осциляторів.
• Натисніть Perturb, щоб штовхнути фази, і подивіться, як швидко система повторно синхронізується.
• Збільшуйте σ (розкид частот), щоб підвищити K_c — потрібно більше зв'язку, щоб подолати безлад.

Застосування в реальному світі

Синхронізація типу Курамото зустрічається в біології та інженерії: клітини серцевого водія в синоатріальному вузлі, нейронні гамма-осциляції в корі, клітини циркадного годинника в супрахіазматичному ядрі, спалахуючі світлячки та стабільність енергосистем змінного струму. Універсальність цієї моделі — один із найбільш вражаючих фактів у науці про складні системи.

Погляд у майбутнє

Кандидати симуляцій для хвилі 60 включають: реакційно-дифузійну систему Тюрінга (відмінну від машини Тюрінга), розв'язувач в'язкопружної рідини, ансамбль подвійного маятника, що показує виникнення хаосу, та візуалізатор мережевої науки для випадкових/безмасштабних/малосвітових графів. Ми також розглядаємо додавання функції "колекцій", щоб дозволити відвідувачам додавати в закладки та ділитися впорядкованими наборами симуляцій.