An efficient algorithm for the discrete-time algebraic Riccati equation

The authors develop an algorithm to solve the standard discrete-time algebraic Riccati equation by using a skew-Hamiltonian transformation and the square-root method. The algorithm is structure-preserving and efficient because the Hamiltonian structure is fully exploited and only orthogonal transformations are used. The efficiency and stability of the algorithm are analyzed, numerical examples are included.