When studying the ε-pseudospectrum of a matrix, one is often interested in computing the extremal points having maximum real part or modulus. This is a crucial step, for example,when computing the distance to instability of a stable system. Using the property that the pseudospectrum is determined via perturbations by rank-1 matrices, we derive differential equations onthe manifold of normalized rank-1 matrices whose solutions tend to the critical rank-1 perturbationsassociated with the extremal points of (locally) maximum real part and modulus. This approach alsoallows us to track the boundary contour of the pseudospectrum in a neighborhood of the extremalpoints. The technique we propose is related to an idea recently developed by Guglielmi and Overton,who derived discrete dynamical systems instead of the continuous ones we present. The methodturns out to be fast in comparison with those previously proposed in the literature and appears tobe promising in dealing with large sparse problems.

Differential equations for roaming pseudospectra: paths to extremal points and boundary tracking.

GUGLIELMI, NICOLA;
2011

Abstract

When studying the ε-pseudospectrum of a matrix, one is often interested in computing the extremal points having maximum real part or modulus. This is a crucial step, for example,when computing the distance to instability of a stable system. Using the property that the pseudospectrum is determined via perturbations by rank-1 matrices, we derive differential equations onthe manifold of normalized rank-1 matrices whose solutions tend to the critical rank-1 perturbationsassociated with the extremal points of (locally) maximum real part and modulus. This approach alsoallows us to track the boundary contour of the pseudospectrum in a neighborhood of the extremalpoints. The technique we propose is related to an idea recently developed by Guglielmi and Overton,who derived discrete dynamical systems instead of the continuous ones we present. The methodturns out to be fast in comparison with those previously proposed in the literature and appears tobe promising in dealing with large sparse problems.
File in questo prodotto:
File Dimensione Formato  
2011_SIAMJNUMERANAL_49_Guglielmi_Lubich.pdf

non disponibili

Tipologia: Versione Editoriale (PDF)
Licenza: Non pubblico
Dimensione 620.68 kB
Formato Adobe PDF
620.68 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: http://hdl.handle.net/20.500.12571/1552
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 31
  • ???jsp.display-item.citation.isi??? 28
social impact