A differential algorithm for the Lyapunov spectrum
arXiv:1008.3368 · doi:10.1016/j.physd.2011.04.007
Abstract
We present a new algorithm for computing the Lyapunov exponents spectrum based on a matrix differential equation. The approach belongs to the so called continuous type, where the rate of expansion of perturbations is obtained for all times, and the exponents are reached as the limit at infinity. It does not involve exponentially divergent quantities so there is no need of rescaling or realigning of the solution. We show the algorithm's advantages and drawbacks using mainly the example of a particle moving between two contracting walls.
11 pages with 8 figures, 3rd version with corrected typos, a new numerical example and slightly expanded conclusions/introduction
References in corpus (4)
Cited by in corpus (5)
- Free motion around black holes with discs or rings: between integrability and chaos - II
- Free motion around black holes with discs or rings: between integrability and chaos - III
- Oblique magnetic fields and the role of frame dragging near rotating black hole
- Dynamics and non-integrability of the variable-length double pendulum: exploring chaos and periodicity via the Lyapunov refined maps
- Transition from regular to chaotic motion in black hole magnetospheres