In this talk, we will show the existence of so-called the real cyclic symplectic Schur decomposition for a cyclic sequence of symplectic matrix pairs. An efficient and robust algorithm for computing this cyclic Schur decomposition, and consequently the stable invariant subspace, is proposed and analyzed. We apply this algorithm to the solution of periodic Riccati equations arising in periodic discrete time optimal control problems.