Structure-preserving algorithms for periodic discrete-time algebraic Riccati equations

In this paper we investigate structure-preserving algorithms for computing the symmetric positive semi-definite solutions to the periodic discrete-time algebraic Riccati equations (P-DAREs). Using a structure-preserving swap and collapse procedure, a single symplectic matrix pair in standard symplectic form is obtained. The P-DAREs can then be solved via a single DARE, using a structure-preserving doubling algorithm. We develop the structure-preserving doubling algorithm from a new point of view and show its quadratic convergence under assumptions which are weaker than stabilizability and detectability. With several numerical results, the algorithm is shown to be efficient, out-performing other algorithms on a large set of benchmark problems.