Простий механізм з складною відповіддю
Підвісьте невеликий магнітний кульок на підвісному маятнику над плоским диском, який містить три (або більше) фіксованих магнетів, розташованих симетрично. Відпустіть кульок з положення спокою десь над диском. Гравітація тягне його назад до центру, магнети тягнуть убік, і невелика тертя виводить енергію з системи при кожному русі. Зрештою, кулька стабілізується в нерухомому положенні безпосередньо над одним магнетом. Питання, яке ставить симуляція, є найпростішим із можливих: який саме магнет?
Для точок відпуску далеко від центру або майже прямо над одним магнетом, відповідь очевидна і маятник там же рухається безпосередньо. Ближче до середини, де три магнітні сили практично збалансовані, результат стає надзвичайно чутливим: незначне зміщення точки відпуску на міліметр може перевернути остаточний магнет, над яким він приземлиться.
Рівняння та необхідність використання реального інтегратора
Рух важеля з маятником у горизонтальному напрямку підкоряється другому закону Ньютона, при якому сумуються три сили: сила відновлення, що тягне його назад до центру (від простого коливання маятника), коефіцієнт тертя та опору повітря, пропорційний швидкості, а також інверсійне квадратичне притягання до кожного магніту:
m·a = -k·x - b·v + Σ_i C · (m_i - x) / |m_i - x|^3 ↑ ↑ ↑ сила відновлення коефіцієнт тертя та опору повітря притягання до кожного магніту m_i
Це гладке, детерміноване РУН (диференціальне рівняння звичайного типу) без будь-якої явної випадковості — проте його розв'язки є хаотичними в строгому сенсі чутливості до початкових умов. Стандартним тут є інтегратор із фіксованим кроком, такий як RK4, оскільки траєкторії біля магніту різко згинаються, і неточний крок може призвести до того, що той чи інший магніт переможе, точно такого роду помилка, яку хаотична система підсилює, а не усереднює.
m·a = -k·x - b·v + Σ_i C · (m_i - x) / |m_i - x|^3
↑ ↑ ↑
restoring damping attraction toward each magnet m_i
Розфарбовування області
Замість того, щоб опускати один важель, забарвте кожну точку на пластині, де магнети, на яких падає перо від маятника, що випущено з спокою там, — магнет 1 червоним, магнет 2 зеленим, магнет 3 блакитним. Удалені від центру, ви отримуєте великі, чисті, одноколірні області: області притягання кожного магнети. Поруч із межами між областями три кольори переплітаються на кожному рівні масштабу, який ви збільшуєте, подібно до фракталу — межа не має найменшої ознаки і незалежно від того, як сильно ви збільшуєте регіон межі, ви продовжуєте знаходити більше структури, що переплітається, а не три чисті смуги.
Чому межа є фракталом, а не просто безладною
Фрактальна структура є прямим відбитком хаосу. Дві точки випуску на протилежних сторонах межі слідують траєкторіям, які спочатку здаються однаковими, оскільки вони відчувають майже однакові сили — але незначна різниця між ними експоненційно зростає з часом (позитивний Ліapunovівський показник), поки врешті-решт дві траєкторії не розходяться до різних магнітів. Оскільки це експоненційне розходження відбувається на кожному рівні, зменшення проміжку між двома початковими точками у десять разів купує лише фіксовану кількість циклів перед тим, як дві точки знову розійдуться — що є точною самоподібністю, яка визначає фрактальну межу.
Чому це має значення за межами іграшки
Магнітний підвіс є улюбленим прикладом для навчання саме тому, що він повністю детермінований, легко збирається на столі, і водночас генерує справді непередбачувану довгострокову поведінку — яскраву спростування інтуїції, що "відсутність випадковості в рівняннях" повинно означати "передбачуваний результат". Та сама структура «басейна притягування» з фрактальними межами зустрічається у набагато більш значущих багатостабільних системах: енергетичних мережах, які можуть встановити різні стабільні режими роботи після аварії, екосистемах із кількома стабільними рівняннями популяцій, і методі Ньютона для пошуку коренів на комплексному плані, фрактальними межами якого він також характеризується з тієї ж причини.
Frequently asked questions
Чи завжди маятник зрештою зупиняється над магнітом?
Якщо в моделі є будь-який опір тертя, то так — енергія розсіюється з кожним рухом, тому вантаж втрачає достатньо кінетичної енергії, щоб бути захопленим потенціальним мінімумом найближчого магніту, коли він вичерпає енергію, необхідну для того, щоб втекти.
Чому я не можу передбачити, який саме магніт обере маятник, лише на основі точки відпуску?
Це пов'язано з хаотичною природою системи: вона має чутливу залежність від початкових умов, тому дві точки відпуску, розділені лише за кілька міліметрів, можуть бути віднесені до різних басейнів після достатньої кількості рухів. Межа між цими базанами є фракталом без найменшого масштабу, тому жодна точність вимірювання не дозволяє передбачити результат поблизу цієї межі.
Чи є фрактальне межу басейну-приваблення однаковим речі з множиною Мандельброта?
Ні, хоча обидва вони є фракталами, народженими від ітерації правила. Множина Мандельброта виникає в результаті ітерації фіксованої комплексної карти та забарвлення точок залежно від того, чи втекло орбіта. Фрактал басейну-приваблення виникає з безперервної динамічної системи з кількома стабільними рівняннями рівноваги; його фрактурна структура відображає чутливість довгострокової долі до початкових умов, а не відхід від карти.
Спробуйте наживо
Усе, що вище, працює прямо у вашому браузері — відкрийте Magnetic Pendulum і змінюйте параметри під час роботи. Нічого не встановлюється, нічого не завантажується на сервер, уся модель живе в одній вкладці.
▶ Відкрити симуляцію Magnetic Pendulum