Abstract:In the last decade matrix polynomials have been investigated with the primary focus on adequate linearizations and good scaling techniques for computing their eigenvalues and eigenvectors. In this article we propose a new method for computing a factored Schur form of the associated companion pencil. The algorithm has a quadratic cost in the degree of the polynomial and a cubic one in the size of the coefficient matrices. Also the eigenvectors can be computed at the same cost. The algorithm is a variant of Francis's implicitly shifted QR algorithm applied on the companion pencil. A preprocessing unitary equivalence is executed on the matrix polynomial to simultaneously bring the leading matrix coefficient and the constant matrix term to triangular form before forming the companion pencil. The resulting structure allows us to stably factor each matrix of the pencil as a product of $k$ matrices of unitary-plus-rank-one form, admitting cheap and numerically reliable storage. The problem is then solved as a product core chasing eigenvalue problem. A backward error analysis is included, implying normwise backward stability after a proper scaling. Computing the eigenvectors via reordering the Schur form is discussed as well. Numerical experiments illustrate stability and efficiency of the proposed methods.

Structured Multi-Matrix Variate, Matrix Polynomial Equations: Solution Techniques

Solving Multi-linear Systems with M-Tensors.

Solving Multi-linear Systems with $$\mathcal {M}$$m-Tensors

Solving Multi-linear Systems with ℳ -Tensors

Parametric Solutions to Sylvester-conjugate Matrix Equations

The Multivariate Eigenvalues of Symmetric Tensors

Bivariate Polynomial Matrix and Smith Form

New Results on Equivalence of Multivariate Polynomial Matrices

Fast and backward stable computation of the eigenvalues and eigenvectors of matrix polynomials

Polynomial-Matrix-Based Approaches

Operational Versus Lie-algebraic Methods and the Theory of Multi-Variable Hermite Polynomials

Tensor Absolute Value Equations

Further results on equivalence of multivariate polynomial matrices

A General Approach to Sylvester-Polynomial-Conjugate Matrix Equations

An Explicit Solution to Polynomial Matrix Right Coprime Factorization with Application in Eigenstructure Assignment

Existence of a unique solution to parametrized systems of generalized polynomial equations

Solutions to a Family of Matrix Equations by Using the Kronecker Matrix Polynomials

Using matrix sparsification to solve tropical linear vector equations

Weyl Tensors, Strongly Regular Graphs, Multiplicative Characters, and a Quadratic Matrix Equation

Lie-algebraic approach to the theory of polynomial solutions. II. Differential equations in one real and one Grassmann variables and 2x2 matrix differential equations

Unified Parametrization for the Solutions to the Polynomial Diophantine Matrix Equation and the Generalized Sylvester Matrix Equation