PL 
EnglishDomyślny
Українська🦋 Atraktor Lorenza
Poznaj słynny dziwny atraktor i efekt motyla. Zmieniaj σ, ρ, β i patrz, jak chaotyczne trajektorie się rozbiegają.
O tej symulacji
Ta symulacja numerycznie całkuje układ Lorenza — dx/dt = σ(y − x), dy/dt = x(ρ − z) − y, dz/dt = xy − βz — metodą Rungego-Kutty czwartego rzędu (RK4) ze stałym krokiem dt ≈ 0,005. Każda klatka przesuwa stan o kilka kroków RK4 (ustawianych przez Speed) i dodaje nowy punkt (x, y, z) do przewijanego bufora do 20 000 punktów, renderowanego jako zanikający, zmieniający kolor ślad 3D za pomocą Three.js. Przy klasycznych parametrach σ=10, ρ=28, β=8/3 trajektoria nigdy się nie powtarza i nigdy się nie ustala, krążąc wokół dwóch niestabilnych punktów równowagi i rysując słynny dwupłatowy kształt motyla — podręcznikowy dziwny atraktor.
Najczęściej zadawane pytania
Co właściwie oznaczają σ, ρ i β?
Pochodzą z uproszczonego modelu konwekcji atmosferycznej Lorenza: σ (liczba Prandtla) to stosunek lepkości płynu do przewodności cieplnej, ρ (liczba Rayleigha) mierzy przyłożoną różnicę temperatur napędzającą konwekcję względem progu potrzebnego do jej rozpoczęcia, a β odnosi się do fizycznego współczynnika kształtu warstwy konwekcyjnej. „Klasyczne” wartości chaosu używane domyślnie to σ=10, ρ=28, β=8/3, zaczerpnięte bezpośrednio z oryginalnej pracy Lorenza z 1963 roku.
Dlaczego nazywa się to „dziwnym atraktorem”?
Atraktor to zbiór stanów, do których system dąży w czasie; jest „dziwny”, ponieważ ma fraktalny (niecałkowity) wymiar, a trajektoria na nim nigdy się nie powtarza ani nie przecina samej siebie, mimo że pozostaje na zawsze ograniczona do skończonego regionu. W tej symulacji widać, jak ślad okrąża oba płaty w nieskończoność, nigdy dokładnie nie powtarzając swojej ścieżki — to definiujące zachowanie dziwnego atraktora.
Jak ta symulacja demonstruje efekt motyla?
Symulacja całkuje te same trzy sprzężone równania różniczkowe w każdej klatce metodą RK4 z ustalonego punktu startowego (0,1, 0, 0). Ponieważ układ Lorenza ma dodatni wykładnik Lapunowa (λ₁ ≈ 0,91 przy klasycznych parametrach), dwie trajektorie zaczynające się w nieskończenie małej odległości od siebie rozbiegają się eksponencjalnie, mniej więcej jak e^(0,91t) — podwajając swoją separację co ≈0,76 jednostki czasu — mimo że obie pozostają ograniczone do tego samego atraktora w kształcie motyla.
Dlaczego całkować metodą RK4 zamiast prostych kroków Eulera?
Układ Lorenza jest chaotyczny i bardzo wrażliwy na błąd numeryczny, więc metoda niskiego rzędu jak Eulera szybko akumulowałaby błąd i zniekształcała długoterminowy kształt trajektorii. Metoda Rungego-Kutty czwartego rzędu znosi wyrazy błędu do O(dt⁴) w każdym kroku, więc przy bardzo małym kroku dt ≈ 0,005 symulowany atraktor ściśle odpowiada prawdziwemu rozwiązaniu ciągłemu, pozostając wystarczająco szybka, by działać w czasie rzeczywistym w przeglądarce.
Co się stanie, jeśli obniżę ρ poniżej około 24,74?
Przy ρ ≈ 24,74 (dla σ=10, β=8/3) układ przechodzi przez subkrytyczną bifurkację Hopfa. Poniżej tego progu dwa symetryczne punkty równowagi stają się stabilne, a zamiast chaotycznego motyla trajektoria zwija się do wewnątrz i osiada na jednym z nich — bez chaosu, bez atraktora, tylko proste zbieganie. Powyżej tego progu punkty równowagi są niestabilne, a trajektorie są odpychane na dziwny atraktor.