Analysis of the rational Krylov subspace projection method for large-scale algebraic Riccati equations

Abstract

In the numerical solution of the algebraic Riccati equation $A^{\ast }X+XA-XB{B}^{\ast }X+C^{\ast }C=0$, where $A$ is large, sparse and stable, and $B$, $C$ have low rank, projection methods have recently emerged as a possible alternative to the more established Newton-Kleinman iteration. In spite of convincing numerical experiments, a systematic matrix analysis of this class of methods is still lacking. We derive new relations for the approximate solution, the residual and the error matrices, giving new insights into the role of the matrix $A-B{B}^{\ast }X$ and of its approximations in the numerical procedure. The new results provide theoretical ground for recently proposed modifications of projection methods onto rational Krylov subspaces.