Block Kronecker linearizations of matrix polynomials and their backward errors

• Published in: Numerische Mathematik 2018-10, Vol.140 (2), p.373-426
• Main Authors: Dopico, Froilán M., Lawrence, Piers W., Pérez, Javier, Dooren, Paul Van
• Format: Article
• Language: English
• Subjects: Algorithms; Error analysis; Mathematical and Computational Engineering; Mathematical and Computational Physics; Mathematical Methods in Physics; Mathematics; Mathematics and Statistics; Numerical Analysis; Numerical and Computational Physics; Pencils; Polynomials; Simulation; Stability analysis; Theoretical
• ISSN: 0029-599X; 0945-3245
• DOI: 10.1007/s00211-018-0969-z

We introduce a new family of strong linearizations of matrix polynomials—which we call “block Kronecker pencils”—and perform a backward stability analysis of complete polynomial eigenproblems. These problems are solved by applying any backward stable algorithm to a block Kronecker pencil, such as the staircase algorithm for singular pencils or the QZ algorithm for regular pencils. This stability analysis allows us to identify those block Kronecker pencils that yield a computed complete eigenstructure which is exactly that of a slightly perturbed matrix polynomial. The global backward error analysis in this work presents for the first time the following key properties: it is a rigorous analysis valid for finite perturbations (i.e., it is not a first order analysis), it provides precise bounds, it is valid simultaneously for a large class of linearizations, and it establishes a framework that may be generalized to other classes of linearizations. These features are related to the fact that block Kronecker pencils are a particular case of the new family of “strong block minimal bases pencils”, which are robust under certain perturbations and, so, include certain perturbations of block Kronecker pencils.