Lyapunov dimension

In the mathematics of dynamical systems, the concept of Lyapunov dimension was suggested by Kaplan and Yorke[1] for estimating the Hausdorff dimension of attractors. Further the concept has been developed and rigorously justified in a number of papers, and nowadays various different approaches to the definition of Lyapunov dimension are used. Remark that the attractors with noninteger Hausdorff dimension are called strange attractors.[2] Since the direct numerical computation of the Hausdorff dimension of attractors is often a problem of high numerical complexity, estimations via the Lyapunov dimension became widely spread. The Lyapunov dimension was named[3] after the Russian mathematician Aleksandr Lyapunov because of the close connection with the Lyapunov exponents.

Definitions

Consider a dynamical system ${\displaystyle {\big (}\{\varphi ^{t}\}_{t\geq 0},(U\subseteq \mathbb {R} ^{n},\|\cdot \|){\big )}}$, where ${\displaystyle \varphi ^{t}}$ is the shift operator along the solutions: ${\displaystyle \varphi ^{t}(u_{0})=u(t,u_{0})}$, of ODE ${\displaystyle {\dot {u}}=f({u})}$, ${\displaystyle t\leq 0}$, or difference equation ${\displaystyle {u}(t+1)=f({u}(t))}$, ${\displaystyle t=0,1,...}$, with continuously differentiable vector-function ${\displaystyle f}$. Then ${\displaystyle D\varphi ^{t}(u)}$ is the fundamental matrix of solutions of linearized system and denote by ${\displaystyle \sigma _{i}(t,u)=\sigma _{i}(D\varphi ^{t}(u)),\ i=1...n}$, singular values with respect to their algebraic multiplicity, ordered by decreasing for any ${\displaystyle u}$ and ${\displaystyle t}$.

Definition via finite-time Lyapunov dimension

The concept of finite-time Lyapunov dimension and related definition of the Lyapunov dimension, developed in the works by N. Kuznetsov,[4][5] is convenient for the numerical experiments where only finite time can be observed. Consider an analog of the Kaplan–Yorke formula for the finite-time Lyapunov exponents:

${\displaystyle d_{\rm {KY}}(\{{\rm {LE}}_{i}(t,u)\}_{i=1}^{n})=j(t,u)+{\frac {{\rm {LE}}_{1}(t,u)+\cdots +{\rm {LE}}_{j(t,u)}(t,u)}{|{\rm {LE}}_{j(t,u)+1}(t,u)|}},}$

${\displaystyle j(t,u)=\max\{m:\sum _{i=1}^{m}{\rm {LE}}_{i}(t,u)\geq 0\},}$

with respect to the ordered set of finite-time Lyapunov exponents ${\displaystyle \{{\rm {LE}}_{i}(t,u)\}_{i=1}^{n}=\{{\frac {1}{t}}\ln \sigma _{i}(t,u)\}_{i=1}^{n}}$ at the point ${\displaystyle u}$. The finite-time Lyapunov dimension of dynamical system with respect to invariant set ${\displaystyle K}$ is defined as follows

${\displaystyle \dim _{\rm {L}}(t,K)=\sup \limits _{u\in K}d_{\rm {KY}}(\{{\rm {LE}}_{i}(t,u)\}_{i=1}^{n}).}$

In this approach the use of the analog of Kaplan–Yorke formula is rigorously justified by the Douady–Oesterlè theorem,[6] which proves that for any fixed ${\displaystyle t>0}$ the finite-time Lyapunov dimension for a closed bounded invariant set ${\displaystyle K}$ is an upper estimate of the Hausdorff dimension:

${\displaystyle \dim _{\rm {H}}K\leq \dim _{\rm {L}}(t,K).}$

Looking for best such estimation ${\displaystyle \inf _{t>0}\dim _{\rm {L}}(t,K)=\liminf _{t\to +\infty }\sup \limits _{u\in K}\dim _{\rm {L}}(t,u)}$, the Lyapunov dimension is defined as follows:[4][5]

${\displaystyle \dim _{\rm {L}}K=\liminf _{t\to +\infty }\sup \limits _{u\in K}\dim _{\rm {L}}(t,u).}$

The possibilities of changing the order of the time limit and the supremum over set is discussed, e.g., in.[7][8]

Note that the above defined Lyapunov dimension is invariant under Lipschitz diffeomorphisms.[4][9]

Exact Lyapunov dimension

Let the Jacobian matrix ${\displaystyle Df(u_{\text{eq}})}$ at one of the equilibria have simple real eigenvalues: ${\displaystyle \{\lambda _{i}(u_{\text{eq}})\}_{i=1}^{n},\lambda _{i}(u_{\text{eq}})\geq \lambda _{i+1}(u_{\text{eq}})}$, then

${\displaystyle \dim _{\rm {L}}u_{\text{eq}}=d_{\rm {KY}}(\{\lambda _{i}(u_{\text{eq}})\}_{i=1}^{n}).}$

If the supremum of local Lyapunov dimensions on the global attractor, which involves all equilibria, is achieved at an equilibrium point, then this allows one to get analytical formula of the exact Lyapunov dimension of the global attractor (see corresponding Eden’s conjecture).

Definition via statistical physics approach and ergodicity

Following the statistical physics approach and assuming the ergodicity the Lyapunov dimension of attractor is estimated[1] by limit value of the local Lyapunov dimension ${\displaystyle \lim _{t\to +\infty }\dim _{\rm {L}}(t,u_{0})}$ of a typical trajectory, which belongs to the attractor. In this case ${\displaystyle \{\lim \limits _{t\to +\infty }{\rm {LE}}_{i}(t,u_{0})\}_{i}^{n}=\{{\rm {LE}}_{i}(u_{0})\}_{1}^{n}}$ and ${\displaystyle \dim _{\rm {L}}u_{0}=d_{\rm {KY}}(\{{\rm {LE}}_{i}(u_{0})\}_{i=1}^{n})=j(u_{0})+{\frac {{\rm {LE}}_{1}(u_{0})+\cdots +{\rm {LE}}_{j(u_{0})}(u_{0})}{|{\rm {LE}}_{j(u_{0})+1}(u_{0})|}}}$. From a practical point of view, the rigorous use of ergodic Oseledec theorem, verification that the considered trajectory ${\displaystyle u(t,u_{0})}$ is a typical trajectory, and the use of corresponding Kaplan–Yorke formula is a challenging task (see, e.g. discussions in[10]). The exact limit values of finite-time Lyapunov exponents, if they exist and are the same for all ${\displaystyle u_{0}\in U}$, are called the absolute ones[3] ${\displaystyle \{\lim \limits _{t\to +\infty }{\rm {LE}}_{i}(t,u_{0})\}_{i}^{n}=\{{\rm {LE}}_{i}(u_{0})\}_{1}^{n}\equiv \{{\rm {LE}}_{i}\}_{1}^{n}}$ and used in the Kaplan–Yorke formula. Examples of the rigorous use of the ergodic theory for the computation of the Lyapunov exponents and dimension can be found in.[11][12][13]

References

1. Kaplan J., Yorke J. (1979). "Functional Differential Equations and Approximations of Fixed Points". Chaotic behavior of multidimensional difference equations. Springer. pp. 204–227.
2. Ruelle D.; Takens F. (1971). "On the nature of turbulence". Communications in Mathematical Physics. 20 (3): 167–192. Bibcode:1971CMaPh..20..167R. doi:10.1007/bf01646553.
3. Frederickson, F.; Kaplan, J.; Yorke, E.; Yorke, J. (1983). "The Liapunov dimension of strange attractors". Journal of Differential Equations. 49 (2): 185–207. Bibcode:1983JDE....49..185F. doi:10.1016/0022-0396(83)90011-6.
4. Kuznetsov, N.V. (2016). "The Lyapunov dimension and its estimation via the Leonov method". Physics Letters A. 380 (25–26): 2142–2149. arXiv:1602.05410. Bibcode:2016PhLA..380.2142K. doi:10.1016/j.physleta.2016.04.036.
5. Kuznetsov, N.V.; Leonov, G.A.; Mokaev, T.N.; Prasad, A.; Shrimali, M.D. (2018). "Finite-time Lyapunov dimension and hidden attractor of the Rabinovich system". Nonlinear Dynamics. 92 (2): 267–285. arXiv:1504.04723. doi:10.1007/s11071-018-4054-z.
6. Douady, A.; Oesterle, J. (1980). "Dimension de Hausdorff des attracteurs". Comptes Rendus de l'Académie des Sciences, Série A. 290 (24): 1135–1138.
7. Constantin, P.; Foias, C.; Temam, R. (1985). "Attractors representing turbulent flows". Memoirs of the American Mathematical Society. 53 (314): 1–67. doi:10.1090/memo/0314.
8. Eden, A.; Foias, C.; Temam, R. (1991). "Local and global Lyapunov exponents". Journal of Dynamics and Differential Equations. 3 (1): 133–177. Bibcode:1991JDDE....3..133E. doi:10.1007/bf01049491.
9. Kuznetsov, N.; Alexeeva, T.; Leonov, G. (2016). "Invariance of Lyapunov exponents and Lyapunov dimension for regular and irregular linearizations". Nonlinear Dynamics. 85 (1): 195–201. arXiv:1410.2016. doi:10.1007/s11071-016-2678-4.
10. P. Cvitanovic; R. Artuso; R. Mainieri; G. Tanner & G. Vattay (2017). Chaos: Classical and Quantum (PDF). Niels Bohr Institute.
11. Ledrappier, F. (1981). "Some relations between dimension and Lyapounov exponents". Communications in Mathematical Physics. 81 (2): 229–238. Bibcode:1981CMaPh..81..229L. doi:10.1007/bf01208896.
12. Benedicks, M.; Young, L.-S. (1993). "Sinai–Bowen–Ruelle measures for certain Henon maps". Inventiones Mathematicae. 112 (1): 541–576. doi:10.1007/bf01232446.
13. Kuznetsov, Nikolay; Reitmann, Volker (2020). Attractor Dimension Estimates for Dynamical Systems: Theory and Computation. Cham: Springer.