Computing unstructured and structured polynomial pseudospectrum approximations

https://doi.org/10.1016/j.cam.2018.09.033Get rights and content
Under an Elsevier user license
open archive

Abstract

In many applications it is important to understand the sensitivity of eigenvalues of a matrix polynomial to perturbations of the polynomial. The sensitivity commonly is described by condition numbers or pseudospectra. However, the determination of pseudospectra of matrix polynomials is very demanding computationally. This paper describes a new approach to computing approximations of pseudospectra of matrix polynomials by using rank-one or projected rank-one perturbations. These perturbations are inspired by Wilkinson’s analysis of eigenvalue sensitivity. This approach allows the approximation of both structured and unstructured pseudospectra. Computed examples show the method to perform much better than a method based on random rank-one perturbations both for the approximation of structured and unstructured (i.e., standard) polynomial pseudospectra.

Keywords

Matrix polynomials
Pseudospectrum
Structured pseudospectrum
Eigenvalue sensitivity
Distance from defectivity
Numerical methods

Cited by (0)