Design of linear-phase IIR digital filters using singular perturbational model reduction

A new technique for the design of a passband linear-phase infinite impulse response (IIR) digital filter is presented. This technique is based on computing orthonormal bases, which span the left and right invariant subspaces associated with the large and small eigenvalues of the cross-Gramian matrix W_c0, which tends to be symmetric and easily constructed. The orthonormal bases are computed by using the ordered real Schur form decomposition, which has been proved to be robust and numerically stable. The proposed technique is illustrated by several examples