Fast Approximation of the H∞ Norm via Optimization over Spectral Value Sets

The H∞ norm of a transfer matrix of a control system is the reciprocal of the largest value of ε such that the associated ε-spectral value set is contained in the stability region for the dynamical system (the left half-plane in the continuous-time case and the unit disk in the discrete-time case). We extend an algorithm recently introduced by Guglielmi and Overton [GO11] for approximating the maximal real part or modulus of points in a matrix pseudospectrum to spectral value sets, characterizing its fixed points. We then introduce a Newton-bisection method to approximate the H∞ norm, for which each step requires optimization of the real part or the modulus over an ε-spectral value set. The algorithm is much faster than the standard Boyd-Balakrishnan-BruinsmaSteinbuch algorithm to compute the H∞ norm when the system matrices are large and sparse and the number of inputs and outputs is small. The main work required by the algorithm is the computation of the spectral abscissa or radius of a sequence of matrices that are rank-one perturbations of a sparse matrix.