Тріплі Віяла: Хаос, Який Можна Спостерігати

З одного суглоба до трьох

Одинокий підвішений маятник є однією з найбільш передбачуваних систем у фізиці: для малих коливань він точно періодичний, і навіть при великих амплітудах його рух залишається обмеженим простим, добре зрозумілим орбітою в фазовому просторі. Додайте другий жорсткий зв’язок і система, відома як подвійний маятник, вже стає хаотичною, коли вона достатньо енергійно відхиляється. Додайте третій зв’язок — три жорсткі стійки, з'єднані безперешкодними осями, кожна з яких несе власну масу, усі вільно коливаються під дією гравітації в вертикальній площині — і хаос не тільки продовжується, він драматично посилюється.

Причина полягає у розмірності. Фазовий простір подвійного маятника описується двома кутами та їх двома кутовими швидкостями, тобто чотиривимірним фазовим простором. Триразовому маятнику потрібно три кути та три кутові швидкості — шість вимірів — і кожен додатковий напрямок надає енергії ще один канал для перемішування між зв’язками. Більше взаємозв’язків, більше способів того, щоб невелике збурення росло, і система, яка швидко плутає спостерігача щодо її майбутнього, ніж у свого дворазового родича.

Визначення руху: щільна матриця з кутовою залежністю

Не існує аналітичного розв’язку для траєкторії потрібного pendula. Він побудований на основі лагранжевої механіки: запишіть кінетичну енергію кожної маси як функцію від усіх трьох кутів та їхніх швидкостей (положення кожної маси залежить від всієї ланцюга над нею, тому швидкості з’являються помноженими на cos і sin терміни різниць кутів), запишіть потенційну енергію від висоти кожної маси та застосуйте рівняння Ейлера-Лагранжа. Результат – три закопчені нелінійні другорядні диференціальні рівняння, які найкраще записано у матричній формі.

M(θ) · θ'' = F(θ, θ') θ'' = M(θ)⁻¹ · F(θ, θ') // M є щільною 3x3 матрицею, яка змінюється з поточними кутами — // прискорення кожного суглоба залежить від стану всіх трьох ланок На кожному кроці симуляції переоцінюється ця матриця мас на поточному конфігурації та інвертується для отримання трьох кутових прискорень. У цьому немає нічого екзотичного в алгебрі, але щільність зв’язку — кожен елемент обернення залежить від кожного кута — є тим, що робить поведінку системи такою непередбачуваною інтуїцією.

M(θ) · θ'' = F(θ, θ')

θ'' = M(θ)⁻¹ · F(θ, θ')

// M is a dense 3×3 matrix that changes with the current angles —
// every joint's acceleration depends on the state of all three links

Чому три удари по дві: більший ляпрок

Хаос — це не просто "складне рух" — у нього є число, пов’язане з ним. Найбільший ляпрок (λ) вимірює, як швидко дві початково близькі траєкторії відокремлюються: мала початкова різниця δ(0) приблизно зростає як δ(0)·e^(λt). Для подвійної маятникової системи з енергетично-вибуховим рухом λ зазвичай становить кілька секунд, що дає межу передбачуваності в кілька секунд. Збільшена взаємодія у потрійному маятнику приблизно подвоює найбільший ляпрок для порівняної кількості енергії, тому така ж початкова невизначеність розростається до абсолютно нерозпізнаваного відмінності у вільному русі за дуже короткий час.

t_horizon ≈ (1 / λ) · ln(Δ_max / Δ0)

// Δ0 = initial uncertainty, Δ_max = the largest error you can tolerate
// doubling λ roughly halves how long your prediction stays useful

Чому інтегратор повинен бути надійною в цьому випадку

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

Чим демонструє хаотичність триплексного маятника п'ять майже ідентичних запусків

Найбільш очевидний спосіб побачити чутливу залежність – не формула, а зображення: запустіть п’ять копій триплексного маятника, відрізняючись у своєму першому куті на величину до 0,0001 радіана, і відстежуйте кожен з них окремим кольором. На мить усі п'ять траєкторій лежать один над одним, оскільки розбіжність починається майже з нуля. Точно тоді, коли експоненційний ріст передбачає закон, вони помітно відокремлюються, і протягом кількох штрихів маятника п’ять маятників повністю відрізняються конфігурації – незважаючи на те, що кожен з них слідує за точно однаковими детермінованими рівняннями із майже ідентичного старту. Ця проміжна ділянка між повним детермінізмом та практичною невловимістю є повною суттю теорії хаосу, а триплексний маятник робить це достатньо наглядно для спостереження в реальному часі.

Frequently asked questions

Спробуйте наживо

Усе, що вище, працює прямо у вашому браузері — відкрийте Triple Pendulum і змінюйте параметри під час роботи. Нічого не встановлюється, нічого не завантажується на сервер, уся модель живе в одній вкладці.