Thumbnail Image

On the nearest singular matrix pencil

Guglielmi, Nicola; Lubich, Christian; Mehrmann, Volker

Preprint-Reihe des Instituts für Mathematik, Technische Universität Berlin

Given a regular matrix pencil A + μE, we consider the problem of determining the nearest singular matrix pencil with respect to the Frobenius norm. We present new approaches based on the solution of matrix differential equations for determining the nearest singular pencil A + ΔA + μ(E + ΔE), one approach for general singular pencils and another one such that A+ ΔA and E + ΔE have a common left/right null vector. For the latter case the nearest singular pencil is shown to differ from the original pencil by rank-one matrices ΔA and ΔE. In both cases we consider also the situation where only A is perturbed. The nearest singular pencil is approached by a two-level iteration, where a gradient flow is driven to a stationary point in the inner iteration and the outer level uses a fast iteration for the distance parameter. This approach extends also to structured matrices A and E.