Prosty mechanizm z skomplikowaną odpowiedzią
Zawiesić mały, magnetyczny zawiesił na pendulu powyżej płaskiej płyty z trzema (lub więcej) magnesami, ułożonymi symetrycznie. Uwolnić zawiesił z położenia spoczynkowego gdzieś nad płytą. Grawitacja przyciąga go do środka, a magnesy odchylają na boki, a niewielki tarcie wypiera energię z systemu podczas każdego ruchu. W końcu zawiesił się nieruchomo, bezpośrednio nad jednym magnesem. Pytanie, które zadaje symulacja, jest najprostsze, jakie można sobie wyobrazić: jaki magnes?
Dla punktów startowych daleko od środka lub prawie idealnie nad jednym magnesem, odpowiedź jest oczywista i zawiesił się tam bezpośrednio. W pobliżu środka, gdzie trzy siły magnetyczne są prawie zrównoważone, wynik staje się niezwykle wrażliwy: delikatne przesunięcie punktu startowego o milimetr może odwrócić ostateczny stan magnesu.
Równania i dlaczego potrzebny jest rzeczywisty licznik
Ruch poziomy bobina poddany jest prawu II Newtona, z uwzględnieniem trzech sił: siły sprężystej przyciągającej ją z powrotem do środka (od ruchu wahadła), siły tłumiącej proporcjonalnej do prędkości (tarcie i opór powietrza) oraz siły przyciągającej w kierunku każdego magnesu:
m·a = -k·x - b·v + Σ_i C · (m_i - x) / |m_i - x|^3 ↑ ↑ ↑ siła sprężysta siła tłumiąca siła przyciągająca w kierunku każdego magnesu m_i Jest to gładkie, deterministyczne równanie różniczkowe oparte na metodach numerycznych, bez żadnej jawnej losowości. Standardowo stosuje się stałkownik RK4 ze względu na to, że trajektorie w pobliżu magnesu wykazują ostry zakrzywienie i niedokładny krok może spowodować odwrócenie, który magnes ostatecznie wygra, dokładnie taki błąd, który system chaotyczny wzmacnia zamiast średnio rozpraszać.
m·a = -k·x - b·v + Σ_i C · (m_i - x) / |m_i - x|^3
↑ ↑ ↑
restoring damping attraction toward each magnet m_i
Malowanie mapy basenu
Zamiast upuszczać jeden bob, koloruj każdy punkt na talerzu, przez który magnes zawiesisty uwalniony z drżenia ostatecznie osiada – magnes 1 czerwony, magnes 2 zielony, magnes 3 niebieski. Z dala od środka otrzymujemy duże, czyste, jednobarwne obszary: baseny sił przyciągania każdego magnesu. W pobliżu granic między basenami trzy kolory przeplatają się na każdym poziomie powiększenia, w sposób podobny do fraktala – granica nie posiada najmniejszej cechy i niezależnie od tego, jak bardzo powiększysz obszar graniczny, nadal znajdujesz więcej struktury z powiązanymi kolorami zamiast trzech czystych pasm.
Dlaczego granica ma strukturę fraktalną, a nie jest po prostu chaotyczna
Struktura fraktalna jest bezpośrednim odciskem chaosu. Dwa punkty uwalniania znajdujące się po przeciwnych stronach granicy podążają za trajektoriami, które początkowo wyglądają identycznie, ponieważ doświadczają prawie tych samych sił – ale niewielka różnica między nimi rośnie wykładniczo z czasem (pozytywny wyznacznik Lyapunov), aż w końcu dwie trajektorie rozchodzą się ku różnym magnesom. Ponieważ to wykładnicze oddzielenie zachodzi na każdym poziomie, zmniejszenie odstępu między dwoma punktami początkowymi o czynnik dziesiętny daje jedynie określoną dodatkową liczbę cykli przed rozdzieleniem tych dwóch punktów ponownie – co dokładnie definiuje samosejność, która charakteryzuje fraktalną granicę.
Dlaczego to ma znaczenie poza zabawką
Magnetyczna penduluma jest ulubionym przykładem dydaktycznym, ponieważ jest w pełni deterministyczna, łatwo się ją zbudować na biurku i nadal generuje naprawdę nieprzewidywalne długoterminowe zachowania — żywy odwet przeciwko intuicji, że "brak losowości w równaniach" oznacza "przewidywalny wynik". Ta sama struktura obszaru przyciągania z fraktalnymi granicami pojawia się w znacznie bardziej istotnych systemach wielostanowiskowych: sieciach energetycznych, które mogą ustalić różne stabilne stany operacyjne po awarii, ekosystemach z wieloma stabilnymi równowagami populacyjnymi oraz metodzie Newtona znajdowania pierwiastków na płaszczyźnie zespolonej, której granice obszaru przyciągania są fraktalne z dokładnie tego samego powodu.
Frequently asked questions
Czy wahadło zawsze ostatecznie zatrzymuje się nad magnesem?
Przy obecności jakiejkolwiek tarcie w modelu, tak — energia jest rozpraszana przy każdym ruchu, więc ciężar ostatecznie traci wystarczającą energię kinetyczną, aby zostać uwięzionym w potencjalnym rowku najbliższego magnesu, gdy wyczerpie się jej energia.
Dlaczego nie mogę przewidzieć, który magnes wahadło wybierze, tylko na podstawie punktu początkowego?
Ponieważ układ jest chaotyczny: wykazuje wrażliwą zależność od warunków początkowych, więc dwa punkty uwalniania oddalone o ułamkę milimetra mogą zostać sklasyfikowane w różnych basenach po wystarczającej liczbie ruchów. Granica między basenami jest fraktalną strukturą bez najmniejszej skali, więc żaden poziom precyzji pomiarowej nie pozwala przewidzieć wyniku w jej pobliżu.
Czy granica basenu fraktala jest tym samym co zbiór Mandelbrota?
Nie, choć oba są fraktalami powstałymi przez iterowanie reguły. Zbiór Mandelbrota pochodzi z iteracji stałej mapy zespolonej i kolorowania punktów w zależności od tego, czy orbitę opuszcza (ucieka). Basen atraktoru fraktalnego pochodzi z ciągłego systemu dynamicznego z kilkoma stabilnymi równowagami; jego struktura fraktalna odzwierciedla wrażliwość długoterminowego losu na warunek początkowy, a nie ucieczkę z mapy.
Wypróbuj na żywo
Wszystko powyżej działa bezpośrednio w Twojej przeglądarce — otwórz Magnetic Pendulum i zmieniaj parametry podczas działania. Nic nie jest instalowane ani przesyłane na serwer, cały model działa w jednej karcie.
▶ Otwórz symulację Magnetic Pendulum