Characterization and construction of the nearest defective matrix via coalescence of pseudospectral components

Rafikul Alam, Shreemayee Bora, Ralph Byers, Michael L. Overton

Research output: Contribution to journalArticle

Abstract

Let A be a matrix with distinct eigenvalues and let w(A) be the distance from A to the set of defective matrices (using either the 2-norm or the Frobenius norm). Define Λ, the -pseudospectrum of A, to be the set of points in the complex plane which are eigenvalues of matrices A+E with E<, and let c(A) be the supremum of all with the property that Λ has n distinct components. Demmel and Wilkinson independently observed in the 1980s that w(A)≥c(A), and equality was established for the 2-norm by Alam and Bora (2005). We give new results on the geometry of the pseudospectrum near points where first coalescence of the components occurs, characterizing such points as the lowest generalized saddle point of the smallest singular value of A-zI over z∈C. One consequence is that w(A)=c(A) for the Frobenius norm too, and another is the perhaps surprising result that the minimal distance is attained by a defective matrix in all cases. Our results suggest a new computational approach to approximating the nearest defective matrix by a variant of Newton's method that is applicable to both generic and nongeneric cases. Construction of the nearest defective matrix involves some subtle numerical issues which we explain, and we present a simple backward error analysis showing that a certain singular vector residual measures how close the computed matrix is to a truly defective matrix. Finally, we present a result giving lower bounds on the angles of wedges contained in the pseudospectrum and emanating from generic coalescence points. Several conjectures and questions remain open.

Original languageEnglish (US)
Pages (from-to)494-513
Number of pages20
JournalLinear Algebra and Its Applications
Volume435
Issue number3
DOIs
StatePublished - Aug 1 2011

Fingerprint

Coalescence
Pseudospectra
Frobenius norm
Backward Error Analysis
Eigenvalue
Distinct
Norm
Singular Vectors
Newton-Raphson method
Singular Values
Saddlepoint
Wedge
Supremum
Newton Methods
Argand diagram
Error analysis
Set of points
Lowest
Equality
Lower bound

Keywords

  • Multiple eigen values
  • Pseudospectrum
  • Saddle point

ASJC Scopus subject areas

  • Algebra and Number Theory
  • Discrete Mathematics and Combinatorics
  • Geometry and Topology
  • Numerical Analysis

Cite this

Characterization and construction of the nearest defective matrix via coalescence of pseudospectral components. / Alam, Rafikul; Bora, Shreemayee; Byers, Ralph; Overton, Michael L.

In: Linear Algebra and Its Applications, Vol. 435, No. 3, 01.08.2011, p. 494-513.

Research output: Contribution to journalArticle

@article{e87e1242bad14195bba6a288713cd6ad,
title = "Characterization and construction of the nearest defective matrix via coalescence of pseudospectral components",
abstract = "Let A be a matrix with distinct eigenvalues and let w(A) be the distance from A to the set of defective matrices (using either the 2-norm or the Frobenius norm). Define Λ, the -pseudospectrum of A, to be the set of points in the complex plane which are eigenvalues of matrices A+E with E<, and let c(A) be the supremum of all with the property that Λ has n distinct components. Demmel and Wilkinson independently observed in the 1980s that w(A)≥c(A), and equality was established for the 2-norm by Alam and Bora (2005). We give new results on the geometry of the pseudospectrum near points where first coalescence of the components occurs, characterizing such points as the lowest generalized saddle point of the smallest singular value of A-zI over z∈C. One consequence is that w(A)=c(A) for the Frobenius norm too, and another is the perhaps surprising result that the minimal distance is attained by a defective matrix in all cases. Our results suggest a new computational approach to approximating the nearest defective matrix by a variant of Newton's method that is applicable to both generic and nongeneric cases. Construction of the nearest defective matrix involves some subtle numerical issues which we explain, and we present a simple backward error analysis showing that a certain singular vector residual measures how close the computed matrix is to a truly defective matrix. Finally, we present a result giving lower bounds on the angles of wedges contained in the pseudospectrum and emanating from generic coalescence points. Several conjectures and questions remain open.",
keywords = "Multiple eigen values, Pseudospectrum, Saddle point",
author = "Rafikul Alam and Shreemayee Bora and Ralph Byers and Overton, {Michael L.}",
year = "2011",
month = "8",
day = "1",
doi = "10.1016/j.laa.2010.09.022",
language = "English (US)",
volume = "435",
pages = "494--513",
journal = "Linear Algebra and Its Applications",
issn = "0024-3795",
publisher = "Elsevier Inc.",
number = "3",

}

TY - JOUR

T1 - Characterization and construction of the nearest defective matrix via coalescence of pseudospectral components

AU - Alam, Rafikul

AU - Bora, Shreemayee

AU - Byers, Ralph

AU - Overton, Michael L.

PY - 2011/8/1

Y1 - 2011/8/1

N2 - Let A be a matrix with distinct eigenvalues and let w(A) be the distance from A to the set of defective matrices (using either the 2-norm or the Frobenius norm). Define Λ, the -pseudospectrum of A, to be the set of points in the complex plane which are eigenvalues of matrices A+E with E<, and let c(A) be the supremum of all with the property that Λ has n distinct components. Demmel and Wilkinson independently observed in the 1980s that w(A)≥c(A), and equality was established for the 2-norm by Alam and Bora (2005). We give new results on the geometry of the pseudospectrum near points where first coalescence of the components occurs, characterizing such points as the lowest generalized saddle point of the smallest singular value of A-zI over z∈C. One consequence is that w(A)=c(A) for the Frobenius norm too, and another is the perhaps surprising result that the minimal distance is attained by a defective matrix in all cases. Our results suggest a new computational approach to approximating the nearest defective matrix by a variant of Newton's method that is applicable to both generic and nongeneric cases. Construction of the nearest defective matrix involves some subtle numerical issues which we explain, and we present a simple backward error analysis showing that a certain singular vector residual measures how close the computed matrix is to a truly defective matrix. Finally, we present a result giving lower bounds on the angles of wedges contained in the pseudospectrum and emanating from generic coalescence points. Several conjectures and questions remain open.

AB - Let A be a matrix with distinct eigenvalues and let w(A) be the distance from A to the set of defective matrices (using either the 2-norm or the Frobenius norm). Define Λ, the -pseudospectrum of A, to be the set of points in the complex plane which are eigenvalues of matrices A+E with E<, and let c(A) be the supremum of all with the property that Λ has n distinct components. Demmel and Wilkinson independently observed in the 1980s that w(A)≥c(A), and equality was established for the 2-norm by Alam and Bora (2005). We give new results on the geometry of the pseudospectrum near points where first coalescence of the components occurs, characterizing such points as the lowest generalized saddle point of the smallest singular value of A-zI over z∈C. One consequence is that w(A)=c(A) for the Frobenius norm too, and another is the perhaps surprising result that the minimal distance is attained by a defective matrix in all cases. Our results suggest a new computational approach to approximating the nearest defective matrix by a variant of Newton's method that is applicable to both generic and nongeneric cases. Construction of the nearest defective matrix involves some subtle numerical issues which we explain, and we present a simple backward error analysis showing that a certain singular vector residual measures how close the computed matrix is to a truly defective matrix. Finally, we present a result giving lower bounds on the angles of wedges contained in the pseudospectrum and emanating from generic coalescence points. Several conjectures and questions remain open.

KW - Multiple eigen values

KW - Pseudospectrum

KW - Saddle point

UR - http://www.scopus.com/inward/record.url?scp=79955652357&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=79955652357&partnerID=8YFLogxK

U2 - 10.1016/j.laa.2010.09.022

DO - 10.1016/j.laa.2010.09.022

M3 - Article

VL - 435

SP - 494

EP - 513

JO - Linear Algebra and Its Applications

JF - Linear Algebra and Its Applications

SN - 0024-3795

IS - 3

ER -