Throughout, the eigenvalues of a Hermitian A ∈ M n A \in M_n A ∈ M n are ordered λ 1 ≤ λ 2 ≤ ⋯ ≤ λ n \lambda_1 \le \lambda_2 \le \cdots \le \lambda_n λ 1 ≤ λ 2 ≤ ⋯ ≤ λ n .
Field of Values
For A ∈ M n A \in M_n A ∈ M n ,
F ( A ) = { x ∗ A x x ∗ x : x ∈ C n , x ≠ 0 } = { x ∗ A x : x ∈ C n , ∥ 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\} F ( A ) = { x ∗ x x ∗ A x : x ∈ C n , x = 0 } = { x ∗ A x : x ∈ C n , ∥ x ∥ 2 = 1 }
If A A A is normal, F ( A ) = C o ( σ ( A ) ) F(A) = \mathrm{Co}(\sigma(A)) F ( A ) = Co ( σ ( A )) , the convex hull of the eigenvalues.
Proof Let S \mathbb{S} S be the unit sphere and A = U D U ∗ A = UDU^* A = U D U ∗ . Since U ∗ U^* U ∗ maps S \mathbb{S} S onto S \mathbb{S} S ,
F ( A ) = { x ∗ U D U ∗ x : x ∈ S } = { y ∗ D y : y ∈ S } = { ∑ i = 1 n ∣ y i ∣ 2 d i i : ∑ i = 1 n ∣ y i ∣ 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\} F ( A ) = { x ∗ U D U ∗ x : x ∈ S } = { y ∗ Dy : y ∈ S } = { i = 1 ∑ n ∣ y i ∣ 2 d ii : i = 1 ∑ n ∣ y i ∣ 2 = 1 }
which are exactly the convex combinations of the eigenvalues.
Rayleigh-Ritz
A ∈ M n A \in M_n A ∈ M n Hermitian is PSD if x ∗ A x ≥ 0 x^*Ax \ge 0 x ∗ A x ≥ 0 for all x ∈ C n x \in \mathbb{C}^n x ∈ C n .
For A ∈ M n A \in M_n A ∈ M n Hermitian,
min x ≠ 0 x ∗ A x x ∗ x = λ 1 ( A ) , max x ≠ 0 x ∗ A x x ∗ 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) x = 0 min x ∗ x x ∗ A x = λ 1 ( A ) , x = 0 max x ∗ x x ∗ A x = λ n ( A )
More generally, let A A A have orthonormal eigenvectors μ 1 , … , μ n \mu_1, \ldots, \mu_n μ 1 , … , μ n with eigenvalues λ 1 ≤ ⋯ ≤ λ n \lambda_1 \le \cdots \le \lambda_n λ 1 ≤ ⋯ ≤ λ n . For all k = 1 , … , n k = 1, \ldots, n k = 1 , … , n
min x ≠ 0 x ⊥ μ 1 , … , μ k − 1 x ∗ A x x ∗ x = λ k , max x ≠ 0 x ⊥ μ k + 1 , … , μ n x ∗ A x x ∗ 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 x = 0 x ⊥ μ 1 , … , μ k − 1 min x ∗ x x ∗ A x = λ k , x = 0 x ⊥ μ k + 1 , … , μ n max x ∗ x x ∗ A x = λ k
Proof Let A = U D U ∗ A = UDU^* A = U D U ∗ with U = [ μ 1 ∣ ⋯ ∣ μ n ] U = [\mu_1 | \cdots | \mu_n] U = [ μ 1 ∣ ⋯ ∣ μ n ] and D = d i a g ( λ 1 , … , λ n ) D = \mathrm{diag}(\lambda_1, \ldots, \lambda_n) D = diag ( λ 1 , … , λ n ) . x ⊥ μ 1 , … , μ k − 1 x \perp \mu_1, \ldots, \mu_{k-1} x ⊥ μ 1 , … , μ k − 1 iff x ∈ S p a n { μ k , … , μ n } x \in \mathrm{Span}\{\mu_k, \ldots, \mu_n\} x ∈ Span { μ k , … , μ n } , i.e. x = ∑ i = k n δ i μ i x = \sum_{i=k}^{n} \delta_i \mu_i x = ∑ i = k n δ i μ i with δ i \delta_i δ i not all 0. Then
x ∗ A x x ∗ x = ∑ i = k n δ ˉ i δ i λ i ∑ i = k n δ ˉ i δ i = ∑ j = k n α j λ j , α j = ∣ δ j ∣ 2 ∑ i = k n ∣ δ 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} x ∗ x x ∗ A x = ∑ i = k n δ ˉ i δ i ∑ i = k n δ ˉ i δ i λ i = j = k ∑ n α j λ j , α j = ∑ i = k n ∣ δ i ∣ 2 ∣ δ j ∣ 2
with α j ≥ 0 \alpha_j \ge 0 α j ≥ 0 and ∑ α j = 1 \sum \alpha_j = 1 ∑ α j = 1 . This is a convex combination of λ k , … , λ n \lambda_k, \ldots, \lambda_n λ k , … , λ n , so the minimum is min C o ( λ k , … , λ n ) = λ k \min \mathrm{Co}(\lambda_k, \ldots, \lambda_n) = \lambda_k min Co ( λ k , … , λ n ) = λ k , attained at x = μ k x = \mu_k x = μ k . The max is the same argument.
Courant-Fisher
For A ∈ M n A \in M_n A ∈ M n Hermitian and all k = 1 , … , n k = 1, \ldots, n k = 1 , … , n ,
max y 1 , … , y k − 1 ∈ C n ( min x ≠ 0 x ⊥ y 1 , … , y k − 1 x ∗ A x x ∗ x ) = λ k , min y k + 1 , … , y n ∈ C n ( max x ≠ 0 x ⊥ y k + 1 , … , y n x ∗ A x x ∗ 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 y 1 , … , y k − 1 ∈ C n max x = 0 x ⊥ y 1 , … , y k − 1 min x ∗ x x ∗ A x = λ k , y k + 1 , … , y n ∈ C n min x = 0 x ⊥ y k + 1 , … , y n max x ∗ x x ∗ A x = λ k
Proof
For any y 1 , … , y k − 1 y_1, \ldots, y_{k-1} y 1 , … , y k − 1 , S p a n { μ 1 , … , μ k } \mathrm{Span}\{\mu_1, \ldots, \mu_k\} Span { μ 1 , … , μ k } has dimension k k k and S p a n { y 1 , … , y k − 1 } ⊥ \mathrm{Span}\{y_1, \ldots, y_{k-1}\}^\perp Span { y 1 , … , y k − 1 } ⊥ has dimension ≥ n − k + 1 \ge n - k + 1 ≥ n − k + 1 , so their intersection is not { 0 ⃗ } \{\vec{0}\} { 0 } .
Let ω ≠ 0 \omega \neq 0 ω = 0 be in the intersection. ω ∗ A ω ω ∗ ω \frac{\omega^*A\omega}{\omega^*\omega} ω ∗ ω ω ∗ A ω is a convex combination of λ 1 , … , λ k \lambda_1, \ldots, \lambda_k λ 1 , … , λ k , so it is ≤ λ k \le \lambda_k ≤ λ k .
So the inner min is ≤ λ k \le \lambda_k ≤ λ k for every choice of y y y 's, with equality at y i = μ i y_i = \mu_i y i = μ i by Rayleigh-Ritz.
Weyl's Inequality
For any A , B ∈ M n A, B \in M_n A , B ∈ M n Hermitian and all k = 1 , … , n k = 1, \ldots, n k = 1 , … , 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) λ k ( A ) + λ 1 ( B ) ≤ λ k ( A + B ) ≤ λ k ( A ) + λ n ( B )
Proof Using the min-max form of Courant-Fisher and Rayleigh-Ritz on B B B ,
λ k ( A + B ) = min y max x ⊥ y ( x ∗ A x x ∗ x + x ∗ B x x ∗ x ) ≤ min y max x ⊥ y ( x ∗ A x x ∗ 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) λ k ( A + B ) = y min x ⊥ y max ( x ∗ x x ∗ A x + x ∗ x x ∗ B x ) ≤ y min x ⊥ y max ( x ∗ x x ∗ A x + λ n ( B ) ) = λ k ( A ) + λ n ( B )
The lower bound is the same with x ∗ B x x ∗ x ≥ λ 1 ( B ) \frac{x^*Bx}{x^*x} \ge \lambda_1(B) x ∗ x x ∗ B x ≥ λ 1 ( B ) .
If B B B is also PSD, then λ 1 ( B ) ≥ 0 \lambda_1(B) \ge 0 λ 1 ( B ) ≥ 0 , so for all k k k
λ k ( A ) ≤ λ k ( A ) + λ 1 ( B ) ≤ λ k ( A + B ) \lambda_k(A) \le \lambda_k(A) + \lambda_1(B) \le \lambda_k(A + B) λ k ( A ) ≤ λ k ( A ) + λ 1 ( B ) ≤ λ k ( A + B )
Interlacing Theorem I
For A ∈ M n A \in M_n A ∈ M n Hermitian, a ∈ R a \in \mathbb{R} a ∈ R , z ∈ C n z \in \mathbb{C}^n z ∈ C n and all k = 1 , … , n − 1 k = 1, \ldots, n - 1 k = 1 , … , n − 1 ,
λ k ( A + a z z ∗ ) ≤ λ k + 1 ( A ) \lambda_k(A + azz^*) \le \lambda_{k+1}(A) λ k ( A + a z z ∗ ) ≤ λ k + 1 ( A )
Proof Adding the constraint x ⊥ z x \perp z x ⊥ z can only decrease the inner min, and on x ⊥ z x \perp z x ⊥ z the term a z z ∗ azz^* a z z ∗ vanishes,
λ k ( A + a z z ∗ ) = max y 1 , … , y k − 1 min x ≠ 0 x ⊥ y 1 , … , y k − 1 x ∗ ( A + a z z ∗ ) x x ∗ x ≤ max y 1 , … , y k − 1 min x ≠ 0 x ⊥ y 1 , … , y k − 1 , z x ∗ A x x ∗ x ≤ max y 1 , … , y k min x ≠ 0 x ⊥ y 1 , … , y k x ∗ A x x ∗ 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} λ k ( A + a z z ∗ ) = y 1 , … , y k − 1 max x = 0 x ⊥ y 1 , … , y k − 1 min x ∗ x x ∗ ( A + a z z ∗ ) x ≤ y 1 , … , y k − 1 max x = 0 x ⊥ y 1 , … , y k − 1 , z min x ∗ x x ∗ A x ≤ y 1 , … , y k max x = 0 x ⊥ y 1 , … , y k min x ∗ x x ∗ A x = λ k + 1 ( A )
Corollary For A , B ∈ M n A, B \in M_n A , B ∈ M n Hermitian with r a n k ( B ) = r \mathrm{rank}(B) = r rank ( B ) = r , and all k k k for which the indices make sense,
λ k ( A + B ) ≤ λ k + r ( A ) \lambda_k(A + B) \le \lambda_{k+r}(A) λ k ( A + B ) ≤ λ k + r ( A )
Write B = U d i a g ( d 11 , … , d r r , 0 , … , 0 ) U ∗ B = U \mathrm{diag}(d_{11}, \ldots, d_{rr}, 0, \ldots, 0) U^* B = U diag ( d 11 , … , d rr , 0 , … , 0 ) U ∗ , so B = ∑ i = 1 r d i i U i U i ∗ B = \sum_{i=1}^{r} d_{ii} U_i U_i^* B = ∑ i = 1 r d ii U i U i ∗ is a sum of r r r terms of the form a z z ∗ azz^* a z z ∗ . The rest follows by induction.
Applying this to A + B A + B A + B and − B -B − B (also rank r r r ) 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) λ k − r ( A ) ≤ λ k ( A + B ) ≤ λ k + r ( A ) , λ k − r ( A + B ) ≤ λ k ( A ) ≤ λ k + r ( A + B )
Interlacing Theorem II
Let A ∈ M n A \in M_n A ∈ M n be Hermitian and B B B an r × r r \times r r × r principal submatrix of A A A , which is Hermitian by default. For all k = 1 , … , r k = 1, \ldots, r k = 1 , … , r ,
λ k ( A ) ≤ λ k ( B ) ≤ λ k + n − r ( A ) \lambda_k(A) \le \lambda_k(B) \le \lambda_{k+n-r}(A) λ k ( A ) ≤ λ k ( B ) ≤ λ k + n − r ( A )
Proof Let the deleted rows and columns be indexed by i 1 , … , i n − r i_1, \ldots, i_{n-r} i 1 , … , i n − r . Adding constraints can only increase the inner min, so
λ k ( A ) = max y 1 , … , y k − 1 ∈ C n min x ≠ 0 x ⊥ y 1 , … , y k − 1 x ∗ A x x ∗ x ≤ max y 1 , … , y k − 1 ∈ C n min x ≠ 0 x ⊥ y 1 , … , y k − 1 x ⊥ e i 1 , … , e i n − r x ∗ A x x ∗ x = max v 1 , … , v k − 1 ∈ C r min z ≠ 0 z ⊥ v 1 , … , v k − 1 z ∗ B z z ∗ 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} λ k ( A ) = y 1 , … , y k − 1 ∈ C n max x = 0 x ⊥ y 1 , … , y k − 1 min x ∗ x x ∗ A x ≤ y 1 , … , y k − 1 ∈ C n max x = 0 x ⊥ y 1 , … , y k − 1 x ⊥ e i 1 , … , e i n − r min x ∗ x x ∗ A x = v 1 , … , v k − 1 ∈ C r max z = 0 z ⊥ v 1 , … , v k − 1 min z ∗ z z ∗ B z = λ k ( B )
The upper bound is the same argument with the min-max form.
Corollary Taking r = 1 r = 1 r = 1 , for all k k k
λ 1 ( A ) ≤ a k k ≤ λ n ( A ) \lambda_1(A) \le a_{kk} \le \lambda_n(A) λ 1 ( A ) ≤ a kk ≤ λ n ( A )