Computing Lyapunov exponents based on the solution expression of the variational system

A simple discrete QR algorithm based on a solution expression of the variational equation of a dynamical system is presented for computing the Lyapunov exponents of n-dimensional continuous dynamical systems. The developed numerical scheme of study is based on a time integration using a constant time-step fourth-order Adams–Bashforth method. Numerical results are presented for a Lorenz system with known Lyapunov exponents, and higher dimensional dynamical systems. The algorithm proposed to compute the Lyapunov exponents is found to be robust, computationally efficient, and stable for a sufficiently small step-size h.

[1]  V. I. Oseledec.  A multiplicative ergodic theorem: Lyapunov characteristic num-bers for dynamical systems , 1968 .

[2]  F. Takens,et al.  On the nature of turbulence , 1971 .

[3]  J. Yorke,et al.  Chaotic behavior of multidimensional difference equations , 1979 .

[4]  W. G. Price,et al.  Transition to chaos in a fluid motion system , 2005 .

[5]  D. Ruelle,et al.  Ergodic theory of chaos and strange attractors , 1985 .

[6]  Cleve B. Moler,et al.  Nineteen Dubious Ways to Compute the Exponential of a Matrix, Twenty-Five Years Later , 1978, SIAM Rev..

[7]  F. Christiansen,et al.  Computing Lyapunov spectra with continuous Gram - Schmidt orthonormalization , 1996, chao-dyn/9611014.

[8]  Otto E. Rössler,et al.  Bifurcation theory and applications in scientific disciplines , 1979 .

[9]  J. Douglas Faires,et al.  Numerical Analysis , 1981 .

[10]  Jia Lu,et al.  Computing Lyapunov exponents of continuous dynamical systems: method of Lyapunov vectors , 2005 .

[11]  A. Wolf,et al.  Determining Lyapunov exponents from a time series , 1985 .

[12]  Hermann Haken,et al.  At least one Lyapunov exponent vanishes if the trajectory of an attractor does not contain a fixed point , 1983 .

[13]  Remarks on the Time Dependent Periodic¶Navier–Stokes Flows on a Two-Dimensional Torus , 1999 .

[14]  Ulrich Parlitz,et al.  Comparison of Different Methods for Computing Lyapunov Exponents , 1990 .

[15]  E. Lorenz.  Deterministic nonperiodic flow , 1963 .

[16]  G. Benettin,et al.  Lyapunov Characteristic Exponents for smooth dynamical systems and for hamiltonian systems; a method for computing all of them. Part 1: Theory , 1980 .