Computation of transfer function matrices of periodic systems

We present a numerical approach to evaluate the transfer function matrices of a periodic system corresponding to lifted state-space representations as constant systems. The proposed pole-zero method determines each entry of the transfer function matrix in a minimal zeros-poles-gain representation. A basic computational ingredient for this method is the extended periodic real Schur form of a periodic matrix, which underlies the computation of minimal realizations and system poles. To compute zeros and gains, fast algorithms are proposed, which are specially tailored to particular single-input single-output periodic systems. The new method relies exclusively on reliable numerical computations and is well suited for robust software implementations.