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 rigurous 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. We hope that this robustness property will allow us to extend the results in this paper to other contexts.
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 ℓ 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.
The standard way of solving the polynomial eigenvalue problem associated with a matrix polynomial is to embed the matrix polynomial into a matrix pencil, transforming the problem into an equivalent generalized eigenvalue problem. Such pencils are known as linearizations. Many of the families of linearizations for matrix polynomials available in the literature are extensions of the so-called family of Fiedler pencils. These families are known as generalized Fiedler pencils, Fiedler pencils with repetition and generalized Fiedler pencils with repetition, or Fiedler-like pencils for simplicity. The goal of this work is to unify the Fiedler-like pencils approach with the more recent one based on strong block minimal bases pencils introduced in [15]. To this end, we introduce a family of pencils that we have named extended block Kronecker pencils, whose members are, under some generic nonsingularity conditions, strong block minimal bases pencils, and show that, with the exception of the non proper generalized Fiedler pencils, all Fiedler-like pencils belong to this family modulo permutations. As a consequence of this result, we obtain a much simpler theory for Fiedler-like pencils than the one available so far. Moreover, we expect this unification to allow for further developments in the theory of Fiedler-like pencils such as global or local backward error analyses and eigenvalue conditioning analyses of polynomial eigenvalue problems solved via Fiedler-like linearizations.Key words. Fiedler pencils, generalized Fiedler pencils, Fiedler pencils with repetition, generalized Fiedler pencils with repetition, matrix polynomials, strong linearizations, block minimal bases pencils, block Kronecker pencils, extended block Kronecker pencils, minimal basis, dual minimal bases AMS subject classifications. 65F15, 15A18, 15A22, 15A54 1. Introduction. Matrix polynomials and their associated polynomial eigenvalue problems appear in many areas of applied mathematics, and they have received in the last years considerable attention. For example, they are ubiquitous in a wide range of problems in engineering, mechanic, control theory, computer-aided graphic design, etc. For detailed discussions of different applications of matrix polynomials, we refer the reader to the classical references [23,28,42], the modern surveys [2, Chapter 12] and [37,43] (and their references therein), and the references [32,33,34]. For those readers not familiar with the theory of matrix polynomials and polynomial eigenvalue problems, those topics are briefly reviewed in Section 2.The standard way of solving the polynomial eigenvalue problem associated with a matrix polynomial is to linearize the polynomial into a matrix pencil (i.e., matrix polynomials of grade 1), known as linearization [13,22,23]. The linearization process transforms the polynomial eigenvalue problem into an equivalent generalized eigenvalue problem, which, then, can be solved using mature and well-understood eigensolvers such as the QZ algorithm or the
Abstract. The need to solve polynomial eigenvalue problems for matrix polynomials expressed in nonmonomial bases has become a very important problem. Among the most important bases in numerical applications are the Chebyshev polynomials of the first and second kind. In this work, we introduce a new approach for constructing strong linearizations for matrix polynomials expressed in Chebyshev bases, generalizing the classical colleague pencil, and expanding the arena in which to look for linearizations of matrix polynomials expressed in Chebyshev bases. We show that any of these linearizations is a strong linearization regardless whether the matrix polynomial is regular or singular. In addition, we show how to recover eigenvectors, minimal indices, and minimal bases of the polynomial from those of any of the new linearizations. As an example, we also construct strong linearizations for matrix polynomials of odd degree that are symmetric whenever the matrix polynomials are symmetric.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.