Krylov subspace methods for discrete-time algebraic Riccati equations

We apply the Krylov subspace methods to large-scale discrete-time algebraic Riccati equations. The solvability of the projected algebraic Riccati equation is not assumed but is shown to be inherited from the original equation. Solvability in terms of stabilizability, detectability, stability radius of the associated Hamiltonian matrix and perturbation theory are considered. We pay particular attention to the stabilizing and the positive semi-definite properties of approximate solutions. Illustrative numerical examples are presented.