[Distinguished Lecture] Computing Closest Stable Non-negative Matrices
Speaker(s): Prof. Yurii Nesterov (CORE, UCL Belgium)
Time: 15:00-16:00 September 6, 2017
Venue: Lecture Hall, Jiayibing Building, Jingchunyuan 82, BICMR
Yurii Nesterov is a Russian mathematician, an internationally recognized expert in convex optimization, especially in the development of efficient algorithms and numerical optimization analysis. He is currently a professor at the Université catholique de Louvain (UCL). Nesterov is most famous for his work in convex optimization, including his 2004 book, considered a canonical reference on the subject. His main novel contribution is an accelerated version of gradient descent that converges considerably faster than ordinary gradient descent. His work with Arkadi Nemirovski in the 1994 book is the first to point out that interior point method can solve convex optimization problems, and the first to make a systematic study of semidefinite programming (SDP). Also in this book, they introduced the self-concordant functions which are useful in the analysis of Newton's method.
Abstract:
Problem of finding the closest stable matrix for a dynamical system has many applications. It is well studied both for continuous and discrete-time systems, and the corresponding optimization problems are formulated for various matrix norms. As a rule, non-convexity of these formulations does not allow finding their global solutions. In this paper, we analyze positive discrete-time systems. They also suffer from non-convexity of the stability region, and the problem in the Frobenius norm or in Euclidean norm remains hard for them. However, it turns out that for certain polyhedral norms, the situation is much better. We show, that for the distances measured in the max-norm, we can find exact solution of the corresponding nonconvex projection problems in polynomial time. For the distance measured in the operator $\ell_{\infty}$-norm or $\ell_{1}$-norm, the exact solution is also efficiently found. To this end, we develop a modification of the recently introduced spectral simplex method. On the other hand, for all these three norms, we obtain exact descriptions of the region of stability around a given stable matrix. In the case of max-norm, this can be seen as an analogue of Kharitonov's theorem for non-negative matrices.
This is a joint work with V. Protasov (Moscow State University).