Extraction of Dynamical Equations from Chaotic Data

G. Rowlands and J. C. Sprott
Department of Physics, University of Warwick, Coventry CV47AL, England
Department of Physics, University of Wisconsin, Madison, WI 53706, USA
Received 26 September 1991
Revised manuscript received 29 January 1992
Accepted 7 February 1992

ABSTRACT

A method is described for extracting from a chaotic time series a system of equations whose solution reproduces the general features of the original data even when these are contaminated with noise. The equations facilitate calculation of fractal dimension, Lyapunov exponents and short-term predictions. The method is applied to data derived from numerical solutions of the logistic equation, the Henon equations with added noise, the Lorenz equations and the Rossler equations.

Ref: G. Rowlands and J. C. Sprott, Physica D 58, 251-259 (1992)

The complete paper is available in PDF format.

Return to Sprott's Books and Publications.

Fig. 1. Phase-space plots of the logistic equation (a) from original data and (b) after singular value decomposition.
[Figure 1a] (a)
[Figure 1b] (b)

Fig. 2. Phase-space plots of the Henon map (a) from original data and (b) after singular value decomposition.
[Figure 2a] (a)
[Figure 2b] (b)

Fig. 3. Phase-space plots of the Henon map (a) with noise added and (b) from a solution of the model equations fit to the noisy data. A comparison of (b) with fig. 2a shows that the method has completely removed the noise and restored the original.
[Figure 3a] (a)
[Figure 3b] (b)

Fig. 4. Three-dimensional phase-space plots of the Lorenz attractor showing that the topology of the attractor is preserved. (a) Original input data, (b) result of singular value decomposition, and (c) solution of the model equations.
[Figure 4a] (a)
[Figure 4b] (b)
[Figure 4c] (c)

Fig. 5. Three-dimensional phase-space plots of the Rossler attractor showing that the topology of the attractor is preserved. (a) Original input data, (b) result of singular value decomposition, and (c) solution of the model equations.
[Figure 5a] (a)
[Figure 5b] (b)
[Figure 5c] (c)

Fig. 6. (a) Variation of the nine largest coefficients of the model equations with the parameter r in the Lorenz equations, (b) along with a least squares fit of each coefficient to a cubic polynomial in r.
[Figure 6a] (a)
[Figure 6b] (b)

Fig. 7. Phase-space portraits for the Lorenz attractor with r=57 (a) obtained directly from singular value decomposition, and (b) resulting from solution of the model equations with coefficients calculated from least squares fits to cubic polynomials in r. [Figure 7]