English

Structural backward stability in rational eigenvalue problems solved via block Kronecker linearizations

Numerical Analysis 2021-03-31 v1 Numerical Analysis Spectral Theory

Abstract

We study the backward stability of running a backward stable eigenstructure solver on a pencil S(λ)S(\lambda) that is a strong linearization of a rational matrix R(λ)R(\lambda) expressed in the form R(λ)=D(λ)+C(λIA)1BR(\lambda)=D(\lambda)+ C(\lambda I_\ell-A)^{-1}B, where D(λ)D(\lambda) is a polynomial matrix and C(λIA)1BC(\lambda I_\ell-A)^{-1}B is a minimal state-space realization. We consider the family of block Kronecker linearizations of R(λ)R(\lambda), which are highly structured pencils. Backward stable eigenstructure solvers applied to S(λ)S(\lambda) will compute the exact eigenstructure of a perturbed pencil S^(λ):=S(λ)+ΔS(λ)\widehat S(\lambda):=S(\lambda)+\Delta_S(\lambda) and the special structure of S(λ)S(\lambda) will be lost. In order to link this perturbed pencil with a nearby rational matrix, we construct a strictly equivalent pencil S~(λ)\widetilde S(\lambda) to S^(λ)\widehat S(\lambda) that restores the original structure, and hence is a block Kronecker linearization of a perturbed rational matrix R~(λ)=D~(λ)+C~(λIA~)1B~\widetilde R(\lambda) = \widetilde D(\lambda)+ \widetilde C(\lambda I_\ell- \widetilde A)^{-1} \widetilde B, where D~(λ)\widetilde D(\lambda) is a polynomial matrix with the same degree as D(λ)D(\lambda). Moreover, we bound appropriate norms of D~(λ)D(λ)\widetilde D(\lambda)- D(\lambda), C~C\widetilde C - C, A~A\widetilde A - A and B~B\widetilde B - B in terms of an appropriate norm of ΔS(λ)\Delta_S(\lambda). 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 QZQZ algorithms compute the exact eigenstructure of a rational matrix R~(λ)\widetilde R(\lambda) that can be expressed in exactly the same form as R(λ)R(\lambda) with the parameters defining the representation very near to those of R(λ)R(\lambda). 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}
}
R2 v1 2026-06-24T00:41:42.630Z