The H∞ norm of a transfer matrix function for a control system is the reciprocal ofthe largest value of ε such that the associated ε-spectral value set is contained in the stability regionfor the dynamical system (the left half-plane in the continuous-time case and the unit disk in thediscrete-time case). After deriving some fundamental properties of spectral value sets, particularlythe intricate relationship between the singular vectors of the transfer matrix and the eigenvectors ofthe corresponding perturbed system matrix, we extend an algorithm recently introduced by Guglielmiand Overton [SIAM J. Matrix Anal. Appl., 32 (2011), pp. 1166–1192] for approximating the maximalreal part or modulus of points in a matrix pseudospectrum to spectral value sets, characterizing itsfixed points. We then introduce a Newton-bisection method to approximate the H∞ norm, forwhich each step requires optimization of the real part or the modulus over an ε-spectral value set.Although the algorithm is guaranteed only to find lower bounds on the H∞ norm, it typically findsgood approximations in cases where we can test this. It is much faster than the standard Boyd–Balakrishnan–Bruinsma–Steinbuch algorithm to compute the H∞ norm when the system matricesare large and sparse and the number of inputs and outputs is small. The main work required by thealgorithm is the computation of the spectral abscissa or radius of a sequence of matrices that arerank-one perturbations of a sparse matrix.

Fast approximation of the H∞ norm via optimization over spectral value sets.

GUGLIELMI N;
2013

Abstract

The H∞ norm of a transfer matrix function for a control system is the reciprocal ofthe largest value of ε such that the associated ε-spectral value set is contained in the stability regionfor the dynamical system (the left half-plane in the continuous-time case and the unit disk in thediscrete-time case). After deriving some fundamental properties of spectral value sets, particularlythe intricate relationship between the singular vectors of the transfer matrix and the eigenvectors ofthe corresponding perturbed system matrix, we extend an algorithm recently introduced by Guglielmiand Overton [SIAM J. Matrix Anal. Appl., 32 (2011), pp. 1166–1192] for approximating the maximalreal part or modulus of points in a matrix pseudospectrum to spectral value sets, characterizing itsfixed points. We then introduce a Newton-bisection method to approximate the H∞ norm, forwhich each step requires optimization of the real part or the modulus over an ε-spectral value set.Although the algorithm is guaranteed only to find lower bounds on the H∞ norm, it typically findsgood approximations in cases where we can test this. It is much faster than the standard Boyd–Balakrishnan–Bruinsma–Steinbuch algorithm to compute the H∞ norm when the system matricesare large and sparse and the number of inputs and outputs is small. The main work required by thealgorithm is the computation of the spectral abscissa or radius of a sequence of matrices that arerank-one perturbations of a sparse matrix.
File in questo prodotto:
File Dimensione Formato  
Postprint_2013_SIAMJMatrixAnalAppl_34_Guglielmi.pdf

non disponibili

Tipologia: Documento in Post-print
Licenza: Non pubblico
Dimensione 386.74 kB
Formato Adobe PDF
386.74 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/1737
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 35
  • ???jsp.display-item.citation.isi??? 29
social impact