17.2  The Min–Max Theorem

Section 1 ended with a description of λk(A) that works and cannot be used: a maximum of the Rayleigh quotient over the orthogonal complement of A’s own top eigenvectors. The defect is that the constraint set is built from the answer. This section repairs it by a move that looks like an extravagance and is exactly the right one: keep the maximum over a subspace, stop insisting that the subspace be the right one, and let it range over all subspaces of the appropriate dimension. The minimum of what comes out is again λk(A), and no eigenvector appears anywhere in the statement. That is the Courant–Fischer theorem, and everything in the rest of the chapter is a consequence of it.

Throughout, F is R or C and every matrix is Hermitian. The eigenvalues of a Hermitian A∈Mn(F) are real (Theorem 12.2.1) and are indexed decreasingly, λ1(A) ≥ λ2(A) ≥ … ≥ λn(A), each repeated according to multiplicity, so that the list has exactly n entries. This is the book’s convention, fixed in Chapter 13 §05, and it is used without exception below. Many sources index the theorem of this section the other way round; the remark after the theorem gives that form, once.

17.2.1. Two subspaces that have to meet

One dimension count does all the work in this chapter, and it is worth a label of its own, because four later proofs cite it.

Two subspaces of an n-dimensional space cannot both be large and disjoint: if their dimensions add to more than n, they share a non-zero vector.

Lemma 17.2.1 (Large Subspaces Intersect)

Let V be a vector space with dim⁡V=n over any field, and let U and W be subspaces of V with dim⁡U+dim⁡W>n. Then U∩W≠{0}. More precisely, dim⁡(U∩W)≥dim⁡U+dim⁡W−n.

Proof. By Theorem 1.7.3, dim⁡(U+W)=dim⁡U+dim⁡W−dim⁡(U∩W). Both U and W are finite-dimensional, being subspaces of the n-dimensional V (Theorem 1.6.7); U+W is a subspace of V (Theorem 1.7.1), so the same theorem gives dim⁡(U+W)≤n, and rearranging gives dim⁡(U∩W) ≥ dim⁡U+dim⁡W−n, which is >0 by hypothesis. A subspace of positive dimension contains a non-zero vector. This proves the lemma.

Remark. The bound is sharp, and the hypothesis cannot be weakened to dim⁡U+dim⁡W=n: two distinct lines in R2 meet only at 0, and so do a line and a plane in R3 chosen in general position. What makes the count useful is that it needs nothing about the two subspaces beyond their dimensions — in every application below, one of them is a span of eigenvectors that we choose and the other is an arbitrary competitor we know nothing about.

17.2.2. The quotient on a subspace

Before quantifying over subspaces, we must know that the inner maximum and minimum exist. They do, and for an algebraic reason: restricting the quadratic form to a subspace produces a smaller Hermitian matrix, to which Proposition 17.1.1 applies verbatim.

Lemma 17.2.2 (The Rayleigh Quotient on a Subspace)

Let A∈Mn(F) be Hermitian and let W⊆Fn be a subspace with dim⁡W=m≥1. Let Q∈Mn×m(F) be a matrix whose columns form an orthonormal basis of W, and put B=Q∗AQ∈Mm(F). Then B is Hermitian and {RA(x):0≠x∈W}=[λm(B), λ1(B)]. In particular max0≠x∈WRA(x) and min0≠x∈WRA(x) exist and are attained.

Idea. Coordinates on W. Writing a vector of W as Qc turns ⟨Ax,x⟩ into ⟨Bc,c⟩ and, because the columns of Q are orthonormal, turns ⟨x,x⟩ into ⟨c,c⟩ with no distortion. So the Rayleigh quotient of A on W is the Rayleigh quotient of B on Fm, and Proposition 17.1.1 answers the question there.

Proof. Such a Q exists: W has an orthonormal basis by Theorem 11.2.2, and we take its members as columns. Orthonormality of the columns says exactly Q∗Q=Im. Then B∗=Q∗A∗Q=Q∗AQ=B, so B is Hermitian and Proposition 17.1.1 applies to it.

The map c↦Qc sends Fm onto W, since the columns of Q span W, and it is injective, since Qc=0 gives c=Q∗Qc=0. So it is a bijection Fm→W carrying 0 to 0, hence a bijection from the non-zero vectors of Fm onto those of W. For c≠0, writing the standard inner product as ⟨u,v⟩=v∗u, ⟨AQc,Qc⟩=c∗Q∗AQc=⟨Bc,c⟩,⟨Qc,Qc⟩=c∗Q∗Qc=⟨c,c⟩, so RA(Qc)=RB(c). Therefore the two value sets {RA(x):0≠x∈W} and {RB(c):0≠c∈Fm} are equal, and the second is [λm(B),λ1(B)] with both ends attained, by Proposition 17.1.1. This proves the lemma.

Remark. Which analysis was invoked: none. The existence of a maximum of a function over an infinite set is normally an analytic fact — the set of unit vectors of W is compact, by fact (A3) of Chapter 16’s introduction, and RA is continuous on it, so the extreme value theorem, fact (A4), applies. That route is available and correct. It is not needed. The spectral theorem, already proved in Chapter 12, produces the maximizing vector by name, and Lemma 17.2.2 inherits that. Nothing in this section, or in §01, uses an item of Chapter 16’s imported-analysis list. The analysis has not vanished, only moved: the spectral theorem rests on the fundamental theorem of algebra (Theorem 6.5.1), whose proof in Chapter 6 §05 uses the extreme value theorem on a closed disc.

17.2.3. Every eigenvalue is a min–max

Theorem 17.2.3 (Courant–Fischer Min–Max Theorem)

Let F be R or C, let A∈Mn(F) be Hermitian and let 1≤k≤n. Then λk(A)=mindim⁡W=n−k+1  max0≠x∈W RA(x)=maxdim⁡W=k  min0≠x∈W RA(x), where in each line W runs over all subspaces of Fn of the stated dimension. All four extrema are attained.

Idea. Read the second line first, because it is the one with a slogan: λk is the best value you can guarantee on a k-dimensional subspace. On a k-dimensional subspace the quotient has a smallest value, which is what that subspace guarantees; the span of the top k eigenvectors guarantees λk, and no other k-dimensional subspace guarantees more.

The proof is three steps, and the first two have entirely different characters. Step 1 exhibits a witness subspace and computes on it, giving one inequality. Step 2 takes an arbitrary competitor subspace and shows it cannot beat the witness; the only thing known about a competitor is its dimension, so the only available tool is Lemma 17.2.1, which forces the competitor to meet a span of eigenvectors from the other end of the list. On that shared non-zero vector the quotient is pinned, and the inequality follows. Step 3 obtains the second line from the first by applying it to −A, which reverses the order of the eigenvalues and turns every maximum into a minimum.

Proof. By Corollary 12.4.3 over C, or Corollary 12.5.2 over R, fix an orthonormal basis (q1,…,qn) of Fn with Aqi=λi(A)qi for every i. Abbreviate λi=λi(A) and put Uj=span(q1,…,qj),Lj=span(qj,…,qn), so that dim⁡Uj=j and dim⁡Lj=n−j+1, the lists being orthonormal and hence independent. By Lemma 17.2.2 the inner extrema below all exist.

Claim. For every j: RA(x)≥λj for all 0≠x∈Uj, and RA(x)≤λj for all 0≠x∈Lj. Both bounds are attained at qj.

Proof. Let 0≠x∈Uj and write ci=⟨x,qi⟩, so x=∑i≤jciqi by Theorem 11.2.1 (a) and the definition of Uj. Then Ax=∑i≤jλiciqi, and Theorem 11.2.1 (b), (c) give RA(x)=∑i≤jλi|ci|2∑i≤j|ci|2, a weighted average of λ1,…,λj with non-negative weights summing to 1. Since the list is indexed decreasingly, the smallest of those eigenvalues is λj, so the average is ≥λj. The statement for Lj is Proposition 17.1.3: Lj is precisely the set of vectors orthogonal to q1,…,qj−1, and the maximum of RA there is λj. In both cases RA(qj)=λj.

Step 1. The witness for the first line. For each subspace W of dimension n−k+1 write f(W)=max0≠x∈WRA(x), which exists by Lemma 17.2.2. The outer minimum of f over all such W is not yet known to exist — there are infinitely many W — so we show that λk is a value of f and a lower bound for all its values; it is then the least value, which is what the minimum means. For the first half take W=Lk, of dimension n−k+1: by the claim, f(Lk)=λk.

Step 2. Every competitor is at least as large. Let W⊆Fn be any subspace with dim⁡W=n−k+1. Then dim⁡W+dim⁡Uk=(n−k+1)+k=n+1>n, so by Lemma 17.2.1 there is a vector x≠0 lying in W∩Uk. Being in Uk, it satisfies RA(x)≥λk by the claim; being in W, it is one of the vectors the inner maximum ranges over. Therefore max0≠y∈WRA(y) ≥ RA(x) ≥ λk. This holds for every competitor W, so f(W)≥λk for all of them. With Step 1, λk is a value of f and a lower bound for every value, so it is the least value: the minimum exists, equals λk, and is attained at W=Lk. This proves the first line.

Step 3. The second line, from the first applied to −A. The matrix −A is Hermitian, its eigenvalues are the numbers −λi, and sorting those decreasingly reverses the list: λj(−A)=−λn+1−j(A)(1≤j≤n). Also R−A(x)=−RA(x) for every x≠0, directly from the definition. Since negation turns a maximum of a set of reals into the minimum of the negated set, and a minimum into a maximum, max0≠x∈WR−A(x)=−min0≠x∈WRA(x),mindim⁡W=m (−f(W))=−maxdim⁡W=mf(W) for any real-valued f on subspaces, where each side exists exactly when the other does. Apply the first line, already proved, to −A with index j: −λn+1−j(A)=λj(−A)=mindim⁡W=n−j+1 max0≠x∈WR−A(x)=−maxdim⁡W=n−j+1 min0≠x∈WRA(x). Now substitute j=n+1−k, so that n−j+1=k and n+1−j=k. Multiplying by −1 gives λk(A)=maxdim⁡W=k min0≠x∈WRA(x), with the maximum attained, the attainment transported from Step 1 by the same negation; explicitly it occurs at W=Uk, which is the subspace Ln+1−k of the eigenbasis of −A taken in its own decreasing order. This proves the theorem.

The theorem removes the eigenvectors from the description of λk(A). What is left mentions only A, through RA, and the subspaces of Fn, which are the same for every matrix. Two Hermitian matrices of the same size can therefore be compared term by term, and that is the whole reason this chapter exists.

Remark. The increasing-index form. If the eigenvalues are listed increasingly as μ1≤μ2≤⋯≤μn, then μj=λn+1−j(A), and substituting k=n+1−j in Theorem 17.2.3 turns it into μj=mindim⁡W=j max0≠x∈WRA(x)=maxdim⁡W=n−j+1 min0≠x∈WRA(x). The pairing of minmax with maxmin is unchanged; only the dimensions move. Many books state the theorem this way, so a reader comparing sources should check which convention is in force before comparing formulas. This book uses the decreasing one everywhere, and this is the only place the other appears.

Warning. Both the order of the two extrema and the dimension attached to each are load-bearing. Neither maxdim⁡W=n−k+1 min0≠x∈WRA(x)normindim⁡W=k max0≠x∈WRA(x) equals λk(A). Each of them is, by Theorem 17.2.3 itself, equal to λn+1−k(A): the first is the second line of the theorem at index n−k+1, and the second is the first line at the index j with n−j+1=k, namely j=n−k+1. So swapping either the order or the dimension returns the eigenvalue counted from the other end of the list.

For a witness take A=diag(2,0) and k=1, so that λ1(A)=2. The first false expression maximizes over the single subspace W=F2 the minimum of RA, which is 0. The second minimizes over lines the maximum on the line, which is the value of RA on that line, and the smallest such value is again 0. Both return λ2(A)=0 instead of 2. The mistake is invisible when n=2k−1, since then n+1−k=k, and this is exactly why it survives in the wild.

Quick check. In the first line of Theorem 17.2.3, why is the dimension n−k+1 rather than k? Check your reason against k=1.

Solution. The competitor W must be forced to meet Uk=span(q1,…,qk), which has dimension k; by Lemma 17.2.1 that needs dim⁡W>n−k, and n−k+1 is the smallest such dimension. It must not be larger either: competitors of dimension n−k+2 would be forced to meet Uk−1 as well, and by the theorem itself the min–max over them is λk−1(A), the wrong eigenvalue.

At k=1 the dimension is n−1+1=n, so the only competitor is W=Fn, and the statement reads λ1(A)=maxx≠0RA(x), which is Proposition 17.1.1. With the dimension k=1 instead, we would be computing min over lines of the single value of RA on the line, which is λn(A) — the wrong end.

Example (Both forms on one matrix)

Let A=(220232024)∈M3(R), whose eigenvalues are 6,3,0, with orthonormal eigenvectors q1=13(1,2,2),q2=13(2,1,−2),q3=13(2,−2,1). Verify both lines of Theorem 17.2.3 for k=2 by computing the witness and two competitors in each case.

Solution. Throughout we use Lemma 17.2.2: for a subspace W with orthonormal basis the columns of Q, the values of RA on W fill the interval between the smallest and largest eigenvalues of the 2×2 matrix B=Q⊤AQ.

Max–min form, dim⁡W=k=2. The witness is U2=span(q1,q2); since q1,q2 are eigenvectors, B=diag(6,3) and the minimum of RA on U2 is 3=λ2(A). Two competitors:

  • W=span(e1,e2), where B is the top-left 2×2 block (2223), with eigenvalues 12(5±17). Its minimum is 12(5−17)≈0.438, well below 3.
  • W=span(e1,e3), where B=diag(2,4) because a13=0. Its minimum is 2, again below 3.

Min–max form, dim⁡W=n−k+1=2 as well, since n=3. The witness is L2=span(q2,q3), where B=diag(3,0) and the maximum is 3=λ2(A). The same two competitors now go the other way:

  • on span(e1,e2) the maximum is 12(5+17)≈4.562>3;
  • on span(e1,e3) the maximum is 4>3.

So each competitor overshoots in the min–max form and undershoots in the max–min form, and the two witnesses pin λ2(A)=3 from both sides. Notice that the competitor span(e1,e2) is the coordinate plane obtained by deleting the third row and column, and its two eigenvalues 0.438… and 4.562… are caught between consecutive eigenvalues 0≤0.438≤3≤4.562≤6 of A. That is no accident; it is the interlacing theorem of §04.

17.2.4. The two debts this pays at once

Corollary 17.2.4 (The Loewner Order Compares Every Eigenvalue)

Let A,B∈Mn(F) be Hermitian with A⪰B. Then λi(A) ≥ λi(B)for every i=1,…,n.

Idea. A⪰B says the Rayleigh quotient of A sits above that of B at every vector. Courant–Fischer builds each λi out of the Rayleigh quotient by a minimum and then a maximum, and both operations preserve “everywhere above”. So the domination passes straight through to every eigenvalue.

Proof. By Definition 13.5.1, A⪰B means A−B⪰0, so ⟨(A−B)x,x⟩≥0 for every x. Dividing by ⟨x,x⟩>0 for x≠0, RA(x) ≥ RB(x)for every x≠0.

Fix i and let W be any subspace of Fn with dim⁡W=i. Let x1∈W∖{0} attain min0≠x∈WRA(x), which exists by Lemma 17.2.2. Then min0≠x∈WRA(x)=RA(x1) ≥ RB(x1) ≥ min0≠x∈WRB(x). So the function W↦min0≠x∈WRA(x) dominates the corresponding function for B at every W, and therefore so do their maxima over the subspaces of dimension i. By Theorem 17.2.3 applied to A and to B, those maxima are λi(A) and λi(B). This proves the corollary.

Example (Checking the corollary on a pair)

Let A=(4113) and B=(2112). Then A−B=diag(2,1), which is positive semidefinite, so A⪰B. The eigenvalues of B are 3 and 1, from trace 4 and determinant 3; those of A are 7±52≈4.618, 2.382, from trace 7 and determinant 11. So λ1(A)≈4.618≥3=λ1(B) and λ2(A)≈2.382≥1=λ2(B), as the corollary demands. Note that the two gaps, 1+52≈1.618 and 5−52≈1.382, are different: the corollary orders the eigenvalues, and says nothing about moving them by the same amount.

This settles an account opened in Chapter 13 §05. Discussing Theorem 13.5.3 (b), which proved only λ1(A)≥λ1(B) and λn(A)≥λn(B), that section said: “The full statement is that λi(A)≥λi(B) for every i, not only for i=1 and i=n; but the intermediate eigenvalues have no description as a plain maximum over the whole unit sphere, and the min–max description that does the job is the Courant–Fischer theorem of Chapter 17.” Corollary 17.2.4 is that full statement, and the diagnosis was exact: what was missing was not an extra trick but a description of λi that two different matrices could be measured against. The proof above is short, and all of its weight is carried by Theorem 17.2.3.

Corollary 17.2.5 (Courant–Fischer for a Self-Adjoint Operator)

Let V be a finite-dimensional inner product space over F=R or F=C, with dim⁡V=n≥1, and let T∈L(V) be self-adjoint, with eigenvalues λ1(T)≥⋯≥λn(T) repeated according to multiplicity. Then for every k with 1≤k≤n, λk(T)=mindim⁡U=n−k+1 max0≠v∈U⟨Tv,v⟩⟨v,v⟩=maxdim⁡U=k min0≠v∈U⟨Tv,v⟩⟨v,v⟩, where U runs over the subspaces of V of the stated dimension.

Idea. Transport the theorem along coordinates. An orthonormal basis of eigenvectors turns V into Fn without changing any inner product, turns T into a real diagonal matrix, and matches subspaces with subspaces of the same dimension. There is nothing to prove beyond checking that each of those three things is true.

Proof. By Proposition 11.6.7 (c) a self-adjoint operator is normal, so Theorem 12.4.2 over C, and Theorem 12.5.1 over R, give an orthonormal basis E=(e1,…,en) of V consisting of eigenvectors of T; relabel it so that Tei=λi(T)ei with the eigenvalues in decreasing order, which is possible because they are real by Theorem 12.2.1. In particular the eigenvalue list of T, repeated according to multiplicity, has exactly n entries.

Let Φ:V→Fn, Φ(v)=[v]E, be the coordinate map, an isomorphism by Corollary 2.5.5. By Theorem 11.2.1 (b), ⟨Φu,Φv⟩=⟨u,v⟩ for all u,v∈V, so Φ is an isometric isomorphism in the sense of Corollary 11.2.4. Put A=[T]EE=diag(λ1(T),…,λn(T)), a real diagonal matrix, hence Hermitian, with λi(A)=λi(T) for every i. By Theorem 4.1.1, Φ(Tv)=AΦ(v).

Consequently, for v≠0, ⟨Tv,v⟩⟨v,v⟩=⟨AΦv,Φv⟩⟨Φv,Φv⟩=RA(Φv), and Φv≠0 because Φ is injective. Finally U↦Φ(U) is a bijection from the subspaces of V to the subspaces of Fn preserving dimension, by Theorem 2.5.2, and it matches the non-zero vectors of U with those of Φ(U). So each of the two displayed expressions equals the corresponding expression for A in Theorem 17.2.3, which is λk(A)=λk(T). This proves the corollary.

This is the statement Chapter 12 §02 promised when it wrote that “the maximization picture returns in Chapter 17, where the Rayleigh quotient and the Courant–Fischer theorem describe every eigenvalue of a self-adjoint operator by optimization”. The picture has returned and the description is complete: every eigenvalue, not merely the extreme two, and stated for an operator on an abstract inner product space, which is the setting Chapter 12 worked in. The remark in §01 records what it cost — nothing. Chapter 12 declined the maximization route because it would have needed the extreme value theorem to know that a maximum exists; here the spectral theorem supplies the maximizing vectors explicitly, so the min–max description is a consequence of Chapter 12 rather than an alternative to it, and no analysis has been imported.

What has been bought is comparison. Two Hermitian matrices are now comparable eigenvalue by eigenvalue whenever their Rayleigh quotients are comparable vector by vector, and Corollary 17.2.4 is the first instance. Section 3 takes the same idea and applies it to A and A+E, where the quotients differ by at most ‖E‖2, and gets the perturbation bound that Chapter 16 could not state.

17.2.5. Exercises

A. Check your understanding

Exercise (A1)
  • (a) State both lines of Theorem 17.2.3, being careful about which dimension goes with which line.
  • (b) State Lemma 17.2.1, and say at which point of the proof of Theorem 17.2.3 it is used and what it is used on.
  • (c) Determine whether the following is correct, and justify your answer: λk(A)=mindim⁡W=k max0≠x∈WRA(x).
  • (d) Write the first line of Theorem 17.2.3 in the increasing-index convention μ1≤⋯≤μn.
  • (e) Explain in one or two sentences why Theorem 13.5.3 (b) could not be extended to intermediate eigenvalues by the argument Chapter 13 used, and what Theorem 17.2.3 supplies instead.

Solution.

  • (a) For Hermitian A∈Mn(F) and 1≤k≤n: λk(A)=minmax over subspaces of dimension n−k+1 of the Rayleigh quotient, and =maxmin over subspaces of dimension k. The large dimension goes with the minimum of maxima; the small one with the maximum of minima.

  • (b) If U,W are subspaces of an n-dimensional space with dim⁡U+dim⁡W>n, then U∩W≠{0}. It is used in Step 2, applied to an arbitrary competitor W of dimension n−k+1 and the span Uk of the top k eigenvectors, whose dimensions add to n+1.

  • (c) Incorrect. By Theorem 17.2.3 itself, that expression is the first line evaluated at the index j with n−j+1=k, so it equals λn−k+1(A). For A=diag(2,0) and k=1 it returns 0, not λ1(A)=2.

  • (d) μj=mindim⁡W=jmax0≠x∈WRA(x), obtained by setting k=n+1−j.

  • (e) Chapter 13’s argument evaluated the quadratic form at a vector that is extremal for the other matrix, which works only when the extremal vector is characterized by a plain maximum over the whole sphere — true for λ1 and λn and for no other index, since the range of RA is an interval (Proposition 17.1.1). Theorem 17.2.3 supplies a description of every λi as an optimization over a family of subspaces that does not depend on the matrix, so the same subspace can be fed to both matrices; that is exactly what the proof of Corollary 17.2.4 does.

B. Practice

Exercise (B1: Both forms on a diagonal matrix)

Let A=diag(5,2,−1)∈M3(R) and k=2.

  • (a) Give the witness subspace for the max–min form, compute the minimum of RA on it, and check the value against Theorem 17.2.3.
  • (b) Do the same for the min–max form.
  • (c) Put u=12(0,1,1). Compute the minimum of RA over the non-zero vectors of span(e1,u), and over those of span(e2,e3), and verify that neither subspace beats the witness in (a).

Solution.

  • (a) The eigenvalues are λ1=5, λ2=2, λ3=−1 with eigenvectors e1,e2,e3. The witness is U2=span(e1,e2), of dimension k=2. There B=diag(5,2), so by Lemma 17.2.2 the values of RA fill [2,5] and the minimum is 2=λ2(A).

  • (b) Here n−k+1=2 as well, and the witness is L2=span(e2,e3), where B=diag(2,−1); the values fill [−1,2] and the maximum is 2=λ2(A).

  • (c) For W=span(e1,u) the two spanning vectors are orthonormal, and B=(50012(2−1))=diag(5,12), the off-diagonal entry vanishing because Ae1=5e1 is orthogonal to (0,1,1). So the minimum is 12<2. For W=span(e2,e3) we get B=diag(2,−1) and minimum −1<2. Neither competitor reaches λ2(A)=2, as Theorem 17.2.3 requires.

Exercise (B2: Shifting by a multiple of the identity)

Let A∈Mn(F) be Hermitian and let c∈R. Prove that λk(A+cI)=λk(A)+c for every k, using the first line of Theorem 17.2.3 for both matrices.

Solution. The matrix A+cI is Hermitian, since c is real. For every x≠0, RA+cI(x)=⟨Ax,x⟩+c⟨x,x⟩⟨x,x⟩=RA(x)+c. Adding a constant to a real-valued function adds it to every maximum and every minimum of that function over any non-empty set. So for each subspace W, max0≠x∈WRA+cI(x)=(max0≠x∈WRA(x))+c, and taking the minimum over all W of dimension n−k+1 adds c again. By Theorem 17.2.3 applied to A+cI on the left and to A on the right, λk(A+cI)=λk(A)+c.

Exercise (B3: Comparing two matrices)

Let A=(3113),B=(2001). Verify that A⪰B, compute both eigenvalue lists, and check Corollary 17.2.4. Check whether λ1(A)−λ1(B)≤λ1(A−B) holds for this pair. (Whether it holds for every pair is a question for Section 3.)

Solution. A−B=(1112) is symmetric with leading principal minors 1>0 and det=2−1=1>0, so A−B≻0 by Theorem 13.1.5; in particular A⪰B.

The eigenvalues: A has pA(x)=(x−3)2−1, so λ1(A)=4, λ2(A)=2; and B is diagonal with λ1(B)=2, λ2(B)=1. Indeed 4≥2 and 2≥1, as Corollary 17.2.4 requires.

For the last question, A−B has trace 3 and determinant 1, so its eigenvalues are 12(3±5) and λ1(A−B)=12(3+5)≈2.618. The inequality asks whether 4−2≤2.618, which holds. (It is not an accident: it is the simplest of the perturbation inequalities of §03, which carry Weyl’s name. Note that Corollary 17.2.4 alone gives no upper bound on the gap λi(A)−λi(B).)

C. Going deeper

Exercise (C1: The definite generalized eigenvalue problem)

Let A,B∈Mn(F) be Hermitian with A≻0, and recall Definition 13.4.1: λ is a generalized eigenvalue of the pair (A,B) when Bx=λAx for some x≠0. Define the generalized Rayleigh quotient RA,B(x)=⟨Bx,x⟩⟨Ax,x⟩(x≠0).

  • (a) Prove that RA,B is well defined and real-valued.
  • (b) Put C=A−1/2BA−1/2. Prove that C is Hermitian and that det(B−λA)=det(A)det(C−λI), so that the generalized eigenvalues of (A,B), with multiplicity, are exactly the eigenvalues of C. Write them d1≥⋯≥dn.
  • (c) Prove that for every k, dk=mindim⁡W=n−k+1 max0≠x∈WRA,B(x).

Hint for (c): the substitution x=A−1/2y.

Solution.

  • (a) By Theorem 13.2.1 and Theorem 13.1.5, A≻0 has a positive definite square root A1/2, which is invertible; write A−1/2 for its inverse, again Hermitian and positive definite. For x≠0, ⟨Ax,x⟩>0 by definiteness, so the denominator is a non-zero real; the numerator is real by Proposition 11.6.7 (a). So the quotient is a well-defined real number.

  • (b) C∗=A−1/2B∗A−1/2=C. Since A1/2CA1/2=B and A1/2A1/2=A, B−λA=A1/2(C−λI)A1/2, so taking determinants and using Theorem 7.5.4 gives det(B−λA)=det(A1/2)2det(C−λI)=det(A)det(C−λI). Since detA≠0, the two polynomials in λ have the same roots with the same multiplicities. By Theorem 13.4.2 (b) the generalized eigenvalues of (A,B) are the roots of det(B−λA), so they are the eigenvalues of C, all real because C is Hermitian.

  • (c) Put S=A−1/2, which is Hermitian and invertible. For y≠0 set x=Sy≠0. Then ⟨Bx,x⟩=⟨SBSy,y⟩=⟨Cy,y⟩,⟨Ax,x⟩=⟨SASy,y⟩=⟨y,y⟩, using SAS=A−1/2AA−1/2=I. Hence RA,B(Sy)=RC(y). The map y↦Sy is an isomorphism of Fn, so W↦S−1W is a bijection of the subspaces of dimension n−k+1 onto themselves, matching non-zero vectors with non-zero vectors. Therefore max0≠x∈WRA,B(x)=max0≠y∈S−1WRC(y), and minimizing both sides over all W of dimension n−k+1 gives, by Theorem 17.2.3 applied to C, the value λk(C)=dk.

Exercise (C2: How much weaker the eigenvalue comparison is)
  • (a) Give Hermitian A,B∈M2(R) with λi(A)≥λi(B) for i=1,2 but A⪰̸B. So the converse of Corollary 17.2.4 is false.
  • (b) Prove that if A and B are Hermitian with λi(A)≥λi(B) for every i and trA=trB, then λi(A)=λi(B) for every i.
  • (c) Prove that if A⪰B and trA=trB, then A=B.

Solution.

  • (a) Take A=diag(2,0) and B=diag(0,1). Then λ1(A)=2≥1=λ1(B) and λ2(A)=0≥0=λ2(B). But A−B=diag(2,−1) has the eigenvalue −1<0, so A−B⪰̸0 by Theorem 13.1.4 and A⪰̸B. Ordering the two spectra says nothing about how the eigenvectors are placed relative to each other, and the Loewner order sees both.

  • (b) By Theorem 9.3.3 the trace is the sum of the eigenvalues with multiplicity, so 0=trA−trB=∑i=1n(λi(A)−λi(B)). Every summand is ≥0 by hypothesis, and a sum of non-negative reals is zero only if each is zero. Hence λi(A)=λi(B) for every i.

  • (c) Put C=A−B, Hermitian with C⪰0, so every eigenvalue of C is ≥0 by Theorem 13.1.4. Also trC=trA−trB=0, which by Theorem 9.3.3 is the sum of those eigenvalues; a sum of non-negative reals vanishing forces each to vanish, so every eigenvalue of C is 0. By Corollary 12.4.3 (or Corollary 12.5.2), C=UDU∗ with D the diagonal matrix of eigenvalues, which is 0; hence C=0 and A=B. (One may also finish with (b) and Corollary 17.2.4, which give the equality of the two spectra but not yet of the matrices; the trace argument on C is what closes the gap.)

Exercise (C3: Does enlarging the field change the answer?)

Let A∈Mn(R) be symmetric. It may be read as an element of Mn(C), where it is Hermitian, and then Theorem 17.2.3 minimizes over a much larger collection of subspaces.

  • (a) Prove that the eigenvalue lists of A over R and over C, with multiplicity, are the same.
  • (b) Deduce that for every k, minW⊆Rndim⁡W=n−k+1max0≠x∈WRA(x)=minW⊆Cndim⁡W=n−k+1max0≠x∈WRA(x).
  • (c) Show that the second minimum really does range over strictly more subspaces, by exhibiting a one-dimensional subspace of C2 that contains no non-zero real vector.

Solution.

  • (a) By Corollary 12.5.2 there is an orthogonal Q∈Mn(R) and a real diagonal D with A=QDQ⊤, and the diagonal of D is the eigenvalue list of A over R with multiplicity. Read over C, the same Q is unitary, since Q∗=Q―⊤=Q⊤=Q−1, and the same equation A=QDQ∗ holds. So A is unitarily similar over C to the same D, and its complex eigenvalue list with multiplicity is again the diagonal of D.

  • (b) Apply Theorem 17.2.3 twice: once with F=R, where the left-hand side equals λk(A) computed over R, and once with F=C, where the right-hand side equals λk(A) computed over C. By (a) these are the same number.

  • (c) Take W=spanC{(1,i)}⊆C2. A non-zero element of W is c(1,i)=(c,ci) with c≠0; for it to be real we would need c∈R and ci∈R, and c≠0 makes the second impossible. So W is not the complex span of any real subspace.

The lesson is worth stating. Enlarging the field adds competitors to the minimization, which can only push the minimum down, and adds vectors to each inner maximization, which can only push those maxima up. The two effects cancel exactly, because both the witness subspace and the bound of Step 2 in the proof of Theorem 17.2.3 survive the enlargement unchanged: the witness is spanned by real eigenvectors, and the dimension count of Lemma 17.2.1 holds over any field.

Found a typo or a mistake? Report it, or read this page’s source

The text of this book is licensed under CC BY-NC-SA 4.0.