In this paper, we propose structured doubling algorithms for the computation of the weakly stabilizing Hermitian solutions of the continuous- and discrete-time algebraic Riccati equations, respectively. Assume that the partial multiplicities of purely imaginary and unimodular eigenvalues (if any) of the associated Hamiltonian and symplectic pencil, respectively, are all even and the C/DARE and the dual C/DARE have weakly stabilizing Hermitian solutions with property (P). Under these assumptions, we prove that if these structured doubling algorithms do not break down, then they converge to the desired Hermitian solutions globally and linearly. Numerical experiments show that the structured doubling algorithms perform efficiently and reliably.
ASJC Scopus subject areas
- Algebra and Number Theory
- Numerical Analysis
- Geometry and Topology
- Discrete Mathematics and Combinatorics