Title
On computing the distance to stability for matrices using linear dissipative Hamiltonian systems.
Abstract
In this paper, we consider the problem of computing the nearest stable matrix to an unstable one. We propose new algorithms to solve this problem based on a reformulation using linear dissipative Hamiltonian systems: we show that a matrix A is stable if and only if it can be written as A=(J−R)Q, where J=−JT, R⪰0 and Q≻0 (that is, R is positive semidefinite and Q is positive definite). This reformulation results in an equivalent optimization problem with a simple convex feasible set. We propose three strategies to solve the problem in variables (J,R,Q): (i) a block coordinate descent method, (ii) a projected gradient descent method, and (iii) a fast gradient method inspired from smooth convex optimization. These methods require O(n3) operations per iteration, where n is the size of A. We show the effectiveness of the fast gradient method compared to the other approaches and to several state-of-the-art algorithms.
Year
DOI
Venue
2017
10.1016/j.automatica.2017.07.047
Automatica
Keywords
Field
DocType
Dissipative Hamiltonian systems,Distance to stability,Convex optimization
Gradient method,Discrete mathematics,Gradient descent,Mathematical optimization,Matrix (mathematics),Proximal Gradient Methods,Random coordinate descent,Coordinate descent,Convex optimization,Linear matrix inequality,Mathematics
Journal
Volume
Issue
ISSN
85
1
0005-1098
Citations 
PageRank 
References 
8
0.84
8
Authors
2
Name
Order
Citations
PageRank
Nicolas Gillis150339.77
Punit Sharma2173.17