On the order reduction of MIMO large-scale systems using block-roots of matrix polynomials