The calculation of Lyapunov spectra
Physica D
Unitary integrators and applications to continuous orthonormalization techniques
SIAM Journal on Numerical Analysis
Matrix computations (3rd ed.)
A New Software Package for Linear Differential-Algebraic Equations
SIAM Journal on Scientific Computing
On the Compuation of Lyapunov Exponents for Continuous Dynamical Systems
SIAM Journal on Numerical Analysis
On Smooth Decompositions of Matrices
SIAM Journal on Matrix Analysis and Applications
Computer Methods for Ordinary Differential Equations and Differential-Algebraic Equations
Computer Methods for Ordinary Differential Equations and Differential-Algebraic Equations
Smoothness and Periodicity of Some Matrix Decompositions
SIAM Journal on Matrix Analysis and Applications
Lyapunov Spectral Intervals: Theory and Computation
SIAM Journal on Numerical Analysis
On the Error in QR Integration
SIAM Journal on Numerical Analysis
Differential-Algebraic Systems: Analytical Aspects and Circuit Applications
Differential-Algebraic Systems: Analytical Aspects and Circuit Applications
Computation of a few Lyapunov exponents for continuous and discrete dynamical systems
Applied Numerical Mathematics
Adjoint pairs of differential-algebraic equations and Hamiltonian systems
Applied Numerical Mathematics
On the Error in the Product QR Decomposition
SIAM Journal on Matrix Analysis and Applications
SIAM Journal on Numerical Analysis
Journal of Computational and Applied Mathematics
Hi-index | 0.00 |
In this paper, we propose and investigate numerical methods based on QR factorization for computing all or some Lyapunov or Sacker---Sell spectral intervals for linear differential-algebraic equations. Furthermore, a perturbation and error analysis for these methods is presented. We investigate how errors in the data and in the numerical integration affect the accuracy of the approximate spectral intervals. Although we need to integrate numerically some differential-algebraic systems on usually very long time-intervals, under certain assumptions, it is shown that the error of the computed spectral intervals can be controlled by the local error of numerical integration and the error in solving the algebraic constraint. Some numerical examples are presented to illustrate the theoretical results.