Structural backward stability in rational eigenvalue problems solved via block Kronecker linearizations
Abstract
We study the backward stability of running a backward stable eigenstructure solver on a pencil that is a strong linearization of a rational matrix expressed in the form , where is a polynomial matrix and is a minimal state-space realization. We consider the family of block Kronecker linearizations of , which are highly structured pencils. Backward stable eigenstructure solvers applied to will compute the exact eigenstructure of a perturbed pencil and the special structure of will be lost. In order to link this perturbed pencil with a nearby rational matrix, we construct a strictly equivalent pencil to that restores the original structure, and hence is a block Kronecker linearization of a perturbed rational matrix , where is a polynomial matrix with the same degree as . Moreover, we bound appropriate norms of , , and in terms of an appropriate norm of . These bounds may be inadmissibly large, but we also introduce a scaling that allows us to make them satisfactorily tiny. Thus, for this scaled representation, we prove that the staircase and the algorithms compute the exact eigenstructure of a rational matrix that can be expressed in exactly the same form as with the parameters defining the representation very near to those of . This shows that this approach is backward stable in a structured sense.
Cite
@article{arxiv.2103.16395,
title = {Structural backward stability in rational eigenvalue problems solved via block Kronecker linearizations},
author = {Froilán M. Dopico and María C. Quintana and Paul Van Dooren},
journal= {arXiv preprint arXiv:2103.16395},
year = {2021}
}