Optimization and pseudospectra, with applications to robust stability

J. V. Burke, A. S. Lewis, M. L. Overton

Research output: Contribution to journalArticle

Abstract

The ε-pseudospectrum of a matrix A is the subset of the complex plane consisting of all eigenvalues of all complex matrices within a distance ε of A, We are interested in two aspects of "optimization and pseudoapectra." The first concerns maximizing the function "real part" over an ε-pseudospectrum of a fixed matrix: this defines a function known as the ε-pseudospectral abscissa of a matrix. We present a bisection algorithm to compute this function. Our second interest is in minimizing the ε-pseudospectral abscissa over a set of feasible matrices. A prerequisite for local optimization of this function is an understanding of its variational properties, the study of which is the main focus of the paper. We show that, in a neighborhood of any nonderogatory matrix, the ε-pseudospectral abscissa is a nonsmooth but locally Lipschitz and subdifferentially regular function for sufficiently small ε; in fact, it can be expressed locally as the maximum of a finite number of smooth functions. Along the way we obtain an eigenvalue perturbation result: near a nonderogatory matrix, the eigenvalues satisfy a Hölder continuity property on matrix space - a property that is well known when only a single perturbation parameter is considered. The pseudospectral abscissa is a powerful modeling tool: not only is it a robust measure of stability, but it also reveals the transient (as opposed to asymptotic) behavior of associated dynamical systems.

Original languageEnglish (US)
Pages (from-to)80-104
Number of pages25
JournalSIAM Journal on Matrix Analysis and Applications
Volume25
Issue number1
DOIs
StatePublished - 2004

Fingerprint

Pseudospectra
Robust Stability
Abscissa
Optimization
Eigenvalue
Local Optimization
Bisection
Robust stability
Parameter Perturbation
Smooth function
Argand diagram
Lipschitz
Dynamical systems
Dynamical system
Asymptotic Behavior
Perturbation
Subset

Keywords

  • Distance to instability
  • Eigenvalue optimization
  • H norm
  • Nonsmooth analysis
  • Pseudospectrum
  • Robust control
  • Robust optimization
  • Spectral abscissa
  • Stability radius
  • Subdifferential regularity

ASJC Scopus subject areas

  • Algebra and Number Theory
  • Applied Mathematics
  • Analysis

Cite this

Optimization and pseudospectra, with applications to robust stability. / Burke, J. V.; Lewis, A. S.; Overton, M. L.

In: SIAM Journal on Matrix Analysis and Applications, Vol. 25, No. 1, 2004, p. 80-104.

Research output: Contribution to journalArticle

@article{0d7e0ddfa12f487d878fcce1fca0b1ff,
title = "Optimization and pseudospectra, with applications to robust stability",
abstract = "The ε-pseudospectrum of a matrix A is the subset of the complex plane consisting of all eigenvalues of all complex matrices within a distance ε of A, We are interested in two aspects of {"}optimization and pseudoapectra.{"} The first concerns maximizing the function {"}real part{"} over an ε-pseudospectrum of a fixed matrix: this defines a function known as the ε-pseudospectral abscissa of a matrix. We present a bisection algorithm to compute this function. Our second interest is in minimizing the ε-pseudospectral abscissa over a set of feasible matrices. A prerequisite for local optimization of this function is an understanding of its variational properties, the study of which is the main focus of the paper. We show that, in a neighborhood of any nonderogatory matrix, the ε-pseudospectral abscissa is a nonsmooth but locally Lipschitz and subdifferentially regular function for sufficiently small ε; in fact, it can be expressed locally as the maximum of a finite number of smooth functions. Along the way we obtain an eigenvalue perturbation result: near a nonderogatory matrix, the eigenvalues satisfy a H{\"o}lder continuity property on matrix space - a property that is well known when only a single perturbation parameter is considered. The pseudospectral abscissa is a powerful modeling tool: not only is it a robust measure of stability, but it also reveals the transient (as opposed to asymptotic) behavior of associated dynamical systems.",
keywords = "Distance to instability, Eigenvalue optimization, H norm, Nonsmooth analysis, Pseudospectrum, Robust control, Robust optimization, Spectral abscissa, Stability radius, Subdifferential regularity",
author = "Burke, {J. V.} and Lewis, {A. S.} and Overton, {M. L.}",
year = "2004",
doi = "10.1137/S0895479802402818",
language = "English (US)",
volume = "25",
pages = "80--104",
journal = "SIAM Journal on Matrix Analysis and Applications",
issn = "0895-4798",
publisher = "Society for Industrial and Applied Mathematics Publications",
number = "1",

}

TY - JOUR

T1 - Optimization and pseudospectra, with applications to robust stability

AU - Burke, J. V.

AU - Lewis, A. S.

AU - Overton, M. L.

PY - 2004

Y1 - 2004

N2 - The ε-pseudospectrum of a matrix A is the subset of the complex plane consisting of all eigenvalues of all complex matrices within a distance ε of A, We are interested in two aspects of "optimization and pseudoapectra." The first concerns maximizing the function "real part" over an ε-pseudospectrum of a fixed matrix: this defines a function known as the ε-pseudospectral abscissa of a matrix. We present a bisection algorithm to compute this function. Our second interest is in minimizing the ε-pseudospectral abscissa over a set of feasible matrices. A prerequisite for local optimization of this function is an understanding of its variational properties, the study of which is the main focus of the paper. We show that, in a neighborhood of any nonderogatory matrix, the ε-pseudospectral abscissa is a nonsmooth but locally Lipschitz and subdifferentially regular function for sufficiently small ε; in fact, it can be expressed locally as the maximum of a finite number of smooth functions. Along the way we obtain an eigenvalue perturbation result: near a nonderogatory matrix, the eigenvalues satisfy a Hölder continuity property on matrix space - a property that is well known when only a single perturbation parameter is considered. The pseudospectral abscissa is a powerful modeling tool: not only is it a robust measure of stability, but it also reveals the transient (as opposed to asymptotic) behavior of associated dynamical systems.

AB - The ε-pseudospectrum of a matrix A is the subset of the complex plane consisting of all eigenvalues of all complex matrices within a distance ε of A, We are interested in two aspects of "optimization and pseudoapectra." The first concerns maximizing the function "real part" over an ε-pseudospectrum of a fixed matrix: this defines a function known as the ε-pseudospectral abscissa of a matrix. We present a bisection algorithm to compute this function. Our second interest is in minimizing the ε-pseudospectral abscissa over a set of feasible matrices. A prerequisite for local optimization of this function is an understanding of its variational properties, the study of which is the main focus of the paper. We show that, in a neighborhood of any nonderogatory matrix, the ε-pseudospectral abscissa is a nonsmooth but locally Lipschitz and subdifferentially regular function for sufficiently small ε; in fact, it can be expressed locally as the maximum of a finite number of smooth functions. Along the way we obtain an eigenvalue perturbation result: near a nonderogatory matrix, the eigenvalues satisfy a Hölder continuity property on matrix space - a property that is well known when only a single perturbation parameter is considered. The pseudospectral abscissa is a powerful modeling tool: not only is it a robust measure of stability, but it also reveals the transient (as opposed to asymptotic) behavior of associated dynamical systems.

KW - Distance to instability

KW - Eigenvalue optimization

KW - H norm

KW - Nonsmooth analysis

KW - Pseudospectrum

KW - Robust control

KW - Robust optimization

KW - Spectral abscissa

KW - Stability radius

KW - Subdifferential regularity

UR - http://www.scopus.com/inward/record.url?scp=1342310021&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=1342310021&partnerID=8YFLogxK

U2 - 10.1137/S0895479802402818

DO - 10.1137/S0895479802402818

M3 - Article

VL - 25

SP - 80

EP - 104

JO - SIAM Journal on Matrix Analysis and Applications

JF - SIAM Journal on Matrix Analysis and Applications

SN - 0895-4798

IS - 1

ER -