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 Computation of ${\mathcal L}_\infty$-Norms by a Greedy Subspace Method

A Subspace Framework for ${\mathcal L}_\infty$ Model Reduction

${L}^{\infty}$-norm computation for linear time-invariant systems depending on parameters

Semismooth Newton Algorithm for Efficient Projections Onto ℓ1, ∞-Norm Ball

Smooth Approximation of L_infinity-Norm for Multi-view Geometry

Norm attaining operators which satisfy a Bollobás type theorem

Semismooth Newton Algorithm for Efficient Projections onto $\ell_1, ∞$-norm Ball

Low-Rank Approximation via Generalized Reweighted Iterative Nuclear and Frobenius Norms

l1-norm regularized l1-norm best-fit lines

Large-Scale Minimization of the Pseudospectral Abscissa

A Joint Bidiagonalization Based Algorithm for Large Scale Linear Discrete Ill-posed Problems in General-Form Regularization

Projection constants for spaces of Dirichlet polynomials

Norm-induced Cuts: Optimization with Lipschitzian Black-box Functions

LatLRR for subspace clustering via reweighted Frobenius norm minimization

Sum-of-squares hierarchies for polynomial optimization and the Christoffel-Darboux kernel

Dominant subspace and low-rank approximations from block Krylov subspaces without a prescribed gap

Multivariate approximation of functions on irregular domains by weighted least-squares methods

Non-Euclidean Monotone Operator Theory and Applications

Lp quasi-norm minimization: algorithm and applications

Efficient Proximal Mapping Computation for Unitarily Invariant Low-Rank Inducing Norms