Weyl's inequality

From HandWiki

In linear algebra, Weyl's inequality is a theorem about the changes to eigenvalues of an Hermitian matrix that is perturbed. It can be used to estimate the eigenvalues of a perturbed Hermitian matrix.

Weyl's inequality about perturbation

Let [math]\displaystyle{ A }[/math] be Hermitian on inner product space [math]\displaystyle{ V }[/math] with dimention [math]\displaystyle{ n }[/math], with spectrum ordered in descending order [math]\displaystyle{ \lambda_1 \geq ... \geq \lambda_n }[/math]. Note that these eigenvalues can be ordered, because they are real (as eigenvalues of Hermitian matrices).[1]

Weyl inequality — [math]\displaystyle{ \lambda_{i+j-1}(A+B) \leq \lambda_i(A)+\lambda_j(B) \leq \lambda_{i+j-n}(A+B) }[/math]

Weyl's inequality states that the spectrum of Hermitian matrices is stable under perturbation. Specifically, we have:[1]

Corollary (Spectral stability) — [math]\displaystyle{ \lambda_k(A+B) - \lambda_k(A) \in [\lambda_n(B), \lambda_1(B)] }[/math] [math]\displaystyle{ |\lambda_k(A+B) - \lambda_k(A)| \leq \|B\|_{op} }[/math] where
[math]\displaystyle{ \|B\|_{op} = \max(|\lambda_1(B)|, |\lambda_n(B)|) }[/math] is the operator norm.

In jargon, it says that [math]\displaystyle{ \lambda_k }[/math] is Lipschitz-continuous on the space of Hermitian matrices with operator norm.

Weyl's inequality between eigenvalues and singular values

Let [math]\displaystyle{ A \in \mathbb{C}^{n \times n} }[/math] have singular values [math]\displaystyle{ \sigma_1(A) \geq \cdots \geq \sigma_n(A) \geq 0 }[/math] and eigenvalues ordered so that [math]\displaystyle{ |\lambda_1(A)| \geq \cdots \geq |\lambda_n(A)| }[/math]. Then

[math]\displaystyle{ |\lambda_1(A) \cdots \lambda_k(A)| \leq \sigma_1(A) \cdots \sigma_k(A) }[/math]

For [math]\displaystyle{ k = 1, \ldots, n }[/math], with equality for [math]\displaystyle{ k=n }[/math]. [2]

Applications

Estimating perturbations of the spectrum

Assume that [math]\displaystyle{ R }[/math] is small in the sense that its spectral norm satisfies [math]\displaystyle{ \|R\|_2 \le \epsilon }[/math] for some small [math]\displaystyle{ \epsilon\gt 0 }[/math]. Then it follows that all the eigenvalues of [math]\displaystyle{ R }[/math] are bounded in absolute value by [math]\displaystyle{ \epsilon }[/math]. Applying Weyl's inequality, it follows that the spectra of the Hermitian matrices M and N are close in the sense that[3]

[math]\displaystyle{ |\mu_i - \nu_i| \le \epsilon \qquad \forall i=1,\ldots,n. }[/math]

Note, however, that this eigenvalue perturbation bound is generally false for non-Hermitian matrices (or more accurately, for non-normal matrices). For a counterexample, let [math]\displaystyle{ t\gt 0 }[/math] be arbitrarily small, and consider

[math]\displaystyle{ M = \begin{bmatrix} 0 & 0 \\ 1/t^2 & 0 \end{bmatrix}, \qquad N = M + R = \begin{bmatrix} 0 & 1 \\ 1/t^2 & 0 \end{bmatrix}, \qquad R = \begin{bmatrix} 0 & 1 \\ 0 & 0 \end{bmatrix}. }[/math]

whose eigenvalues [math]\displaystyle{ \mu_1 = \mu_2 = 0 }[/math] and [math]\displaystyle{ \nu_1 = +1/t, \nu_2 = -1/t }[/math] do not satisfy [math]\displaystyle{ |\mu_i - \nu_i| \le \|R\|_2 = 1 }[/math].

Weyl's inequality for singular values

Let [math]\displaystyle{ M }[/math] be a [math]\displaystyle{ p \times n }[/math] matrix with [math]\displaystyle{ 1 \le p \le n }[/math]. Its singular values [math]\displaystyle{ \sigma_k(M) }[/math] are the [math]\displaystyle{ p }[/math] positive eigenvalues of the [math]\displaystyle{ (p+n) \times (p+n) }[/math] Hermitian augmented matrix

[math]\displaystyle{ \begin{bmatrix} 0 & M \\ M^* & 0 \end{bmatrix}. }[/math]

Therefore, Weyl's eigenvalue perturbation inequality for Hermitian matrices extends naturally to perturbation of singular values.[1] This result gives the bound for the perturbation in the singular values of a matrix [math]\displaystyle{ M }[/math] due to an additive perturbation [math]\displaystyle{ \Delta }[/math]:

[math]\displaystyle{ |\sigma_k(M+\Delta) - \sigma_k(M)| \le \sigma_1(\Delta) }[/math]

where we note that the largest singular value [math]\displaystyle{ \sigma_1(\Delta) }[/math] coincides with the spectral norm [math]\displaystyle{ \|\Delta\|_2 }[/math].

Notes

  1. 1.0 1.1 1.2 Tao, Terence (2010-01-13). "254A, Notes 3a: Eigenvalues and sums of Hermitian matrices". https://terrytao.wordpress.com/2010/01/12/254a-notes-3a-eigenvalues-and-sums-of-hermitian-matrices/. Retrieved 25 May 2015. 
  2. Roger A. Horn, and Charles R. Johnson Topics in Matrix Analysis. Cambridge, 1st Edition, 1991. p.171
  3. Weyl, Hermann. "Das asymptotische Verteilungsgesetz der Eigenwerte linearer partieller Differentialgleichungen (mit einer Anwendung auf die Theorie der Hohlraumstrahlung)." Mathematische Annalen 71, no. 4 (1912): 441-479.

References

  • Matrix Theory, Joel N. Franklin, (Dover Publications, 1993) ISBN:0-486-41179-6
  • "Das asymptotische Verteilungsgesetz der Eigenwerte linearer partieller Differentialgleichungen", H. Weyl, Math. Ann., 71 (1912), 441–479