We propose a method for computing the distance of a stable polynomial to the set of unstable ones (both in the Hurwitz and in the Schur case). The method is based on the reformulation of the problem as the structured distance to instability of a companion matrix associated to a polynomial. We first introduce the structured ε-pseudospectrum of a companion matrix and write a system of ordinary differential equations which maximize the real part (or the absolute value) of elements of the structured ε-pseudospectrum and then exploit the knowledge of the derivative of the maximizers with respect to ε to devise a quadratically convergent iteration. Furthermore we use a variant of the same ODEs to compute the boundary of structured pseudospectra and compare them to unstructured ones. An extension to constrained perturbations is also considered.

An iterative method for computing robustness of polynomial stability

GUGLIELMI, NICOLA;
2016-01-01

Abstract

We propose a method for computing the distance of a stable polynomial to the set of unstable ones (both in the Hurwitz and in the Schur case). The method is based on the reformulation of the problem as the structured distance to instability of a companion matrix associated to a polynomial. We first introduce the structured ε-pseudospectrum of a companion matrix and write a system of ordinary differential equations which maximize the real part (or the absolute value) of elements of the structured ε-pseudospectrum and then exploit the knowledge of the derivative of the maximizers with respect to ε to devise a quadratically convergent iteration. Furthermore we use a variant of the same ODEs to compute the boundary of structured pseudospectra and compare them to unstructured ones. An extension to constrained perturbations is also considered.
2016
Companion matrices; Perturbations of polynomials; Polynomial stability; Pseudozero set; Stability radii; Structured pseudospectra; Computational Mathematics; Applied Mathematics
File in questo prodotto:
File Dimensione Formato  
2016_JComputApplMath_292_638_Guglielmi.pdf

non disponibili

Tipologia: Versione Editoriale (PDF)
Licenza: Non pubblico
Dimensione 620.18 kB
Formato Adobe PDF
620.18 kB Adobe PDF   Visualizza/Apri   Richiedi una copia

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/20.500.12571/914
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 1
  • ???jsp.display-item.citation.isi??? 1
social impact