Mehmet Emre Tiryaki
Back to Notes

Eigenvalues of Hermitian Matrices

Throughout, the eigenvalues of a Hermitian A∈MnA \in M_n are ordered λ1≤λ2≤⋯≤λn\lambda_1 \le \lambda_2 \le \cdots \le \lambda_n.

Field of Values

For A∈MnA \in M_n,

F(A)={x∗Axx∗x:x∈Cn,x≠0}={x∗Ax:x∈Cn,∥x∥2=1}F(A) = \left\{\frac{x^*Ax}{x^*x} : x \in \mathbb{C}^n, x \neq 0\right\} = \{x^*Ax : x \in \mathbb{C}^n, \|x\|_2 = 1\}

If AA is normal, F(A)=Co(σ(A))F(A) = \mathrm{Co}(\sigma(A)), the convex hull of the eigenvalues.

Proof Let S\mathbb{S} be the unit sphere and A=UDU∗A = UDU^*. Since U∗U^* maps S\mathbb{S} onto S\mathbb{S},

F(A)={x∗UDU∗x:x∈S}={y∗Dy:y∈S}={∑i=1n∣yi∣2dii:∑i=1n∣yi∣2=1}F(A) = \{x^*UDU^*x : x \in \mathbb{S}\} = \{y^*Dy : y \in \mathbb{S}\} = \Big\{\sum_{i=1}^{n} |y_i|^2 d_{ii} : \sum_{i=1}^{n} |y_i|^2 = 1\Big\}

which are exactly the convex combinations of the eigenvalues.

Rayleigh-Ritz

A∈MnA \in M_n Hermitian is PSD if x∗Ax≥0x^*Ax \ge 0 for all x∈Cnx \in \mathbb{C}^n.

For A∈MnA \in M_n Hermitian,

min⁡x≠0x∗Axx∗x=λ1(A),max⁡x≠0x∗Axx∗x=λn(A)\min_{x \neq 0} \frac{x^*Ax}{x^*x} = \lambda_1(A), \qquad \max_{x \neq 0} \frac{x^*Ax}{x^*x} = \lambda_n(A)

More generally, let AA have orthonormal eigenvectors μ1,…,μn\mu_1, \ldots, \mu_n with eigenvalues λ1≤⋯≤λn\lambda_1 \le \cdots \le \lambda_n. For all k=1,…,nk = 1, \ldots, n

min⁡x≠0x⊥μ1,…,μk−1x∗Axx∗x=λk,max⁡x≠0x⊥μk+1,…,μnx∗Axx∗x=λk\min_{\substack{x \neq 0 \\ x \perp \mu_1, \ldots, \mu_{k-1}}} \frac{x^*Ax}{x^*x} = \lambda_k, \qquad \max_{\substack{x \neq 0 \\ x \perp \mu_{k+1}, \ldots, \mu_n}} \frac{x^*Ax}{x^*x} = \lambda_k

Proof Let A=UDU∗A = UDU^* with U=[μ1∣⋯∣μn]U = [\mu_1 | \cdots | \mu_n] and D=diag(λ1,…,λn)D = \mathrm{diag}(\lambda_1, \ldots, \lambda_n). x⊥μ1,…,μk−1x \perp \mu_1, \ldots, \mu_{k-1} iff x∈Span{μk,…,μn}x \in \mathrm{Span}\{\mu_k, \ldots, \mu_n\}, i.e. x=∑i=knδiμix = \sum_{i=k}^{n} \delta_i \mu_i with δi\delta_i not all 0. Then

x∗Axx∗x=∑i=knδˉiδiλi∑i=knδˉiδi=∑j=knαjλj,αj=∣δj∣2∑i=kn∣δi∣2\frac{x^*Ax}{x^*x} = \frac{\sum_{i=k}^{n} \bar{\delta}_i \delta_i \lambda_i}{\sum_{i=k}^{n} \bar{\delta}_i \delta_i} = \sum_{j=k}^{n} \alpha_j \lambda_j, \qquad \alpha_j = \frac{|\delta_j|^2}{\sum_{i=k}^{n} |\delta_i|^2}

with αj≥0\alpha_j \ge 0 and ∑αj=1\sum \alpha_j = 1. This is a convex combination of λk,…,λn\lambda_k, \ldots, \lambda_n, so the minimum is min⁡Co(λk,…,λn)=λk\min \mathrm{Co}(\lambda_k, \ldots, \lambda_n) = \lambda_k, attained at x=μkx = \mu_k. The max is the same argument.

Courant-Fisher

For A∈MnA \in M_n Hermitian and all k=1,…,nk = 1, \ldots, n,

max⁡y1,…,yk−1∈Cn(min⁡x≠0x⊥y1,…,yk−1x∗Axx∗x)=λk,min⁡yk+1,…,yn∈Cn(max⁡x≠0x⊥yk+1,…,ynx∗Axx∗x)=λk\max_{y_1, \ldots, y_{k-1} \in \mathbb{C}^n} \left( \min_{\substack{x \neq 0 \\ x \perp y_1, \ldots, y_{k-1}}} \frac{x^*Ax}{x^*x} \right) = \lambda_k, \qquad \min_{y_{k+1}, \ldots, y_n \in \mathbb{C}^n} \left( \max_{\substack{x \neq 0 \\ x \perp y_{k+1}, \ldots, y_n}} \frac{x^*Ax}{x^*x} \right) = \lambda_k

Proof

  1. For any y1,…,yk−1y_1, \ldots, y_{k-1}, Span{μ1,…,μk}\mathrm{Span}\{\mu_1, \ldots, \mu_k\} has dimension kk and Span{y1,…,yk−1}⊥\mathrm{Span}\{y_1, \ldots, y_{k-1}\}^\perp has dimension ≥n−k+1\ge n - k + 1, so their intersection is not {0⃗}\{\vec{0}\}.
  2. Let ω≠0\omega \neq 0 be in the intersection. ω∗Aωω∗ω\frac{\omega^*A\omega}{\omega^*\omega} is a convex combination of λ1,…,λk\lambda_1, \ldots, \lambda_k, so it is ≤λk\le \lambda_k.
  3. So the inner min is ≤λk\le \lambda_k for every choice of yy's, with equality at yi=μiy_i = \mu_i by Rayleigh-Ritz.

Weyl's Inequality

For any A,B∈MnA, B \in M_n Hermitian and all k=1,…,nk = 1, \ldots, n,

λk(A)+λ1(B)≤λk(A+B)≤λk(A)+λn(B)\lambda_k(A) + \lambda_1(B) \le \lambda_k(A + B) \le \lambda_k(A) + \lambda_n(B)

Proof Using the min-max form of Courant-Fisher and Rayleigh-Ritz on BB,

λk(A+B)=min⁡ymax⁡x⊥y(x∗Axx∗x+x∗Bxx∗x)≤min⁡ymax⁡x⊥y(x∗Axx∗x+λn(B))=λk(A)+λn(B)\lambda_k(A + B) = \min_{y} \max_{x \perp y} \left(\frac{x^*Ax}{x^*x} + \frac{x^*Bx}{x^*x}\right) \le \min_{y} \max_{x \perp y} \left(\frac{x^*Ax}{x^*x} + \lambda_n(B)\right) = \lambda_k(A) + \lambda_n(B)

The lower bound is the same with x∗Bxx∗x≥λ1(B)\frac{x^*Bx}{x^*x} \ge \lambda_1(B).

If BB is also PSD, then λ1(B)≥0\lambda_1(B) \ge 0, so for all kk

λk(A)≤λk(A)+λ1(B)≤λk(A+B)\lambda_k(A) \le \lambda_k(A) + \lambda_1(B) \le \lambda_k(A + B)

Interlacing Theorem I

For A∈MnA \in M_n Hermitian, a∈Ra \in \mathbb{R}, z∈Cnz \in \mathbb{C}^n and all k=1,…,n−1k = 1, \ldots, n - 1,

λk(A+azz∗)≤λk+1(A)\lambda_k(A + azz^*) \le \lambda_{k+1}(A)

Proof Adding the constraint x⊥zx \perp z can only decrease the inner min, and on x⊥zx \perp z the term azz∗azz^* vanishes,

λk(A+azz∗)=max⁡y1,…,yk−1min⁡x≠0x⊥y1,…,yk−1x∗(A+azz∗)xx∗x≤max⁡y1,…,yk−1min⁡x≠0x⊥y1,…,yk−1,zx∗Axx∗x≤max⁡y1,…,ykmin⁡x≠0x⊥y1,…,ykx∗Axx∗x=λk+1(A)\begin{aligned} \lambda_k(A + azz^*) &= \max_{y_1, \ldots, y_{k-1}} \min_{\substack{x \neq 0 \\ x \perp y_1, \ldots, y_{k-1}}} \frac{x^*(A + azz^*)x}{x^*x} \le \max_{y_1, \ldots, y_{k-1}} \min_{\substack{x \neq 0 \\ x \perp y_1, \ldots, y_{k-1}, z}} \frac{x^*Ax}{x^*x} \\ &\le \max_{y_1, \ldots, y_k} \min_{\substack{x \neq 0 \\ x \perp y_1, \ldots, y_k}} \frac{x^*Ax}{x^*x} = \lambda_{k+1}(A) \end{aligned}

Corollary For A,B∈MnA, B \in M_n Hermitian with rank(B)=r\mathrm{rank}(B) = r, and all kk for which the indices make sense,

λk(A+B)≤λk+r(A)\lambda_k(A + B) \le \lambda_{k+r}(A)

Write B=Udiag(d11,…,drr,0,…,0)U∗B = U \mathrm{diag}(d_{11}, \ldots, d_{rr}, 0, \ldots, 0) U^*, so B=∑i=1rdiiUiUi∗B = \sum_{i=1}^{r} d_{ii} U_i U_i^* is a sum of rr terms of the form azz∗azz^*. The rest follows by induction.

Applying this to A+BA + B and −B-B (also rank rr) gives both directions,

λk−r(A)≤λk(A+B)≤λk+r(A),λk−r(A+B)≤λk(A)≤λk+r(A+B)\lambda_{k-r}(A) \le \lambda_k(A + B) \le \lambda_{k+r}(A), \qquad \lambda_{k-r}(A + B) \le \lambda_k(A) \le \lambda_{k+r}(A + B)

Interlacing Theorem II

Let A∈MnA \in M_n be Hermitian and BB an r×rr \times r principal submatrix of AA, which is Hermitian by default. For all k=1,…,rk = 1, \ldots, r,

λk(A)≤λk(B)≤λk+n−r(A)\lambda_k(A) \le \lambda_k(B) \le \lambda_{k+n-r}(A)

Proof Let the deleted rows and columns be indexed by i1,…,in−ri_1, \ldots, i_{n-r}. Adding constraints can only increase the inner min, so

λk(A)=max⁡y1,…,yk−1∈Cnmin⁡x≠0x⊥y1,…,yk−1x∗Axx∗x≤max⁡y1,…,yk−1∈Cnmin⁡x≠0x⊥y1,…,yk−1x⊥ei1,…,ein−rx∗Axx∗x=max⁡v1,…,vk−1∈Crmin⁡z≠0z⊥v1,…,vk−1z∗Bzz∗z=λk(B)\begin{aligned} \lambda_k(A) &= \max_{y_1, \ldots, y_{k-1} \in \mathbb{C}^n} \min_{\substack{x \neq 0 \\ x \perp y_1, \ldots, y_{k-1}}} \frac{x^*Ax}{x^*x} \le \max_{y_1, \ldots, y_{k-1} \in \mathbb{C}^n} \min_{\substack{x \neq 0 \\ x \perp y_1, \ldots, y_{k-1} \\ x \perp e_{i_1}, \ldots, e_{i_{n-r}}}} \frac{x^*Ax}{x^*x} \\ &= \max_{v_1, \ldots, v_{k-1} \in \mathbb{C}^r} \min_{\substack{z \neq 0 \\ z \perp v_1, \ldots, v_{k-1}}} \frac{z^*Bz}{z^*z} = \lambda_k(B) \end{aligned}

The upper bound is the same argument with the min-max form.

Corollary Taking r=1r = 1, for all kk

λ1(A)≤akk≤λn(A)\lambda_1(A) \le a_{kk} \le \lambda_n(A)