Numerical Methods That Preserve a Lyapunov Function for Ordinary Differential Equations

Yadira Hernández-Solano and Miguel Atencia ()
Additional contact information
Yadira Hernández-Solano: Departamento de Matemática Aplicada, Universidad de Málaga, 29071 Málaga, Spain
Miguel Atencia: Departamento de Matemática Aplicada, Universidad de Málaga, 29071 Málaga, Spain

Mathematics, 2022, vol. 11, issue 1, 1-21

Abstract: The paper studies numerical methods that preserve a Lyapunov function of a dynamical system, i.e., numerical approximations whose energy decreases, just like in the original differential equation. With this aim, a discrete gradient method is implemented for the numerical integration of a system of ordinary differential equations. In principle, this procedure yields first-order methods, but the analysis paves the way for the design of higher-order methods. As a case in point, the proposed method is applied to the Duffing equation without external forcing, considering that, in this case, preserving the Lyapunov function is more important than the accuracy of particular trajectories. Results are validated by means of numerical experiments, where the discrete gradient method is compared to standard Runge–Kutta methods. As predicted by the theory, discrete gradient methods preserve the Lyapunov function, whereas conventional methods fail to do so, since either periodic solutions appear or the energy does not decrease. Moreover, the discrete gradient method outperforms conventional schemes when these do preserve the Lyapunov function, in terms of computational cost; thus, the proposed method is promising.

Keywords: geometric numerical integration; dynamical systems; Lyapunov function; stability; numerical methods; ordinary differential equations; discrete gradient (search for similar items in EconPapers)
JEL-codes: C (search for similar items in EconPapers)
Date: 2022
References: View references in EconPapers View complete reference list from CitEc
Citations:

Downloads: (external link)
https://www.mdpi.com/2227-7390/11/1/71/pdf (application/pdf)
https://www.mdpi.com/2227-7390/11/1/71/ (text/html)

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:gam:jmathe:v:11:y:2022:i:1:p:71-:d:1014282

Access Statistics for this article

Mathematics is currently edited by Ms. Emma He

More articles in Mathematics from MDPI
Bibliographic data for series maintained by MDPI Indexing Manager ().