A new look at partial fraction expansion of transfer function matrices from a computational viewpoint (Q2367556)

A computational procedure is given for obtaining the matrix residues of the partial fraction expansion (PFE) of the transfer function matrix \(H(s) = C(sI_ n - A)^{-1}B\) for the dynamical system described by \(dx(t)/dt = Ax(t) + Bu(t)\), \(y(t) = Cx(t)\). No knowledge is necessary of the minimal polynomial or the particular structure of the Jordan canonical form of the system matrix \(A\), which is a desirable advantage of the procedure. It gives all matrix residues at the same time and no derivatives are required for repeated eigenvalues. The procedure is particularly suitable for digital computer implementation with MATLAB. Numerical examples are also given to illustrate the performance of the procedure.