Finite Difference Models of Dynamical Systems with Quadratic Right-Hand Side

Mikhail Malykh (), Mark Gambaryan, Oleg Kroytor and Alexander Zorin
Additional contact information
Mikhail Malykh: Department of Computational Mathematics and Artificial Intelligence, RUDN University, 117198 Moscow, Russia
Mark Gambaryan: Department of Computational Mathematics and Artificial Intelligence, RUDN University, 117198 Moscow, Russia
Oleg Kroytor: Department of Computational Mathematics and Artificial Intelligence, RUDN University, 117198 Moscow, Russia
Alexander Zorin: Department of Computational Mathematics and Artificial Intelligence, RUDN University, 117198 Moscow, Russia

Mathematics, 2024, vol. 12, issue 1, 1-17

Abstract: Difference schemes that approximate dynamic systems are considered discrete models of the same phenomena that are described by continuous dynamic systems. Difference schemes with t -symmetry and midpoint and trapezoid schemes are considered. It is shown that these schemes are dual to each other, and, from this fact, we derive theorems on the inheritance of quadratic integrals by these schemes (Cooper’s theorem and its dual theorem on the trapezoidal scheme). Using examples of nonlinear oscillators, it is shown that these schemes poses challenges for theoretical research and practical application due to the problem of extra roots: these schemes do not allow one to unambiguously determine the final values from the initial values and vice versa. Therefore, we consider difference schemes in which the transitions from layer to layer in time are carried out using birational transformations (Cremona transformations). Such schemes are called reversible. It is shown that reversible schemes with t -symmetry can be easily constructed for any dynamical system with a quadratic right-hand side. As an example of such a dynamic system, a top fixed at its center of gravity is considered in detail. In this case, the discrete theory repeats the continuous theory completely: (1) the points of the approximate solution lie on some elliptic curve, which at Δ t → 0 turns into an integral curve; (2) the difference scheme can be represented using quadrature; and (3) the approximate solution can be represented using an elliptic function of a discrete argument. The last section considers the general case. The integral curves are replaced with closures of the orbits of the corresponding Cremona transformation as sets in the projective space over R . The problem of the dimension of this set is discussed.

Keywords: finite difference method; dynamical systems; Cremona transformations (search for similar items in EconPapers)
JEL-codes: C (search for similar items in EconPapers)
Date: 2024
References: View complete reference list from CitEc
Citations:

Downloads: (external link)
https://www.mdpi.com/2227-7390/12/1/167/pdf (application/pdf)
https://www.mdpi.com/2227-7390/12/1/167/ (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:12:y:2024:i:1:p:167-:d:1313449

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 ().