Abstract:SIAM Journal on Matrix Analysis and Applications, Volume 45, Issue 4, Page 2104-2134, December 2024. This work concerns the minimization of the pseudospectral abscissa of a matrix-valued function dependent on parameters analytically. The problem is motivated by robust stability and transient behavior considerations for a linear control system that has optimization parameters. We describe a subspace procedure to cope with the setting when the matrix-valued function is of large size. The proposed subspace procedure solves a sequence of reduced problems obtained by restricting the matrix-valued function to small subspaces, whose dimensions increase gradually. It possesses desirable features such as a superlinear convergence exhibited by the decay in the errors of the minimizers of the reduced problems. In mathematical terms, the problem we consider is a large-scale nonconvex minimax eigenvalue optimization problem such that the eigenvalue function appears in the constraint of the inner maximization problem. Devising and analyzing a subspace framework for the minimax eigenvalue optimization problem at hand with the eigenvalue function in the constraint require special treatment that makes use of a Lagrangian and dual variables. There are notable advantages in minimizing the pseudospectral abscissa over maximizing the distance to instability or minimizing the [math] norm; the optimized pseudospectral abscissa provides quantitative information about the worst-case transient growth, and the initial guesses for the parameter values to optimize the pseudospectral abscissa can be arbitrary, unlike the case to optimize the distance to instability and [math] norm that would normally require initial guesses yielding asymptotically stable systems. Reproducibility of computational results. This paper has been awarded the "SIAM Reproducibility Badge: Code and data available" as a recognition that the authors have followed reproducibility principles valued by SIMAX and the scientific computing community. Code and data that allow readers to reproduce the results in this paper are available at https://zenodo.org/records/6992092.

Large-Scale Minimization of the Pseudospectral Abscissa

Large-Scale Minimization of the Pseudospectral Abscissa

Minimization of the Pseudospectral Abscissa of a Quadratic Matrix Polynomial

An Effective Pseudospectral Optimization Approach with Sparse Variable Time Nodes for Maximum Production of Chemical Engineering Problems

Lp quasi-norm minimization: algorithm and applications

Scalable Algorithms for Tractable Schatten Quasi-Norm Minimization

An Alternating Minimization Method for Robust Principal Component Analysis

Uniform Approximation of Eigenproblems of a Large-Scale Parameter-Dependent Hermitian Matrix

A continuous approach for computing the pseudospectra of linear operators

A new constrained optimization model for solving the nonsymmetric stochastic inverse eigenvalue problem

Alternating Direction and Taylor Expansion Minimization Algorithms for Unconstrained Nuclear Norm Optimization.

Low-complexity subspace-descent over symmetric positive definite manifold

l1-norm regularized l1-norm best-fit lines

Large-Scale Computation of ${\mathcal L}_\infty$-Norms by a Greedy Subspace Method

The Minimization of Piecewise Functions: Pseudo Stationarity

The Projection Methods for Computing the Pseudospectra of Large Scale Matrices

A Gradient Descent Based Algorithm for ℓ Minimization.

Solving Minimal Residual Methods in with Large Exponents p

A Hermite‐Lobatto Pseudospectral Method for Optimal Control

An efficient active-set method with applications to sparse approximations and risk minimization

Approximation Accuracy of the Krylov Subspaces for Linear Discrete Ill-Posed Problems