Lyapunov Exponents without Rescaling and Reorthogonalization
نویسندگان
چکیده
Govindan Rangarajan,1,* Salman Habib,2,† and Robert D. Ryne3,‡ 1Department of Mathematics and Center for Theoretical Studies, Indian Institute of Science, Bangalore 560 012, India 2T-8, Theoretical Division, MS B285, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 3LANSCE-1, LANSCE Division, MS H817, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 (Received 27 May 1997; revised manuscript received 9 March 1998)
منابع مشابه
Scaling and interleaving of subsystem Lyapunov exponents for spatio-temporal systems.
The computation of the entire Lyapunov spectrum for extended dynamical systems is a very time consuming task. If the system is in a chaotic spatio-temporal regime it is possible to approximately reconstruct the Lyapunov spectrum from the spectrum of a subsystem by a suitable rescaling in a very cost effective way. We compute the Lyapunov spectrum for the subsystem by truncating the original Jac...
متن کاملAlternative algorithm for the computation of Lyapunov spectra of dynamical systems.
Recently a new method for the computation of Lyapunov exponents that does not require rescaling and reorthogonalization was proposed [Rangarajan, Habib, and Ryne, Phys. Rev. Lett. 80, 3747 (1998)]. In this paper we make a detailed numerical comparison of the new method and a standard algorithm, as regards accuracy and efficiency, by applying them to some typical two-, three-, and four-dimension...
متن کاملComment on "ruling out chaos in compact binary systems".
Intuitively, black hole binaries are obvious candidates for chaotic dynamics. Firstly, they are highly non-linear systems and chaos is an expression of extreme nonlinearity. Secondly, there are isolated unstable orbits around a Schwarzschild black hole, and unstable orbits are a red flag for the onset of chaos. Surely the orbits only become more unstable and more complex when there are two blac...
متن کاملA comparative study of computation of Lyapunov spectra with different algorithms
In this paper we make a detailed numerical comparison between three algorithms for the computation of the full Lyapunov spectrum as well as the associated eigen-vectors of general dynamical systems. They are : (a) the standard method, (b) a differential formulation of the standard method, and (c) a new algorithm which does not require rescaling and reorthogonalization. We also bring out the rel...
متن کاملJacobian Free Computation of Lyapunov Exponents
The purpose of this paper is to present new algorithms to approximate Lyapunov exponents of nonlinear differential equations, without using Jacobian matrices. We first derive first order methods for both continuous and discrete QR approaches, and then second order methods. Numerical testing is given, showing considerable savings with respect to existing implementations.
متن کامل