Inversion components of block Hankel-like matrices (Q1208260)

The inversion problem for square matrices having the structure of a block Hankel-like (i.e. Hankel-striped, Hankel-layered, vector Hankel matrices) is studied. The author shows that the components that both determine nonsingularity and construct the inverse of such matrices are closely related to multidimensional generalizations of Padé-Hermite and simultaneous Padé approximants. Matrix Padé-Hermite and matrix simultaneous Padé systems are shown to provide a second set of inverse components for block Hankel-like matrices. A recurrence relation is presented that allows for efficient computation of matrix Padé-Hermite and matrix simultaneous Padé systems. Thus the inverse components can be computed via either the matrix Euclidean algorithm or a matrix Berlekamp-Massey algorithm applied to an associated matrix power series. An alternative algorithm based on this recurrence relation is also presented. This algorithm has the advantage that no extra conditions are required on the input matrix.