Abstract:SIAM Journal on Matrix Analysis and Applications, Volume 45, Issue 4, Page 2007-2038, December 2024. Given a square pencil [math], where [math] and [math] are [math] complex (resp., real) matrices, we consider the problem of finding the singular complex (resp., real) pencil nearest to it in the Frobenius distance. This problem is known to be very difficult, and the few algorithms available in the literature can only deal efficiently with pencils of very small size. We show that the problem is equivalent to minimizing a certain objective function [math] over the Riemannian manifold [math] (resp., [math] if the nearest real singular pencil is sought), where [math] denotes the special unitary group (resp., [math] denotes the special orthogonal group). This novel perspective is based on the generalized Schur form of pencils, and yields competitive numerical methods, by pairing it with algorithms capable of doing optimization on Riemannian manifolds. We propose one algorithm that directly minimizes the (almost everywhere, but not everywhere, differentiable) function [math], as well as a smoothed alternative and a third algorithm that deals with smooth functions and can also solve the problem of finding a nearest singular pencil with a specified minimal index. We provide numerical experiments that show that the resulting methods allow us to deal with pencils of much larger size than alternative techniques, yielding candidate minimizers of comparable or better quality. In the course of our analysis, we also obtain a number of new theoretical results related to the generalized Schur form of a (regular or singular) square pencil and to the minimal index of a singular square pencil whose nullity is 1.

Weierstrass Structure and Eigenvalue Placement of Regular Matrix Pencils under Low Rank Perturbations

Generic Change of the Partial Multiplicities of Regular Matrix Pencils under Low-Rank Perturbations

Eigenvalue characterization of some structured matrix pencils under linear perturbation

On the change of the Weyr characteristics of matrix pencils after rank-one perturbations

The shift-and-invert Arnoldi method for singular matrix pencils

Generic Eigenstructures of Hermitian Pencils

Analysis of eigenvalue condition numbers for a class of randomized numerical methods for singular matrix pencils

Solving singular generalized eigenvalue problems. Part III: structure preservation

Deflating invariant subspaces for rank structured pencils

Generic symmetric matrix pencils with bounded rank

On the minimal ranks of matrix pencils and the existence of a best approximate block-term tensor decomposition

Global properties of eigenvalues of parametric rank one perturbations for unstructured and structured matrices II

Partial isospectrality of a matrix pencil and circularity of the c -numerical range

Eigenvalues of graph Laplacians via rank-one perturbations

A Riemannian Optimization Method to Compute the Nearest Singular Pencil

Structured Perturbations of Group Inverse and Singular Linear System with Index One

Computing a canonical form of a matrix pencil

First order structure-preserving perturbation theory for eigenvalues of symplectic matrices: Part II

A Schatten-$q$ Low-rank Matrix Perturbation Analysis via Perturbation Projection Error Bound

Nearest linearly structured polynomial matrix with some prescribed distinct eigenvalues