Nearest Ω -stable matrix via Riemannian optimization

Vanni Noferini*, Federico Poloni

*Corresponding author for this work

Research output: Contribution to journalArticleScientificpeer-review

4 Downloads (Pure)

Abstract

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.

Original languageEnglish
Pages (from-to)817–851
JournalNumerische Mathematik
DOIs
Publication statusE-pub ahead of print - 2021
MoE publication typeA1 Journal article-refereed

Fingerprint

Dive into the research topics of 'Nearest Ω -stable matrix via Riemannian optimization'. Together they form a unique fingerprint.

Cite this