Books
Theme
Text size
Width

[ and ] go to the previous and next page.

Not every operator can be diagonalized, even over ℂ. The shift 𝑁=(0100) has the single eigenvalue 0 and ker𝑁=span(𝐞1), so ℂ2 has no basis of eigenvectors of 𝑁 (Section 7.4). But 𝑁 is already upper triangular, and an upper-triangular matrix shows most of what a diagonal one does: its eigenvalues, its determinant, and its trace can all be read off the diagonal. This section shows that over ℂ every operator has an upper-triangular matrix (Theorem 7.5.13), and that on an inner product space the basis can be taken orthonormal (Schur’s theorem, Theorem 7.5.18). The cost of settling for triangular is the entries above the diagonal; in exchange, the result holds for every operator.

Throughout, 𝑉 is a finite-dimensional vector space over a field 𝐹 with dim𝑉=𝑛≥1, and 𝑇∈L(𝑉). We abbreviate [𝑇]ℬ︀ℬ︀ to [𝑇]ℬ︀.

7.5.1Upper-triangular matrices and invariant flags

Recall (Definition 3.3.31) that 𝐴∈𝐹𝑛×𝑛 is upper triangular if 𝑎𝑖𝑗=0 whenever 𝑖>𝑗. For the matrix of an operator this condition has a clean meaning. Column 𝑗 of [𝑇]ℬ︀ holds the coordinates of 𝑇𝐛𝑗, so the zeros below the diagonal in column 𝑗 say exactly that 𝑇𝐛𝑗 is a combination of 𝐛1,…,𝐛𝑗 alone: each basis vector is sent into the span of itself and the vectors before it.

Definition 7.5.1 (Triangularizable). Let ℬ︀ be a basis of 𝑉. The operator 𝑇 is upper triangular with respect to ℬ︀ if [𝑇]ℬ︀ is upper triangular, and 𝑇 is triangularizable if it is upper triangular with respect to some basis of 𝑉. A matrix 𝐴∈𝐹𝑛×𝑛 is triangularizable if it is similar to an upper-triangular matrix, that is (Theorem 4.5.25), if 𝑇𝐴 is triangularizable.

Proposition 7.5.2 (Upper-triangular matrices and invariant subspaces). Let ℬ︀=𝐛1,…,𝐛𝑛 be a basis of 𝑉, and let 𝑉𝑗=span(𝐛1,…,𝐛𝑗) for 0≤𝑗≤𝑛, so that 𝑉0={𝟎}. The following are equivalent.

  1. (I)[𝑇]ℬ︀ is upper triangular.
  2. (II)𝑇𝐛𝑗∈𝑉𝑗 for each 𝑗=1,…,𝑛.
  3. (III)𝑉𝑗 is 𝑇-invariant for each 𝑗=1,…,𝑛.

Proof.routine

Uses. Definition 2.3.1 (Matrix of a linear map), Definition 1.6.14 (Coordinates), Proposition 7.2.4 (Invariance is checked on a spanning list), Definition 7.2.1 (Invariant subspace)

Techniques. Proving TFAE (Technique 2.7.1), Check it on a basis (Technique 2.2.7)

Hint 1. Column 𝑗 of [𝑇]ℬ︀ is [𝑇𝐛𝑗]ℬ︀. When does a vector lie in 𝑉𝑗, in terms of its coordinates?
Hint 2. For (II) ⇒ (III), 𝑉𝑗 is spanned by 𝐛1,…,𝐛𝑗, and 𝑇𝐛𝑖∈𝑉𝑖⊆𝑉𝑗 for 𝑖≤𝑗. Use Proposition 7.2.4.
Sketch. A vector lies in 𝑉𝑗 exactly when its coordinates after position 𝑗 vanish, so (I) ⇔ (II) column by column. Since 𝑉𝑖⊆𝑉𝑗 for 𝑖≤𝑗, (II) gives 𝑇𝐛𝑖∈𝑉𝑗 for all 𝑖≤𝑗, and invariance follows from Proposition 7.2.4. Conversely, 𝐛𝑗∈𝑉𝑗, so invariance of 𝑉𝑗 gives (II).

We prove (I) ⇔ (II) and (II) ⇔ (III).

(I) ⇔ (II). A vector of 𝑉 lies in 𝑉𝑗 if and only if its ℬ︀-coordinates in positions 𝑗+1,…,𝑛 are 0, by uniqueness of coordinates (Definition 1.6.14). Column 𝑗 of [𝑇]ℬ︀ is [𝑇𝐛𝑗]ℬ︀ (Definition 2.3.1), so 𝑇𝐛𝑗∈𝑉𝑗 if and only if column 𝑗 has zeros in rows 𝑗+1,…,𝑛. These conditions for 𝑗=1,…,𝑛 together say that 𝑎𝑖𝑗=0 whenever 𝑖>𝑗.

(II) ⇒ (III). Fix 𝑗. For 𝑖≤𝑗 we have 𝑇𝐛𝑖∈𝑉𝑖⊆𝑉𝑗. So 𝑇 maps each vector of the spanning list 𝐛1,…,𝐛𝑗 of 𝑉𝑗 into 𝑉𝑗, and 𝑉𝑗 is 𝑇-invariant by Proposition 7.2.4.

(III) ⇒ (II). 𝐛𝑗∈𝑉𝑗, and 𝑉𝑗 is 𝑇-invariant, so 𝑇𝐛𝑗∈𝑉𝑗. □

The subspaces 𝑉0⊊𝑉1⊊⋯⊊𝑉𝑛 in the proposition are nested and grow by one dimension at a time. Such chains deserve a name.

Definition 7.5.3 (Flag). A flag in 𝑉 is a list of subspaces 𝑉0⊆𝑉1⊆⋯⊆𝑉𝑛 with dim𝑉𝑗=𝑗 for each 𝑗 (so 𝑉0={𝟎} and 𝑉𝑛=𝑉). The flag of a basis 𝐛1,…,𝐛𝑛 is given by 𝑉𝑗=span(𝐛1,…,𝐛𝑗). A flag is 𝑇-invariant if every 𝑉𝑗 is 𝑇-invariant.

Every flag is the flag of some basis: choose 𝐛𝑗∈𝑉𝑗 with 𝐛𝑗∉𝑉𝑗−1 for each 𝑗. By induction on 𝑗, the list 𝐛1,…,𝐛𝑗 is linearly independent (Proposition 1.5.14) and lies in 𝑉𝑗, so it is a basis of 𝑉𝑗 (Theorem 4.1.10). In ℝ3, a flag is a line through 𝟎 lying in a plane through 𝟎. So Proposition 7.5.2 says that 𝑇 is triangularizable if and only if 𝑉 has a 𝑇-invariant flag, and that then every basis of that flag triangularizes 𝑇. Only the flag matters, not the particular vectors. In particular, the first vector 𝐛1 of a triangularizing basis spans an invariant line, so it is an eigenvector of 𝑇, with eigenvalue the (1,1) entry (Proposition 7.2.7). A triangularizable operator has at least one eigenvalue.

Technique 7.5.4 (Build a flag of invariant subspaces). To find a basis in which 𝑇 is upper triangular, look for nested invariant subspaces of every dimension instead of for matrix entries (Proposition 7.5.2). The flag can be built from the bottom (first an eigenvector, then a vector sent into the span of it and itself, and so on) or from the top (a large invariant subspace, then invariant subspaces inside it). Any basis adapted to the flag will do, so the vectors can afterwards be changed freely within the flag, for instance orthonormalized.

Reach for it whenever a statement is about triangular matrices of an operator, both to produce them and to use them: “[𝑇]ℬ︀ is upper triangular” is the statement “𝑇 maps span(𝐛1,…,𝐛𝑗) into itself for every 𝑗”, which survives any operation that preserves these spans. The canonical example is the second proof of Theorem 7.5.13. Later uses: Schur’s theorem (Theorem 7.5.18), where Gram–Schmidt keeps the flag, and simultaneous triangularization of commuting operators (Section 7.6).

Example 7.5.5 (Differentiation and shifts).
  1. (a)Let 𝐷∈L(𝐹[𝑥]≤3) be differentiation, and use the basis 1,𝑥,𝑥2,𝑥3. Differentiation lowers degrees, so 𝐷 maps 𝐹[𝑥]≤𝑘 into 𝐹[𝑥]≤𝑘−1⊆𝐹[𝑥]≤𝑘. The flag 𝐹[𝑥]≤0⊆𝐹[𝑥]≤1⊆𝐹[𝑥]≤2⊆𝐹[𝑥]≤3 (plus {𝟎} at the bottom) is 𝐷-invariant, and indeed

    [𝐷]1,𝑥,𝑥2,𝑥3=(0100002000030000).
  2. (b)The shift 𝑁∈L(𝐹𝑛) with 𝑁𝐞1=𝟎 and 𝑁𝐞𝑗=𝐞𝑗−1 for 𝑗≥2 (Example 7.2.20) is upper triangular in the standard basis: the flag of coordinate subspaces span(𝐞1,…,𝐞𝑗) is invariant.
  3. (c)Every upper-triangular 𝐴∈𝐹𝑛×𝑛 preserves the coordinate subspaces span(𝐞1,…,𝐞𝑗), the vectors whose last 𝑛−𝑗 coordinates vanish.
Non-example 7.5.6 (A rotation). Let 𝑅=(0−110)∈ℝ2×2, the rotation by a right angle. If 𝑇𝑅 were upper triangular with respect to a basis 𝐛1,𝐛2 of ℝ2, then 𝐛1 would be an eigenvector. But 𝑅 has no real eigenvalues (Example 7.2.18): it moves every line through 𝟎. So 𝑅 is not triangularizable over ℝ. Over ℂ it is, and even diagonalizable: 𝑅(1,−𝑖)=(𝑖,1)=𝑖(1,−𝑖) and 𝑅(1,𝑖)=−𝑖(1,𝑖), and these two vectors form a basis of ℂ2.
𝐴=(2101),𝐴𝐞1=2𝐞1,𝐴𝐞2=𝐞1+𝐞2
The upper-triangular matrix 𝐴=(2101). The 𝑥-axis is invariant, since 𝐴𝐞1=2𝐞1, and {𝟎}⊆span(𝐞1)⊆ℝ2 is an invariant flag. In the interactive version, drag 𝐴𝐞1 off the axis: the matrix is no longer upper triangular in the standard basis, but it is still triangularizable over ℝ as long as some line through 𝟎 is mapped into itself. A rotation by an angle that is not a multiple of 𝜋 leaves no line in place.
Exercise 7.5.7.routine Let ℬ︀=𝐛1,…,𝐛𝑛 be a basis of 𝑉 and ℬ︀′=𝐛𝑛,…,𝐛1 the same basis in reverse order. Show that [𝑇]ℬ︀ is upper triangular if and only if [𝑇]ℬ︀′ is lower triangular. Write down [𝐷]ℬ︀′ for the differentiation operator of Example 7.5.5(a) and ℬ︀′=𝑥3,𝑥2,𝑥,1.

Solution.

Uses. Definition 2.3.1 (Matrix of a linear map), Proposition 7.5.2 (Upper-triangular matrices and invariant subspaces)

Hint 1. The (𝑖,𝑗) entry of [𝑇]ℬ︀′ is the coefficient of 𝐛𝑛+1−𝑖 in 𝑇𝐛𝑛+1−𝑗. Relate it to an entry of [𝑇]ℬ︀.
Hint 2. Reversing both the rows and the columns sends position (𝑖,𝑗) to (𝑛+1−𝑖,𝑛+1−𝑗), and 𝑖>𝑗 if and only if 𝑛+1−𝑖<𝑛+1−𝑗.
Sketch. The (𝑖,𝑗) entry of [𝑇]ℬ︀′ is the (𝑛+1−𝑖,𝑛+1−𝑗) entry of [𝑇]ℬ︀. Below-diagonal positions of one matrix correspond to above-diagonal positions of the other. For 𝐷, [𝐷]ℬ︀′ has 3,2,1 just below the diagonal.

Write 𝐛𝑗′=𝐛𝑛+1−𝑗. The (𝑖,𝑗) entry of [𝑇]ℬ︀′ is the coefficient of 𝐛𝑖′=𝐛𝑛+1−𝑖 in 𝑇𝐛𝑗′=𝑇𝐛𝑛+1−𝑗 (Definition 2.3.1), which is the (𝑛+1−𝑖,𝑛+1−𝑗) entry of [𝑇]ℬ︀. Now 𝑖<𝑗 if and only if 𝑛+1−𝑖>𝑛+1−𝑗. So the entries of [𝑇]ℬ︀′ above the diagonal are exactly the entries of [𝑇]ℬ︀ below the diagonal, and one set vanishes if and only if the other does.

For 𝐷, 𝐷(𝑥3)=3𝑥2, 𝐷(𝑥2)=2𝑥, 𝐷(𝑥)=1, and 𝐷(1)=𝟎, so

[𝐷]𝑥3,𝑥2,𝑥,1=(0000300002000010).

This is [𝐷]1,𝑥,𝑥2,𝑥3 with the order of the rows and of the columns reversed.

Reading the diagonal

The diagonal of a triangular matrix carries the eigenvalues. We give two proofs: one computes a determinant, the other sees where a zero on the diagonal loses a dimension.

Proposition 7.5.8 (The diagonal of a triangular matrix). Suppose [𝑇]ℬ︀ is upper triangular with diagonal entries 𝜆1,…,𝜆𝑛.

  1. (a)𝑇 is invertible if and only if 𝜆𝑗≠0 for every 𝑗.
  2. (b)The eigenvalues of 𝑇 are exactly 𝜆1,…,𝜆𝑛.

Proof (via determinants).routine

Uses. Definition 5.4.18 (Determinant of an operator), Proposition 5.3.15 (Triangular matrices), Proposition 5.4.20 (Properties of the determinant of an operator), Theorem 7.2.9 (Eigenvalues and invertibility), Proposition 2.4.12 (Matrices of sums and scalar multiples)

Hint 1. What is det𝑇, computed from [𝑇]ℬ︀?
Hint 2. For (b), [𝑇−𝜆id𝑉]ℬ︀=[𝑇]ℬ︀−𝜆𝐼𝑛 is again upper triangular. Apply (a) to 𝑇−𝜆id𝑉.
Sketch. det𝑇=det[𝑇]ℬ︀=𝜆1⋯𝜆𝑛, which is nonzero exactly when every 𝜆𝑗 is. For (b), 𝑇−𝜆id𝑉 has the upper-triangular matrix [𝑇]ℬ︀−𝜆𝐼𝑛 with diagonal 𝜆𝑗−𝜆, and 𝜆 is an eigenvalue exactly when 𝑇−𝜆id𝑉 is not invertible.

(a) By Definition 5.4.18 and Proposition 5.3.15, det𝑇=det[𝑇]ℬ︀=𝜆1𝜆2⋯𝜆𝑛. By Proposition 5.4.20(c), 𝑇 is invertible if and only if this product is nonzero, that is, if and only if every 𝜆𝑗 is nonzero.

(b) Let 𝜆∈𝐹. By Proposition 2.4.12, [𝑇−𝜆id𝑉]ℬ︀=[𝑇]ℬ︀−𝜆𝐼𝑛, which is upper triangular with diagonal entries 𝜆1−𝜆,…,𝜆𝑛−𝜆. By Theorem 7.2.9, 𝜆 is an eigenvalue of 𝑇 if and only if 𝑇−𝜆id𝑉 is not invertible, which by (a) happens if and only if 𝜆𝑗−𝜆=0 for some 𝑗. □

Intuition (Where the dimension goes). Suppose 𝜆𝑗=0. Then 𝑇𝐛𝑗 has no 𝐛𝑗-component, so 𝑇 maps 𝑉𝑗=span(𝐛1,…,𝐛𝑗) into the smaller space 𝑉𝑗−1. A map from a 𝑗-dimensional space into a (𝑗−1)-dimensional one cannot be injective. Conversely, if every 𝜆𝑗≠0, we can solve for 𝐛𝑗 in 𝑇𝐛𝑗=𝜆𝑗𝐛𝑗+(⋯) and see, one index at a time, that every basis vector lies in the image. This argument uses no determinants, and it is the version that survives where determinants are unavailable.

Proof (without determinants).standard

Uses. Proposition 7.5.2 (Upper-triangular matrices and invariant subspaces), Theorem 4.2.11 (Rank–nullity), Corollary 4.2.18 (Injective if and only if surjective), Theorem 7.2.9 (Eigenvalues and invertibility), Proposition 2.4.12 (Matrices of sums and scalar multiples)

Techniques. Dimension counting (Technique 4.2.24), Strong induction and well-ordering (Technique 5.2.12)

Hint 1. Let 𝑉𝑗=span(𝐛1,…,𝐛𝑗). If 𝜆𝑗=0, where does 𝑇 send 𝑉𝑗?
Hint 2. If 𝜆𝑗=0, then 𝑇 maps the 𝑗-dimensional 𝑉𝑗 into 𝑉𝑗−1, so 𝑇↾𝑉𝑗 has a nonzero kernel by rank–nullity.
Hint 3. If all 𝜆𝑗≠0, show 𝐛𝑗∈im𝑇 by strong induction on 𝑗, solving 𝑇𝐛𝑗=𝜆𝑗𝐛𝑗+∑𝑖<𝑗𝑎𝑖𝑗𝐛𝑖 for 𝐛𝑗.
Sketch. If 𝜆𝑗=0, then 𝑇𝐛𝑖∈𝑉𝑖⊆𝑉𝑗−1 for 𝑖<𝑗 and 𝑇𝐛𝑗=∑𝑖<𝑗𝑎𝑖𝑗𝐛𝑖∈𝑉𝑗−1, so 𝑇↾𝑉𝑗:𝑉𝑗→𝑉𝑗−1 is not injective. If all 𝜆𝑗≠0, then 𝐛𝑗=𝜆𝑗−1(𝑇𝐛𝑗−∑𝑖<𝑗𝑎𝑖𝑗𝐛𝑖)∈im𝑇 by induction, so 𝑇 is surjective, hence invertible. (b) follows from (a) applied to 𝑇−𝜆id𝑉.

Write [𝑇]ℬ︀=(𝑎𝑖𝑗), so 𝑎𝑗𝑗=𝜆𝑗, and 𝑉𝑗=span(𝐛1,…,𝐛𝑗). By Proposition 7.5.2, 𝑇𝐛𝑗=𝜆𝑗𝐛𝑗+∑𝑖<𝑗𝑎𝑖𝑗𝐛𝑖 for each 𝑗.

(a) Suppose 𝜆𝑗=0 for some 𝑗. For 𝑖<𝑗, 𝑇𝐛𝑖∈𝑉𝑖⊆𝑉𝑗−1, and 𝑇𝐛𝑗=∑𝑖<𝑗𝑎𝑖𝑗𝐛𝑖∈𝑉𝑗−1. So 𝑇 maps the spanning list 𝐛1,…,𝐛𝑗 of 𝑉𝑗 into the subspace 𝑉𝑗−1, and im(𝑇↾𝑉𝑗)⊆𝑉𝑗−1. By rank–nullity (Theorem 4.2.11) applied to 𝑇↾𝑉𝑗:𝑉𝑗→𝑉,

dimker(𝑇↾𝑉𝑗)=𝑗−dimim(𝑇↾𝑉𝑗)≥𝑗−(𝑗−1)=1.

So some nonzero vector of 𝑉𝑗 is sent to 𝟎, and 𝑇 is not invertible.

Conversely, suppose every 𝜆𝑗≠0. We show 𝐛𝑗∈im𝑇 by strong induction on 𝑗. Assume 𝐛𝑖∈im𝑇 for all 𝑖<𝑗 (for 𝑗=1 there is nothing to assume). Then

𝐛𝑗=𝜆𝑗−1(𝑇𝐛𝑗−∑𝑖<𝑗𝑎𝑖𝑗𝐛𝑖)

is a linear combination of 𝑇𝐛𝑗 and the 𝐛𝑖 with 𝑖<𝑗, all of which lie in the subspace im𝑇. So 𝐛𝑗∈im𝑇. Hence im𝑇 contains a basis of 𝑉, and 𝑇 is surjective, so it is invertible by Corollary 4.2.18.

(b) As in the first proof, [𝑇−𝜆id𝑉]ℬ︀ is upper triangular with diagonal entries 𝜆𝑗−𝜆, and 𝜆 is an eigenvalue if and only if 𝑇−𝜆id𝑉 is not invertible (Theorem 7.2.9), that is, by (a), if and only if 𝜆=𝜆𝑗 for some 𝑗. □

The diagonal also records how often each eigenvalue occurs. Write 𝑝𝑇 for the characteristic polynomial of 𝑇 (Definition 7.3.4), so that 𝑝𝑇=det(𝑥𝐼𝑛−𝐴) for every matrix 𝐴 of 𝑇. Here the determinant of a matrix with entries in 𝐹[𝑥] is defined by the Leibniz formula, and Section 7.3 shows that the results of Chapter 5 proved by expanding that formula, such as Proposition 5.3.15 and Theorem 5.4.26, hold for such matrices.

Proposition 7.5.9 (Multiplicities on the diagonal). Suppose [𝑇]ℬ︀ is upper triangular with diagonal entries 𝜆1,…,𝜆𝑛. Then

𝑝𝑇=(𝑥−𝜆1)(𝑥−𝜆2)⋯(𝑥−𝜆𝑛).

Consequently each eigenvalue 𝜆 of 𝑇 occurs among 𝜆1,…,𝜆𝑛 exactly as many times as its algebraic multiplicity (Definition 7.4.11). In particular, this number is the same for every basis in which 𝑇 is upper triangular.

Proof.routine

Uses. Definition 7.3.4 (Characteristic polynomial), Proposition 5.3.15 (Triangular matrices), Proposition 7.1.34 (Roots of a split polynomial), Definition 7.4.11 (Algebraic multiplicity)

Hint 1. 𝑥𝐼𝑛−[𝑇]ℬ︀ is upper triangular. What is its determinant?
Hint 2. The algebraic multiplicity of 𝜆 is its multiplicity as a root of 𝑝𝑇. Use Proposition 7.1.34.
Sketch. 𝑥𝐼𝑛−[𝑇]ℬ︀ is upper triangular with diagonal entries 𝑥−𝜆𝑗, so its determinant is their product. A root of a split polynomial occurs in its list of linear factors exactly as often as its multiplicity.
The matrix 𝑥𝐼𝑛−[𝑇]ℬ︀ has entries in 𝐹[𝑥]; it is upper triangular, with diagonal entries 𝑥−𝜆1,…,𝑥−𝜆𝑛. By Proposition 5.3.15 (whose proof only uses the Leibniz formula), 𝑝𝑇=det(𝑥𝐼𝑛−[𝑇]ℬ︀)=(𝑥−𝜆1)⋯(𝑥−𝜆𝑛). By Proposition 7.1.34, each root of 𝑝𝑇 occurs in the list 𝜆1,…,𝜆𝑛 exactly as many times as its multiplicity as a root of 𝑝𝑇, which is its algebraic multiplicity as an eigenvalue (Definition 7.4.11). The polynomial 𝑝𝑇 does not depend on the basis, so neither do these counts. □

When 𝑝𝑇 splits over 𝐹, say 𝑝𝑇=(𝑥−𝜆1)⋯(𝑥−𝜆𝑛) (Definition 7.1.33), we call 𝜆1,…,𝜆𝑛 the eigenvalues of 𝑇 counted with multiplicity. Each eigenvalue occurs in this list as often as its algebraic multiplicity, and the list is determined by 𝑇 up to order (Proposition 7.1.34). Over ℂ, 𝑝𝑇 always splits (Corollary 7.1.35).

Exercise 7.5.10.routine Let 𝐛1,𝐛2,𝐛3 be a basis of a complex vector space 𝑉, and let 𝑇∈L(𝑉) satisfy

𝑇𝐛1=3𝐛1,𝑇𝐛2=𝐛1−𝐛2,𝑇𝐛3=𝐛1+𝐛2+3𝐛3.

Find the eigenvalues of 𝑇 with their algebraic multiplicities, decide whether 𝑇 is invertible, and find an eigenvector for each eigenvalue. Is 𝑇 diagonalizable?

Solution.

Uses. Proposition 7.5.8 (The diagonal of a triangular matrix), Proposition 7.5.9 (Multiplicities on the diagonal), Theorem 7.4.2 (Diagonalizable means an eigenbasis)

Hint 1. Each 𝑇𝐛𝑗 lies in span(𝐛1,…,𝐛𝑗), so [𝑇]ℬ︀ is upper triangular. Read off the diagonal.
Hint 2. For the eigenvectors, solve (𝑇−𝜆id)(𝑥𝐛1+𝑦𝐛2+𝑧𝐛3)=𝟎 coefficient by coefficient.
Sketch. The diagonal is 3,−1,3, so the eigenvalues are 3 (multiplicity 2) and −1 (multiplicity 1), and 𝑇 is invertible. 𝐸3=span(𝐛1) and 𝐸−1=span(𝐛1−4𝐛2). Only two independent eigenvectors exist, so 𝑇 is not diagonalizable.

The matrix

[𝑇]ℬ︀=(3110−11003)

is upper triangular. By Proposition 7.5.8 and Proposition 7.5.9, the eigenvalues are 3, with algebraic multiplicity 2, and −1, with algebraic multiplicity 1. No diagonal entry is 0, so 𝑇 is invertible.

For 𝜆=3: (𝑇−3id)(𝑥𝐛1+𝑦𝐛2+𝑧𝐛3)=𝑦(𝐛1−4𝐛2)+𝑧(𝐛1+𝐛2)=(𝑦+𝑧)𝐛1+(𝑧−4𝑦)𝐛2. This vanishes if and only if 𝑦+𝑧=0 and 𝑧=4𝑦, that is, 𝑦=𝑧=0. So 𝐸3=span(𝐛1).

For 𝜆=−1: (𝑇+id)(𝑥𝐛1+𝑦𝐛2+𝑧𝐛3)=(4𝑥+𝑦+𝑧)𝐛1+𝑧𝐛2+4𝑧𝐛3, which vanishes if and only if 𝑧=0 and 𝑦=−4𝑥. So 𝐸−1=span(𝐛1−4𝐛2). Check: 𝑇(𝐛1−4𝐛2)=3𝐛1−4𝐛1+4𝐛2=−(𝐛1−4𝐛2).

Every eigenvector lies in 𝐸3∪𝐸−1, a union of two lines, so there is no basis of eigenvectors, and 𝑇 is not diagonalizable (Theorem 7.4.2).

7.5.2Triangularization over ℂ

Intuition (How to find an invariant flag). By Proposition 7.5.2 we need an invariant flag, and its bottom 𝑉1 is an invariant line, spanned by an eigenvector. So we need an eigenvalue, and over ℂ there always is one (Theorem 7.3.27). Over ℝ the rotation of Non-example 7.5.6 shows that the theorem can fail at this very first step.

Having found an eigenvector 𝐛1, we would like to finish by induction on dim𝑉. The obvious move is to restrict 𝑇 to a complement 𝑊 of span(𝐛1) and use the induction hypothesis there. It fails: 𝑇 need not map 𝑊 into 𝑊, so 𝑇↾𝑊 is a map 𝑊→𝑉, not an operator on 𝑊, and the induction hypothesis says nothing about it. (Exercise A.2.13 makes exactly this mistake.) The fix is to throw away the part of 𝑇𝐰 that lies in span(𝐛1), by projecting 𝑇𝐰 back into 𝑊 along span(𝐛1). That produces an honest operator on 𝑊. The part thrown away is harmless, because span(𝐛1) sits at the bottom of the flag, and every 𝑉𝑗 will contain it anyway.

Lemma 7.5.11 (Compressing an operator to a complement). Let 𝑈 be a 𝑇-invariant subspace of 𝑉, and let 𝑉=𝑈⊕𝑊. The compression of 𝑇 to 𝑊 is the operator

𝑆=Pr𝑊,𝑈∘𝑇↾𝑊∈L(𝑊),

where Pr𝑊,𝑈 is the projection onto 𝑊 along 𝑈 (Definition 4.4.24).

  1. (a)𝑇𝐰−𝑆𝐰∈𝑈 for every 𝐰∈𝑊.
  2. (b)Let 𝒞︀ be a basis of 𝑈 and 𝒟︀ a basis of 𝑊, both nonempty. Then ℬ︀=𝒞︀,𝒟︀ is a basis of 𝑉, and

    [𝑇]ℬ︀=([𝑇↾𝑈]𝒞︀𝑋0[𝑆]𝒟︀)

    for some matrix 𝑋, where 𝑇↾𝑈 is regarded as an operator on 𝑈.

Proof.routine

Uses. Definition 4.4.24 (Projection along a complement), Proposition 4.4.26 (Properties of projections), Theorem 4.4.13 (Bases and dimension of a sum), Proposition 4.5.35 (Preserved subspaces give block triangular matrices), Definition 2.3.1 (Matrix of a linear map)

Techniques. Choose a basis adapted to a preserved subspace (Technique 4.5.36)

Hint 1. For (a), Pr𝑈,𝑊+Pr𝑊,𝑈=id𝑉. Apply both sides to 𝑇𝐰.
Hint 2. For (b), Proposition 4.5.35 gives the zero block and the top left block. For a column belonging to 𝐝𝑗∈𝒟︀, write 𝑇𝐝𝑗=(𝑇𝐝𝑗−𝑆𝐝𝑗)+𝑆𝐝𝑗 and read off the 𝒟︀-coordinates.
Sketch. 𝑇𝐰=Pr𝑈,𝑊(𝑇𝐰)+Pr𝑊,𝑈(𝑇𝐰), and the first term lies in 𝑈; this is (a). In (b), ℬ︀ is a basis because the sum is direct, Proposition 4.5.35 gives the block form with [𝑇↾𝑈]𝒞︀ in the corner, and by (a) the 𝒟︀-coordinates of 𝑇𝐝𝑗 are those of 𝑆𝐝𝑗, which form column 𝑗 of [𝑆]𝒟︀.

(a) By Proposition 4.4.26, Pr𝑈,𝑊+Pr𝑊,𝑈=id𝑉 and imPr𝑈,𝑊=𝑈. So for 𝐰∈𝑊,

𝑇𝐰−𝑆𝐰=𝑇𝐰−Pr𝑊,𝑈(𝑇𝐰)=Pr𝑈,𝑊(𝑇𝐰)∈𝑈.

(b) Since 𝑉=𝑈⊕𝑊, the list ℬ︀=𝒞︀,𝒟︀ is a basis of 𝑉 (Theorem 4.4.13). Let 𝑘 be the length of 𝒞︀. Since 𝑈=span(𝒞︀) is 𝑇-invariant, Proposition 4.5.35 shows that [𝑇]ℬ︀ has a zero block below the first 𝑘 columns and that its top left block is [𝑇↾𝑈]𝒞︀. Let 𝒟︀=𝐝1,…,𝐝𝑚 and [𝑆]𝒟︀=(𝑠𝑖𝑗), so 𝑆𝐝𝑗=∑𝑖𝑠𝑖𝑗𝐝𝑖 (Definition 2.3.1). By (a), 𝑇𝐝𝑗=𝐮+∑𝑖𝑠𝑖𝑗𝐝𝑖 for some 𝐮∈𝑈=span(𝒞︀). So the last 𝑚 coordinates of 𝑇𝐝𝑗 in ℬ︀ are 𝑠1𝑗,…,𝑠𝑚𝑗, and the bottom right block of [𝑇]ℬ︀ is [𝑆]𝒟︀. □

The compression depends on the choice of 𝑊, but its matrix is always the bottom right block of [𝑇]ℬ︀. This is the move the triangularization theorem needs, and it is an instance of a general technique.

Technique 7.5.12 (Induction on dimension via an invariant subspace). To prove that every operator on every finite-dimensional space (of some kind) has a property, induct on dim𝑉. In the inductive step, find a 𝑇-invariant subspace 𝑈 with {𝟎}≠𝑈≠𝑉, apply the induction hypothesis to an operator on a smaller space built from 𝑈, and glue the pieces with the block triangular form of Proposition 4.5.35. The smaller operator can be

  • the restriction 𝑇↾𝑈, an operator on 𝑈 because 𝑈 is invariant;
  • the compression of 𝑇 to a complement 𝑊 of 𝑈 (Lemma 7.5.11), or to 𝑈⟂ when there is an inner product;
  • the quotient operator 𝑇/𝑈 on 𝑉/𝑈 (Proposition 4.6.32), if quotients are available.

The hypothesis applies only to operators, so whatever is passed down must map the smaller space into itself; a restriction to a subspace that is not invariant does not qualify (Exercise A.2.13). The invariant subspace usually comes from an eigenvector, which is why such inductions tend to need ℂ, or from a kernel or image of a polynomial in 𝑇 (Technique 7.2.31).

Reach for it for existence statements about bases adapted to an operator: triangular, block triangular, orthonormal and triangular, diagonal. The canonical example is Theorem 7.5.13. Later uses: Schur’s theorem (Theorem 7.5.18), the spectral theorems (Section 8.2), and the structure of nilpotent operators (Section 10.3).

Theorem 7.5.13 (Upper-triangular form). Let 𝑉 be a finite-dimensional complex vector space with dim𝑉≥1, and let 𝑇∈L(𝑉). Then there is a basis ℬ︀ of 𝑉 such that [𝑇]ℬ︀ is upper triangular. Its diagonal entries are the eigenvalues of 𝑇, each occurring as many times as its algebraic multiplicity.

Proof (compress to a complement).standard

Uses. Theorem 7.3.27 (Existence of eigenvalues over ℂ), Proposition 7.2.7 (Invariant lines are spanned by eigenvectors), Proposition 4.4.16 (Complements exist), Lemma 7.5.11 (Compressing an operator to a complement), Proposition 7.5.2 (Upper-triangular matrices and invariant subspaces), Theorem 4.4.13 (Bases and dimension of a sum), Proposition 7.5.9 (Multiplicities on the diagonal)

Techniques. Induction on dimension via an invariant subspace (Technique 7.5.12), Build a flag of invariant subspaces (Technique 7.5.4)

Hint 1. Induct on 𝑛=dim𝑉. Over ℂ, 𝑇 has an eigenvector 𝐛1, and span(𝐛1) is invariant. It should be the bottom of the flag.
Hint 2. Choose a complement 𝑊 of span(𝐛1). The restriction of 𝑇 to 𝑊 is not an operator on 𝑊; what is?
Hint 3. Apply the induction hypothesis to the compression 𝑆 of 𝑇 to 𝑊, getting 𝐛2,…,𝐛𝑛. Then use 𝑇𝐰−𝑆𝐰∈span(𝐛1) to check that 𝑇𝐛𝑗∈span(𝐛1,…,𝐛𝑗).
Sketch. Induct on 𝑛; 𝑛=1 is trivial. Take an eigenvector 𝐛1 (Theorem 7.3.27), a complement 𝑊 of 𝑈=span(𝐛1), and the compression 𝑆∈L(𝑊). By induction, 𝑊 has a basis 𝐛2,…,𝐛𝑛 with 𝑆𝐛𝑗∈span(𝐛2,…,𝐛𝑗). Then 𝑇𝐛1=𝜆1𝐛1 and 𝑇𝐛𝑗=(𝑇𝐛𝑗−𝑆𝐛𝑗)+𝑆𝐛𝑗∈span(𝐛1,…,𝐛𝑗), so 𝐛1,…,𝐛𝑛 triangularizes 𝑇. The diagonal is described by Proposition 7.5.9.

The statement about the diagonal is Proposition 7.5.9, so we prove existence, by induction on 𝑛=dim𝑉. If 𝑛=1, every 1×1 matrix is upper triangular, and any basis works.

Let 𝑛≥2, and assume the theorem for complex vector spaces of dimension 𝑛−1. By Theorem 7.3.27, 𝑇 has an eigenvalue 𝜆1; let 𝐛1 be an eigenvector for it. Then 𝑈=span(𝐛1) is 𝑇-invariant (Proposition 7.2.7). By Proposition 4.4.16, 𝑈 has a complement 𝑊, and dim𝑊=𝑛−1. Let 𝑆∈L(𝑊) be the compression of 𝑇 to 𝑊 (Lemma 7.5.11). By the induction hypothesis and Proposition 7.5.2, 𝑊 has a basis 𝐛2,…,𝐛𝑛 with

𝑆𝐛𝑗∈span(𝐛2,…,𝐛𝑗)for𝑗=2,…,𝑛.

The list ℬ︀=𝐛1,𝐛2,…,𝐛𝑛 is a basis of 𝑉=𝑈⊕𝑊 (Theorem 4.4.13). We check condition (II) of Proposition 7.5.2. First, 𝑇𝐛1=𝜆1𝐛1∈span(𝐛1). For 𝑗≥2, by Lemma 7.5.11(a),

𝑇𝐛𝑗=(𝑇𝐛𝑗−𝑆𝐛𝑗)⏟∈span(𝐛1)+𝑆𝐛𝑗⏟∈span(𝐛2,…,𝐛𝑗)∈span(𝐛1,…,𝐛𝑗).

So [𝑇]ℬ︀ is upper triangular. In block form (Lemma 7.5.11(b)), [𝑇]ℬ︀=(𝜆1𝑋0[𝑆]𝐛2,…,𝐛𝑛). □

Intuition (A second route, from the top). The first proof builds the flag from the bottom: an invariant line first, then whatever the smaller problem returns. The second proof (LADR’s) works from the top. Fix an eigenvalue 𝜆, and consider 𝑈=im(𝑇−𝜆id𝑉). Since 𝑇𝐯=𝜆𝐯+(𝑇−𝜆id𝑉)𝐯, every subspace that contains 𝑈 is 𝑇-invariant: 𝑇 moves each vector only within its own span plus 𝑈. And 𝑈≠𝑉, because 𝑇−𝜆id𝑉 is not injective. So above 𝑈 the flag can be chosen arbitrarily, and the only real work happens inside 𝑈, where the restriction of 𝑇 is an honest operator on a smaller space. This route needs no complement and no projection, only restriction.

Proof (via the invariant subspace im(𝑇−𝜆id𝑉)).standard

Uses. Theorem 7.3.27 (Existence of eigenvalues over ℂ), Theorem 7.2.9 (Eigenvalues and invertibility), Corollary 4.2.18 (Injective if and only if surjective), Proposition 7.2.30 (Kernels and images of commuting operators are invariant), Theorem 4.1.16 (Extending to a basis), Proposition 7.5.2 (Upper-triangular matrices and invariant subspaces)

Techniques. Induction on dimension via an invariant subspace (Technique 7.5.12), Build a flag of invariant subspaces (Technique 7.5.4), Strong induction and well-ordering (Technique 5.2.12)

Hint 1. Let 𝜆 be an eigenvalue and 𝑈=im(𝑇−𝜆id𝑉). Show that 𝑈 is 𝑇-invariant and dim𝑈<𝑛.
Hint 2. Apply the induction hypothesis to 𝑇↾𝑈 (strong induction, since dim𝑈 can be anything below 𝑛), and extend the resulting basis of 𝑈 to a basis of 𝑉.
Hint 3. For an added vector 𝐯, write 𝑇𝐯=𝜆𝐯+(𝑇−𝜆id𝑉)𝐯.
Sketch. By strong induction on 𝑛. Let 𝜆 be an eigenvalue and 𝑈=im(𝑇−𝜆id𝑉); it is 𝑇-invariant because 𝑇 commutes with 𝑇−𝜆id𝑉, and 𝑈≠𝑉 because 𝑇−𝜆id𝑉 is not injective, hence not surjective. Triangularize 𝑇↾𝑈 by 𝐮1,…,𝐮𝑚 and extend by 𝐯1,…,𝐯𝑘. Then 𝑇𝐯𝑗=𝜆𝐯𝑗+(𝑇−𝜆id𝑉)𝐯𝑗∈span(𝐮1,…,𝐮𝑚,𝐯𝑗).

We use strong induction on 𝑛=dim𝑉. For 𝑛=1 there is nothing to prove. Let 𝑛≥2, and assume the theorem for all complex spaces of dimension less than 𝑛.

By Theorem 7.3.27, 𝑇 has an eigenvalue 𝜆. Let 𝑈=im(𝑇−𝜆id𝑉). Since 𝑇 commutes with 𝑇−𝜆id𝑉, the image 𝑈 is 𝑇-invariant (Proposition 7.2.30). Since 𝜆 is an eigenvalue, 𝑇−𝜆id𝑉 is not injective (Theorem 7.2.9), hence not surjective (Corollary 4.2.18), so 𝑚=dim𝑈<𝑛.

If 𝑚=0, then 𝑇=𝜆id𝑉, and every basis triangularizes 𝑇. Otherwise 𝑇↾𝑈 is an operator on 𝑈, and by the induction hypothesis and Proposition 7.5.2, 𝑈 has a basis 𝐮1,…,𝐮𝑚 with 𝑇𝐮𝑗∈span(𝐮1,…,𝐮𝑗) for each 𝑗. Extend it to a basis 𝐮1,…,𝐮𝑚,𝐯1,…,𝐯𝑘 of 𝑉 (Theorem 4.1.16). For each 𝑗≤𝑘,

𝑇𝐯𝑗=𝜆𝐯𝑗+(𝑇−𝜆id𝑉)𝐯𝑗⏟∈𝑈∈span(𝐮1,…,𝐮𝑚,𝐯1,…,𝐯𝑗).

So every vector of the basis is sent into the span of itself and the vectors before it, and 𝑇 is upper triangular with respect to this basis by Proposition 7.5.2. □

The two proofs pass different objects down the induction. The first passes a compression to a space of dimension exactly 𝑛−1, and so needs a complement and a projection; the second passes a restriction to a space of unknown smaller dimension, and so needs strong induction. The second proof also says something about the diagonal: the last 𝑘 basis vectors have diagonal entry 𝜆, and 𝑘=𝑛−dimim(𝑇−𝜆id𝑉)=dimker(𝑇−𝜆id𝑉) by rank–nullity. With Proposition 7.5.9, the algebraic multiplicity of 𝜆 is at least dim𝐸𝜆, its geometric multiplicity. So this argument recovers Theorem 7.4.15 for complex operators.

Corollary 7.5.14 (Triangularizing a matrix). Every 𝐴∈ℂ𝑛×𝑛 is similar to an upper-triangular matrix: 𝐴=𝑆𝑅𝑆−1 with 𝑆∈ℂ𝑛×𝑛 invertible and 𝑅 upper triangular. The columns of 𝑆 form a basis of ℂ𝑛 whose flag is 𝑇𝐴-invariant, and the diagonal of 𝑅 lists the eigenvalues of 𝐴 counted with multiplicity.

Proof.routine

Uses. Theorem 7.5.13 (Upper-triangular form), Theorem 4.5.17 (Change of basis)

Hint. Apply Theorem 7.5.13 to 𝑇𝐴 and translate the change of basis into matrices.
Sketch. If ℬ︀ triangularizes 𝑇𝐴 and 𝑆 has the vectors of ℬ︀ as columns, then 𝑆=[id]ℰ︀ℬ︀ and [𝑇𝐴]ℬ︀=𝑆−1𝐴𝑆 by Theorem 4.5.17.
By Theorem 7.5.13 there is a basis ℬ︀=𝐛1,…,𝐛𝑛 of ℂ𝑛 with 𝑅=[𝑇𝐴]ℬ︀ upper triangular. Let 𝑆 be the matrix with columns 𝐛1,…,𝐛𝑛, that is, 𝑆=[id]ℰ︀ℬ︀ for the standard basis ℰ︀. It is invertible, and Theorem 4.5.17 gives 𝑅=𝑆−1[𝑇𝐴]ℰ︀𝑆=𝑆−1𝐴𝑆, so 𝐴=𝑆𝑅𝑆−1. The rest is Theorem 7.5.13. □

Example 7.5.15 (Triangularizing a 3×3 matrix). Let

𝐴=(20−102−2−110)∈ℂ3×3.

We follow the first proof of Theorem 7.5.13. Its characteristic polynomial is (𝑥−2)(𝑥−1)2, and 𝐛1=(1,1,0) is an eigenvector: 𝐴𝐛1=(2,2,0)=2𝐛1. A convenient complement of span(𝐛1) is 𝑊=span(𝐞2,𝐞3), since 𝐛1 has a nonzero first coordinate. Every vector decomposes as

(𝑥1,𝑥2,𝑥3)=𝑥1𝐛1+(0,𝑥2−𝑥1,𝑥3),

so the projection onto 𝑊 along span(𝐛1) is (𝑥1,𝑥2,𝑥3)↦(0,𝑥2−𝑥1,𝑥3). The compression 𝑆 sends

𝐞2↦Pr(𝐴𝐞2)=Pr(0,2,1)=2𝐞2+𝐞3,𝐞3↦Pr(𝐴𝐞3)=Pr(−1,−2,0)=−𝐞2,

so [𝑆]𝐞2,𝐞3=(2−110). This 2×2 problem has the eigenvalue 1 with eigenvector (1,1) in these coordinates, that is, 𝐛2=𝐞2+𝐞3, and we complete it with 𝐛3=𝐞3. Now

𝐴𝐛1=2𝐛1,𝐴𝐛2=(−1,0,1)=−𝐛1+𝐛2,𝐴𝐛3=(−1,−2,0)=−𝐛1−𝐛2+𝐛3,

so with 𝑃=(100110011), whose columns are 𝐛1,𝐛2,𝐛3,

𝑃−1𝐴𝑃=(2−1−101−1001).

The diagonal 2,1,1 matches 𝑝𝐴=(𝑥−2)(𝑥−1)2. The matrix 𝐴 is not diagonalizable: 𝐴−𝐼3 has rank 2, so 𝐸1 is a line, while 1 has algebraic multiplicity 2. The entries above the diagonal cannot all be removed.

Remark 7.5.16 (The quotient version).Beyond 21-242 With quotient spaces (Section 4.6), the first proof needs no choice of complement. Let 𝑈=span(𝐛1). The operator 𝑇/𝑈 on 𝑉/𝑈 (Proposition 4.6.32) is triangularized by some 𝐛2+𝑈,…,𝐛𝑛+𝑈, by induction, since dim𝑉/𝑈=𝑛−1. Then 𝐛1,𝐛2,…,𝐛𝑛 is a basis of 𝑉 (Proposition 4.6.28), and Proposition 4.6.32(b) gives [𝑇]=(𝜆1𝑋0[𝑇/𝑈]), which is upper triangular. The compression is 𝑇/𝑈 in disguise: the isomorphism 𝑊→𝑉/𝑈, 𝐰↦𝐰+𝑈 (Proposition 4.6.29), carries 𝑆 to 𝑇/𝑈, because 𝑆𝐰+𝑈=𝑇𝐰+𝑈 by Lemma 7.5.11(a). So different complements give different compressions, but they are all similar to one operator.

Only one step of the first proof used ℂ: the existence of an eigenvalue, for 𝑇 and then for each compression. Over a general field we can ask for exactly that, and the characteristic polynomial tells us when it holds.

Theorem 7.5.17 (Triangularizable means the characteristic polynomial splits).Beyond 21-242 Let 𝑉 be a finite-dimensional vector space over any field 𝐹, with dim𝑉≥1, and let 𝑇∈L(𝑉). Then 𝑇 is triangularizable if and only if 𝑝𝑇 splits over 𝐹.

Proof.hard

Uses. Proposition 7.5.9 (Multiplicities on the diagonal), Theorem 7.3.6 (Eigenvalues are the roots of the characteristic polynomial), Lemma 7.5.11 (Compressing an operator to a complement), Theorem 5.4.26 (Block-triangular matrices), Proposition 7.1.34 (Roots of a split polynomial), Corollary 7.1.7 (No zero divisors, cancellation, and units), Proposition 4.4.16 (Complements exist)

Techniques. Induction on dimension via an invariant subspace (Technique 7.5.12)

Hint 1. One direction is Proposition 7.5.9. For the other, run the first proof of Theorem 7.5.13. What has to be checked at each step?
Hint 2. The eigenvalue for the first step is a root of 𝑝𝑇. For the induction, show that 𝑝𝑇=(𝑥−𝜆1)𝑝𝑆 for the compression 𝑆, using Lemma 7.5.11(b) and block triangular determinants.
Hint 3. If (𝑥−𝜆1)𝑞=(𝑥−𝜇1)⋯(𝑥−𝜇𝑛), then 𝜆1 is one of the 𝜇𝑖, and cancelling shows that 𝑞 splits.
Sketch. If 𝑇 is upper triangular, 𝑝𝑇=∏(𝑥−𝜆𝑗) splits. Conversely, by induction on 𝑛: a root 𝜆1 of 𝑝𝑇 is an eigenvalue; take an eigenvector, a complement 𝑊 and the compression 𝑆. The block form of Lemma 7.5.11(b) gives 𝑝𝑇=(𝑥−𝜆1)𝑝𝑆, so 𝑝𝑆 splits after cancelling the factor 𝑥−𝜆1. Triangularize 𝑆 by induction and finish as in the first proof of Theorem 7.5.13.

If [𝑇]ℬ︀ is upper triangular, then 𝑝𝑇 splits by Proposition 7.5.9.

Conversely, suppose 𝑝𝑇 splits. We induct on 𝑛=dim𝑉; for 𝑛=1 every basis works. Let 𝑛≥2 and 𝑝𝑇=(𝑥−𝜇1)⋯(𝑥−𝜇𝑛). Then 𝜆1=𝜇1 is a root of 𝑝𝑇, hence an eigenvalue of 𝑇 (Theorem 7.3.6); let 𝐛1 be an eigenvector, 𝑈=span(𝐛1), 𝑊 a complement of 𝑈 (Proposition 4.4.16), and 𝑆 the compression of 𝑇 to 𝑊. For any basis 𝒟︀ of 𝑊 and ℬ︀=𝐛1,𝒟︀, Lemma 7.5.11(b) gives

𝑥𝐼𝑛−[𝑇]ℬ︀=(𝑥−𝜆1−𝑋0𝑥𝐼𝑛−1−[𝑆]𝒟︀),

so by Theorem 5.4.26 (over 𝐹[𝑥]), 𝑝𝑇=(𝑥−𝜆1)𝑝𝑆. Hence (𝑥−𝜇1)𝑝𝑆=(𝑥−𝜇1)(𝑥−𝜇2)⋯(𝑥−𝜇𝑛), and cancelling the nonzero factor 𝑥−𝜇1 (Corollary 7.1.7) gives 𝑝𝑆=(𝑥−𝜇2)⋯(𝑥−𝜇𝑛), which splits. By the induction hypothesis, 𝑊 has a basis 𝐛2,…,𝐛𝑛 in which 𝑆 is upper triangular, and the computation in the first proof of Theorem 7.5.13, which used only Lemma 7.5.11(a), shows that 𝑇 is upper triangular with respect to 𝐛1,…,𝐛𝑛. □

For example, a real matrix is triangularizable over ℝ exactly when all its complex eigenvalues are real, and (0210)∈ℚ2×2 is triangularizable over ℝ but not over ℚ, since 𝑥2−2 splits over ℝ but not over ℚ.

7.5.3Schur’s theorem

Now let 𝑉 be a finite-dimensional inner product space over 𝔽=ℝ or ℂ (Definition 6.1.2). Orthonormal bases are the best bases to compute in: coordinates are inner products (Proposition 6.3.18), the (𝑖,𝑗) entry of [𝑇]𝒪︀ is ⟨𝑇𝐮𝑗,𝐮𝑖⟩, and the change-of-basis matrix to an orthonormal basis of 𝔽𝑛 is unitary, so its inverse is its conjugate transpose (Definition 6.3.33). Can a triangularizing basis be chosen orthonormal?

By Proposition 7.5.2, being upper triangular depends only on the flag of the basis. And Gram–Schmidt replaces a basis by an orthonormal one with the same flag: span(𝐮1,…,𝐮𝑗)=span(𝐛1,…,𝐛𝑗) for every 𝑗 (Theorem 6.4.15(d)). So the answer is yes, at no extra cost.

Theorem 7.5.18 (Schur’s theorem). Let 𝑉 be a finite-dimensional complex inner product space with dim𝑉≥1, and let 𝑇∈L(𝑉).

  1. (a)There is an orthonormal basis 𝒪︀ of 𝑉 such that [𝑇]𝒪︀ is upper triangular.
  2. (b)Every 𝐴∈ℂ𝑛×𝑛 can be written as 𝐴=𝑈𝑅𝑈∗ with 𝑈 unitary and 𝑅 upper triangular. The diagonal entries of 𝑅 are the eigenvalues of 𝐴, counted with multiplicity.

Proof (Gram–Schmidt on a triangularizing basis).routine

Uses. Theorem 7.5.13 (Upper-triangular form), Proposition 7.5.2 (Upper-triangular matrices and invariant subspaces), Theorem 6.4.15 (Gram–Schmidt is correct), Theorem 4.5.17 (Change of basis), Definition 6.3.33 (Unitary and orthogonal matrices), Theorem 6.3.29 (Orthonormal columns)

Techniques. Build a flag of invariant subspaces (Technique 7.5.4)

Hint 1. Start from a basis ℬ︀ with [𝑇]ℬ︀ upper triangular. What does Proposition 7.5.2 say about the subspaces span(𝐛1,…,𝐛𝑗)?
Hint 2. Gram–Schmidt does not change these subspaces (Theorem 6.4.15(d)).
Sketch. Triangularize 𝑇 in some basis ℬ︀ (Theorem 7.5.13), so the flag of ℬ︀ is invariant. Gram–Schmidt turns ℬ︀ into an orthonormal basis 𝒪︀ with the same flag, so 𝑇 is upper triangular with respect to 𝒪︀. For (b), apply (a) to 𝑇𝐴 on ℂ𝑛 with the dot product; the matrix 𝑈 whose columns are the vectors of 𝒪︀ is unitary.

(a) By Theorem 7.5.13, there is a basis ℬ︀=𝐛1,…,𝐛𝑛 of 𝑉 with [𝑇]ℬ︀ upper triangular, and by Proposition 7.5.2 every 𝑉𝑗=span(𝐛1,…,𝐛𝑗) is 𝑇-invariant. Apply Gram–Schmidt to ℬ︀. Since ℬ︀ is linearly independent, Theorem 6.4.15(c),(d) gives an orthonormal basis 𝒪︀=𝐮1,…,𝐮𝑛 of 𝑉 with span(𝐮1,…,𝐮𝑗)=𝑉𝑗 for every 𝑗. These subspaces are 𝑇-invariant, so [𝑇]𝒪︀ is upper triangular by Proposition 7.5.2.

(b) Apply (a) to 𝑇𝐴 on ℂ𝑛 with the dot product, and let 𝑈 be the matrix whose columns are the vectors of 𝒪︀. Its columns are orthonormal, so 𝑈 is unitary (Theorem 6.3.29, Definition 6.3.33) and 𝑈−1=𝑈∗. By Theorem 4.5.17, 𝑅=[𝑇𝐴]𝒪︀=𝑈−1𝐴𝑈=𝑈∗𝐴𝑈 is upper triangular, and 𝐴=𝑈𝑅𝑈∗. The diagonal of 𝑅 is described by Theorem 7.5.13. □

Intuition (The lecture’s proof). The second proof does the algebra and the geometry at once, without quoting Theorem 7.5.13. It is the first proof of Theorem 7.5.13 with the complement chosen orthogonal: take a unit eigenvector 𝐮1, use 𝑊=span(𝐮1)⟂ as the complement, and compress 𝑇 to 𝑊 with the orthogonal projection. The compression is the operator on a smaller space to which the induction hypothesis applies (Technique 7.5.12). This is the template for the spectral theorem in Section 8.2: for a self-adjoint 𝑇, the complement 𝑊 turns out to be invariant itself, the compression is just the restriction, and the entries above the diagonal vanish too.

Proof (induction with an orthogonal complement).standard

Uses. Theorem 7.3.27 (Existence of eigenvalues over ℂ), Theorem 6.4.24 (Orthogonal decomposition), Corollary 6.4.25 (Dimension of the orthogonal complement), Theorem 4.4.13 (Bases and dimension of a sum), Definition 6.3.15 (Orthonormal list, orthonormal basis), Proposition 7.5.2 (Upper-triangular matrices and invariant subspaces), Definition 6.1.2 (Inner product)

Techniques. Induction on dimension via an invariant subspace (Technique 7.5.12)

Hint 1. Induct on 𝑛=dim𝑉. Take a unit eigenvector 𝐮1 and 𝑊=span(𝐮1)⟂. Which operator on 𝑊 should the induction hypothesis be applied to?
Hint 2. Define 𝑆𝐰=𝑇𝐰−⟨𝑇𝐰,𝐮1⟩𝐮1. Check that 𝑆𝐰∈𝑊.
Hint 3. With an orthonormal basis 𝐮2,…,𝐮𝑛 of 𝑊 triangularizing 𝑆, write 𝑇𝐮𝑗=⟨𝑇𝐮𝑗,𝐮1⟩𝐮1+𝑆𝐮𝑗.
Sketch. Induct on 𝑛. Take a unit eigenvector 𝐮1 and 𝑊=span(𝐮1)⟂, of dimension 𝑛−1. The map 𝑆𝐰=𝑇𝐰−⟨𝑇𝐰,𝐮1⟩𝐮1 is an operator on 𝑊 (the compression of 𝑇). By induction, 𝑊 has an orthonormal basis 𝐮2,…,𝐮𝑛 triangularizing 𝑆. Then 𝐮1,…,𝐮𝑛 is orthonormal, and 𝑇𝐮𝑗=⟨𝑇𝐮𝑗,𝐮1⟩𝐮1+𝑆𝐮𝑗∈span(𝐮1,…,𝐮𝑗).

We prove (a) by induction on 𝑛=dim𝑉; part (b) follows from (a) as in the first proof. If 𝑛=1, any unit vector forms an orthonormal basis, and 1×1 matrices are upper triangular.

Let 𝑛≥2, and assume (a) for complex inner product spaces of dimension 𝑛−1. By Theorem 7.3.27, 𝑇 has an eigenvector; dividing it by its norm gives an eigenvector 𝐮1 with ‖𝐮1‖=1 and 𝑇𝐮1=𝜆1𝐮1. Let 𝑊=span(𝐮1)⟂, the set of vectors orthogonal to 𝐮1. By Theorem 6.4.24 and Corollary 6.4.25, 𝑉=span(𝐮1)⊕𝑊 and dim𝑊=𝑛−1; with the inner product of 𝑉 restricted to it, 𝑊 is a complex inner product space. Define

𝑆:𝑊→𝑊,𝑆𝐰=𝑇𝐰−⟨𝑇𝐰,𝐮1⟩𝐮1.

This lands in 𝑊, since ⟨𝑆𝐰,𝐮1⟩=⟨𝑇𝐰,𝐮1⟩−⟨𝑇𝐰,𝐮1⟩⟨𝐮1,𝐮1⟩=0, and it is linear because 𝑇 is linear and the inner product is linear in its first slot (Definition 6.1.2). (𝑆𝐰 is the orthogonal projection of 𝑇𝐰 onto 𝑊, so 𝑆 is the compression of 𝑇 to 𝑊.)

By the induction hypothesis and Proposition 7.5.2, 𝑊 has an orthonormal basis 𝐮2,…,𝐮𝑛 with 𝑆𝐮𝑗∈span(𝐮2,…,𝐮𝑗) for 2≤𝑗≤𝑛. Each 𝐮𝑗 with 𝑗≥2 lies in 𝑊, so it is orthogonal to 𝐮1, and 𝒪︀=𝐮1,…,𝐮𝑛 is orthonormal (Definition 6.3.15); it is a basis of 𝑉=span(𝐮1)⊕𝑊 by Theorem 4.4.13. Finally, 𝑇𝐮1=𝜆1𝐮1, and for 𝑗≥2,

𝑇𝐮𝑗=⟨𝑇𝐮𝑗,𝐮1⟩𝐮1+𝑆𝐮𝑗∈span(𝐮1,𝐮2,…,𝐮𝑗).

By Proposition 7.5.2, [𝑇]𝒪︀ is upper triangular. □

The first proof separates the two ingredients: the algebra (an invariant flag exists over ℂ) and the geometry (Gram–Schmidt keeps flags). It shows that any triangularizing basis can be orthonormalized, so for instance the basis of Example 7.5.15 turns into a Schur basis of that matrix. The second proof is self-contained and computes the entries as it goes: [𝑇]𝒪︀ has (𝑖,𝑗) entry ⟨𝑇𝐮𝑗,𝐮𝑖⟩=⟨𝑆𝐮𝑗,𝐮𝑖⟩ for 𝑖,𝑗≥2.

Example 7.5.19 (A Schur form). Let 𝐴=(12−25)∈ℂ2×2. Then det(𝐴−𝜆𝐼2)=(1−𝜆)(5−𝜆)+4=(𝜆−3)2, so 3 is the only eigenvalue, and 𝐴−3𝐼2=(−22−22) gives the unit eigenvector 𝐮1=12(1,1). Complete it to the orthonormal basis 𝐮1,𝐮2 with 𝐮2=12(−1,1). Since 𝐴𝐮2=12(1,7),

⟨𝐴𝐮2,𝐮1⟩=12(1+7)=4,⟨𝐴𝐮2,𝐮2⟩=12(−1+7)=3,

and so

𝐴=𝑈𝑅𝑈∗,𝑈=12(1−111),𝑅=(3403).

Here 𝐸3 is a line, so 𝐴 is not diagonalizable, and no unitary change of basis can remove the entry 4. In fact no change of basis at all can make it 0, although a non-orthonormal one can make it 1.

Remark 7.5.20 (What is unique). A Schur form is far from unique. The eigenvalues can be made to appear on the diagonal in any order (Exercise 7.5.38), and even for a fixed order, replacing 𝑈 by 𝑈𝐷 with 𝐷=diag(𝑒𝑖𝜃1,…,𝑒𝑖𝜃𝑛) replaces 𝑅 by 𝐷∗𝑅𝐷, which multiplies 𝑟𝑖𝑗 by 𝑒𝑖(𝜃𝑗−𝜃𝑖). What does not change is the list of diagonal entries, up to order (Proposition 7.5.9), and the total size ∑𝑖,𝑗|𝑟𝑖𝑗|2=tr(𝐴∗𝐴) (Exercise 7.5.41).

For a real matrix with nonreal eigenvalues, 𝑈 and 𝑅 cannot both be real, since a real upper-triangular 𝑅 has real eigenvalues. If all eigenvalues of 𝐴∈ℝ𝑛×𝑛 are real, then 𝐴=𝑄𝑅𝑄𝖳 with 𝑄 orthogonal and 𝑅 real (Exercise 7.5.42). For general real matrices, see Section 7.5.5.

Exercise 7.5.21.routine Find a unitary 𝑈 and an upper-triangular 𝑅 with (3−111)=𝑈𝑅𝑈∗.

Solution.

Uses. Theorem 7.5.18 (Schur’s theorem), Theorem 7.2.9 (Eigenvalues and invertibility), Proposition 6.3.18 (Coordinates in an orthonormal basis)

Hint 1. Find an eigenvalue and a unit eigenvector 𝐮1; complete it to an orthonormal basis of ℂ2.
Hint 2. The entries of 𝑅 are 𝑟𝑖𝑗=⟨𝐴𝐮𝑗,𝐮𝑖⟩.
Sketch. det(𝐴−𝜆𝐼2)=(𝜆−2)2; 𝐮1=12(1,1), 𝐮2=12(−1,1); 𝐴𝐮2=12(−4,0), so 𝑅=(2−202).

det(𝐴−𝜆𝐼2)=(3−𝜆)(1−𝜆)+1=𝜆2−4𝜆+4=(𝜆−2)2, so 2 is the only eigenvalue. From 𝐴−2𝐼2=(1−11−1), a unit eigenvector is 𝐮1=12(1,1); complete it with 𝐮2=12(−1,1). Then 𝐴𝐮1=2𝐮1, and 𝐴𝐮2=12(−4,0), so by Proposition 6.3.18

𝑟12=⟨𝐴𝐮2,𝐮1⟩=12(−4)=−2,𝑟22=⟨𝐴𝐮2,𝐮2⟩=12(4)=2.

Hence 𝐴=𝑈𝑅𝑈∗ with 𝑈=12(1−111) and 𝑅=(2−202). Check: 𝐴𝑈=12(2−420)=𝑈𝑅.

7.5.4Trace and determinant from eigenvalues

Every operator on a complex space has a triangular matrix whose diagonal lists its eigenvalues. Trace and determinant are both read off the diagonal of a triangular matrix, and both are independent of the basis. Putting these facts together gives formulas for them that mention no basis at all.

Theorem 7.5.22 (Trace and determinant from eigenvalues). Let 𝑉 be a finite-dimensional vector space over 𝐹 with dim𝑉=𝑛≥1, let 𝑇∈L(𝑉), and suppose that 𝑝𝑇 splits over 𝐹, as it always does when 𝐹=ℂ. Let 𝜆1,…,𝜆𝑛 be the eigenvalues of 𝑇 counted with multiplicity. Then

tr𝑇=𝜆1+𝜆2+⋯+𝜆𝑛anddet𝑇=𝜆1𝜆2⋯𝜆𝑛.

In particular, for 𝐴∈ℂ𝑛×𝑛, and so for every real matrix regarded as a complex one, tr𝐴 is the sum and det𝐴 the product of the complex eigenvalues of 𝐴, counted with multiplicity.

Proof (triangularize).routine

Uses. Theorem 7.5.13 (Upper-triangular form), Theorem 7.5.17 (Triangularizable means the characteristic polynomial splits), Definition 4.5.30 (Trace of an operator), Definition 5.4.18 (Determinant of an operator), Definition 2.4.26 (Trace), Proposition 5.3.15 (Triangular matrices), Proposition 7.5.9 (Multiplicities on the diagonal)

Techniques. Prove it for triangular matrices first (Technique 7.5.23), Choose a convenient basis (Technique 4.5.33)

Hint 1. Both sides do not depend on the basis. Which basis makes the left sides easy to compute?
Hint 2. In a triangularizing basis ℬ︀, what are tr[𝑇]ℬ︀ and det[𝑇]ℬ︀, and what is the diagonal?
Sketch. Choose ℬ︀ with 𝑅=[𝑇]ℬ︀ upper triangular (Theorem 7.5.13 over ℂ, Theorem 7.5.17 in general). Then tr𝑇=tr𝑅 and det𝑇=det𝑅 are the sum and the product of the diagonal entries, and the diagonal entries are the eigenvalues counted with multiplicity (Proposition 7.5.9).

By Theorem 7.5.13 when 𝐹=ℂ, and by Theorem 7.5.17 in general, there is a basis ℬ︀ with 𝑅=[𝑇]ℬ︀ upper triangular. By Definition 4.5.30 and Definition 5.4.18, tr𝑇=tr𝑅 and det𝑇=det𝑅. By Definition 2.4.26 and Proposition 5.3.15,

tr𝑅=𝑟11+⋯+𝑟𝑛𝑛,det𝑅=𝑟11𝑟22⋯𝑟𝑛𝑛.

By Proposition 7.5.9, 𝑝𝑇=(𝑥−𝑟11)⋯(𝑥−𝑟𝑛𝑛), so 𝑟11,…,𝑟𝑛𝑛 are the eigenvalues of 𝑇 counted with multiplicity, in some order. Sums and products do not depend on the order. □

Technique 7.5.23 (Prove it for triangular matrices first). If a statement about an operator 𝑇 (or a square matrix 𝐴) does not change under change of basis (or similarity), it is enough to prove it for an upper-triangular matrix of 𝑇, which exists over ℂ (Theorem 7.5.13), or over any field once the characteristic polynomial splits. Triangular matrices are easy to compute with: products, powers, inverses, and polynomials of upper-triangular matrices are upper triangular, and their diagonals are computed entry by entry (Proposition 2.4.36, Exercise 7.5.34), while the diagonal itself is the list of eigenvalues.

Reach for it for identities that involve eigenvalues together with similarity invariants: trace, determinant, the characteristic polynomial, tr(𝐴𝑘), invertibility, nilpotence. It is a special case of Technique 4.5.33. Over ℝ, either check that the characteristic polynomial splits, or regard the real matrix as a complex one first. The canonical example is Theorem 7.5.22. Later it proves the Cayley–Hamilton theorem (Section 7.6).

The second proof avoids triangularization. It compares two descriptions of the coefficients of the characteristic polynomial.

Proof (read off the characteristic polynomial).standard

Hint 1. Section 7.3 computes two coefficients of 𝑝𝐴=det(𝑥𝐼𝑛−𝐴) directly from the entries of 𝐴. Which ones?
Hint 2. Compare them with the coefficients of (𝑥−𝜆1)⋯(𝑥−𝜆𝑛) from Proposition 7.1.36.
Sketch. For any matrix 𝐴 of 𝑇, Section 7.3 shows 𝑝𝑇=𝑥𝑛−(tr𝐴)𝑥𝑛−1+⋯+(−1)𝑛det𝐴: the constant term is det(−𝐴), and in the Leibniz expansion only the identity permutation reaches degree 𝑛−1. On the other hand 𝑝𝑇=∏(𝑥−𝜆𝑖) has these coefficients equal to −∑𝜆𝑖 and (−1)𝑛∏𝜆𝑖.
Let 𝐴=[𝑇]ℬ︀ for some basis ℬ︀, so 𝑝𝑇=det(𝑥𝐼𝑛−𝐴), tr𝑇=tr𝐴, and det𝑇=det𝐴. By Section 7.3, the coefficient of 𝑥𝑛−1 in 𝑝𝑇 is −tr𝐴 and its constant term is (−1)𝑛det𝐴. (The constant term is det(−𝐴), by setting 𝑥=0; for the other coefficient, a permutation other than the identity moves at least two indices, so its term in the Leibniz expansion has at most 𝑛−2 factors containing 𝑥, and only (𝑥−𝑎11)⋯(𝑥−𝑎𝑛𝑛) contributes.) On the other hand 𝑝𝑇=(𝑥−𝜆1)⋯(𝑥−𝜆𝑛), whose coefficient of 𝑥𝑛−1 is −(𝜆1+⋯+𝜆𝑛) and whose constant term is (−1)𝑛𝜆1⋯𝜆𝑛 (Proposition 7.1.36). Comparing coefficients, tr𝐴=𝜆1+⋯+𝜆𝑛 and det𝐴=𝜆1⋯𝜆𝑛. □

The first proof is the more flexible one: the same reduction gives the eigenvalues of 𝐴𝑘 and of 𝐴−1 (Exercise 7.5.34). The second proof needs no triangular form, only two coefficients of the characteristic polynomial, which are −tr𝑇 and (−1)𝑛det𝑇 over every field, whether or not 𝑝𝑇 splits.

Example 7.5.24 (Rotations and a 3×3 check). The rotation 𝑅𝜃=(cos𝜃−sin𝜃sin𝜃cos𝜃)∈ℝ2×2 has the complex eigenvalues cos𝜃±𝑖sin𝜃 (Exercise 7.2.42). Their sum is 2cos𝜃=tr𝑅𝜃 and their product is cos2𝜃+sin2𝜃=1=det𝑅𝜃. For a real matrix the nonreal eigenvalues come in conjugate pairs (Proposition 7.1.38), and each pair 𝜆,𝜆 contributes 2Re𝜆 to the trace and |𝜆|2>0 to the determinant.

For the matrix 𝐴 of Example 7.5.15, with eigenvalues 2,1,1, the theorem predicts tr𝐴=4 and det𝐴=2. Directly: tr𝐴=2+2+0=4, and expanding along the first row, det𝐴=2(0+2)−0+(−1)(0+2)=2.

The trace of an operator

The trace of an operator was defined through a matrix (Definition 4.5.30). Over ℂ, Theorem 7.5.22 gives it a meaning without any basis: it is the sum of the eigenvalues. The matrix identities of Section 2.4.7 also have versions for maps, and for those, maps between different spaces are allowed.

Proposition 7.5.25 (Cyclicity of the trace for linear maps). Let 𝑉 and 𝑊 be finite-dimensional, both of dimension at least 1, and let 𝑆∈L(𝑉,𝑊) and 𝑇∈L(𝑊,𝑉). Then

tr(𝑇𝑆)=tr(𝑆𝑇),

where 𝑇𝑆∈L(𝑉) and 𝑆𝑇∈L(𝑊). In particular tr(𝑆𝑇)=tr(𝑇𝑆) for all 𝑆,𝑇∈L(𝑉).

Proof.routine

Uses. Definition 4.5.30 (Trace of an operator), Theorem 2.4.2 (The matrix of a composite), Proposition 2.4.28 (Cyclicity of the trace)

Techniques. Translate between maps and matrices (Technique 2.3.27), Cyclicity of trace (Technique 2.4.29)

Hint 1. Choose bases ℬ︀ of 𝑉 and 𝒞︀ of 𝑊, and write the matrices of 𝑇𝑆 and 𝑆𝑇 as products.
Hint 2. Proposition 2.4.28 holds for rectangular matrices.
Sketch. With 𝑋=[𝑆]𝒞︀ℬ︀ and 𝑌=[𝑇]ℬ︀𝒞︀, we have [𝑇𝑆]ℬ︀=𝑌𝑋 and [𝑆𝑇]𝒞︀=𝑋𝑌, and tr(𝑌𝑋)=tr(𝑋𝑌).

Let ℬ︀ and 𝒞︀ be bases of 𝑉 and 𝑊, with dim𝑉=𝑛 and dim𝑊=𝑚, and let 𝑋=[𝑆]𝒞︀ℬ︀∈𝐹𝑚×𝑛 and 𝑌=[𝑇]ℬ︀𝒞︀∈𝐹𝑛×𝑚. By Theorem 2.4.2, [𝑇𝑆]ℬ︀=𝑌𝑋 and [𝑆𝑇]𝒞︀=𝑋𝑌. So by Definition 4.5.30 and Proposition 2.4.28,

tr(𝑇𝑆)=tr(𝑌𝑋)=tr(𝑋𝑌)=tr(𝑆𝑇).

□

Corollary 7.5.26 (𝑆𝑇−𝑇𝑆≠id𝑉). Let 𝑉 be finite-dimensional with 𝑛=dim𝑉≥1, and suppose that 𝑛∗1≠0 in 𝐹 (for instance, 𝐹=ℝ or ℂ). Then 𝑆𝑇−𝑇𝑆≠id𝑉 for all 𝑆,𝑇∈L(𝑉).

Proof.routine

Uses. Proposition 7.5.25 (Cyclicity of the trace for linear maps), Definition 4.5.30 (Trace of an operator), Proposition 2.4.12 (Matrices of sums and scalar multiples), Proposition 2.4.27 (The trace is linear)

Techniques. Cyclicity of trace (Technique 2.4.29)

Hint. Take the trace of both sides.
Sketch. tr(𝑆𝑇−𝑇𝑆)=tr(𝑆𝑇)−tr(𝑇𝑆)=0, while trid𝑉=tr𝐼𝑛=𝑛∗1≠0.
The trace of operators is linear, by Definition 4.5.30, Proposition 2.4.12, and Proposition 2.4.27. So by Proposition 7.5.25, tr(𝑆𝑇−𝑇𝑆)=tr(𝑆𝑇)−tr(𝑇𝑆)=0. But trid𝑉=tr𝐼𝑛=𝑛∗1≠0. Hence 𝑆𝑇−𝑇𝑆≠id𝑉. □

This is Proposition 2.4.30 restated for maps, and the hypotheses matter in the same ways. On the infinite-dimensional space 𝐹[𝑥], differentiation and multiplication by 𝑥 satisfy 𝐷𝑀−𝑀𝐷=id (Exercise 2.4.51), and over 𝔽𝑝 there are 𝑝×𝑝 matrices with 𝐴𝐵−𝐵𝐴=𝐼𝑝 (Exercise 2.4.52). The trace argument shows more than the corollary states: 𝑆𝑇−𝑇𝑆 always has trace 0. Over ℂ the converse also holds, a theorem of Shoda: every matrix of trace 0 is of the form 𝐴𝐵−𝐵𝐴. The practice bank has a proof.

7.5.5Real operators: block upper-triangular form Beyond 21-242

The real block form and the real Schur form are not covered in lecture.

Over ℝ an operator need not have any eigenvalue, so it need not be triangularizable. But it always comes close: the bottom of an invariant flag may have to be a plane instead of a line, and the matrix becomes upper triangular with some 2×2 blocks on the diagonal. The key fact, proved in Section 7.3, is that small invariant subspaces always exist: every operator on a nonzero finite-dimensional real vector space has an invariant subspace of dimension 1 or 2. One proof complexifies a matrix 𝐴 of 𝑇: a complex eigenvector 𝐱+𝑖𝐲 (𝐱,𝐲∈ℝ𝑛) for an eigenvalue 𝑎+𝑏𝑖 gives 𝐴𝐱=𝑎𝐱−𝑏𝐲 and 𝐴𝐲=𝑏𝐱+𝑎𝐲, by comparing real and imaginary parts, so the vectors with coordinates 𝐱 and 𝐲 span an invariant subspace. The induction of Theorem 7.5.13 then goes through with this subspace in place of an eigenvector.

Theorem 7.5.27 (Real block upper-triangular form). Let 𝑉 be a real vector space with 1≤dim𝑉<∞, and let 𝑇∈L(𝑉).

  1. (a)There is a basis ℬ︀ of 𝑉 such that [𝑇]ℬ︀ is block upper triangular, with every diagonal block of size 1×1 or 2×2, and with no 2×2 diagonal block having a real eigenvalue.
  2. (b)If 𝑉 is a real inner product space, ℬ︀ can be chosen orthonormal. In particular, every 𝐴∈ℝ𝑛×𝑛 can be written as 𝐴=𝑄𝑅𝑄𝖳 with 𝑄 orthogonal and 𝑅 block upper triangular as in (a) (the real Schur form).

Proof.hard

Uses. an invariant subspace of dimension 1 or 2 exists (Section 7.3), Proposition 7.2.7 (Invariant lines are spanned by eigenvectors), Proposition 4.4.16 (Complements exist), Lemma 7.5.11 (Compressing an operator to a complement), Theorem 6.4.24 (Orthogonal decomposition), Theorem 4.4.13 (Bases and dimension of a sum), Theorem 4.5.17 (Change of basis)

Techniques. Induction on dimension via an invariant subspace (Technique 7.5.12)

Hint 1. Strong induction on dim𝑉. Choose an invariant subspace 𝑈 of dimension 1 if 𝑇 has a real eigenvalue, and of dimension 2 otherwise.
Hint 2. Compress 𝑇 to a complement 𝑊 of 𝑈, apply the induction hypothesis to the compression, and read off the block form from Lemma 7.5.11(b).
Hint 3. For (b), use the orthogonal complement and orthonormal bases throughout.
Sketch. If 𝑇 has a real eigenvalue, let 𝑈 be an eigenvector’s span; otherwise an invariant subspace of dimension 1 or 2 (Section 7.3) must be a plane 𝑈, and 𝑇↾𝑈 has no real eigenvalue. If 𝑈=𝑉 we are done. Otherwise compress 𝑇 to a complement 𝑊 (for (b), to 𝑈⟂), use the induction hypothesis on the compression, and assemble with Lemma 7.5.11(b): the first diagonal block is [𝑇↾𝑈], the rest are the blocks of the compression.

(a) We use strong induction on 𝑛=dim𝑉. Choose a 𝑇-invariant subspace 𝑈 as follows. If 𝑇 has a real eigenvalue, let 𝑈 be the span of an eigenvector, of dimension 1 (Proposition 7.2.7). Otherwise, by Section 7.3 there is an invariant subspace 𝑈 of dimension 1 or 2, and dimension 1 is impossible since it would be spanned by an eigenvector; so dim𝑈=2. In that case 𝑇↾𝑈 has no real eigenvalue, since an eigenvector of 𝑇↾𝑈 would be an eigenvector of 𝑇.

If 𝑈=𝑉, any basis of 𝑉 works: [𝑇]ℬ︀ is a single block of size 1, or of size 2 with no real eigenvalue. Otherwise let 𝑊 be a complement of 𝑈 (Proposition 4.4.16), so 1≤dim𝑊<𝑛, and let 𝑆∈L(𝑊) be the compression of 𝑇 to 𝑊. By the induction hypothesis, 𝑊 has a basis 𝒟︀ with [𝑆]𝒟︀ block upper triangular as in (a). Let 𝒞︀ be any basis of 𝑈, and ℬ︀=𝒞︀,𝒟︀. By Lemma 7.5.11(b),

[𝑇]ℬ︀=([𝑇↾𝑈]𝒞︀𝑋0[𝑆]𝒟︀),

which is block upper triangular, with first diagonal block [𝑇↾𝑈]𝒞︀ and the remaining diagonal blocks those of [𝑆]𝒟︀. All of them have the required form.

(b) Run the same induction with 𝑊=𝑈⟂, which is a complement of 𝑈 by Theorem 6.4.24, and with an orthonormal basis 𝒞︀ of 𝑈 and, by the induction hypothesis applied to the real inner product space 𝑊, an orthonormal basis 𝒟︀ of 𝑊. Then 𝒞︀,𝒟︀ is orthonormal, since vectors of 𝑈 and of 𝑈⟂ are orthogonal. For 𝐴∈ℝ𝑛×𝑛, apply this to 𝑇𝐴 on ℝ𝑛 with the dot product, and let 𝑄 have the resulting orthonormal basis as columns. Then 𝑄 is orthogonal, and 𝑅=𝑄𝖳𝐴𝑄=[𝑇𝐴]ℬ︀ by Theorem 4.5.17. □

Example 7.5.28 (An invariant plane). Let

𝐴=(101010−22−1)∈ℝ3×3.

Its characteristic polynomial is (𝑥−1)(𝑥2+1), so 1 is its only real eigenvalue, with eigenvector 𝐯=(1,1,0). Over ℂ, 𝐳=(1,0,−1+𝑖) is an eigenvector for 𝑖, as one checks: 𝐴𝐳=(𝑖,0,−1−𝑖)=𝑖𝐳. As above, split 𝐳=𝐱+𝑖𝐲 with 𝐱=(1,0,−1) and 𝐲=(0,0,1). Then

𝐴𝐱=(0,0,−1)=−𝐲,𝐴𝐲=(1,0,−1)=𝐱,

so span(𝐱,𝐲) is an invariant plane on which 𝐴 acts by (01−10), a rotation by a right angle, with eigenvalues ±𝑖. In the basis 𝐱,𝐲,𝐯,

[𝑇𝐴]=(010−100001).

Here the plane happens to have an invariant complement, span(𝐯), so the matrix is even block diagonal. As a check, tr𝐴=1=𝑖+(−𝑖)+1 and det𝐴=1=𝑖⋅(−𝑖)⋅1.

Exercises

Exercise 7.5.29.routine Find an invertible 𝑆 and an upper-triangular 𝑅 with diagonal 1,1,2 such that 𝐴=𝑆𝑅𝑆−1, where

𝐴=(2−202−1−2−103)∈ℂ3×3.

Solution.

Uses. Theorem 7.5.13 (Upper-triangular form), Lemma 7.5.11 (Compressing an operator to a complement), Corollary 7.5.14 (Triangularizing a matrix)

Hint 1. (2,1,1) is an eigenvector. Use the complement span(𝐞2,𝐞3) and compute the compression.
Hint 2. In the basis (2,1,1),𝐞2,𝐞3, the matrix is (1−1000−2013). Triangularize the bottom right 2×2 block, starting with an eigenvector for 1.
Sketch. 𝐴(2,1,1)=(2,1,1). With 𝑀 having columns (2,1,1),𝐞2,𝐞3, 𝑀−1𝐴𝑀=(1−1000−2013). The block (0−213) has eigenvector (−2,1) for 1, giving 𝐛2=(0,−2,1). With 𝑆 having columns (2,1,1),(0,−2,1),𝐞3, 𝑆−1𝐴𝑆=(120011002).

First, 𝐴(2,1,1)=(4−2,4−1−2,−2+3)=(2,1,1), so 𝐛1=(2,1,1) is an eigenvector for 1. Since its first coordinate is nonzero, 𝐛1,𝐞2,𝐞3 is a basis, and 𝑊=span(𝐞2,𝐞3) is a complement of span(𝐛1). Writing (𝑥1,𝑥2,𝑥3)=(𝑥1/2)𝐛1+(0,𝑥2−𝑥1/2,𝑥3−𝑥1/2), the projection onto 𝑊 along span(𝐛1) is (𝑥1,𝑥2,𝑥3)↦(0,𝑥2−𝑥1/2,𝑥3−𝑥1/2). Since 𝐴𝐞2=(−2,−1,0) and 𝐴𝐞3=(0,−2,3), the compression sends

𝐞2↦(0,−1+1,0+1)=𝐞3,𝐞3↦(0,−2,3)=−2𝐞2+3𝐞3,

so its matrix is (0−213), with characteristic polynomial 𝑥2−3𝑥+2=(𝑥−1)(𝑥−2). An eigenvector for 1 is (−2,1), that is, 𝐛2=−2𝐞2+𝐞3=(0,−2,1). Complete with 𝐛3=𝐞3. Then

𝐴𝐛2=(4,0,3)=2𝐛1+𝐛2,𝐴𝐛3=(0,−2,3)=𝐛2+2𝐛3.

So with

𝑆=(2001−20111),𝑅=(120011002),

we have 𝐴𝑆=𝑆𝑅 column by column, that is, 𝐴=𝑆𝑅𝑆−1. The diagonal 1,1,2 agrees with 𝑝𝐴=(𝑥−1)2(𝑥−2): the first entry comes from 𝐛1 and the other two from the compression. (Starting instead with the eigenvector (1,0,1) for 2 gives a triangular form with diagonal 2,1,1.)

Exercise 7.5.30.standard
  1. (a)Let 𝑉 be a complex vector space with dim𝑉=𝑛≥1, and 𝑇∈L(𝑉). Show that for each 𝑘=0,1,…,𝑛, there is a 𝑇-invariant subspace of 𝑉 of dimension 𝑘.
  2. (b)Let 𝑇∈L(ℝ4) be 𝑇(𝑥1,𝑥2,𝑥3,𝑥4)=(−𝑥2,𝑥1,−𝑥4,𝑥3). Show that 𝑇 has no invariant subspace of dimension 1 or 3.

Solution.

Uses. Theorem 7.5.13 (Upper-triangular form), Proposition 7.5.2 (Upper-triangular matrices and invariant subspaces), Proposition 7.2.7 (Invariant lines are spanned by eigenvectors), Theorem 7.3.6 (Eigenvalues are the roots of the characteristic polynomial), Corollary 7.1.43 (Odd degree gives a real root)

Techniques. Build a flag of invariant subspaces (Technique 7.5.4)

Hint 1. For (a), a triangularizing basis comes with an invariant flag.
Hint 2. For (b), compute 𝑇2. Why does it follow that 𝑇 has no real eigenvalue?
Hint 3. If 𝑈 is invariant of dimension 3, then 𝑇↾𝑈 is an operator on a 3-dimensional real space. Its characteristic polynomial has odd degree.
Sketch. (a) The flag span(𝐛1,…,𝐛𝑘) of a triangularizing basis. (b) 𝑇2=−id, so 𝑇𝐯=𝜆𝐯 forces (𝜆2+1)𝐯=𝟎 and 𝑇 has no eigenvector, hence no invariant line. An invariant 𝑈 of dimension 3 would give an operator 𝑇↾𝑈 whose characteristic polynomial has odd degree, hence a real root, hence an eigenvector of 𝑇.

(a) By Theorem 7.5.13 there is a basis 𝐛1,…,𝐛𝑛 with [𝑇] upper triangular, and by Proposition 7.5.2 the subspaces span(𝐛1,…,𝐛𝑘), of dimension 𝑘, are 𝑇-invariant (𝑘=0 gives {𝟎}).

(b) 𝑇2(𝑥1,𝑥2,𝑥3,𝑥4)=𝑇(−𝑥2,𝑥1,−𝑥4,𝑥3)=(−𝑥1,−𝑥2,−𝑥3,−𝑥4), so 𝑇2=−id. If 𝑇𝐯=𝜆𝐯 with 𝜆∈ℝ, then −𝐯=𝑇2𝐯=𝜆2𝐯, so (𝜆2+1)𝐯=𝟎 and 𝐯=𝟎, as 𝜆2+1>0. So 𝑇 has no eigenvector, and no invariant subspace of dimension 1 (Proposition 7.2.7).

Suppose 𝑈 were invariant with dim𝑈=3. Then 𝑇↾𝑈 is an operator on 𝑈, and its characteristic polynomial has degree 3, so it has a real root (Corollary 7.1.43), which is an eigenvalue of 𝑇↾𝑈 (Theorem 7.3.6). An eigenvector of 𝑇↾𝑈 is an eigenvector of 𝑇, which we just showed does not exist. So no such 𝑈 exists. (Invariant subspaces of dimension 2 do exist, such as span(𝐞1,𝐞2).)

Exercise 7.5.31.standard How many of the 16 matrices in 𝔽22×2 are triangularizable over 𝔽2? List the ones that are not.

Solution.

Uses. Proposition 7.5.2 (Upper-triangular matrices and invariant subspaces), Theorem 4.1.16 (Extending to a basis), Theorem 7.2.9 (Eigenvalues and invertibility), Theorem 5.1.13 (The 2×2 determinant)

Techniques. Try small cases and small fields (Technique 1.2.38)

Hint 1. In dimension 2, triangularizable means: there is an eigenvector. Why is that enough?
Hint 2. 𝜆 is an eigenvalue if and only if det(𝐴−𝜆𝐼2)=0. There are only two candidates, 𝜆=0 and 𝜆=1. Express both determinants in terms of tr𝐴 and det𝐴.
Sketch. A 2×2 matrix is triangularizable if and only if it has an eigenvector 𝐛1: extend it to a basis, and 𝑇𝐛2∈span(𝐛1,𝐛2)=𝔽22 automatically. Since det(𝐴−𝜆𝐼2)=𝜆2−(tr𝐴)𝜆+det𝐴, the candidates fail exactly when det𝐴=1 and tr𝐴=1, which forces 𝑏=𝑐=1 and {𝑎,𝑑}={0,1}. So 14 are triangularizable, and the exceptions are (1110) and (0111).

In dimension 2, 𝐴 is triangularizable if and only if it has an eigenvector. If 𝑇𝐴 is upper triangular with respect to 𝐛1,𝐛2, then 𝐛1 is an eigenvector. Conversely, if 𝐛1 is an eigenvector, extend it to a basis 𝐛1,𝐛2 (Theorem 4.1.16); then 𝑇𝐛1∈span(𝐛1) and 𝑇𝐛2∈span(𝐛1,𝐛2), so 𝑇𝐴 is upper triangular by Proposition 7.5.2.

For 𝐴=(𝑎𝑏𝑐𝑑), Theorem 5.1.13 gives det(𝐴−𝜆𝐼2)=(𝑎−𝜆)(𝑑−𝜆)−𝑏𝑐=𝜆2−(𝑎+𝑑)𝜆+(𝑎𝑑−𝑏𝑐). By Theorem 7.2.9, 𝐴 has no eigenvalue in 𝔽2 if and only if this is nonzero for 𝜆=0 and 𝜆=1, that is,

det𝐴=1and1+tr𝐴+det𝐴=1,

using −1=1 in 𝔽2. Substituting det𝐴=1 into the second condition gives tr𝐴=1, since 1+1=0. From 𝑎+𝑑=1, one of 𝑎,𝑑 is 0, so 𝑎𝑑=0 and det𝐴=−𝑏𝑐=𝑏𝑐=1, giving 𝑏=𝑐=1. So the matrices with no eigenvalue are exactly

(1110)and(0111),

whose characteristic polynomial is 𝑥2+𝑥+1, the one irreducible quadratic over 𝔽2. The other 14 matrices are triangularizable.

Exercise 7.5.32.routine
  1. (a)Find 𝐴∈ℚ2×2 that is triangularizable over ℝ but not over ℚ.
  2. (b)Find 𝐴∈ℝ3×3 that has a real eigenvalue but is not triangularizable over ℝ.

Solution.

Uses. Proposition 7.2.7 (Invariant lines are spanned by eigenvectors), Proposition 7.5.9 (Multiplicities on the diagonal), Theorem 7.2.9 (Eigenvalues and invertibility)

Hint 1. A triangularizable operator has an eigenvector. Look for a rational matrix whose eigenvalues are irrational.
Hint 2. For (b), if 𝐴 is upper triangular in a real basis, its characteristic polynomial is a product of real linear factors (Proposition 7.5.9). Put a rotation block next to a 1×1 block.
Sketch. (a) 𝐴=(0210) has eigenvalues ±2: no eigenvector in ℚ2, but an eigenvector (2,1) in ℝ2, which a second vector completes to a triangularizing basis. (b) 𝐴=(0−10100001) has eigenvalue 1, but 𝑝𝐴=(𝑥−1)(𝑥2+1) is not a product of real linear factors.

(a) Let 𝐴=(0210). Then det(𝐴−𝜆𝐼2)=𝜆2−2, which has no rational root, so 𝐴 has no eigenvalue in ℚ (Theorem 7.2.9). The first vector of a triangularizing basis would be an eigenvector, so 𝐴 is not triangularizable over ℚ. Over ℝ, 𝐴(2,1)=(2,2)=2(2,1), and in the basis (2,1),𝐞2 the first vector spans an invariant line, so the flag {𝟎}⊆span((2,1))⊆ℝ2 is invariant and 𝐴 is triangularizable over ℝ.

(b) Let 𝐴=(0−10100001). Then 𝐴𝐞3=𝐞3, so 1 is an eigenvalue. By cofactor expansion along the last row, 𝑝𝐴=(𝑥−1)(𝑥2+1). If 𝐴=𝑆𝑅𝑆−1 with 𝑆,𝑅 real and 𝑅 upper triangular, then 𝑝𝐴=(𝑥−𝑟11)(𝑥−𝑟22)(𝑥−𝑟33) with real 𝑟𝑗𝑗 (Proposition 7.5.9). Evaluating at 𝑖 gives 0=𝑝𝐴(𝑖)=(𝑖−𝑟11)(𝑖−𝑟22)(𝑖−𝑟33), a product of nonzero complex numbers, which is impossible. So 𝐴 is not triangularizable over ℝ.

Exercise 7.5.33.standard Prove or disprove.

  1. (a)If 𝐴∈𝐹𝑛×𝑛 is triangularizable, then so is 𝐴𝖳.
  2. (b)If 𝐴,𝐵∈ℝ2×2 are triangularizable over ℝ, then so is 𝐴+𝐵.
  3. (c)If 𝑇 is upper triangular with respect to both 𝐛1,…,𝐛𝑛 and 𝐜1,…,𝐜𝑛, then 𝐜1 is a multiple of 𝐛1.
  4. (d)Two upper-triangular matrices with the same diagonal are similar.

Solution.

Hint 1. For (a), transpose 𝐴=𝑆𝑅𝑆−1. The result has a lower-triangular middle factor; reverse the order of a basis.
Hint 2. For (b), add an upper-triangular and a lower-triangular nilpotent matrix.
Hint 3. For (c) and (d), look at diagonal matrices, and at the identity.
Sketch. (a) True: 𝐴𝖳=(𝑆𝖳)−1𝑅𝖳𝑆𝖳, and 𝑅𝖳 is lower triangular, so with the reversal permutation matrix 𝑃, 𝑃𝑅𝖳𝑃 is upper triangular and similar to 𝐴𝖳. (b) False: (0100)+(00−10) is a rotation. (c) False: diag(1,2) with the bases 𝐞1,𝐞2 and 𝐞2,𝐞1. (d) False: 𝐼2 and (1101).

(a) True. Write 𝐴=𝑆𝑅𝑆−1 with 𝑅 upper triangular. Then 𝐴𝖳=(𝑆−1)𝖳𝑅𝖳𝑆𝖳, and 𝑅𝖳 is lower triangular. Let 𝑃 be the permutation matrix with 𝑃𝐞𝑗=𝐞𝑛+1−𝑗, so 𝑃2=𝐼𝑛. The (𝑖,𝑗) entry of 𝑃𝑅𝖳𝑃 is the (𝑛+1−𝑖,𝑛+1−𝑗) entry of 𝑅𝖳, which is 𝑟𝑛+1−𝑗,𝑛+1−𝑖; it vanishes when 𝑛+1−𝑗>𝑛+1−𝑖, that is, when 𝑖>𝑗. So 𝑅′=𝑃𝑅𝖳𝑃 is upper triangular (this is Exercise 7.5.7 in matrix form), and with 𝑀=(𝑆−1)𝖳𝑃,

𝐴𝖳=(𝑆−1)𝖳𝑃𝑅′𝑃𝑆𝖳=𝑀𝑅′𝑀−1,

since 𝑀−1=𝑃−1𝑆𝖳=𝑃𝑆𝖳.

(b) False. 𝐴=(0100) is upper triangular, and 𝐵=(00−10) is lower triangular, hence triangularizable by (a). But 𝐴+𝐵=(01−10) is the rotation by −𝜋/2, which has no real eigenvalue and is not triangularizable over ℝ (as in Non-example 7.5.6).

(c) False. 𝑇=𝑇𝐷 with 𝐷=diag(1,2) is upper triangular (even diagonal) with respect to 𝐞1,𝐞2 and with respect to 𝐞2,𝐞1, and 𝐞2 is not a multiple of 𝐞1.

(d) False. 𝐼2 and 𝐽=(1101) have the same diagonal, but 𝑆−1𝐼2𝑆=𝐼2≠𝐽 for every invertible 𝑆, so 𝐼2 is similar only to itself.

Exercise 7.5.34.standard Let 𝐴∈ℂ𝑛×𝑛 have eigenvalues 𝜆1,…,𝜆𝑛, counted with multiplicity.

  1. (a)Show that for every 𝑘≥1, the eigenvalues of 𝐴𝑘 counted with multiplicity are 𝜆1𝑘,…,𝜆𝑛𝑘. Deduce that tr(𝐴𝑘)=𝜆1𝑘+⋯+𝜆𝑛𝑘.
  2. (b)Suppose 𝐴 is invertible. Show that the eigenvalues of 𝐴−1 counted with multiplicity are 𝜆1−1,…,𝜆𝑛−1, and that tr(𝐴−1)=𝜆1−1+⋯+𝜆𝑛−1.

Solution.

Uses. Corollary 7.5.14 (Triangularizing a matrix), Proposition 2.4.36 (Products of upper triangular matrices), Proposition 7.5.2 (Upper-triangular matrices and invariant subspaces), Proposition 7.5.9 (Multiplicities on the diagonal), Theorem 7.5.22 (Trace and determinant from eigenvalues), Proposition 7.5.8 (The diagonal of a triangular matrix)

Techniques. Build a flag of invariant subspaces (Technique 7.5.4), Prove it for triangular matrices first (Technique 7.5.23)

Hint 1. Write 𝐴=𝑆𝑅𝑆−1 with 𝑅 upper triangular. Then 𝐴𝑘=𝑆𝑅𝑘𝑆−1 and 𝐴−1=𝑆𝑅−1𝑆−1.
Hint 2. Powers of 𝑅 are upper triangular with the powers of the diagonal entries on the diagonal (Proposition 2.4.36). For 𝑅−1, use flags: 𝑇𝑅 maps each span(𝐞1,…,𝐞𝑗) into itself and is injective.
Sketch. (a) 𝑅𝑘 is upper triangular with diagonal 𝑟𝑗𝑗𝑘 by Proposition 2.4.36 and induction, and 𝐴𝑘=𝑆𝑅𝑘𝑆−1, so Proposition 7.5.9 applies. (b) 𝑇𝑅 maps 𝑉𝑗=span(𝐞1,…,𝐞𝑗) injectively into itself, hence onto it, so 𝑇𝑅−1 also preserves every 𝑉𝑗 and 𝑅−1 is upper triangular; 𝑅𝑅−1=𝐼𝑛 and Proposition 2.4.36 give diagonal entries 𝑟𝑗𝑗−1. Then take traces.

By Corollary 7.5.14, 𝐴=𝑆𝑅𝑆−1 with 𝑅 upper triangular and diagonal entries 𝑟𝑗𝑗=𝜆𝑗, after reordering the 𝜆𝑗.

(a) 𝐴𝑘=𝑆𝑅𝑘𝑆−1, by induction on 𝑘, since 𝑆−1𝑆=𝐼𝑛 cancels between consecutive factors. By Proposition 2.4.36 and induction on 𝑘, 𝑅𝑘 is upper triangular with diagonal entries 𝜆1𝑘,…,𝜆𝑛𝑘. So 𝑇𝐴𝑘 is upper triangular with respect to the columns of 𝑆, with this diagonal, and by Proposition 7.5.9 these are the eigenvalues of 𝐴𝑘 counted with multiplicity. By Theorem 7.5.22, tr(𝐴𝑘)=∑𝑗𝜆𝑗𝑘.

(b) 𝐴 is invertible, so 𝑅=𝑆−1𝐴𝑆 is too, and 𝐴−1=𝑆𝑅−1𝑆−1. Let 𝑉𝑗=span(𝐞1,…,𝐞𝑗). By Proposition 7.5.2, 𝑇𝑅 maps 𝑉𝑗 into 𝑉𝑗, and it is injective, so 𝑇𝑅↾𝑉𝑗 is an injective map from 𝑉𝑗 into itself, hence onto 𝑉𝑗 (Corollary 4.2.18). So each 𝐯∈𝑉𝑗 is 𝑇𝑅𝐰 for some 𝐰∈𝑉𝑗, and 𝑇𝑅−1𝐯=𝐰∈𝑉𝑗. Thus 𝑇𝑅−1=𝑇𝑅−1 preserves every 𝑉𝑗, and 𝑅−1 is upper triangular by Proposition 7.5.2. By Proposition 2.4.36, the diagonal entries of 𝑅𝑅−1=𝐼𝑛 are 𝑟𝑗𝑗(𝑅−1)𝑗𝑗=1, so (𝑅−1)𝑗𝑗=𝜆𝑗−1 (each 𝜆𝑗≠0 by Proposition 7.5.8). Now argue as in (a).

Exercise 7.5.35.standard A matrix 𝐴∈ℂ3×3 satisfies tr𝐴=2, tr(𝐴2)=6, and det𝐴=−2. Find the eigenvalues of 𝐴, and decide whether 𝐴 is diagonalizable.

Solution.

Uses. Theorem 7.5.22 (Trace and determinant from eigenvalues), Exercise 7.5.34, Proposition 7.5.9 (Multiplicities on the diagonal), Theorem 7.2.21 (Eigenvectors for distinct eigenvalues are independent), Theorem 7.4.2 (Diagonalizable means an eigenbasis)

Hint 1. Let 𝜆1,𝜆2,𝜆3 be the eigenvalues. You know ∑𝜆𝑖, ∑𝜆𝑖2, and ∏𝜆𝑖.
Hint 2. (𝜆1+𝜆2+𝜆3)2=∑𝜆𝑖2+2(𝜆1𝜆2+𝜆1𝜆3+𝜆2𝜆3). This gives every coefficient of 𝑝𝐴=(𝑥−𝜆1)(𝑥−𝜆2)(𝑥−𝜆3).
Sketch. 𝑒1=∑𝜆𝑖=2, ∑𝜆𝑖2=6 gives 𝑒2=∑𝑖<𝑗𝜆𝑖𝜆𝑗=(4−6)/2=−1, and 𝑒3=∏𝜆𝑖=−2. So 𝑝𝐴=𝑥3−2𝑥2−𝑥+2=(𝑥−1)(𝑥+1)(𝑥−2), with three distinct eigenvalues, so 𝐴 is diagonalizable.

Let 𝜆1,𝜆2,𝜆3 be the eigenvalues of 𝐴 counted with multiplicity. By Theorem 7.5.22 and Exercise 7.5.34(a),

𝜆1+𝜆2+𝜆3=2,𝜆12+𝜆22+𝜆32=6,𝜆1𝜆2𝜆3=−2.

Let 𝑒2=𝜆1𝜆2+𝜆1𝜆3+𝜆2𝜆3. Expanding the square, (𝜆1+𝜆2+𝜆3)2=∑𝜆𝑖2+2𝑒2, so 4=6+2𝑒2 and 𝑒2=−1. Expanding the product,

𝑝𝐴=(𝑥−𝜆1)(𝑥−𝜆2)(𝑥−𝜆3)=𝑥3−2𝑥2+𝑒2𝑥−𝜆1𝜆2𝜆3=𝑥3−2𝑥2−𝑥+2.

This factors as 𝑥2(𝑥−2)−(𝑥−2)=(𝑥−2)(𝑥−1)(𝑥+1), so the eigenvalues are 1, −1, and 2. (Check: 1+1+4=6.) Eigenvectors for the three distinct eigenvalues are linearly independent (Theorem 7.2.21), so they form a basis of ℂ3, and 𝐴 is diagonalizable (Theorem 7.4.2).

Exercise 7.5.36.standard
  1. (a)Let 𝐴∈ℂ𝑛×𝑛. Show that 0 is the only eigenvalue of 𝐴 if and only if 𝐴𝑛=0.
  2. (b)Find 𝐴∈ℝ3×3 whose only real eigenvalue is 0 but with 𝐴𝑘≠0 for every 𝑘≥1.

Solution.

Uses. Corollary 7.5.14 (Triangularizing a matrix), Proposition 7.5.8 (The diagonal of a triangular matrix), Example 2.8.3 (Strictly upper triangular matrices are nilpotent, three ways), Theorem 7.3.27 (Existence of eigenvalues over ℂ)

Techniques. Prove it for triangular matrices first (Technique 7.5.23)

Hint 1. If 0 is the only eigenvalue, a triangular form of 𝐴 has zeros on the diagonal.
Hint 2. Strictly upper-triangular 𝑛×𝑛 matrices satisfy 𝑁𝑛=0 (Example 2.8.3). For the converse, apply 𝐴𝑛 to an eigenvector.
Hint 3. For (b), the real eigenvalues do not see a rotation block.
Sketch. (a) 𝐴=𝑆𝑅𝑆−1 with 𝑅 upper triangular and diagonal the eigenvalues, all 0; so 𝑅 is strictly upper triangular, 𝑅𝑛=0, and 𝐴𝑛=𝑆𝑅𝑛𝑆−1=0. Conversely 𝐴𝐯=𝜆𝐯 gives 𝟎=𝐴𝑛𝐯=𝜆𝑛𝐯, so 𝜆=0, and an eigenvalue exists over ℂ. (b) 𝐴=(00000−1010).

(a) Suppose 0 is the only eigenvalue. By Corollary 7.5.14, 𝐴=𝑆𝑅𝑆−1 with 𝑅 upper triangular, and the diagonal entries of 𝑅 are eigenvalues of 𝐴 (Proposition 7.5.8), hence all 0. So 𝑅 is strictly upper triangular, and 𝑅𝑛=0 by Example 2.8.3. Then 𝐴𝑛=𝑆𝑅𝑛𝑆−1=0.

Conversely, suppose 𝐴𝑛=0, and let 𝐴𝐯=𝜆𝐯 with 𝐯≠𝟎. Then 𝐴𝑛𝐯=𝜆𝑛𝐯 (induction), so 𝜆𝑛𝐯=𝟎, 𝜆𝑛=0, and 𝜆=0. Since 𝐴 has at least one eigenvalue (Theorem 7.3.27), 0 is the only one.

(b) Let 𝐴=(00000−1010), which is block diagonal with blocks (0) and the rotation 𝑅=(0−110). A real eigenvector (𝑥,𝑦,𝑧) with eigenvalue 𝜆 satisfies 0=𝜆𝑥, −𝑧=𝜆𝑦, 𝑦=𝜆𝑧, so 𝑦=−𝜆2𝑦 and 𝑧=−𝜆2𝑧, which forces 𝑦=𝑧=0, and then 𝑥≠0 and 𝜆=0. So 0 is the only real eigenvalue. But 𝐴𝑘 is block diagonal with blocks (0) and 𝑅𝑘, and 𝑅𝑘≠0 since 𝑅 is invertible. Over ℂ, 𝐴 also has the eigenvalues ±𝑖, which is why (a) does not apply.

Exercise 7.5.37.standard Find the error in the following argument.

Claim. If 𝐴,𝐵∈ℂ𝑛×𝑛 have eigenvalues 𝜆1,…,𝜆𝑛 and 𝜇1,…,𝜇𝑛 (counted with multiplicity), then the eigenvalues of 𝐴𝐵 are 𝜆1𝜇1,…,𝜆𝑛𝜇𝑛, in some order.

Proof. By Schur’s theorem, 𝐴=𝑈𝑅𝑈∗ and 𝐵=𝑈𝑅′𝑈∗ with 𝑈 unitary and 𝑅,𝑅′ upper triangular, with diagonals 𝜆1,…,𝜆𝑛 and 𝜇1,…,𝜇𝑛. Then 𝐴𝐵=𝑈𝑅𝑅′𝑈∗, and 𝑅𝑅′ is upper triangular with diagonal 𝜆1𝜇1,…,𝜆𝑛𝜇𝑛.

Solution.

Uses. Theorem 7.5.18 (Schur’s theorem), Theorem 5.4.10 (Multiplicativity), Theorem 7.5.22 (Trace and determinant from eigenvalues)

Techniques. Finding counterexamples and minimal examples (Technique 1.4.18)

Hint 1. Schur’s theorem gives a unitary 𝑈 for each matrix separately. Does it give one 𝑈 that works for both?
Hint 2. Try 𝐴=(0100) and 𝐵=(0010).
Sketch. The proof uses the same 𝑈 for 𝐴 and 𝐵, but Schur’s theorem only provides a 𝑈 depending on the matrix. The claim is false: 𝐴=(0100) and 𝐵=(0010) have only the eigenvalue 0, but 𝐴𝐵=(1000) has the eigenvalue 1.

Schur’s theorem applied to 𝐴 gives some unitary 𝑈, and applied to 𝐵 gives some unitary 𝑈′, which in general differs from 𝑈. The argument silently takes 𝑈′=𝑈. A common triangularizing basis exists only in special situations, for instance for commuting matrices (Section 7.6).

The claim itself is false. Take 𝐴=(0100) and 𝐵=(0010); both have 0 as their only eigenvalue, with multiplicity 2. But 𝐴𝐵=(1000) has the eigenvalues 1 and 0, while the claim predicts 0,0. What survives is the product of all of them: det(𝐴𝐵)=det𝐴det𝐵 (Theorem 5.4.10), so by Theorem 7.5.22 the product of the eigenvalues of 𝐴𝐵 is ∏𝜆𝑖∏𝜇𝑖.

Exercise 7.5.38.hard Let 𝑉 be a complex vector space with dim𝑉=𝑛≥1, and let 𝑇∈L(𝑉) have eigenvalues 𝜆1,…,𝜆𝑛, counted with multiplicity and listed in any order chosen in advance. Show that there is a basis ℬ︀ of 𝑉 such that [𝑇]ℬ︀ is upper triangular with diagonal entries 𝜆1,𝜆2,…,𝜆𝑛 in this order. If 𝑉 is an inner product space, show that ℬ︀ can be taken orthonormal.

Solution.

Uses. Theorem 7.3.6 (Eigenvalues are the roots of the characteristic polynomial), Lemma 7.5.11 (Compressing an operator to a complement), Theorem 5.4.26 (Block-triangular matrices), Corollary 7.1.7 (No zero divisors, cancellation, and units), Proposition 4.4.16 (Complements exist), Theorem 6.4.24 (Orthogonal decomposition)

Techniques. Induction on dimension via an invariant subspace (Technique 7.5.12)

Hint 1. Follow the first proof of Theorem 7.5.13, but choose the first eigenvector for 𝜆1.
Hint 2. For the induction, you need the eigenvalues of the compression 𝑆. Show 𝑝𝑇=(𝑥−𝜆1)𝑝𝑆 from the block form of Lemma 7.5.11(b).
Hint 3. For the orthonormal version, use 𝑊=span(𝐛1)⟂ with 𝐛1 a unit vector, as in the second proof of Theorem 7.5.18.
Sketch. Induct on 𝑛. 𝜆1 is a root of 𝑝𝑇, hence an eigenvalue; take an eigenvector 𝐛1, a complement 𝑊, and the compression 𝑆. The block form [𝑇]=(𝜆1𝑋0[𝑆]) gives 𝑝𝑇=(𝑥−𝜆1)𝑝𝑆, so 𝑝𝑆=(𝑥−𝜆2)⋯(𝑥−𝜆𝑛). By induction 𝑊 has a basis triangularizing 𝑆 with diagonal 𝜆2,…,𝜆𝑛, and the block form gives the diagonal 𝜆1,…,𝜆𝑛 for 𝑇. With a unit 𝐛1 and 𝑊=span(𝐛1)⟂ everything can be made orthonormal.

We induct on 𝑛; for 𝑛=1 there is nothing to show. Let 𝑛≥2. Since 𝑝𝑇=(𝑥−𝜆1)⋯(𝑥−𝜆𝑛), 𝜆1 is a root of 𝑝𝑇, hence an eigenvalue (Theorem 7.3.6). Let 𝐛1 be an eigenvector for 𝜆1, 𝑊 a complement of span(𝐛1) (Proposition 4.4.16), and 𝑆 the compression of 𝑇 to 𝑊. For a basis 𝒟︀ of 𝑊 and ℬ︀=𝐛1,𝒟︀, Lemma 7.5.11(b) gives

[𝑇]ℬ︀=(𝜆1𝑋0[𝑆]𝒟︀),

so, taking det(𝑥𝐼𝑛−⋅) with Theorem 5.4.26, 𝑝𝑇=(𝑥−𝜆1)𝑝𝑆. Cancelling 𝑥−𝜆1 from (𝑥−𝜆1)𝑝𝑆=(𝑥−𝜆1)(𝑥−𝜆2)⋯(𝑥−𝜆𝑛) (Corollary 7.1.7) gives 𝑝𝑆=(𝑥−𝜆2)⋯(𝑥−𝜆𝑛): the eigenvalues of 𝑆 counted with multiplicity are 𝜆2,…,𝜆𝑛. By the induction hypothesis, 𝑊 has a basis 𝒟︀=𝐛2,…,𝐛𝑛 with [𝑆]𝒟︀ upper triangular with diagonal 𝜆2,…,𝜆𝑛 in this order. Then the displayed block matrix is upper triangular with diagonal 𝜆1,𝜆2,…,𝜆𝑛.

If 𝑉 is an inner product space, choose 𝐛1 of norm 1 and 𝑊=span(𝐛1)⟂, a complement by Theorem 6.4.24, and apply the inner product version of the induction hypothesis to 𝑊. Then 𝐛1,𝐛2,…,𝐛𝑛 is orthonormal.

Exercise 7.5.39.standardBeyond 21-242 Let 𝐴∈ℂ𝑛×𝑛, and write 𝐴=𝑆𝑅𝑆−1 as in Corollary 7.5.14. Let 𝑆=𝑄𝑅′ be the QR factorization of 𝑆 (Theorem 6.4.38). Show that 𝐴=𝑄𝑅″𝑄∗ with 𝑅″=𝑅′𝑅𝑅′−1 upper triangular with the same diagonal as 𝑅. Explain why this is the first proof of Theorem 7.5.18 in matrix form.

Solution.

Uses. Theorem 6.4.38 (QR factorization), Corollary 7.5.14 (Triangularizing a matrix), Proposition 2.4.36 (Products of upper triangular matrices), Exercise 7.5.34, Definition 6.3.33 (Unitary and orthogonal matrices)

Hint 1. 𝑆 is invertible, so its columns are independent and Theorem 6.4.38 applies, with 𝑄 square.
Hint 2. Substitute 𝑆=𝑄𝑅′ into 𝐴=𝑆𝑅𝑆−1. Inverses and products of invertible upper-triangular matrices are upper triangular (Exercise 7.5.34(b), Proposition 2.4.36).
Sketch. 𝐴=𝑄𝑅′𝑅𝑅′−1𝑄−1 and 𝑄−1=𝑄∗. The product 𝑅′𝑅𝑅′−1 is upper triangular with diagonal 𝑟𝑗𝑗′𝑟𝑗𝑗𝑟𝑗𝑗′−1=𝑟𝑗𝑗. The columns of 𝑄 are Gram–Schmidt applied to the columns of 𝑆, the triangularizing basis, which is exactly the first proof.

The columns of 𝑆 form a basis of ℂ𝑛, so by Theorem 6.4.38, 𝑆=𝑄𝑅′ with 𝑄∈ℂ𝑛×𝑛 having orthonormal columns, hence unitary (Definition 6.3.33), and 𝑅′ upper triangular with positive diagonal entries, hence invertible. Then

𝐴=𝑆𝑅𝑆−1=𝑄𝑅′𝑅𝑅′−1𝑄−1=𝑄𝑅″𝑄∗.

By Exercise 7.5.34(b) (whose argument applies to any invertible upper-triangular matrix), 𝑅′−1 is upper triangular with diagonal entries 1/𝑟𝑗𝑗′, and by Proposition 2.4.36, 𝑅″=𝑅′𝑅𝑅′−1 is upper triangular with diagonal entries 𝑟𝑗𝑗′𝑟𝑗𝑗(1/𝑟𝑗𝑗′)=𝑟𝑗𝑗.

By Theorem 6.4.38, the columns of 𝑄 are the output of Gram–Schmidt applied to the columns of 𝑆, which form a basis triangularizing 𝑇𝐴. So this is the first proof of Theorem 7.5.18: Gram–Schmidt applied to a triangularizing basis. The triangular factor 𝑅′ records the change of basis between the two bases with the same flag, and conjugating by it keeps 𝑅 triangular.

Exercise 7.5.40.standard
  1. (a)Show that an upper-triangular unitary matrix 𝑅∈ℂ𝑛×𝑛 is diagonal, with |𝑟𝑗𝑗|=1 for every 𝑗.
  2. (b)Deduce that every unitary 𝐴∈ℂ𝑛×𝑛 can be written as 𝐴=𝑈𝐷𝑈∗ with 𝑈 unitary and 𝐷 diagonal.

Solution.

Uses. Definition 6.3.33 (Unitary and orthogonal matrices), Exercise 7.5.34, Theorem 7.5.18 (Schur’s theorem), Exercise 6.3.37, Proposition 6.3.27 (Properties of the conjugate transpose)

Hint 1. 𝑅−1=𝑅∗. One side is upper triangular and the other lower triangular.
Hint 2. Alternatively, look at the columns: the first column (𝑟11,0,…,0) is a unit vector, and every other column is orthogonal to it.
Hint 3. For (b), apply Schur’s theorem and check that 𝑅=𝑈∗𝐴𝑈 is unitary.
Sketch. (a) 𝑅∗=𝑅−1 is upper triangular (Exercise 7.5.34(b)) and lower triangular (as the conjugate transpose of an upper-triangular matrix), so it is diagonal, and so is 𝑅; then 𝑅∗𝑅=𝐼𝑛 gives |𝑟𝑗𝑗|2=1. (b) Schur gives 𝐴=𝑈𝑅𝑈∗, and 𝑅=𝑈∗𝐴𝑈 is a product of unitary matrices.

(a) Since 𝑅 is unitary, 𝑅−1=𝑅∗. By (the argument of) Exercise 7.5.34(b), the inverse of an invertible upper-triangular matrix is upper triangular, so 𝑅∗ is upper triangular. But the (𝑖,𝑗) entry of 𝑅∗ is 𝑟𝑗𝑖, which is 0 for 𝑖<𝑗, so 𝑅∗ is also lower triangular. Hence 𝑅∗ is diagonal, and so is 𝑅=(𝑅∗)∗. Finally the (𝑗,𝑗) entry of 𝑅∗𝑅=𝐼𝑛 is 𝑟𝑗𝑗𝑟𝑗𝑗=|𝑟𝑗𝑗|2=1.

(b) By Theorem 7.5.18, 𝐴=𝑈𝑅𝑈∗ with 𝑈 unitary and 𝑅 upper triangular. Then 𝑅=𝑈∗𝐴𝑈 is a product of unitary matrices, so it is unitary (Exercise 6.3.37). By (a), 𝑅 is diagonal; take 𝐷=𝑅.

Exercise 7.5.41.hard Let 𝐴∈ℂ𝑛×𝑛 have eigenvalues 𝜆1,…,𝜆𝑛, counted with multiplicity. Prove Schur’s inequality

∑𝑗=1𝑛|𝜆𝑗|2≤∑𝑖,𝑗=1𝑛|𝑎𝑖𝑗|2,

with equality if and only if 𝐴=𝑈𝐷𝑈∗ for some unitary 𝑈 and diagonal 𝐷.

Solution.

Uses. Theorem 7.5.18 (Schur’s theorem), Proposition 6.1.16 (Frobenius inner product), Proposition 2.4.28 (Cyclicity of the trace), Proposition 6.3.27 (Properties of the conjugate transpose), Proposition 7.5.8 (The diagonal of a triangular matrix)

Techniques. The trace of 𝐴∗𝐴 (Technique 6.1.18), Cyclicity of trace (Technique 2.4.29)

Hint 1. The right side is tr(𝐴∗𝐴) (Proposition 6.1.16). How does it change when 𝐴 is replaced by 𝑈∗𝐴𝑈?
Hint 2. For a Schur form 𝑅=𝑈∗𝐴𝑈, tr(𝐴∗𝐴)=tr(𝑅∗𝑅)=∑𝑖,𝑗|𝑟𝑖𝑗|2, and the diagonal of 𝑅 lists the eigenvalues.
Sketch. With 𝐴=𝑈𝑅𝑈∗, 𝐴∗𝐴=𝑈𝑅∗𝑅𝑈∗, so tr(𝐴∗𝐴)=tr(𝑅∗𝑅𝑈∗𝑈)=tr(𝑅∗𝑅). Hence ∑|𝑎𝑖𝑗|2=∑|𝑟𝑖𝑗|2=∑𝑗|𝜆𝑗|2+∑𝑖<𝑗|𝑟𝑖𝑗|2. Equality holds exactly when 𝑅 is diagonal.

By Theorem 7.5.18, 𝐴=𝑈𝑅𝑈∗ with 𝑈 unitary and 𝑅 upper triangular with diagonal 𝜆1,…,𝜆𝑛 (in some order). By Proposition 6.3.27, 𝐴∗=𝑈𝑅∗𝑈∗, so 𝐴∗𝐴=𝑈𝑅∗𝑈∗𝑈𝑅𝑈∗=𝑈𝑅∗𝑅𝑈∗. By Proposition 6.1.16 and Proposition 2.4.28,

∑𝑖,𝑗|𝑎𝑖𝑗|2=tr(𝐴∗𝐴)=tr(𝑈(𝑅∗𝑅𝑈∗))=tr(𝑅∗𝑅𝑈∗𝑈)=tr(𝑅∗𝑅)=∑𝑖,𝑗|𝑟𝑖𝑗|2.

Since 𝑅 is upper triangular,

∑𝑖,𝑗|𝑟𝑖𝑗|2=∑𝑗=1𝑛|𝜆𝑗|2+∑𝑖<𝑗|𝑟𝑖𝑗|2≥∑𝑗=1𝑛|𝜆𝑗|2,

which is the inequality.

If equality holds, then 𝑟𝑖𝑗=0 for all 𝑖<𝑗, so 𝑅 is diagonal and 𝐴=𝑈𝑅𝑈∗ has the required form. Conversely, if 𝐴=𝑈𝐷𝑈∗ with 𝐷 diagonal, then 𝐷 is upper triangular, its diagonal lists the eigenvalues of 𝐴 (Proposition 7.5.8, Proposition 7.5.9), and the computation above with 𝑅=𝐷 gives equality. (In Section 8.2 such matrices are identified as the normal ones, 𝐴∗𝐴=𝐴𝐴∗.)

Exercise 7.5.42.hard Let 𝐴∈ℝ𝑛×𝑛, and suppose that 𝑝𝐴 splits over ℝ, that is, all eigenvalues of 𝐴 are real. Show that 𝐴=𝑄𝑅𝑄𝖳 with 𝑄∈ℝ𝑛×𝑛 orthogonal and 𝑅∈ℝ𝑛×𝑛 upper triangular.

Solution.

Uses. Theorem 7.3.6 (Eigenvalues are the roots of the characteristic polynomial), Lemma 7.5.11 (Compressing an operator to a complement), Theorem 5.4.26 (Block-triangular matrices), Corollary 7.1.7 (No zero divisors, cancellation, and units), Theorem 6.4.24 (Orthogonal decomposition), Proposition 7.5.2 (Upper-triangular matrices and invariant subspaces), Theorem 4.5.17 (Change of basis)

Techniques. Induction on dimension via an invariant subspace (Technique 7.5.12)

Hint 1. The second proof of Theorem 7.5.18 only used ℂ to find an eigenvalue. Over ℝ, with a split characteristic polynomial, an eigenvalue exists. What has to be checked for the compression?
Hint 2. Show that the compression 𝑆 to span(𝐮1)⟂ satisfies 𝑝𝐴=(𝑥−𝜆1)𝑝𝑆, so 𝑝𝑆 splits too.
Sketch. Prove the operator version on a real inner product space by induction on the dimension: a real root 𝜆1 of 𝑝𝑇 gives a unit eigenvector 𝐮1; the compression 𝑆 to 𝑊=span(𝐮1)⟂ has 𝑝𝑇=(𝑥−𝜆1)𝑝𝑆 by Lemma 7.5.11(b), so 𝑝𝑆 splits; by induction 𝑊 has an orthonormal basis triangularizing 𝑆, and 𝐮1 followed by it triangularizes 𝑇. Apply this to 𝑇𝐴 on ℝ𝑛.

We prove: if 𝑉 is a real inner product space with 1≤dim𝑉=𝑛<∞ and 𝑇∈L(𝑉) has 𝑝𝑇 splitting over ℝ, then 𝑉 has an orthonormal basis triangularizing 𝑇. Applied to 𝑇𝐴 on ℝ𝑛 with the dot product, with 𝑄 the matrix whose columns are that basis, this gives 𝑄𝖳𝐴𝑄=𝑅 upper triangular by Theorem 4.5.17, since 𝑄−1=𝑄𝖳 for orthogonal 𝑄.

Induct on 𝑛; for 𝑛=1 any unit vector works. Let 𝑛≥2 and 𝑝𝑇=(𝑥−𝜇1)⋯(𝑥−𝜇𝑛) with real 𝜇𝑖. Then 𝜇1 is an eigenvalue (Theorem 7.3.6); let 𝐮1 be an eigenvector of norm 1. Let 𝑊=span(𝐮1)⟂, a complement of span(𝐮1) of dimension 𝑛−1 (Theorem 6.4.24), and let 𝑆𝐰=𝑇𝐰−⟨𝑇𝐰,𝐮1⟩𝐮1, the compression of 𝑇 to 𝑊 (as in the second proof of Theorem 7.5.18). For any basis 𝒟︀ of 𝑊, Lemma 7.5.11(b) and Theorem 5.4.26 give 𝑝𝑇=(𝑥−𝜇1)𝑝𝑆, and cancelling (Corollary 7.1.7) gives 𝑝𝑆=(𝑥−𝜇2)⋯(𝑥−𝜇𝑛), which splits. By the induction hypothesis, 𝑊 has an orthonormal basis 𝐮2,…,𝐮𝑛 with 𝑆𝐮𝑗∈span(𝐮2,…,𝐮𝑗). Exactly as in the second proof of Theorem 7.5.18, 𝐮1,…,𝐮𝑛 is an orthonormal basis of 𝑉 and 𝑇𝐮𝑗=⟨𝑇𝐮𝑗,𝐮1⟩𝐮1+𝑆𝐮𝑗∈span(𝐮1,…,𝐮𝑗), so 𝑇 is upper triangular with respect to it (Proposition 7.5.2).

Exercise 7.5.43.hard Let 𝐴∈ℂ𝑛×𝑛 with tr(𝐴𝑘)=0 for 𝑘=1,2,…,𝑛. Show that 𝐴𝑛=0. Then show that the statement fails over 𝔽𝑝.

Solution.

Uses. Exercise 7.5.34, Exercise 7.5.36, Theorem 5.5.21 (The Vandermonde determinant), Corollary 3.5.25 (The invertible matrix theorem), Theorem 5.4.14 (The transpose)

Hint 1. Group the nonzero eigenvalues: let 𝜇1,…,𝜇𝑟 be the distinct nonzero ones, with algebraic multiplicities 𝑚1,…,𝑚𝑟≥1. By Exercise 7.5.36 it suffices to show 𝑟=0.
Hint 2. By Exercise 7.5.34(a), ∑𝑗𝑚𝑗𝜇𝑗𝑘=0 for 𝑘=1,…,𝑟. This is a square linear system for (𝑚1,…,𝑚𝑟).
Hint 3. Its matrix is a Vandermonde matrix times diag(𝜇1,…,𝜇𝑟). For 𝔽𝑝, try 𝐼𝑝.
Sketch. If there are distinct nonzero eigenvalues 𝜇1,…,𝜇𝑟 with multiplicities 𝑚𝑗, then ∑𝑗𝑚𝑗𝜇𝑗𝑘=tr(𝐴𝑘)=0 for 𝑘=1,…,𝑟. The coefficient matrix (𝜇𝑗𝑘)𝑘,𝑗 has determinant 𝜇1⋯𝜇𝑟∏𝑖<𝑗(𝜇𝑗−𝜇𝑖)≠0, so all 𝑚𝑗=0 in ℂ, which is absurd. So 0 is the only eigenvalue, and 𝐴𝑛=0. Over 𝔽𝑝, tr(𝐼𝑝𝑘)=𝑝∗1=0 for all 𝑘, but 𝐼𝑝𝑝≠0.

Let 𝜇1,…,𝜇𝑟 be the distinct nonzero eigenvalues of 𝐴, with algebraic multiplicities 𝑚1,…,𝑚𝑟≥1, and suppose 𝑟≥1. The eigenvalues of 𝐴 counted with multiplicity are these, each repeated 𝑚𝑗 times, together with 0 repeated 𝑛−∑𝑚𝑗 times. By Exercise 7.5.34(a), for 𝑘=1,…,𝑟 (note 𝑟≤𝑛),

0=tr(𝐴𝑘)=𝑚1𝜇1𝑘+𝑚2𝜇2𝑘+⋯+𝑚𝑟𝜇𝑟𝑘,

where 𝑚𝑗 stands for 𝑚𝑗∗1∈ℂ. So 𝑀𝐜=𝟎, where 𝐜=(𝑚1,…,𝑚𝑟) and 𝑀∈ℂ𝑟×𝑟 has entries 𝑀𝑘𝑗=𝜇𝑗𝑘. Factoring 𝜇𝑗 out of column 𝑗, 𝑀=𝑉′diag(𝜇1,…,𝜇𝑟) with 𝑉𝑘𝑗′=𝜇𝑗𝑘−1, the transpose of the Vandermonde matrix of 𝜇1,…,𝜇𝑟. By Theorem 5.5.21 and Theorem 5.4.14,

det𝑀=𝜇1𝜇2⋯𝜇𝑟∏𝑖<𝑗(𝜇𝑗−𝜇𝑖)≠0,

since the 𝜇𝑗 are nonzero and distinct. So 𝑀 is invertible (Corollary 3.5.25), 𝐜=𝟎, and 𝑚1=0 in ℂ, contradicting 𝑚1≥1. Hence 𝑟=0: 0 is the only eigenvalue of 𝐴, and 𝐴𝑛=0 by Exercise 7.5.36.

Over 𝔽𝑝, take 𝐴=𝐼𝑝∈𝔽𝑝𝑝×𝑝. Then tr(𝐴𝑘)=tr𝐼𝑝=𝑝∗1=0 for every 𝑘, but 𝐴𝑝=𝐼𝑝≠0. The proof breaks at its last step: in 𝔽𝑝, the conclusion 𝑚1∗1=0 is no contradiction when 𝑝 divides 𝑚1, as here, where 𝑚1=𝑝.

Exercise 7.5.44.hard Let 𝐴,𝐵∈ℂ𝑛×𝑛 satisfy 𝐴𝐵−𝐵𝐴=𝐴. Show that 𝐴𝑛=0.

Solution (traces of powers).

Uses. Proposition 2.4.28 (Cyclicity of the trace), Proposition 2.4.27 (The trace is linear), Exercise 7.5.43

Techniques. Cyclicity of trace (Technique 2.4.29)

Hint 1. Compute 𝐴2𝐵−𝐵𝐴2 in terms of 𝐴2. Guess a formula for 𝐴𝑘𝐵−𝐵𝐴𝑘.
Hint 2. Prove 𝐴𝑘𝐵−𝐵𝐴𝑘=𝑘𝐴𝑘 by induction, then take the trace.
Sketch. By induction, 𝐴𝑘+1𝐵−𝐵𝐴𝑘+1=𝐴(𝐴𝑘𝐵−𝐵𝐴𝑘)+(𝐴𝐵−𝐵𝐴)𝐴𝑘=𝑘𝐴𝑘+1+𝐴𝑘+1, so 𝐴𝑘𝐵−𝐵𝐴𝑘=𝑘𝐴𝑘. Taking traces, 𝑘tr(𝐴𝑘)=0, so tr(𝐴𝑘)=0 for all 𝑘≥1, and Exercise 7.5.43 gives 𝐴𝑛=0.

We show 𝐴𝑘𝐵−𝐵𝐴𝑘=𝑘𝐴𝑘 for 𝑘≥1 by induction. For 𝑘=1 this is the hypothesis. If it holds for 𝑘, then

𝐴𝑘+1𝐵−𝐵𝐴𝑘+1=𝐴(𝐴𝑘𝐵−𝐵𝐴𝑘)+(𝐴𝐵−𝐵𝐴)𝐴𝑘=𝐴(𝑘𝐴𝑘)+𝐴𝐴𝑘=(𝑘+1)𝐴𝑘+1.

Taking traces, using Proposition 2.4.27 and Proposition 2.4.28, 𝑘tr(𝐴𝑘)=tr(𝐴𝑘𝐵)−tr(𝐵𝐴𝑘)=0, so tr(𝐴𝑘)=0 for every 𝑘≥1, since 𝑘≠0 in ℂ. By Exercise 7.5.43, 𝐴𝑛=0.

Solution (eigenvalues of a map on matrices).

Uses. Corollary 7.2.22 (At most dim𝑉 eigenvalues), Corollary 4.1.27 (A nilpotent operator dies within dim𝑉 steps)

Techniques. Study a linear map on a space of matrices (Technique 4.2.29), Cyclicity of trace (Technique 2.4.29)

Hint 1. Consider the operator 𝐿(𝑋)=𝐵𝑋−𝑋𝐵 on ℂ𝑛×𝑛. What does it do to 𝐴, to 𝐴2, to 𝐴𝑘?
Hint 2. If no power of 𝐴 were 0, 𝐿 would have infinitely many eigenvalues.
Sketch. The identity 𝐴𝑘𝐵−𝐵𝐴𝑘=𝑘𝐴𝑘 says 𝐿(𝐴𝑘)=−𝑘𝐴𝑘 for 𝐿(𝑋)=𝐵𝑋−𝑋𝐵. If 𝐴𝑘≠0 for every 𝑘, then −1,−2,−3,… are all eigenvalues of 𝐿, an operator on the 𝑛2-dimensional space ℂ𝑛×𝑛, which has at most 𝑛2 eigenvalues. So some 𝐴𝑘=0, and then 𝐴𝑛=0.
As in the first approach, 𝐴𝑘𝐵−𝐵𝐴𝑘=𝑘𝐴𝑘 for all 𝑘≥1. Let 𝐿∈L(ℂ𝑛×𝑛) be 𝐿(𝑋)=𝐵𝑋−𝑋𝐵. Then 𝐿(𝐴𝑘)=−𝑘𝐴𝑘. If 𝐴𝑘≠0 for every 𝑘≥1, then −1,−2,…,−(𝑛2+1) are 𝑛2+1 distinct eigenvalues of 𝐿, which is impossible, since dimℂ𝑛×𝑛=𝑛2 (Corollary 7.2.22). So 𝐴𝑘=0 for some 𝑘, and then 𝐴𝑛=0 (Corollary 4.1.27).

Exercise 7.5.45.hard Let 𝐴∈ℂ𝑛×𝑛 have eigenvalues 𝜆1,…,𝜆𝑛, counted with multiplicity, and let 𝐿𝐴∈L(ℂ𝑛×𝑛) be 𝐿𝐴(𝑋)=𝐴𝑋−𝑋𝐴. Use the matrix units 𝑈𝑖𝑗=𝐞𝑖𝐞𝑗𝖳 (Lemma 2.4.24).

  1. (a)Suppose 𝐴 is upper triangular. List the matrix units row by row, and within each row from right to left: 𝑈1𝑛,𝑈1,𝑛−1,…,𝑈11,𝑈2𝑛,…,𝑈21,…,𝑈𝑛1. Show that 𝐿𝐴 is upper triangular with respect to this basis, with diagonal entry 𝑎𝑖𝑖−𝑎𝑗𝑗 at 𝑈𝑖𝑗.
  2. (b)Deduce that for every 𝐴, the eigenvalues of 𝐿𝐴 counted with multiplicity are the 𝑛2 numbers 𝜆𝑖−𝜆𝑗 (1≤𝑖,𝑗≤𝑛).
  3. (c)Deduce again that tr𝐿𝐴=0, and that 𝐿𝐴 is not invertible.

Solution.

Uses. Lemma 2.4.24 (Matrix units), Proposition 7.5.2 (Upper-triangular matrices and invariant subspaces), Proposition 7.5.9 (Multiplicities on the diagonal), Corollary 7.5.14 (Triangularizing a matrix), Theorem 7.5.22 (Trace and determinant from eigenvalues)

Techniques. Study a linear map on a space of matrices (Technique 4.2.29), Prove it for triangular matrices first (Technique 7.5.23)

Hint 1. Compute 𝐴𝑈𝑖𝑗 (column 𝑖 of 𝐴, moved to column 𝑗) and 𝑈𝑖𝑗𝐴 (row 𝑗 of 𝐴, moved to row 𝑖), using that 𝐴 is upper triangular.
Hint 2. Which matrix units appear in 𝐿𝐴(𝑈𝑖𝑗) besides 𝑈𝑖𝑗? Do they come before 𝑈𝑖𝑗 in the given order?
Hint 3. For (b), if 𝐴=𝑆𝑅𝑆−1, then 𝐿𝐴 is conjugate to 𝐿𝑅 by the invertible operator 𝑋↦𝑆𝑋𝑆−1.
Sketch. (a) 𝐴𝑈𝑖𝑗=∑𝑘≤𝑖𝑎𝑘𝑖𝑈𝑘𝑗 and 𝑈𝑖𝑗𝐴=∑𝑙≥𝑗𝑎𝑗𝑙𝑈𝑖𝑙, so 𝐿𝐴(𝑈𝑖𝑗)=(𝑎𝑖𝑖−𝑎𝑗𝑗)𝑈𝑖𝑗 plus units 𝑈𝑘𝑗 with 𝑘<𝑖 (earlier rows) and 𝑈𝑖𝑙 with 𝑙>𝑗 (earlier in the same row). (b) With 𝐴=𝑆𝑅𝑆−1 and 𝐶(𝑋)=𝑆𝑋𝑆−1, 𝐿𝐴∘𝐶=𝐶∘𝐿𝑅, so 𝐿𝐴 is upper triangular in the basis 𝐶(𝑈𝑖𝑗) with diagonal 𝑟𝑖𝑖−𝑟𝑗𝑗. (c) The sum of all 𝜆𝑖−𝜆𝑗 is 0, and the terms with 𝑖=𝑗 are 0, so det𝐿𝐴=0.

(a) By Lemma 2.4.24, 𝐴𝑈𝑖𝑗=(𝐴𝐞𝑖)𝐞𝑗𝖳=∑𝑘𝑎𝑘𝑖𝑈𝑘𝑗 and 𝑈𝑖𝑗𝐴=𝐞𝑖(𝐞𝑗𝖳𝐴)=∑𝑙𝑎𝑗𝑙𝑈𝑖𝑙. Since 𝐴 is upper triangular, 𝑎𝑘𝑖=0 for 𝑘>𝑖 and 𝑎𝑗𝑙=0 for 𝑙<𝑗, so

𝐿𝐴(𝑈𝑖𝑗)=(𝑎𝑖𝑖−𝑎𝑗𝑗)𝑈𝑖𝑗+∑𝑘<𝑖𝑎𝑘𝑖𝑈𝑘𝑗−∑𝑙>𝑗𝑎𝑗𝑙𝑈𝑖𝑙.

In the given order, every 𝑈𝑘𝑗 with 𝑘<𝑖 lies in an earlier row, and every 𝑈𝑖𝑙 with 𝑙>𝑗 lies earlier in the same row. So 𝐿𝐴(𝑈𝑖𝑗) lies in the span of 𝑈𝑖𝑗 and the units before it, and by Proposition 7.5.2, 𝐿𝐴 is upper triangular in this basis, with diagonal entry 𝑎𝑖𝑖−𝑎𝑗𝑗 at 𝑈𝑖𝑗.

(b) By Corollary 7.5.14, 𝐴=𝑆𝑅𝑆−1 with 𝑅 upper triangular, 𝑟𝑖𝑖=𝜆𝑖 after reordering. Let 𝐶∈L(ℂ𝑛×𝑛) be 𝐶(𝑋)=𝑆𝑋𝑆−1, which is invertible with inverse 𝑋↦𝑆−1𝑋𝑆. For every 𝑋,

𝐿𝐴(𝐶(𝑋))=𝑆𝑅𝑆−1𝑆𝑋𝑆−1−𝑆𝑋𝑆−1𝑆𝑅𝑆−1=𝑆(𝑅𝑋−𝑋𝑅)𝑆−1=𝐶(𝐿𝑅(𝑋)).

So if 𝒰︀ is the ordered basis of (a), then 𝐶(𝒰︀) is a basis (the image of a basis under an isomorphism), and 𝐿𝐴 sends 𝐶(𝑈𝑖𝑗) to 𝐶(𝐿𝑅(𝑈𝑖𝑗)), which has the same coordinates in 𝐶(𝒰︀) as 𝐿𝑅(𝑈𝑖𝑗) has in 𝒰︀. Hence [𝐿𝐴]𝐶(𝒰︀)=[𝐿𝑅]𝒰︀ is upper triangular with diagonal entries 𝜆𝑖−𝜆𝑗, and by Proposition 7.5.9 these are the eigenvalues of 𝐿𝐴 counted with multiplicity.

(c) By Theorem 7.5.22, tr𝐿𝐴=∑𝑖,𝑗(𝜆𝑖−𝜆𝑗)=𝑛∑𝑖𝜆𝑖−𝑛∑𝑗𝜆𝑗=0, and det𝐿𝐴=∏𝑖,𝑗(𝜆𝑖−𝜆𝑗)=0, because the factors with 𝑖=𝑗 vanish. So 𝐿𝐴 is not invertible. (Directly: 𝐿𝐴(𝐼𝑛)=0.)

Exercise 7.5.46.standardBeyond 21-242 Let 𝐴=(02−10101−10)∈ℝ3×3.

  1. (a)Check that (1,0,−𝑖) is an eigenvector of 𝐴 over ℂ, and use it to find a two-dimensional 𝑇𝐴-invariant subspace of ℝ3.
  2. (b)Find a basis of ℝ3 in which 𝑇𝐴 has the form of Theorem 7.5.27(a).
  3. (c)Find an orthogonal 𝑄 such that 𝑄𝖳𝐴𝑄 is a real Schur form.

Solution.

Uses. Theorem 7.5.27 (Real block upper-triangular form), Lemma 7.5.11 (Compressing an operator to a complement)

Hint 1. Split (1,0,−𝑖) into real and imaginary parts 𝐱+𝑖𝐲 and compute 𝐴𝐱 and 𝐴𝐲.
Hint 2. For (b), add a real eigenvector. For (c), the invariant plane already has an orthonormal basis; compress to its orthogonal complement.
Sketch. 𝐴(1,0,−𝑖)=(𝑖,0,1)=𝑖(1,0,−𝑖). With 𝐱=𝐞1, 𝐲=−𝐞3: 𝐴𝐱=−𝐲, 𝐴𝐲=𝐱, so span(𝐞1,𝐞3) is invariant. (b) (3,2,1) is an eigenvector for 1, and in the basis 𝐱,𝐲,(3,2,1) the matrix is (010−100001). (c) 𝑄 with columns 𝐞1,−𝐞3,𝐞2 gives 𝑄𝖳𝐴𝑄=(012−101001).

(a) 𝐴(1,0,−𝑖)=(0+0+𝑖,0,1+0+0)=(𝑖,0,1)=𝑖(1,0,−𝑖), since 𝑖(−𝑖)=1. Write (1,0,−𝑖)=𝐱+𝑖𝐲 with 𝐱=(1,0,0) and 𝐲=(0,0,−1). Then

𝐴𝐱=(0,0,1)=−𝐲,𝐴𝐲=(1,0,0)=𝐱,

as in the complexification argument before Theorem 7.5.27, with 𝑎=0 and 𝑏=1. So span(𝐱,𝐲)=span(𝐞1,𝐞3) is invariant, and 𝐴 acts on it by (01−10) in the basis 𝐱,𝐲, a block with eigenvalues ±𝑖.

(b) For the eigenvalue 1, 𝐴−𝐼3=(−12−10001−1−1); adding the first and last rows gives 𝑦−2𝑧=0, so 𝑦=2𝑧 and 𝑥=𝑦+𝑧=3𝑧, and 𝐯=(3,2,1) is an eigenvector. In the basis 𝐱,𝐲,𝐯,

[𝑇𝐴]=(010−100001),

a 2×2 block without real eigenvalues followed by a 1×1 block.

(c) The plane 𝑈=span(𝐞1,𝐞3) has the orthonormal basis 𝐞1,−𝐞3, and 𝑈⟂=span(𝐞2). Let 𝑄 have columns 𝐞1,−𝐞3,𝐞2:

𝑄=(1000010−10).

Its columns are orthonormal, so 𝑄 is orthogonal. The columns of 𝑄𝖳𝐴𝑄 are the coordinates of 𝐴𝐞1=−(−𝐞3), 𝐴(−𝐞3)=𝐞1, and 𝐴𝐞2=(2,1,−1)=2𝐞1+(−𝐞3)+𝐞2 in this basis:

𝑄𝖳𝐴𝑄=(012−101001).

The compression of 𝑇𝐴 to 𝑈⟂ is multiplication by 1 on span(𝐞2), which is the bottom right entry.

More practice

The practice bank has 22 exercises on this section’s topics.