Computing Lyapunov Exponents using Weighted Birkhoff Averages
arXiv:2409.08496 · doi:10.1088/1751-8121/adfac1
Abstract
The Lyapunov exponents of a dynamical system measure the average rate of exponential stretching along an orbit. Positive exponents are often taken as a defining characteristic of chaotic dynamics. However, the standard orthogonalization-based method for computing Lyapunov exponents converges slowly -- if at all. Many alternatively techniques have been developed to distinguish between regular and chaotic orbits, though most do not compute the exponents. We compute the Lyapunov spectrum in three ways: the standard method, the weighted Birkhoff average (WBA), and the ``mean exponential growth rate for nearby orbits'' (MEGNO). The latter two improve convergence for nonchaotic orbits, but the WBA is fastest. However, for chaotic orbits the three methods convergence at similar, slow rates. Though the original MEGNO method does not compute Lyapunov exponents, we show how to reformulate it as a weighted average that does.
10 pages, 10 figures
References in corpus (13)
- On the Implementation of the 0-1 Test for Chaos
- The Lyapunov Characteristic Exponents and their computation
- Geometrical properties of local dynamics in Hamiltonian systems: the Generalized Alignment Index (GALI) method
- Computing Lyapunov spectra with continuous Gram-Schmidt orthonormalization
- Lyapunov Exponents without Rescaling and Reorthogonalization
- A comparison of different indicators of chaos based on the deviation vectors. Application to symplectic mappings
- Birkhoff Averages and Rotational Invariant Circles for Area-Preserving Maps
- Quantitative Quasiperiodicity
- How chaotic are strange nonchaotic attractors
- Distinguishing between Regular and Chaotic orbits of Flows by the Weighted Birkhoff Average
- Performance analysis of indicators of chaos for nonlinear dynamical systems
- Exponential convergence of weighted Birkhoff average
- Resonance and Weak Chaos in Quasiperiodically-Forced Circle Maps