Trójkatny wahadło: chaos, który możesz obserwować

Od jednego związku do trzech

Wahadło pojedyncze jest jedną z najbardziej przewidywalnych systemów w fizyce: dla małych drgań jest dokładnie okresowy, a nawet przy dużym amplitudzie ruch pozostaje ograniczony do prostej, dobrze zrozumianej orbity w przestrzeni fazowej. Dodanie drugiej twardej ramki daje sytemowi, znanej jako wahadło podwójne, już chaotyczny charakter, gdy drga energicznie. Dodanie trzeciego linka — trzech twardych prętów połączonych bezfrictionowymi osiami obrotu, każdy z własną masą, wszystkie drgające w polu grawitacyjnym w pionowej płaszczyźnie — spowoduje, że chaos nie tylko kontynuuje się, lecz intensyfikuje się znacząco.

Przyczyną jest wymiarowość. Stan wahadła podwójnego opisany jest dwiema kątami i ich dwiema prędkościami katowymi, czterodimensionalną przestrzenią fazową. Wahadło trójkatne potrzebuje trzech kątów i trzech prędkości katowych — sześciu wymiarów — a każda z dodatkowych kierunków daje energii kolejny kanal do przepływu między linkami. Więcej połączeń, więcej sposobów, w jaki mała perturbacja może się rozszerzyć, i system, który szybciej zaciąga niepewność obserwatora co jego przyszłość niż jego dwulinkowy krewny.

Pochodzenie ruchu: gęsta macierz masy zależna od kątów

Nie ma rozwiązania zamkniętego postaci dla ścieżki trójkatnego wahadła. Jest ono budowane na mechanice Lagrange'a: zapisz kinetyczną energię każdego masztu jako funkcję wszystkich trzech katów skrzydłowych i ich prędkości (położenie każdego masztu zależy od całości łańcucha powyżej, więc prędkości pojawiają się mnożone razem z cosinusem i sinusem różnic kątowych), zapisz potencjalną energię na podstawie wysokości każdego masztu, a następnie zastosuj równania Eulera-Lagrange'a. Wynikiem jest trzykrotne, nieliniowe, drugiego rzędu równanie różniczkowe, najlepiej zapisać w formie macierzowej.

M(θ) · θ'' = F(θ, θ') θ'' = M(θ)⁻¹ · F(θ, θ') // M jest gęstą 3×3 macierzą zależną od bieżących katów — // każda prędkość kątowa zależy od stanu wszystkich trzech linków W każdym kroku symulacji ponownie ocenia się tę macierz masy na podstawie bieżacej konfiguracji i odwzorowuje ją, aby uzyskać trzy prędkości kątowe. Algebra nie jest zbyt skomplikowana, ale gęstość sprzężenia — każda wpis w odwrotność zależy od każdego kąta — to dokładnie to, co sprawia, że zachowanie systemu jest tak trudne do przewidzenia na podstawie intuicji.

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

Dlaczego trzy jest lepsze od dwóch: większa wartość wykładnika Lyapunowa

Zamiana skomplikowanej ruchu na chaos polega nie tylko na tym, że jest ono złożone — chaotyczne ruchy mają przypisany im numer. Największy wykładnik Lyapunowa (λ) pomaga miary szybkości rozdzielenia dwóch początkowo bliskich trajektorii: mała początkowa różnica δ(0) rosnie około jako δ(0)·e^(λt). Dla podwójnego wahadła, które ruchem mechanicznym drży energetycznie, λ jest zwykle kilka na sekundę, co daje okres przewidywalności około dwóch sekund. Dodatkowe połączenia w trójwahadle roughly podwójne największego wykładnika Lyapunowa dla porównywalnej ilości energii, więc taka sama początkowa niepewność rozszerza się do pełnego, nienagranego różnicy wzoru na ruch w znacznie krótszym czasie.

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

Dlaczego integrator musi być zaufany tutaj

Przykładowa artykuły na tej stronie, Verlet, Leapfrog & RK4, przeanalizują ogólne trade-off między integratorami numerycznymi. W przypadku chaotycznego systemu jak ten wybór ma specyficzny i nieco przerażający charakter: niestosowny, niskorzędnikowy integrator wprowadza własny błąd numeryczny na każdym kroku, a w chaotycznym systemie ten błąd jest ekspansywnie zwiększany tym samym mechanizmem, który zwiększa rzeczywiste perturbacje fizyczne. Twarda schemat nie tylko robi symulację nieco niepoprawną — po granicy przewidywalności może cicho wygenerować całkowicie inną, plynąco sztuczną chaotyczną trajektorię, która ma nic wspólnego z prawdziwą fizykę. To dlaczego ta symulacja integruje za pomocą klasycznej czwartorzędnikowej metody Runge-Kutta przy małym, ustalonej kroku czasowym, a dlaczego rysowanie całkowitego energetycznego mechanicznego po przebiegu jest tanim i skutecznym sprawdzeniem zdrowia: dla nieszczelnego, konserwatywnego systemu energia powinna być prawie stała, a jakakolwiek ciągła migracja oznacza, że numeryka, a nie fizyka, przejęła kontrolę.

Co pokazują pięć niemal identycznych wystrzałów

Najjawniejszy sposób na zrozumienie zależności czułej jest nie wzór, ale obraz: wystrzelić pięć kopii tego samego wahadła trójkatnego, różniących się tylko bardzo początkowym kątem o 0,0001 radian. Następnie śledzić każdy z nich innym kolorem. Na chwilę wszystkie pięć śladów leżą na siebie, ponieważ odchylenie zaczyna się od niemal niczego. Potem dokładnie tak, jak przewiduje prawo wzrostu wykładniczego, one widocznie oddzielają się i w ciągu kilku drgań trzy wahadły znajdą się w całkowicie innych konfiguracjach — nawet mimo to, że każda z nich śledzi dokładnie te same deterministyczne równania z niemal identycznego początku. Ta przestrzeń między pełną determinizmem a praktycznym nieprzewidywalnością jest całą treścią teorii chaosu, a wahadło trójkatne sprawia, że można to obserwować w czasie rzeczywistym.

Często zadawane pytania

Wypróbuj na żywo

Wszystko powyżej działa bezpośrednio w Twojej przeglądarce — otwórz Triple Pendulum i zmieniaj parametry podczas działania. Nic nie jest instalowane ani przesyłane na serwer, cały model działa w jednej karcie.