Computing a compact local Smith–McMillan form

Abstract

We define a compact local Smith–McMillan form of a rational matrix R ( λ ) as the diagonal matrix whose diagonal elements are the nonzero entries of a local Smith-McMillan form of R ( λ ) . We show that a recursive rank search procedure, applied to a block-Toeplitz matrix built on the Laurent expansion of R ( λ ) around an arbitrary complex point λ 0 , allows us to compute a compact local Smith-McMillan form of that rational matrix R ( λ ) at the point λ 0 , provided we keep track of the transformation matrices used in the rank search. It also allows us to recover the root polynomials of a polynomial matrix and root vectors of a rational matrix, at an expansion point λ 0 . Numerical tests illustrate the promising performance of the resulting algorithm.