On the Closest Stable/Unstable Nonnegative Matrix and Related Stability Radii

We consider the problem of computing the closest stable/unstable non-negative matrix to a given real matrix. This problem is important in the study of linear dynamical systems, numerical methods, etc. The distance between matrices is measured in the Frobenius norm. The problem is addressed for two types of stability: the Schur stability (the matrix is stable if its spectral radius is smaller than one) and the Hurwitz stability (the matrix is stable if its spectral abscissa is negative). We show that the closest unstable matrix can always be explicitly found. For the closest stable matrix, we present an iterative algorithm which converges to a local minimum with a linear rate. It is shown that the total number of local minima can be exponential in the dimension. Numerical results and the complexity estimates are presented.

[1]  Punit Sharma,et al.  Stability Radii for Linear Hamiltonian Systems with Dissipation Under Structure-Preserving Perturbations , 2016, SIAM J. Matrix Anal. Appl..

[2]  Tosio Kato.  Perturbation theory for linear operators , 1966 .

[3]  Robert J. Plemmons,et al.  Nonnegative Matrices in the Mathematical Sciences , 1979, Classics in Applied Mathematics.

[4]  F. R. Gantmakher.  The Theory of Matrices , 1984 .

[5]  Dmitriĭ Olegovich Logofet,et al.  Matrices and Graphs Stability Problems in Mathematical Ecology , 1993 .

[6]  Punit Sharma,et al.  On computing the distance to stability for matrices using linear dissipative Hamiltonian systems , 2016, Autom..

[7]  Paul Van Dooren,et al.  Nearest stable system using successive convex approximations , 2013, Autom..

[8]  Nicola Guglielmi,et al.  Approximating real stability radii , 2015 .

[9]  Michael Margaliot,et al.  A Maximum Principle for the Stability Analysis of Positive Bilinear Control Systems with Applications to Positive Linear Switched Systems , 2012, SIAM J. Control. Optim..

[10]  James Anderson,et al.  Distance to the nearest stable Metzler matrix , 2017, 2017 IEEE 56th Annual Conference on Decision and Control (CDC).

[11]  R. Byers.  A Bisection Method for Measuring the Distance of a Stable Matrix to the Unstable Matrices , 1988 .

[12]  R. Fletcher.  Practical Methods of Optimization , 1988 .

[13]  G. Stewart,et al.  Matrix Perturbation Theory , 1990 .

[14]  Daniel Liberzon,et al.  Switching in Systems and Control , 2003, Systems & Control: Foundations & Applications.