Nearest Ω-stable matrix via Riemannian optimization

02/17/2020
by   Vanni Noferini, et al.
0

We study the problem of finding the nearest Ω-stable matrix to a certain matrix A, i.e., the nearest matrix with all its eigenvalues in a prescribed closed set Ω. Distances are measured in the Frobenius norm. An important special case is finding the nearest Hurwitz or Schur stable matrix, which has applications in systems theory. We describe a reformulation of the task as an optimization problem on the Riemannian manifold of orthogonal (or unitary) matrices. The problem can then be solved using standard methods from the theory of Riemannian optimization. The resulting algorithm is remarkably fast on small-scale and medium-scale matrices, and returns directly a Schur factorization of the minimizer, sidestepping the numerical difficulties associated with eigenvalues with high multiplicity.

READ FULL TEXT
POST COMMENT

Comments

There are no comments yet.

Authors

page 10

08/27/2021

The real symplectic Stiefel and Grassmann manifolds: metrics, geodesics and applications

The real symplectic Stiefel manifold is the manifold of symplectic bases...
09/20/2018

On constructing orthogonal generalized doubly stochastic matrices

A real quadratic matrix is generalized doubly stochastic (g.d.s.) if all...
03/02/2021

A numerical method for computing the Jordan Canonical Form

The Jordan Canonical Form of a matrix is highly sensitive to perturbatio...
09/30/2021

On Riemannian Approach for Constrained Optimization Model in Extreme Classification Problems

We propose a novel Riemannian method for solving the Extreme multi-label...
09/24/2020

Sampling the eigenvalues of random orthogonal and unitary matrices

We develop an efficient algorithm for sampling the eigenvalues of random...
10/15/2020

Solving Trust Region Subproblems Using Riemannian Optimization

The Trust Region Subproblem is a fundamental optimization problem that t...
07/18/2019

Fast approximation of orthogonal matrices and application to PCA

We study the problem of approximating orthogonal matrices so that their ...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.