PL 
EnglishDomyślny
Українська🌀 Atraktor Lorenza — efekt motyla w 3D
Atraktor Lorenza 3D z interaktywnymi suwakami parametrów (sigma, rho, beta). Wypuść dwie bliskie trajektorie, by zobaczyć efekt motyla: wykładniczą rozbieżność mierzoną wykładnikami Lapunowa.
Równania Lorenza
dx/dt = σ(y−x); dy/dt = x(ρ−z)−y; dz/dt = xy−βz. Klasyczne parametry: σ=10, ρ=28, β=8/3. Atraktor Lorenza ma wymiar fraktalny D≈2,06 i dodatni wykładnik Lapunowa λ₁≈0,906, potwierdzający zachowanie chaotyczne. Trajektorie nigdy się nie powtarzają, ale pozostają w ograniczonym obszarze przestrzeni fazowej.
Efekt motyla
Dwie trajektorie startujące w odległości ε=10⁻⁵ od siebie rozbiegają się wykładniczo: |δ(t)| ≈ |δ₀|·e^(λ₁t). Przy λ₁≈0,9 horyzont przewidywalności T_pred ≈ (1/λ₁)·ln(1/ε) ≈ 13 jednostek czasu. W modelowaniu atmosferycznym (gdzie Lorenz to odkrył) ogranicza to prognozowanie pogody do około 2 tygodni, niezależnie od dokładności pomiarów początkowych.
Bifurkacje i reżimy
Dla ρ < 1: stabilny początek (wszystkie trajektorie → 0). Dla 1 < ρ < 24,74: dwa stabilne punkty stałe (C⁺ i C⁻). Przy ρ=24,74: rozpoczyna się kaskada bifurkacji podwajania okresu. Dla ρ > 24,74: chaotyczny dziwny atraktor. Przy ρ≈99,96: okno okresowe wewnątrz chaosu. Przeciągnij suwak ρ, aby odkrywać te reżimy.
O tej symulacji
Ta symulacja numerycznie całkuje układ Lorenza, dx/dt = σ(y−x), dy/dt = x(ρ−z)−y, dz/dt = xy−βz, za pomocą czterorzędowej metody Rungego–Kutty (RK4) ze stałym krokiem dt = 0,005. Dwie trajektorie są przesuwane równolegle z niemal identycznych punktów startowych — (0,1, 0, 0) i (0,1+ε, ε, 0) — dzięki czemu można obserwować ich rozdzielanie się. Każda klatka animacji wykonuje 5 podkroków RK4, a ilekroć obie trajektorie zbliżą się do siebie na mniej niż 10⁻³, kod mierzy ich odległość, dodaje ją do bieżącego szacunku wykładnika Lapunowa, a następnie renormalizuje trajektorię 2 z powrotem do odległości ε od trajektorii 1. Widok 3D to prosta projekcja perspektywiczna z obrotem przez przeciąganie myszą i przybliżaniem przez przewijanie, wyśrodkowana na średniej wysokości atraktora (z−25).
🔬 Co pokazuje
Dwa punkty materialne ewoluują pod tymi samymi równaniami Lorenza z niemal identycznych warunków początkowych. Trajektoria 1 (różowa) i Trajektoria 2 (niebieska) kreślą klasyczny motylkowy dziwny atraktor, przeplatając się między dwoma płatami, nigdy się dokładnie nie powtarzając. Odczyty na żywo pokazują pozycję (x, y, z) każdego punktu, bieżącą odległość między dwiema trajektoriami oraz szacowany największy wykładnik Lapunowa λ₁, obliczony na podstawie tempa wzrostu tej odległości między renormalizacjami.
🎮 Jak korzystać
Przeciągnij σ (sigma), ρ (rho) i β (beta), aby zmienić parametry równań — domyślne klasyczne wartości to σ=10, ρ=28, β=8/3 (≈2,67). Początkowa separacja ε ustawia, jak blisko siebie są dwa punkty startowe (jako wielokrotność 10⁻⁵), a Długość śladu kontroluje, ile ostatnich punktów jest rysowanych. Kliknij Pauza/Uruchom, aby zamrozić lub wznowić całkowanie, Restart, aby zacząć obie trajektorie od nowa, i +Trajektoria, aby ponownie utworzyć trajektorię 2 obok bieżącej pozycji trajektorii 1. Przeciągnij canvas, aby obrócić widok, i przewiń, aby przybliżyć; cztery przyciski presetów przenoszą bezpośrednio do Klasycznego chaosu, Blisko bifurkacji (ρ=24), Okresowego (ρ=99,96) i Chaosu przejściowego.
💡 Czy wiesz, że...?
Edward Lorenz odkrył ten układ w 1963 roku, ponownie uruchamiając model pogody z zaokrągloną wartością początkową — drobne obcięcie dało zupełnie inną prognozę, co później opisał jako „efekt motyla”. Przy klasycznych parametrach (σ=10, ρ=28, β=8/3) atraktor ma fraktalny (niecałkowity) wymiar wynoszący w przybliżeniu 2,06, co oznacza, że jest czymś więcej niż powierzchnią 2D, ale mniej niż pełną objętością 3D.
Najczęściej zadawane pytania
Co kontrolują σ, ρ i β w równaniach Lorenza?
σ (sigma) to liczba Prandtla, kontrolująca, jak szybko zmienna x reaguje na różnice z y; ρ (rho) to (przeskalowana) liczba Rayleigha, napędzająca ogólną energię układu i, po przekroczeniu progu, jego zachowanie chaotyczne; β (beta) to czynnik geometryczny związany z fizycznym pudełkiem, które reprezentował oryginalny model konwekcji. Domyślne wartości symulacji, σ=10, ρ=28, β=8/3, to dokładnie parametry, których Edward Lorenz użył w swojej pracy z 1963 roku, dlatego są oznaczone jako preset „Klasyczny chaos”.
Jak faktycznie obliczana jest trajektoria — czy jest dokładna?
Symulacja numerycznie całkuje trzy sprzężone równania różniczkowe klasyczną czterorzędową metodą Rungego–Kutty (RK4) ze stałym krokiem czasowym dt = 0,005, przesuwając każdą trajektorię o 5 podkroków na renderowaną klatkę. RK4 jest znacznie dokładniejsza niż prosty krok Eulera przy danym kroku czasowym, ale wciąż jest przybliżeniem: jak każdy integrator numeryczny zastosowany do układu chaotycznego, drobne błędy dyskretyzacji ostatecznie rosną wykładniczo, co samo w sobie jest demonstracją czułości, którą ilustruje symulacja.
Dlaczego są dwie trajektorie i czym jest wykładnik Lapunowa?
Trajektoria 2 zaczyna w odległości ε od Trajektorii 1 (ε domyślnie wynosi 5×10⁻⁵, regulowane suwakiem Początkowej separacji). Ponieważ układ Lorenza jest chaotyczny, ta znikoma początkowa luka rośnie wykładniczo w czasie — „efekt motyla”. Ilekroć oba punkty zbliżą się do siebie na mniej niż 10⁻³, kod mierzy rzeczywisty wzrost odległości od ostatniej renormalizacji, dodaje go do bieżącej sumy i ponownie renormalizuje trajektorię 2 do odległości ε; podzielenie bieżącej średniej przez dt daje szacunek na żywo największego wykładnika Lapunowa λ₁ pokazanego w statystykach.
Co się dzieje, gdy zmienię ρ (rho)?
Zmniejszanie ρ wyprowadza układ z chaosu: poniżej ρ≈24,74 punkty stałe obu płatów znów stają się stabilne, a trajektorie zwijają się spiralnie w jeden z nich zamiast wiecznie błądzić, co pokazuje preset „Blisko bifurkacji (ρ=24)”. Zwiększenie ρ do konkretnych wartości może dać okna okresowe wewnątrz w innym razie chaotycznego zachowania — preset „Okresowy (ρ=99,96)” trafia w jedno z takich okien, gdzie trajektoria ustala się w powtarzającej się pętli zamiast chaotycznego błądzenia.
Dlaczego kształt wygląda jak motyl?
Atraktor Lorenza ma dwa „skrzydła” odpowiadające dwóm niestabilnym punktom stałym równań. Trajektoria krąży wokół jednego skrzydła przez pewien czas, po czym nieprzewidywalnie przechodzi na drugie i krąży wokół niego, tworząc dwupłatową, motylkową sylwetkę. Ponieważ przełączanie między skrzydłami zależy czule od dokładnej trajektorii, żadne dwa przebiegi (ani nawet dwie trajektorie renderowane tutaj obok siebie) nie kreślą dokładnie tej samej ścieżki, mimo że obie pozostają w tym samym ogólnym motylkowym atraktorze.
Trójwymiarowy atraktor motyla Lorenza z interaktywnymi suwakami parametrów (sigma, rho, beta). Uruchom dwie bliskie trajektorie, aby zwizualizować efekt motyla: wykładnicze rozbieganie kwantyfikowane wykładnikami Lapunowa.
3D · Renderer Three.js / WebGL · cel 60 FPS · działa w całości po stronie klienta, bez instalacji