În teoria sistemelor dinamice , un exponent Lyapunov al unui sistem dinamic (determinist) într-un punct din spațiul de fază oferă o măsură a modului în care semnificativ orbitele sistemului sunt dependente de datele inițiale, caracterizând prezența dinamicii haotice . Exponenții Lyapunov măsoară în special viteza medie de îndepărtare a două orbite infinit de apropiate pentru perioade suficient de lungi.
Un număr de exponenți Lyapunov egali cu dimensiunea spațiului sunt asociați cu un punct din spațiul de fază; dacă exponentul maxim Lyapunov este {\ displaystyle \ lambda} , și dacă distanța {\ displaystyle \ delta \ mathbf {Z} _ {0}} între orbite este destul de mic, apoi vectorul {\ displaystyle \ delta \ mathbf {Z} _ {0}} are o evoluție în timp (rata de separare a celor două orbite) decât în timp {\ displaystyle t} mare este dat aproximativ de:
- {\ displaystyle | \ delta \ mathbf {Z} (t) | \ approx Ce ^ {\ lambda t} | \ delta \ mathbf {Z} _ {0} |}
De sine {\ displaystyle \ lambda} este pozitiv, atunci sistemul are o dependență sensibilă de datele inițiale (exponențial) și, prin urmare, este un sistem haotic. Momentul în care un sistem devine haotic este dat de reciprocitatea lui {\ displaystyle \ lambda} , și se numește timp caracteristic sau timp Lyapunov al sistemului. Reprezintă limita de predictibilitate a sistemului.
Hărți unidimensionale
Este {\ displaystyle I \ subset \ mathbb {R}} Și {\ displaystyle T: I \ rightarrow I} o funcție diferențiată și ia în considerare sistemul dinamic discret dat de iterația hărții {\ displaystyle T} . Se definește exponentul Lyapunov al punctului {\ displaystyle x_ {0}} , adică a orbitei {\ displaystyle \ {x_ {0}, x_ {1}, ..., x_ {n}, ... \}} , ca:
- {\ displaystyle \ lambda (x_ {0}): = \ lim _ {n \ rightarrow \ infty} {\ frac {1} {n}} \ sum _ {k = 0} ^ {n-1} \ log | T '(x_ {k}) |}
sau echivalent ca:
- {\ displaystyle \ lambda (x_ {0}): = \ lim _ {n \ rightarrow \ infty} {\ frac {1} {n}} \ log \ left | {\ frac {d} {dx}} T ^ {n} (x_ {0}) \ right |}
acolo unde există limita.
Pentru a justifica această definiție, se poate observa în primul rând că derivatul lui {\ displaystyle T} intr-un loc {\ displaystyle x_ {0}} dă rata la care puncte apropiate de {\ displaystyle x_ {0}} s-au îndepărtat după o singură iterație: dacă distanța inițială dintre două puncte apropiate de {\ displaystyle x_ {0}} Și {\ displaystyle \ delta} , după aplicare {\ displaystyle T} aceasta devine {\ displaystyle \ delta | T ^ {\ prime} (x_ {0}) |} , adică {\ displaystyle \ delta e ^ {\ log | T '(x_ {0}) |}} . În plus, produsul {\ displaystyle | T ^ {\ prime} (x_ {0}) | \ cdots | T ^ {\ prime} (x_ {n}) |} dă derivata iterației {\ displaystyle T ^ {n}} în sens {\ displaystyle x_ {0}} , de la care avem viteza cu care se apropie punctele {\ displaystyle x_ {0}} s-au îndepărtat după aceea {\ displaystyle n} iterații. Mai exact, dacă distanța inițială dintre două puncte apropiate de {\ displaystyle x_ {0}} Și {\ displaystyle \ delta} , după aplicare {\ displaystyle T ^ {n}} aceasta devine {\ displaystyle \ delta | T ^ {\ prime} (x_ {1}) | \ cdots | T ^ {\ prime} (x_ {n}) |} , adică:
- {\ displaystyle \ delta e ^ {\ log (| T '(x_ {0}) | \ cdots | T' (x_ {n}) |)} = \ delta e ^ {\ log | T '(x_ {0 }) | + ... + \ log | T '(x_ {n}) |}}
care poate fi scris (ținând cont de discursul inițial și de faptul că timpul pe care îl luăm în considerare este {\ displaystyle n} ) ca:
- {\ displaystyle \ delta e ^ {n {\ frac {\ log | T '(x_ {0}) | + ... + \ log | T' (x_ {n}) |} {n}}}}
Din aceste observații se concluzionează că, dacă există limita {\ displaystyle \ lambda (x_ {0})} pentru cantitate:
- {\ displaystyle {\ frac {\ log | T '(x_ {0}) | + ... + \ log | T' (x_ {n}) |} {n}}}
apoi pentru vremuri {\ displaystyle n} foarte mare, distanța dintre două orbite aproape de {\ displaystyle x_ {0}} a crescut cu un factor de multiplicare aproximativ egal cu {\ displaystyle e ^ {\ lambda (x_ {0}) n}} .
Hărți multidimensionale
Pentru o hartă {\ displaystyle f: \ mathbb {R} ^ {m} \ longrightarrow \ mathbb {R} ^ {m}} diferențiat și orbita sa poate fi definită {\ displaystyle m} exponenții lui Lyapunov {\ displaystyle \ lambda _ {1}, ..., \ lambda _ {m}} care măsoară viteza de separare de orbită în {\ displaystyle m} direcții ortogonale astfel încât de-a lungul direcției {\ displaystyle i} -distanțele dintre punctele apropiate de orbită vor evolua pe măsură ce {\ displaystyle \ delta _ {0} e ^ {\ lambda _ {i} n}} pentru {\ displaystyle n} mare. Prima direcție va fi aceea în care această viteză este maximă, a doua va fi aleasă ca cea a vitezei maxime în setul de direcții ortogonale față de prima și așa mai departe. În direcții care sunt combinații liniare de două direcții asociate cu exponenți Lyapunov diferiți, viteza de separare este stabilită de cel mai mare exponent Lyapunov.
Se definește exponentul Lyapunov asociat cu un punct {\ displaystyle x_ {0}} și o direcție {\ displaystyle v} ca viteza medie de separare a unui punct {\ displaystyle x} aproape {\ displaystyle x_ {0}} astfel încât vectorul de îmbinare {\ displaystyle x-x_ {0}} are direcția {\ displaystyle v} . După {\ displaystyle n} iterații distanța dintre {\ displaystyle F ^ {n} (x)} Și {\ displaystyle F ^ {n} (x_ {0})} care a fost inițial {\ displaystyle \ left \ | x-x_ {0} \ right \ |} a devenit aprox {\ displaystyle \ left \ | DF ^ {n} (x_ {0}) (x-x_ {0}) \ right \ |} , rata medie de creștere pentru fiecare pas este dată de:
- {\ displaystyle \ left ({\ frac {\ left \ | DF ^ {n} (x_ {0}) (x-x_ {0}) \ right \ |} {\ left \ | x-x_ {0} \ right \ |}} \ right) ^ {\ frac {1} {n}} = \ left \ | DF ^ {n} (x_ {0}) v \ right \ | ^ {\ frac {1} {n} }}
unde este {\ displaystyle v} este vectorul unitar al direcției {\ displaystyle (x-x_ {0})} . Dacă luăm în considerare logaritmul:
- {\ displaystyle L = \ log \ left \ | DF ^ {n} (x_ {0}) v \ right \ | ^ {\ frac {1} {n}} = {\ frac {1} {n}} \ log \ left \ | DF ^ {n} (x_ {0}) v \ right \ |}
se poate spune că sistemul a evoluat astfel încât distanța inițială {\ displaystyle \ delta _ {0}} a devenit {\ displaystyle \ delta _ {0} e ^ {nL}} . Cu toate acestea, media a fost făcută pe un număr finit de pași, dacă luăm în considerare întreaga traiectorie putem defini exponentul Ljapunov al {\ displaystyle x_ {0}} in directia {\ displaystyle v} ca rata medie de creștere exponențială după cum urmează:
- {\ displaystyle \ lambda (x_ {0}, v): = \ lim _ {n \ to \ infty} {\ frac {1} {n}} \ log \ left \ | DF ^ {n} (x_ {0 }) v \ right \ |}
Din această definiție deducem că dacă vectorul de legătură are direcția {\ displaystyle v} apoi distanța {\ displaystyle \ left \ | x-x_ {0} \ right \ |} evoluează ca {\ displaystyle Ce ^ {n \ lambda (x_ {0}, v)}} pentru {\ displaystyle n} mare.
Pentru a evalua valoarea valorii {\ displaystyle \ lambda} poate varia dacă luăm în considerare direcții diferite, se arată că {\ displaystyle \ lambda} poate presupune cel mult un număr de valori egal cu dimensiunea {\ displaystyle m} a spațiului și că pentru aproape toate punctele spațiului își asumă aceeași valoare: valoarea maximă.
Exemplu
În cele ce urmează arătăm un caz în care aproximarea liniară a {\ displaystyle F} rămâne mereu la fel. Luați în considerare sistemul dinamic discret dat de iterația hărții {\ displaystyle F (x): = Ax} cu {\ displaystyle A} matrice {\ displaystyle m \ times m} echipat cu {\ displaystyle m} valori proprii {\ displaystyle 0 \ leq \ lambda _ {1} \ leq \ lambda _ {2} \ leq ... \ leq \ lambda _ {m}} . La pasul n-o avem {\ displaystyle F ^ {n} (x) -F ^ {n} (x_ {0}) \ sim A ^ {n} (x-x_ {0})} , de aici și distanța de plecare {\ displaystyle \ delta _ {0} = \ left \ | x-x_ {0} \ right \ |} a devenit {\ displaystyle \ left \ | A ^ {n} (x-x_ {0}) \ right \ |} . Dacă transportatorul {\ displaystyle (x-x_ {0})} se află în spațiul auto asociat cu {\ displaystyle \ lambda _ {k}} asa de:
- {\ displaystyle \ left \ | A ^ {n} (x-x_ {0}) \ right \ | = | \ lambda _ {k} | ^ {n} \ left \ | x-x_ {0} \ right \ | = \ delta _ {0} e ^ {n \ log | \ lambda _ {k} |}}
Dacă transportatorul {\ displaystyle (x-x_ {0})} are o componentă non-nulă în spațiul automat asociată cu {\ displaystyle \ lambda _ {m}} (care este maximul valorilor proprii așa cum le-am numerotat), atunci acesta poate fi exprimat {\ displaystyle (x-x_ {0})} ca o combinație liniară :
- {\ displaystyle x-x_ {0} = a_ {1} v_ {1} + ... + a_ {m} v_ {m}} cu {\ displaystyle a_ {m} \ neq 0}
unde este {\ displaystyle v_ {1}, ..., v_ {m}} este o bază ortonormală a vectorilor proprii (se presupune pentru simplitate că există o astfel de bază). Asa de:
- {\ displaystyle A ^ {n} (x-x_ {0}) = \ lambda _ {1} ^ {n} a_ {1} v_ {1} + \ lambda _ {2} ^ {n} a_ {2} v_ {2} + ... + \ lambda _ {m} ^ {n} a_ {m} v_ {m} =}
- {\ displaystyle = \ lambda _ {m} ^ {n} \ left (\ left ({\ frac {\ lambda _ {1}} {\ lambda _ {m}}} \ right) ^ {n} a_ {1 } v_ {1} + \ left ({\ frac {\ lambda _ {2}} {\ lambda _ {m}}} \ right) ^ {n} a_ {2} v_ {2} + ... + \ left ({\ frac {\ lambda _ {m-1}} {\ lambda _ {m}}} \ right) ^ {n} a_ {m-1} v_ {m-1} + a_ {m} v_ { m} \ dreapta)}
Pentru a vă face o idee despre factorul mediu de expansiune pentru fiecare pas, se poate calcula limita mediei geometrice:
- {\ displaystyle \ lim _ {n \ to \ infty} \ left ({\ frac {\ left \ | A ^ {n} (x-x_ {0}) \ right \ |} {\ left \ | x-x_ {0} \ right \ |}} \ right) ^ {1 \ left / n \ right.}}
care din calculele anterioare se dovedește a fi egală cu {\ displaystyle \ left | \ lambda _ {m} \ right |} . De aici și distanța {\ displaystyle \ left \ | x-x_ {0} \ right \ |} va evolua mult timp ca {\ displaystyle C | \ lambda _ {m} | ^ {n} = Ce ^ {n \ log | \ lambda _ {m} |}} . Aceasta înseamnă toate punctele {\ displaystyle x} aproape de {\ displaystyle x_ {0}} pentru care vectorul de unire {\ displaystyle (x-x_ {0})} are o lungă componentă nulă {\ displaystyle v_ {m}} au o viteză medie asimptotică de separare (sau apropiere) de {\ displaystyle x_ {0}} determinat exclusiv de maximul valorilor proprii ale {\ displaystyle A} .
Calculul exponentului Lyapunov pe baza relațiilor stabilite mai sus prevede de fapt:
- {\ displaystyle \ lim _ {n \ to \ infty} {\ frac {1} {n}} \ log \ left \ | A ^ {n} {\ frac {x-x_ {0}} {\ left \ | x-x_ {0} \ right \ |}} \ right \ | = \ lim _ {n \ to \ infty} {\ frac {1} {n}} \ log \ left \ | A ^ {n} (x -x_ {0}) \ right \ | - {\ frac {1} {n}} \ log \ left \ | x-x_ {0} \ right \ | = \ log | \ lambda _ {m} |}
Cu un discurs analog se poate arăta că dacă vectorul de îmbinare {\ displaystyle x-x_ {0}} este ortogonală la spațiul egal raportat la valoarea proprie maximă, dar are o componentă diferită de zero față de a doua cea mai mare valoare proprie {\ displaystyle \ lambda _ {m-1}} atunci exponentul Lyapunov asociat cu această direcție este {\ displaystyle \ log | \ lambda _ {m-1} |} . Mai general, exponentul lui Lyapunov în {\ displaystyle x_ {0}} de-a lungul direcției {\ displaystyle v} este dat de logaritmul celei mai mari valori proprii {\ displaystyle \ lambda _ {k}} asociat cu un vector propriu cu privire la care {\ displaystyle v} nu este ortogonală.
Pentru a vizualiza intuitiv conceptul, putem considera o sferă infinitesimală în jurul punctului {\ displaystyle p} a unei orbite: aceasta după fiecare iterație a hărții {\ displaystyle F} este deformat într-un elipsoid obținut ca imagine a sferei prin intermediul aplicației liniare date de matricea iacobiană {\ displaystyle D (f ^ {n}) (p)} . Elipsoidul oferă informații despre comportamentul local al hărții, în special despre direcțiile în care se contractă sau extinde cel mai mult spațiul. Pot fi identificate axele principale ale acestui elipsoid care corespund direcțiilor de contracție sau expansiune. Cu toate acestea, la fiecare iterație transformarea liniară este diferită, la fel și vectorii proprii și valorile proprii și, prin urmare, axele și forma elipsoidului. Teorema lui Oseledec asigură că pentru aproape fiecare punct acțiunea transformărilor liniare dată de diferențiale {\ displaystyle DF (x_ {n})} , calculat de-a lungul traiectoriei, în medie tinde asimptotic să fie echivalent cu acțiunea aceleiași matrice cu {\ displaystyle m} valori proprii ale căror logaritmi dau exponenții Lyapunov și ai căror vectori proprii dau direcțiile de expansiune și contracție corespunzătoare axelor unui elipsoid „mediu”.
Bibliografie
- ( EN ) R. Temam, Sisteme dinamice dimensionale infinite în mecanică și fizică , Cambridge: Springer-Verlag, 1988.
- ( EN ) J. Kaplan și J. Yorke , Comportamentul haotic al ecuațiilor diferențiale multidimensionale , în HO Peitgen și HO Walther (eds), Ecuații funcționale diferențiale și Aproximarea punctelor fixe , New York, Springer, 1979, ISBN 3-540 -09518 -7 .
- ( EN ) Cvitanović P., Artuso R., Mainieri R., Tanner G., Vattay G.; Haos: Institutul Niels Bohr clasic și cuantic , Copenhaga 2005.
Elemente conexe
linkuri externe