Stability and Lyapunov stability of dynamical systems: A differential approach and a numerical method

https://doi.org/10.1016/0167-2789(87)90034-0Get rights and content

Abstract

A set of differential equations for the eigenvalues and eigenvectors of the stability matrix of a dynamical system, as well as for the Lyapunov exponents and the corresponding eigenvectors, is derived. The rate of convergence of the Lyapunov eigenvectors is shown to be exponential. The eigenvectors of the stability matrix can be grouped into sets, each spanning a subspace which converges at an exponential rate. It is demonstrated that, generically, the real parts of the eigenvalues of the stability matrix equal the corresponding Lyapunov exponents. This statement has been tested numerically. The values of the Lyapunov exponents, μi, are shown to be related to the corresponding finite time values of the Lyapunov exponents (e.g. those deduced from a finite time numerical simulation), μi(t), by: μi(t) = μi + (bi + ξi(t))/t. The bi's are constants and ξi(t) are “noise” terms of zero mean. This observation leads to a method of extrapolation, which has been used to predict Lyapunov exponents from a finite amount of data. It is shown that the use of the standard (numerical) methods to compute Lyapunov exponents introduces an error ai/t in the value of μi(t), where the ai's are constants. Thus the standard method has a rate of convergence which is the same as that of the exact μi(t)'s. Finally, we have shown how one can compute the eigenvectors associated with each of the eigenvalues of the stability matrix as well as the Lyapunov eigenvectors.

References (25)

  • J.D. Farmer

    Chaotic attractors of an infinite-dimensional dynamical system

    Physica

    (1982)
  • P. Grassberger et al.

    Characterization of strange attractors

    Phys. Rev. Lett.

    (1983)
  • B.J. Bayley

    Ph.D. Thesis

    (1986)
  • M.C. Mackey et al.

    Oscillation and chaos in physiological control systems

    Science

    (1977)
  • V.I. Oseledec

    A multiplicative ergodic theorem. Lyapunov characteristic numbers for dynamical systems

    Trans. Moscow Math. Soc.

    (1968)

    Ann. Math. Study

    (1977)
    M.S. Raghunathatan

    Israel J. Math.

    (1979)
    H. Kesten, in: Random Matrices and Their Applications, J.E. Cohen, H. Kesten and C.M. Newman, eds., in AMS Series in...
  • J.D. Farmer et al.

    The dimension of chaotic attractors

    Physica

    (1983)
  • D.A. Russell et al.

    Phys. Rev. Lett.

    (1980)
  • J. Froyland

    Lyapunov exponents for multidimensional orbits

    Phys. Lett. A

    (1983)
  • J. Kaplan et al.P. Frederickson et al.

    J. Diff. Eq.

    (1983)
    H. Mori

    Prog. Theor. Phys.

    (1980)
  • B.B. Mandelbrot

    Fractals-Form, Chance and Dimension

    (1977)
  • P. Grassberger et al.

    Estimation of the Kolmogorov entropy from a chaotic signal

    Phys. Rev. A

    (1983)
  • G. Bennetin et al.

    Phys. Rev. A

    (1976)
  • Cited by (266)

    • Dynamic snap-through of functionally graded shallow arches under rapid surface heating

      2022, Thin-Walled Structures
      Citation Excerpt :

      As mentioned earlier, there are different criteria which are introduced by researchers to recognize the buckling load of structures. To introduce some of the most important criteria, one may refer to Volmir criterion [32] (as the first criterion of detection for the moment of dynamic buckling), Budiansky–Roth [3] (or Budiansky–Hutchinson [4] that developed for structures with initial imperfection), phase-plane approach criterion [14], Hoff–Simitses criterion [2], Kounadis criterion [33], Lyapunov exponent criterion [34], and Kleiber–Katula–Saran criterion [35]. In this study Budiansky–Roth criterion that is a popular and widely accepted criterion among the researchers is used to estimate the dynamic buckling load.

    View all citing articles on Scopus
    View full text