Jeden parametr, trzy symetryczne równania
Attractor Thomasa, wprowadzony przez fizykochemika Renę Thomasa w 1999 roku, jest jednym z najprostszych układów chaotycznych kiedykolwiek opisanych. Jego trzy sprzężone równania różniczkowe są cykliczne dla x, y, z – obracając trzy zmienne, równania mapują się na siebie – a każda niejednolitość jest jednym sinusem:
x' = sin(y) - b*x y' = sin(z) - b*y z' = sin(x) - b*z Istnieje dokładnie jeden wolny parametr, b, który odgrywa rolę współczynnika tłumienia liniowego lub rozpuszczalnego. Ponieważ układ jest niezmieniony pod wpływem cyklicznego permutowania (x, y, z), wynikający attractor – kiedy istnieje – ma trójkrotne symetrię rotacyjną: wygląda tak samo po przekształceniu osi, co daje Attractorowi Thomasa charakterystyczną strukturę węknianą, siatkowatą zamiast dwukroplistego wygląd attractora Lorenza.
x' = sin(y) - b*x y' = sin(z) - b*y z' = sin(x) - b*z
Dissycacja i dlaczego trajektorie nie rozszerzają się
Czy przepływ takowy kontraktuje objętość przestrzeni fazowej, a więc czy może wspierać ograniczonego atraktora w ogóle, jest regulowane przez rozwazanie (divergencję) pola wektorowego. Oto div(F) = -3b, stała, niezależna od położenia. Dla b > 0 rozwazanie jest jednolitnie ujemne: każda niewielka ilość bliskich warunków początkowych zawsze zmniejsza swoją objętość podczas ewolucji, więc trajektorie nie mogą rozszerzyć się do nieskończoności. Jednak nie mogą również w ogólności skonvergować do pojedynczego punktu, chyba że rzeczywisty punkt równowagi jest tam stabilny. Ta kombinacja – kontrakcja objętości na ogół, ale bez upadku do punktu – jest dokładnie geometrycznym oznakiem atraktora nieregularnego: nieskończenie skrzyżowanego, frakcyjnego obiektu z objętością zerową, na którym trajektorie są przyciągane i spacerują w nieskończoność bez nigdy dokładnie powtarzania się.
Od punktu stałego do chaosu przy zredukowaniu b
Jedyna korekta b przewodni na system przez typowy trasz dla przejścia do chaosu. Dla dużych wartości b (silne tłumienie, około b > 1), źródło i jakiejkolwiek bliska punkt stały są stabilne i wszystko spirali ku spoczynkowi -- nudno, ale przydatny test sensowności. Przy zredukowaniu b, punkty stałe stracają stabilność poprzez bifurację Hopfa i pojawia się cykl graniczny stabilny -- pojedyncza zamknięta okrągła ścieżka. Dalsze zredukowanie b, przekraczając około b ~ 0.32-0.20 (w zależności od dokładnego miejsca obserwacji), spowoduje cascade bifurcacji podwójnej okresowości, w której liczba pętli przed powtarzaniem ścieżki zwiększa się dwukrotnie na każdym kroku, a szybkość zbiegania się do chaosu przyspiesza. Ostatecznie okres staje się efektywnie nieskończony i ruch nigdy nie dokładnie się powtarza -- pełny chaos. Znane jako "klasyczne" rejon chaotyczny znajduje się około b ~ 0.19-0.20, który jest domyślnym ustawieniem dla większości renderów (w tym tego serwisu), ponieważ tworzy on najbardziej zrozumiałą strukturę siatkową.
b ~ 1.0 -> stable fixed point, trajectory spirals to rest b ~ 0.5 -> stable limit cycle, one closed loop b ~ 0.32 -> period-doubling cascade begins b ~ 0.19 -> fully chaotic, fractal strange attractor
Renderowanie: RK4 i wolne akumulowanie
Dzięki temu, że pole wektorowe zawiera tylko sinus i wyrażenie liniowe, jest dobrej jakości do oceny, więc czwarte stopnia krok Rungego-Kutta na małym ustalonej wartości dt (0,01-0,05) śledzi trajektorię dokładnie po tylu tysiącach kroków, ile potrzeba, aby struktura atraktora wypełniła się widocznie. Nie tak jak w mapie, na przykład trąciłki Barnsleya, nie ma tu prostej formy zamkniętej – każdy punkt na atraktorze jest dozwołony tylko przez rzeczywiste integrowanie równania różniczkowego zaczynając od pewnego punktu początkowego. Dlatego te renderowania zawsze pokazują jedno ciągle rosnące ścieżki, a nie rozrzucane pośrednie kropki. Kolor jest zwykle mapowany na chwilową prędkość lub czas od początku, co konwertuje inaczej monokromatyczną siatkę na czytelny wizualnie węzeł pokazujący, które obszary trajektoria odwiedza często.
Często zadawane pytania
Dlaczego atraktor Thomasa wygląda tak inaczej niż atraktor Lorenza?
Równania systemu Thomasa są symetryczne w cykliczny sposób względem x, y i z, dlatego atraktor dziedziczy trójwymiarową simetrię obrotową i tworzy przewinątko, podobne do siatki. System Lorenza nie ma takiej symetrii między swoimi zmiennymi, co powoduje, że jego atraktor tworzy znany kształt skrzydłowej motyla o dwóch lobach.
Co fizycznie reprezentuje parametr b?
Parametr b jest współczynnikiem zanikania liniowego, który jest stosowany równomiernie do wszystkich trzech zmiennych. Ustala rozbieżność przepływu (-3b), która kontroluje, jak silnie sąsiadujące ścieżki są przyciągane do siebie; zmniejszanie b w dół prowadzi system od punktu stałego stabilnego przez cykliczną granicę ograniczoną do pełnej chaoty.
Czy atraktor Thomasa jest czuły na swój punkt początkowy?
Tak, w klasycznym chaotycznym obszarze dla b ~ 0.19 pokazuje ono cechy chaosu: dwie ścieżki zaczynające się na niezauważalnie od siebie różnych odległości oddziałują wykładniczo i wkrótce prowadzą do całkowicie innych ścieżek wewnątrz atraktora, mimo że obie pozostają ograniczone do tej samej ogólnego kształtu frakcyjnego na zawsze.
Wypróbuj na żywo
Wszystko powyżej działa bezpośrednio w Twojej przeglądarce — otwórz Thomas Attractor 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ę Thomas Attractor