Orthosymplectic integration of linear Hamiltonian systems
long-term integrationsLyapunov exponentsnumerical examplesorthosymplectic integrationtimesteppingunstable linear Hamiltonian systems
Linear ordinary differential equations and systems (34A30) Characteristic and Lyapunov exponents of ordinary differential equations (34D08) Dynamical systems and ergodic theory (37-XX) Numerical methods for initial value problems involving ordinary differential equations (65L05) Numerical investigation of stability of solutions to ordinary differential equations (65L07)
A continuous, orthogonal and symplectic factorization procedure for integrating unstable linear Hamiltonian systems is presented. The method is based on the development of an orthogonal, symplectic change of variables to block triangular Hamiltonian form. Integration is thus carried out within the class of linear Hamiltonian systems. Use of an appropriate timestepping strategy ensures that the symplectic pairing of eigenvalues is automatically preserved. For long-term integrations, which are necessary in the calculation of Lyapunov exponents, the properties of such a symplectic framework can be expected to yield improved estimates. The method is illustrated and compared with other techniques by means of numerical experiments on some equations.
- SVD algorithms to approximate spectra of dynamical systems
- Computation of functions of Hamiltonian and skew-symmetric matrices
- Lyapunov exponents and localization in randomly layered media
- A class of ergodic solutions of nonlinear differential equations and numerical treatment
- Applications of the Cayley approach in the numerical solution of matrix differential systems on quadratic groups
- Preserving Poisson structure and orthogonality in numerical integration of differential equations
- Rotation number and exponential dichotomy for linear Hamiltonian systems: from theoretical to numerical results
- Generalized inverse eigenvalue problems for Hermitian and \(J\)-Hamiltonian/skew-Hamiltonian matrices
- A structure preserving flow for computing Hamiltonian matrix exponential
- The generalized inverse eigenvalue problem of Hamiltonian matrices and its approximation
- A structure preserving approximation method for Hamiltonian exponential matrices
- On the Floquet-Magnus expansion: applications in solid-state nuclear magnetic resonance and physics
- Computing the Maslov index of solitary waves. II: Phase space with dimension greater than four
- Preserving geometric properties of the exponential matrix by block Krylov subspace methods
- Smooth singular value decomposition on the symplectic group and Lyapunov exponents approximation
This page was built for publication: Orthosymplectic integration of linear Hamiltonian systems
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1363213)
