Block minimal bases ℓ-ifications of matrix polynomials

Abstract

The standard way of solving a polynomial eigenvalue problem associated with a matrix polynomial starts by embedding the matrix coefficients of the polynomial into a matrix pencil, known as a strong linearization. This process transforms the problem into an equivalent generalized eigenvalue problem. However, there are some situations in which is more convenient to replace linearizations by other low degree matrix polynomials. This has motivated the idea of a strong ℓ-ification of a matrix polynomial, which is a matrix polynomial of degree at most ℓ having the same finite and infinite elementary divisors, and the same numbers of left and right minimal indices as the original matrix polynomial. We present in this work a novel method for constructing strong ℓ-ifications of matrix polynomials of size m×n and grade d when ℓ<d, and ℓ divides nd or md. This method is based on a family called “strong block minimal bases matrix polynomials”, and relies heavily on properties of dual minimal bases. We show how strong block minimal bases ℓ-ifications can be constructed from the coefficients of a given matrix polynomial P(λ). We also show that these ℓ-ifications satisfy many desirable properties for numerical applications: they are strong ℓ-ifications regardless of whether P(λ) is regular or singular, the minimal indices of the ℓ-ifications are related to those of P(λ) via constant uniform shifts, and eigenvectors and minimal bases of P(λ) can be recovered from those of any of the strong block minimal bases ℓ-ifications. In the special case where ℓ divides d, we introduce a subfamily of strong block minimal bases matrix polynomials named “block Kronecker matrix polynomials”, which is shown to be a fruitful source of companion ℓ-ifications.