A numerical evaluation of solvers for the periodic Riccati differential equation

Abstract

Efficient and accurate structure exploiting numerical methods for solving the periodic Riccati differential equation (PRDE) are addressed. Such methods are essential, for example, to design periodic feedback controllers for periodic control systems. Three recently proposed methods for solving the PRDE are presented and evaluated on challenging periodic linear artificial systems with known solutions and applied to the stabilization of periodic motions of mechanical systems. The first two methods are of the type multiple shooting and rely on computing the stable invariant subspace on associated Hamiltonian system. The stable subspace is determined using either algorithms for computing an ordered periodic real Schur form of a cyclic matrix sequence, or a recently proposed method which implicitly constructs a stable deflating subspace from an associated lifted pencil. The third method reformulates the PRDE as a convex optimization problem where the stabilizing solution is approximated by its truncated Fourier series. As known, this reformulations leads to a semi-definite programming problem with linear matrix inequality constraints admitting an effective nuermical realization. The numerical evaluation of the PRDE methods, with focus on the number of states (n) and the length of the period (T) of the periodic systems considered, includes both quantitative and qualitative results.