← ⚛ Fyzika

𝚲

Lorenzův atraktor

Tři vázané ODR, tři parametry a podivný atraktor, který se nikdy neopakuje. Interaktivní simulace deterministického chaosu — efekt motýlích křídel zviditelněný.

V roce 1963 zjednodušil meteorolog MIT Edward Lorenz atmosferickou konvekci na tři vázané diferenciální rovnice. Zjistil, že jejich řešení — přes veškerou zdánlivou jednoduchost — nikdy nesplývají do periodických orbit. Místo toho sledují podivný atraktor: fraktální strukturu nekonečné topologické složitosti, ohraničenou v prostoru, ale aperiodickou v čase.

Interaktivní simulace výše umožňuje prozkoumávat prostor parametrů — σ (Prandtlovo číslo), ρ (Rayleighovo číslo), β (poměr stran) — a přímé pozorování citlivé závislosti prostřednictvím trajektorie Ducha.

Historický kontext — konvekce a nepředvídatelnost

Lorenz studoval zjednodušený model Rayleighovy-Bénardovy konvekce. Chaotické chování objevil náhodou: restartem výpočtu z výtisku (zaokrouhleného na tři desetinná místa místo šesti) vznikla zcela odlišná trajektorie. Rozdíl 0,000127 vzrostl za méně než dva simulované dny natolik, aby pokryl celý atraktor.

Pevné body a bifurkace

Lorenzova soustava má tři pevné body: počátek (0, 0, 0) a dva symetrické body C± = (±√(β(ρ−1)), ±√(β(ρ−1)), ρ−1). Při klasickém ρ = 28 jsou všechny pevné body nestabilní — trajektorie chaoticky bloudí mezi dvěma laloky.

Geometrie podivného atraktoru

Lorenzův atraktor je podivný atraktor s fraktální dimenzí d_KY ≈ 2,062. Největší Ljapunovův exponent λ₁ ≈ 0,906 s⁻¹ kvantifikuje průměrnou rychlost exponenciální divergence. Ljapunovův čas τ_L = 1/λ₁ ≈ 1,1 s je charakteristická doba, po níž perturbace velikosti ε vzroste na řád 1.

∞ Lorenzův atraktor — interaktivní simulace

Systém se vyvíjí podle dx/dt = σ(y−x), dy/dt = x(ρ−z)−y, dz/dt = xy−βz. Klasické parametry σ = 10, ρ = 28, β = 8/3 produkují podivný atraktor. Stín (červený) začíná s Δx = 0,01 — zpočátku nerozlišitelný, exponenciálně divergující.

Lorenzovy rovnice (σ = 10, ρ = 28, β = 8/3):
dx/dt = σ(y − x)
dy/dt = x(ρ − z) − y
dz/dt = xy − βz
Pevné body (netriviální):
C± = (±√(β(ρ−1)), ±√(β(ρ−1)), ρ−1)

Táhněte na plátně pro orbitu. Barva kóduje velikost rychlosti (modrá = pomalá, červená = rychlá). Stínová trajektorie (červená čára) začíná Δx = 0,01 od primárního bodu — sledujte exponenciální divergenci.

λ Podivný atraktor — topologie a chaos

Ljapunovovy exponenty (σ=10, ρ=28, β=8/3):
λ₁ ≈ +0.906, λ₂ ≈ 0, λ₃ ≈ −14.572
Kaplan-Yorkova dimenze:
d_KY = 2 + λ₁/|λ₃| ≈ 2.062
Ljapunovův čas (zdvojení perturbace):
τ_L = 1/λ₁ ≈ 1.1 s

Atraktor je fraktální plocha dimenze ≈ 2,06 — ne 2D plocha a ne 3D objem, ale něco mezi. Sousední trajektorie se oddělují rychlostí e^{λ₁t}, zatímco atraktor jako celek je disipativní (součet Ljapunových exponentů < 0): objem ve fázovém prostoru se smršťuje rychlostí σ + 1 + β ≈ 13,7.