Jak změřit, že je systém chaotický?

Praktický průvodce Ljapunovovým exponentem — co to je, jak ho číst a které tři omyly vedou k tomu, že se čte špatně. Ukázáno na konkrétním příkladu z rotační dynamiky, se skutečnými čísly. Česká verze anglického textu; čísla a tvrzení jsou shodná, stavba vět je česká.

Otázka za otázkou

„Je tenhle systém chaotický?“ se ptáme mnohem častěji, než na to správně odpovídáme, protože slovo chaotický se používá pro tři různé věci, které je třeba držet od sebe:

  1. Vypadá to složitě. Hodně pohybu, okem se to nedá sledovat.
  2. Je to nestabilní. Malé poruchy rostou — aspoň lokálně, poblíž nějaké rovnováhy.
  3. Je to skutečně nepředvídatelné. Chyby v počátečním stavu rostou exponenciálně, takže dlouhodobá předpověď má tvrdý horizont.

Jen to třetí znamená „chaos“ v technickém smyslu, a první dvě se za něj běžně zaměňují. Systém může vypadat divoce a být dokonale předvídatelný; systém může mít nestabilní rovnováhu a přitom být globálně pravidelný. Nástroj, který ty tři věci od sebe odděluje, je Ljapunovův exponent a smyslem tohoto textu je naučit vás ho číst plynule.

Definice, která je použitelná

Vezměte systém a k němu virtuální druhou kopii, která startuje v nekonečně malé vzdálenosti δ0. Nechte oba vyvíjet. Měřte odstup δ(t). Ljapunovův exponent λ je průměrná rychlost exponenciálního růstu:

δ(t) ∼ δ0 eλt.

  • λ > 0 → chaotický. Horizont předpovědi ~ 1/λ: to je časové měřítko, na němž počáteční chyba vzroste e-krát. Za několika jeho násobky je předpověď bezcenná, ať je model jakkoli dobrý.
  • λ = 0 (při omezeném pohybu) → pravidelný / integrabilní. Chyby rostou pomalu (polynomiálně), předpověď zůstává zvládnutelná libovolně dlouho.
  • λ < 0 → systém se usazuje na atraktoru (disipativní).

U konzervativního (celkovou energii zachovávajícího systému bez disipace) mechanického systému přicházejí exponenty v párech  ± λ, jejichž součet je nula — užitečná kontrola, k níž se vrátíme.

To číslo se počítá, ne odhaduje: standardní metoda (Benettin a kol., 1980) nese podél trajektorie malou sadu poruchových vektorů, pravidelně je znovu normuje a sčítá jejich rychlosti změny velikosti. Je robustní a desítky let stará; inženýrská výzva je zde správně číst její výstup, nikoliv ho vyrobit.

Řešený příklad: převalující se těleso v silném couplingu s centrálním tělesem, malá sonda u velké planety

Vezměte tuhé těleso — představte si malou sondu — které se převaluje poblíž mnohem větší hmoty. Tři konfigurace, tři odpovědi. Čísla níže pocházejí z otevřené referenční implementace (odkazy na konci); jsou skutečná, ne ilustrační.

Případ 1 — volné převalování (křídlová matice). Nesymetrické těleso roztočené kolem své prostřední osy, bez vnějšího momentu. To je ta slavná nestabilita překlápění (efekt tenisové rakety). Vypadá chaoticky. Naměřený Ljapunovův exponent: ≈ 0 (zbytek 2.3 × 10 − 4 s − 1, který se s delší integrací dál zmenšuje). Verdikt: není chaotický. Převalování má řešení v uzavřeném tvaru; je plně předvídatelné. To je omyl č. 1 — vypadá to složitě — přistižený při činu.

Případ 2 — mírné kolébání (librace). Těleso se kolébá sem a tam kolem stabilní orientace jako kyvadlo, které nikdy nepřejde přes vrchol. Naměřený exponent: 2.4 × 10 − 6 s − 1, víc než 200× pod frekvencí kolébání. Verdikt: není chaotický. Kmitat není rozbíhat se. To je omyl č. 2 — záměna frekvence kmitání za rychlost nestability.

Případ 3 — převalování + gravitačně-gradientní vazba. Totéž převalování kolem prostřední osy, ale nyní se zapnutou vazbou mezi orientací a polohou (gravitační gradient). Naměřený exponent: 1.5 × 10 − 5 s − 1 > 0. Verdikt: chaotický — slabě. Horizont předpovědi 1/λ je asi 18 hodin; chaos je tenká vrstva, ne divoká bouře, ale je skutečný a kladný. Skutečná nepředvídatelnost se objevila teprve tehdy, když vazba zničila jednu zachovávanou veličinu.

To je diagnostika v akci: totéž těleso, tři režimy, a Ljapunovův exponent je jediná ze tří intuic „je to chaotické?“, která trefí všechny tři případy.

Tři omyly, řečeno rovnou

Polní průvodce k uschování:

Příznak Špatné čtení Správné čtení
Pohyb vypadá divoce/složitě „Je to chaotické“ Složitost ≠ chaos. Křídlová matice je integrabilní. Změřte λ.
Existuje nestabilní rovnováha „Je to chaotické“ Lokální nestabilita ≠ globální chaos. Prostřední osa bez momentu je nestabilní, a přesto je pohyb integrabilní.
Něco rychle kmitá „λ ≈ ta frekvence“ Frekvence kmitání a λ jsou různé pozorovatelné veličiny. Librující systém má λ ≈ 0 při jakékoli rychlosti kolébání.

Jak číst spočtené Ljapunovovo spektrum a neošálit sám sebe

Tři praktická varování, naučená po zlém na příkladu výše:

1. Integrabilní systémy konvergují pomalu — nepřitahujte práh příliš. U pravidelného systému klesá exponent na konečném čase ke své pravé hodnotě 0 jen jako O(log T/T). Uvidíte něco jako 2 × 10 − 4 při délce běhu, u níž jste doufali v 10 − 8. To není selhání vašeho kódu; to je fyzika. Otázku „je to nula?“ posuzujte vůči relevantní dynamické frekvenci systému, ne vůči doslovné nule — a podle toho nastavte bránu prošel/neprošel (my používáme 10 − 3, ne 10 − 8).

2. Pomalá konvergence je obvykle vlastní dynamice, ne artefakt souřadnic. Láká svést pomalou konvergenci na volbu proměnných a předpokládat, že ji vyřeší rafinovanější metoda. Na našem příkladu nevyřešila: přeformulování respektující geometrii dalo tentýž pomalý pokles, čímž se ukázalo, že jev je v dynamice, ne v účetnictví. Ověřte tenhle předpoklad, než podle něj jednáte — rafinovanější metoda možná řeší jiný problém než ten, z něhož jste ji vinili.

3. Zkontrolujte, že součet spektra je ~0 (konzervativní systémy). Pokud spočítáte celé spektrum a jeho součet není nula, pravděpodobně ho znečišťují falešné nefyzikální směry (v našem případě dva umělé módy z příliš velkého souřadnicového prostoru). Náprava je počítat ve skutečně dynamickém podprostoru. Párování  ± λ a nulový součet jsou vestavěné kontroly zadarmo — používejte je.

Co si odnést, mluvíme-li k inženýrům

Jedno provozní pravidlo, pokud si máte odnést jen jedno: než systém nazvete chaotickým, spočítejte λ, a než λ uvěříte, ověřte ho na případu, u něhož odpověď už znáte. Kalibrujte na integrabilní mezi (má dát ~0), na známé kmitavé mezi (má dát ≪ frekvence) a na známé chaotické mezi (má dát stabilní kladnou hodnotu). Ljapunovovský řetězec, který kdy běžel jen na zkoumaném systému, je teploměr, který nikdy neviděl led ani vařící vodu. Tříbodová kalibrace výše je celá disciplína v jednom návyku.

A hlubší poučení, které stojí za to odnést si dál než k sondě: systém může být plně deterministický — řízený přesnými, známými rovnicemi — a přesto být v praxi nepředvídatelný, protože informace potřebná k jeho předpovědi roste rychleji, než ji jakýkoli přístroj dokáže dodat. Znát zákon není totéž co umět ho používat neomezeně dlouho. Ljapunovův exponent je číslo, které přesně říká, kde ta hranice leží.

Dál do hloubky. Tento text je pohled praktika. Odvození z prvních principů — Ljapunovův algoritmus, kritérium integrability a plnou tříbodovou kalibraci s každým číslem — najdete v kanonické kapitole docs/canonical/cs/07-ljapunovovo-spektrum.md v repozitáři https://github.com/lege-artis/jwst-l2-coupled-dynamics (anglický originál: docs/canonical/en/07-lyapunov-spectrum.md). Chcete-li nástroje spustit sami — kterou funkci volat, jak číst výstup, výše uvedené zrady v podobě kódu — viz inženýrský návod docs/engineer/en/07-lyapunov-howto.md v témže repozitáři. To je naše referenční dokumentace k tomuto textu.

Odkazy

  • G. Benettin, L. Galgani, A. Giorgilli & J.-M. Strelcyn, „Lyapunov Characteristic Exponents for smooth dynamical systems and for Hamiltonian systems; a method for computing all of them,“ Meccanica 15, 9–30 (1980). DOI 10.1007/BF02128236
  • C. Skokos, „The Lyapunov Characteristic Exponents and Their Computation,“ in Dynamics of Small Solar System Bodies and Exoplanets, Lecture Notes in Physics 790, 63–135 (2010). DOI 10.1007/978-3-642-04458-8_2
  • L. D. Landau & E. M. Lifšic, Mechanika (Course of Theoretical Physics, sv. 1),
    1. vyd., Pergamon 1976, §37.
  • P. C. Hughes, Spacecraft Attitude Dynamics, Dover 2004.
  • lege-artis/jwst-l2-coupled-dynamics v0.3.1 — github.com/lege-artis/jwst-l2-coupled-dynamics (kód Apache-2.0, dokumentace CC-BY-SA-4.0).

Součást série o aplikované matematice ve výuce: nástroje, které z mlhavé inženýrské otázky udělají číslo, jež lze spočítat a ověřit. Připravujeme navazující text o tom, jak tutéž kalibrační disciplínu uplatnit na statistické testy v posuzování výzkumu.