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 is stable if and only if it can be written as , where , and (that is, is positive semidefinite and 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 : (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 operations per iteration, where is the size of . We show the effectiveness of the fast gradient method compared to the other approaches and to several state-of-the-art algorithms.
Cite
@article{arxiv.1611.00595,
title = {On computing the distance to stability for matrices using linear dissipative Hamiltonian systems},
author = {Nicolas Gillis and Punit Sharma},
journal= {arXiv preprint arXiv:1611.00595},
year = {2017}
}
Comments
21 pages, 5 figures, 4 tables. Some typos and errors fixed, new remark on uniqueness issues