Equivalences for Linearizations of Matrix Polynomials
Abstract
One useful standard method to compute eigenvalues of matrix polynomials of degree at most in (denoted of grade , for short) is to first transform to an equivalent linear matrix polynomial , called a companion pencil, where and are usually of larger dimension than but is now only of grade in . The eigenvalues and eigenvectors of can be computed numerically by, for instance, the QZ algorithm. The eigenvectors of , including those for infinite eigenvalues, can also be recovered from eigenvectors of if is what is called a "strong linearization" of . In this paper we show how to use algorithms for computing the Hermite Normal Form of a companion matrix for a scalar polynomial to direct the discovery of unimodular matrix polynomial cofactors and which, via the equation , explicitly show the equivalence of and . By this method we give new explicit constructions for several linearizations using different polynomial bases. We contrast these new unimodular pairs with those constructed by strict equivalence, some of which are also new to this paper. We discuss the limitations of this experimental, computational discovery method of finding unimodular cofactors.
Cite
@article{arxiv.2102.09726,
title = {Equivalences for Linearizations of Matrix Polynomials},
author = {Robert M. Corless and Leili Rafiee Sevyeri and B. David Saunders},
journal= {arXiv preprint arXiv:2102.09726},
year = {2021}
}