Algorithm for identification of multivariable discrete-time linear systems

The method presented is for the z-transfer function matrix G(z) identification for constant-coefficient linear systems from samples of the impulse response matrix. By the proposed approach the roots of the minimal polynomial of G(z) and the corresponding residues are sequentially determined. The computational effort is here greatly reduced.