Oscylator Duffinga: Chaos w Nieliniowej Sprężynie

Sprężyna przestaje być prosta

Prosty oscylator harmoniczny posiada siłę zwrotności proporcjonalną do przesunięcia, F = -kx, a jego ruch to czysta, nudna funkcja sinusowa na zawsze. Oscylator Duffinga dodaje jeden dodatkowy człon – kubiczny βx³ – aby modelować sprężynę, która staje się twardsza lub bardziej miękka w miarę jak jest rozciągana dalej, oraz liniowe tłumienie δ i okresowy napęd zewnętrzny γ·cos(ωt) w celu utrzymania energii.

ẍ + δ ẋ + α x + β x³ = γ cos(ωt)
   damping     linear     cubic       periodic drive
                spring     spring

Podwójnie głęboki potencjał ukryty

Gdy α jest ujemne a β jest dodatnie, podłoża krzywa energii potencjalnej ma dwa dołki zamiast jednego – podwójnie głęboki potencjał – co oznacza, że ​​system ma dwie stabilne pozycje równowagi zamiast jednej. Słabo napędzony oscylator może spokojnie kołować się w jednym dołku, ale gdy amplituda napędu γ rośnie, zyskuje wystarczającą energię, aby przeskoczyć przez barierę do drugiego dołka w nieprzewidywalny sposób, mimo że równanie go napędzające nie zawiera losowości.

Od czystych cykli do podwójnego okresu

Powoli zwiększaj amplitudę napędową γ z małej wartości i długoterminowe ruch oscylatora przechodzi przez dobrze udokumentowaną sekwencję: najpierw ustabilizuje się w cyklu, który powtarza się co okres napędu; następnie, przy krytycznej wartości γ, ten cykl rozszczepia się tak, że ruch powtarza się co dwa okresy napędu; dalsze zwiększenie γ powoduje ponowne rozszczepienie do cyklu czterokrotnego okresu, a następnie ośmiokrotnego i szesnastokrotnego – kaskady podwójnego okresu, w którym krytyczne wartości γ zgrupowane są geometrycznie, zbliżając się do skończonej wartości, powyżej której ruch staje się chaotyczny: chaos. Ta sama droga do chaosu, z tym samym uniwersalnym stosunkiem między kolejnymi przerwami rozszczepienia (stała Feigenbaum, ≈ 4,669), występuje w ogromnej liczbie niezależnych systemów nieliniowych, od kapania wody po modele populacji.

Czytanie chaosu za pomocą sekcji Poincarégo

Losowa trajektoria narysowana bezpośrednio w przestrzeni położenia-prędkości jest zagmatany zgrabą, z której trudno wyciągnąć jakiekolwiek wnioski. Sekcja Poincarégo to rozwiązanie tego problemu, poprzez próbkowanie stanu tylko raz na cykl jazdy – w zasadzie robi to, jakby robiono zdjęcie w tym samym momencie obrotowym – co upraszcza ciągły chaos do sekwencji punktów. Dla okresowej orbity punkty te nakładają się na siebie, tworząc pojedynczy punkt (lub kilka, dla okręgu o okresie n), a w przypadku chaotycznego ruchu rysują skomplikowaną, samopodobną strukturę zwaną atraktorem strasznego, ujawniając organizację fraktalną, która była niewidoczna w pierwotnym chaosie.

Chaos nie oznacza losowości

Każda wielkość w równaniu Duffinga jest deterministyczna – przy identycznych początkowych warunkach, czyli dokładnej samej pozycji, prędkości i czasie, przyszłość jest całkowicie wyznaczona, bez żadnych elementów losowych. To wrażliwość na pierwotne warunki sprawia, że zachowanie wydaje się chaotyczne: dwie ścieżki, które rozpoczynają się niewiarygodnie blisko siebie, w rozdziale wykładniczym oddalają się od siebie bardzo szybko, więc nawet najdrobniejsza niepewność w stanie początkowym – w tym również błąd zaokrąglania wynikający z użycia liczb zmiennoprzecinkowych w symulacji – rośnie i wkrótce staje się dominującym czynnikiem w prognozie. Dlatego chaos wydaje się niedeterministyczny, choć zasadniczo pozostaje to deterministycznym procesem.

Frequently asked questions

Wypróbuj na żywo

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