Books
Theme
Text size
Width

[ and ] go to the previous and next page.

Let 𝐴=(4−211), and suppose we want 𝐴100. Multiplying out 100 copies of 𝐴 is hopeless by hand, and the entries of the partial products show no obvious pattern. But 𝐴 has two eigenvectors,

𝐴(11)=(22)=2(11),𝐴(21)=(63)=3(21),

and they form a basis of ℝ2. In this basis 𝐴 only scales: it doubles the first coordinate and triples the second, so 𝐴100 multiplies them by 2100 and 3100. All that remains is to convert coordinates back. This section finds which operators admit such a basis, by a test in terms of two kinds of multiplicity, and then uses such bases for matrix powers, linear recurrences, and the long-run behavior of linear dynamical systems.

Throughout, 𝑉 is a finite-dimensional vector space over a field 𝐹 with dim𝑉=𝑛≥1, and 𝑇∈L(𝑉).

7.4.1Diagonalizable operators

Definition 7.4.1 (Diagonalizable). An operator 𝑇∈L(𝑉) is diagonalizable if there is a basis ℬ︀ of 𝑉 for which [𝑇]ℬ︀ℬ︀ is a diagonal matrix. A matrix 𝐴∈𝐹𝑛×𝑛 is diagonalizable if it is similar to a diagonal matrix, that is, if

𝐴=𝑆𝐷𝑆−1

for some invertible 𝑆∈𝐹𝑛×𝑛 and some diagonal 𝐷∈𝐹𝑛×𝑛. Such a factorization is a diagonalization of 𝐴.

The lecture defines a diagonalizable operator as one with a basis of eigenvectors. The next theorem shows that this is the same notion, and that the matrix and operator versions agree. A basis of 𝑉 consisting of eigenvectors of 𝑇 is called an eigenbasis for 𝑇.

Theorem 7.4.2 (Diagonalizable means an eigenbasis). Let ℬ︀=𝐛1,…,𝐛𝑛 be a basis of 𝑉 and 𝜆1,…,𝜆𝑛∈𝐹. Then

[𝑇]ℬ︀ℬ︀=diag(𝜆1,…,𝜆𝑛)⟺𝑇𝐛𝑗=𝜆𝑗𝐛𝑗for𝑗=1,…,𝑛.

In particular, 𝑇 is diagonalizable if and only if 𝑉 has an eigenbasis for 𝑇, and then the diagonal entries of [𝑇]ℬ︀ℬ︀ are the eigenvalues of the vectors of the eigenbasis, in order.

For 𝐴∈𝐹𝑛×𝑛 the following are equivalent:

  1. (I)𝐴 is diagonalizable;
  2. (II)the operator 𝑇𝐴 on 𝐹𝑛 is diagonalizable;
  3. (III)𝐹𝑛 has a basis of eigenvectors of 𝐴.

More precisely, if 𝑆 is invertible with columns 𝐬1,…,𝐬𝑛 and 𝐷=diag(𝜆1,…,𝜆𝑛), then 𝐴=𝑆𝐷𝑆−1 if and only if 𝐴𝐬𝑗=𝜆𝑗𝐬𝑗 for every 𝑗.

Proof.routine

Uses. Definition 2.3.1 (Matrix of a linear map), Definition 1.6.14 (Coordinates), Definition 7.2.6 (Eigenvalue, eigenvector), Example 2.4.10 (Diagonal matrices), Corollary 3.5.25 (The invertible matrix theorem), Proposition 2.3.16 (The matrix of 𝑇𝐴)

Techniques. Translate between maps and matrices (Technique 2.3.27)

Hint 1. What is the 𝑗-th column of [𝑇]ℬ︀ℬ︀? When is it 𝜆𝑗𝐞𝑗?
Hint 2. For matrices, rewrite 𝐴=𝑆𝐷𝑆−1 as 𝐴𝑆=𝑆𝐷 and compare the two sides column by column. Right multiplication by a diagonal matrix scales columns.
Hint 3. Columns of an invertible matrix form a basis of 𝐹𝑛, and conversely (Corollary 3.5.25).
Sketch. The 𝑗-th column of [𝑇]ℬ︀ℬ︀ is [𝑇𝐛𝑗]ℬ︀, and it equals 𝜆𝑗𝐞𝑗 exactly when 𝑇𝐛𝑗=𝜆𝑗𝐛𝑗; basis vectors are nonzero, so these are eigenvectors. For matrices, 𝐴=𝑆𝐷𝑆−1 iff 𝐴𝑆=𝑆𝐷, whose 𝑗-th columns are 𝐴𝐬𝑗 and 𝜆𝑗𝐬𝑗. An invertible 𝑆 is the same as a basis of columns, which gives (I) ⇔ (III), and (II) ⇔ (III) is the first part for 𝑇𝐴.

By Definition 2.3.1, the 𝑗-th column of [𝑇]ℬ︀ℬ︀ is the coordinate vector [𝑇𝐛𝑗]ℬ︀. Since coordinates are unique (Definition 1.6.14), [𝑇𝐛𝑗]ℬ︀=𝜆𝑗𝐞𝑗 if and only if 𝑇𝐛𝑗=𝜆𝑗𝐛𝑗. So [𝑇]ℬ︀ℬ︀=diag(𝜆1,…,𝜆𝑛) exactly when 𝑇𝐛𝑗=𝜆𝑗𝐛𝑗 for all 𝑗. The vectors of a basis are nonzero, so in that case each 𝐛𝑗 is an eigenvector with eigenvalue 𝜆𝑗 (Definition 7.2.6), and ℬ︀ is an eigenbasis. Conversely, if ℬ︀ is an eigenbasis with 𝑇𝐛𝑗=𝜆𝑗𝐛𝑗, then [𝑇]ℬ︀ℬ︀=diag(𝜆1,…,𝜆𝑛). This proves the first two claims.

For matrices, let 𝑆 be invertible with columns 𝐬𝑗 and 𝐷=diag(𝜆1,…,𝜆𝑛). Multiplying on the right by 𝑆 or by 𝑆−1 shows that 𝐴=𝑆𝐷𝑆−1 if and only if 𝐴𝑆=𝑆𝐷. The 𝑗-th column of 𝐴𝑆 is 𝐴𝐬𝑗, and the 𝑗-th column of 𝑆𝐷 is 𝜆𝑗𝐬𝑗, since right multiplication by a diagonal matrix scales the columns (Example 2.4.10). So 𝐴=𝑆𝐷𝑆−1 if and only if 𝐴𝐬𝑗=𝜆𝑗𝐬𝑗 for every 𝑗.

(I) ⇔ (III). If 𝐴=𝑆𝐷𝑆−1, the columns of 𝑆 form a basis of 𝐹𝑛 by Corollary 3.5.25, and by the previous paragraph they are eigenvectors of 𝐴. Conversely, if 𝐬1,…,𝐬𝑛 is a basis of 𝐹𝑛 with 𝐴𝐬𝑗=𝜆𝑗𝐬𝑗, the matrix 𝑆 with these columns is invertible (Corollary 3.5.25), and 𝐴=𝑆𝐷𝑆−1 with 𝐷=diag(𝜆1,…,𝜆𝑛).

(II) ⇔ (III). The eigenvectors of 𝐴 are by definition those of 𝑇𝐴, so this is the first part of the theorem applied to 𝑇=𝑇𝐴. □

Remark 7.4.3 (Two readings of 𝐴=𝑆𝐷𝑆−1). The proof read 𝐴=𝑆𝐷𝑆−1 column by column: 𝐴 scales the columns of 𝑆. There is also the reading of Section 4.5. If ℬ︀̊ is the list of columns of 𝑆, then 𝑆=[id]ℰ︀ℬ︀̊, and Theorem 4.5.17 gives

𝐷=𝑆−1𝐴𝑆=[id]ℬ︀̊ℰ︀[𝑇𝐴]ℰ︀ℰ︀[id]ℰ︀ℬ︀̊=[𝑇𝐴]ℬ︀̊ℬ︀̊.

Read from right to left, 𝐴=𝑆𝐷𝑆−1 says: convert standard coordinates into eigenbasis coordinates (𝑆−1), scale each coordinate by its eigenvalue (𝐷), and convert back (𝑆). Diagonalization is the strongest instance of Technique 4.5.33: the operator is as simple as an operator can be, and the basis that achieves this is made of the operator’s own eigenvectors.

Example 7.4.4 (Diagonalizing a 2×2 matrix). Let 𝐴=(4−211)∈ℝ2×2, the matrix of the introduction. By Theorem 7.2.9, 𝜆 is an eigenvalue exactly when

det(𝐴−𝜆𝐼2)=(4−𝜆)(1−𝜆)+2=𝜆2−5𝜆+6=(𝜆−2)(𝜆−3)

vanishes, so the eigenvalues are 2 and 3. Row reducing,

𝐴−2𝐼2=(2−21−1)⇝(1−100),𝐴−3𝐼2=(1−21−2)⇝(1−200),

so 𝐸2=span((1,1)) and 𝐸3=span((2,1)). The vectors (1,1) and (2,1) are not parallel, so they form an eigenbasis, and by Theorem 7.4.2

𝐴=𝑆𝐷𝑆−1,𝑆=(1211),𝐷=(2003),𝑆−1=(−121−1).

To check the work it is enough to verify 𝐴𝑆=𝑆𝐷: both sides equal (2623). The order of the columns is our choice. Swapping them swaps the diagonal entries of 𝐷, and rescaling a column of 𝑆 by 𝑐≠0 changes nothing in 𝐷. So a diagonalization is far from unique, but 𝐷 is unique up to the order of its entries (Theorem 7.4.16 will show why).

[𝐞1]ℬ︀=(−11),[𝐴𝐞1]ℬ︀=[(41)]ℬ︀=(−23)=(2003)(−11).
The standard grid and the grid of the eigenbasis ℬ︀=(1,1),(2,1) of 𝐴=(4−211). In ℬ︀-coordinates 𝐴 only scales, [𝐴𝐯]ℬ︀=𝐷[𝐯]ℬ︀ with 𝐷=diag(2,3). The vector 𝐯=𝐞1 starts with [𝐯]ℬ︀=(−1,1), so 𝐴𝐞1=−2(1,1)+3(2,1)=(4,1), which is the first column of 𝐴. In the interactive version, type 𝐯=(4,1) and confirm that its ℬ︀-coordinates are (−2,3). Moving 𝐯 along a line of the ℬ︀-grid changes only one coordinate, and 𝐴 moves the point along the corresponding line of the stretched grid.

In Example 7.4.4 the two eigenvalues were different, and that alone guaranteed the eigenbasis.

Corollary 7.4.5 (Distinct eigenvalues give a diagonalization). If 𝑇 has 𝑛=dim𝑉 distinct eigenvalues 𝜆1,…,𝜆𝑛, then 𝑇 is diagonalizable: any choice of eigenvectors 𝐯1,…,𝐯𝑛 for 𝜆1,…,𝜆𝑛 is an eigenbasis, and in it 𝑇 has matrix diag(𝜆1,…,𝜆𝑛). In particular, 𝐴∈𝐹𝑛×𝑛 with 𝑛 distinct eigenvalues is diagonalizable.

Proof.routine

Uses. Theorem 7.2.21 (Eigenvectors for distinct eigenvalues are independent), Theorem 4.1.10 (Independent lists are short), Theorem 7.4.2 (Diagonalizable means an eigenbasis)

Hint 1. Pick one eigenvector for each eigenvalue. What do you know about this list?
Hint 2. It is independent and has length dim𝑉.
Sketch. Eigenvectors for distinct eigenvalues are independent (Theorem 7.2.21), and 𝑛 independent vectors in an 𝑛-dimensional space form a basis (Theorem 4.1.10(b)). Apply Theorem 7.4.2.
The list 𝐯1,…,𝐯𝑛 is independent by Theorem 7.2.21, and it has length 𝑛=dim𝑉, so it is a basis by Theorem 4.1.10(b). It consists of eigenvectors, so Theorem 7.4.2 gives [𝑇]ℬ︀ℬ︀=diag(𝜆1,…,𝜆𝑛) for ℬ︀=𝐯1,…,𝐯𝑛. □

The converse fails: for 𝑛≥2, 𝐼𝑛 is diagonal, with the single eigenvalue 1. Repeated eigenvalues are allowed as long as each one comes with enough eigenvectors.

Example 7.4.6 (A repeated eigenvalue with enough eigenvectors). The matrix 𝐴=2𝐼3+𝐽=(311131113) of Example 7.2.14 has only the eigenvalues 2 and 5, with

𝐸2=span((1,−1,0),(0,1,−1)),𝐸5=span((1,1,1)).

These three vectors form a basis of ℝ3: the sum 𝐸2+𝐸5 is direct (Corollary 7.2.23), so concatenating bases gives an independent list (Theorem 4.4.13), of length 3. Hence

𝐴=𝑆(200020005)𝑆−1,𝑆=(101−1110−11).

Obstructions to diagonalizability

There are two ways an operator can fail to be diagonalizable: it may have too few eigenvalues in 𝐹, or an eigenvalue may have too few eigenvectors. The simplest instance of the second failure is an operator with one eigenvalue.

Proposition 7.4.7 (One eigenvalue). Suppose 𝑇 is diagonalizable and 𝜆 is its only eigenvalue. Then 𝑇=𝜆id𝑉. In particular, if 𝐴∈𝐹𝑛×𝑛 is diagonalizable with only eigenvalue 𝜆, then 𝐴=𝜆𝐼𝑛.

Proof.routine

Uses. Theorem 7.4.2 (Diagonalizable means an eigenbasis), Definition 7.2.6 (Eigenvalue, eigenvector)

Techniques. Check it on a basis (Technique 2.2.7)

Hint 1. Take an eigenbasis. What is the eigenvalue of each of its vectors?
Hint 2. Two linear maps that agree on a basis are equal.
Sketch. Every vector of an eigenbasis is an eigenvector, so its eigenvalue is 𝜆. Then 𝑇 and 𝜆id agree on a basis, so they are equal.
Let 𝐛1,…,𝐛𝑛 be an eigenbasis for 𝑇. The eigenvalue of each 𝐛𝑗 is an eigenvalue of 𝑇, hence equals 𝜆, so 𝑇𝐛𝑗=𝜆𝐛𝑗=(𝜆id𝑉)𝐛𝑗 for every 𝑗. Two linear maps that agree on a basis are equal (Technique 2.2.7), so 𝑇=𝜆id𝑉. For 𝐴, apply this to 𝑇𝐴: then 𝑇𝐴=𝜆id, and its matrix in the standard basis is 𝐴=𝜆𝐼𝑛. □

Example 7.4.8 (Shears and nilpotent matrices). Let 𝐽=(1101), a shear, over any field. Then 𝐽−𝜇𝐼2=(1−𝜇101−𝜇) has determinant (1−𝜇)2, so 1 is the only eigenvalue (Theorem 7.2.9). Since 𝐽≠𝐼2, Proposition 7.4.7 shows that 𝐽 is not diagonalizable. Directly: 𝐸1=ker(0100)=span(𝐞1) is a line, so there is no basis of eigenvectors.

More generally, let 𝑁∈𝐹𝑛×𝑛 be nilpotent and nonzero, say 𝑁𝑘=0. If 𝑁𝐯=𝜆𝐯 with 𝐯≠𝟎, then 𝟎=𝑁𝑘𝐯=𝜆𝑘𝐯, so 𝜆=0; and 0 is an eigenvalue, because 𝑁 is not injective (otherwise 𝑁𝑘 would be injective too, while 𝑁𝑘=0 on 𝐹𝑛≠{𝟎}). So 0 is the only eigenvalue, and 𝑁≠0=0⋅𝐼𝑛 is not diagonalizable. The same argument works for an upper triangular matrix with all diagonal entries equal to 𝜆: its only eigenvalue is 𝜆, because 𝐴−𝜇𝐼𝑛 is triangular with diagonal entries 𝜆−𝜇, hence invertible for 𝜇≠𝜆 (Proposition 5.3.15, Theorem 5.4.7). So it is diagonalizable only if it equals 𝜆𝐼𝑛.

These are the minimal counterexamples of Technique 1.4.18 for diagonalizability. In particular, being invertible (𝐽) and being diagonalizable are unrelated: 𝐽 is invertible but not diagonalizable, and diag(1,0) is diagonal but not invertible.

Example 7.4.9 (Rotations: the field matters). Let 𝐴=(0−110), the rotation by a right angle. By Example 7.2.18 its eigenvalues are the square roots of −1 in 𝐹.

  • Over ℝ there are none, so there are no eigenvectors at all, and 𝐴 is not diagonalizable over ℝ.
  • Over ℂ the eigenvalues 𝑖 and −𝑖 are distinct, with eigenvectors (1,−𝑖) and (1,𝑖), so by Corollary 7.4.5

    𝐴=(11−𝑖𝑖)(𝑖00−𝑖)(11−𝑖𝑖)−1.
  • Over 𝔽5, where 22=32=−1, it is diagonalizable with eigenvalues 2 and 3; over 𝔽3 there are no eigenvalues.
  • Over 𝔽2, 𝐴 is the swap (𝑥,𝑦)↦(𝑦,𝑥), whose only eigenvalue is 1 (Example 7.2.17). Since 𝐴≠𝐼2, it is not diagonalizable over 𝔽2, by Proposition 7.4.7. Over every field with 1+1≠0 the swap is diagonalizable, with eigenbasis (1,1),(1,−1).

So “diagonalizable” is a property of a matrix together with a field. The real matrix 𝐴 fails for lack of eigenvalues, while over 𝔽2 the matrix fails for lack of eigenvectors.

Exercise 7.4.10.routine Without computing any eigenvectors, decide which of the following real matrices are diagonalizable over ℝ, and which over ℂ.

𝐴1=(1502),𝐴2=(2502),𝐴3=(1−21−1),𝐴4=(700070007).

Solution.

Uses. Corollary 7.4.5 (Distinct eigenvalues give a diagonalization), Proposition 7.4.7 (One eigenvalue), Theorem 7.2.9 (Eigenvalues and invertibility), Proposition 5.3.15 (Triangular matrices)

Hint 1. Find the eigenvalues from det(𝐴−𝜆𝐼)=0. Triangular matrices show them on the diagonal.
Hint 2. Distinct eigenvalues settle the question at once. A matrix with a single eigenvalue is diagonalizable only if it is a scalar matrix.
Sketch. 𝐴1 has distinct eigenvalues 1,2: diagonalizable over both fields. 𝐴2 has only the eigenvalue 2 and is not 2𝐼2: diagonalizable over neither. 𝐴3 has det(𝐴3−𝜆𝐼)=𝜆2+1: no real eigenvalues, so not over ℝ; distinct eigenvalues ±𝑖 over ℂ. 𝐴4=7𝐼3 is already diagonal.

The matrices 𝐴1 and 𝐴2 are upper triangular, so 𝐴−𝜆𝐼 is triangular with diagonal entries 𝑎𝑖𝑖−𝜆, and by Proposition 5.3.15 and Theorem 7.2.9 the eigenvalues are the diagonal entries. For 𝐴1 they are 1 and 2, which are distinct, so 𝐴1 is diagonalizable over ℝ, and over ℂ, by Corollary 7.4.5. For 𝐴2 the only eigenvalue is 2, over either field, and 𝐴2≠2𝐼2, so 𝐴2 is diagonalizable over neither field (Proposition 7.4.7).

For 𝐴3, det(𝐴3−𝜆𝐼2)=(1−𝜆)(−1−𝜆)+2=𝜆2+1. Over ℝ it never vanishes, so 𝐴3 has no eigenvalues and no eigenvectors, and is not diagonalizable. Over ℂ it has the two distinct eigenvalues 𝑖 and −𝑖, so 𝐴3 is diagonalizable over ℂ.

Finally 𝐴4=7𝐼3 is diagonal, so it is diagonalizable over both fields (take 𝑆=𝐼3), even though it has only one eigenvalue.

7.4.2Algebraic and geometric multiplicity

The obstructions above suggest counting, for each eigenvalue, how many times it “should” occur on the diagonal and how many independent eigenvectors it actually has. The first count comes from the characteristic polynomial, the second from the eigenspace.

Recall from Section 7.3 that the characteristic polynomial of 𝐴∈𝐹𝑛×𝑛 is the monic polynomial 𝑝𝐴=det(𝑥𝐼𝑛−𝐴)∈𝐹[𝑥] of degree 𝑛 (Definition 7.3.4), and that its roots in 𝐹 are exactly the eigenvalues of 𝐴 (Theorem 7.3.6). Similar matrices have the same characteristic polynomial, so the characteristic polynomial of 𝑇, 𝑝𝑇=𝑝𝐴 for 𝐴=[𝑇]ℬ︀ℬ︀, does not depend on the basis ℬ︀. Determinants of matrices with entries in 𝐹[𝑥] obey the same rules as over 𝐹 (triangular and block triangular matrices, products), as explained in Section 7.3; we use these rules freely.

Definition 7.4.11 (Algebraic multiplicity). Let 𝜆∈𝐹. The algebraic multiplicity of 𝜆 for 𝑇, written am𝑇(𝜆), is the multiplicity of 𝜆 as a root of 𝑝𝑇 (Definition 7.1.24): the largest 𝑚≥0 with (𝑥−𝜆)𝑚 dividing 𝑝𝑇. For a matrix 𝐴, am𝐴(𝜆) is the multiplicity of 𝜆 as a root of 𝑝𝐴.

Definition 7.4.12 (Geometric multiplicity). Let 𝜆∈𝐹. The geometric multiplicity of 𝜆 for 𝑇 is

gm𝑇(𝜆)=dim𝐸𝜆(𝑇)=dimker(𝑇−𝜆id𝑉),

and for a matrix 𝐴, gm𝐴(𝜆)=dimker(𝐴−𝜆𝐼𝑛).

We drop the subscript when the operator is clear. Both multiplicities are defined for every 𝜆∈𝐹, and by Theorem 7.3.6 and Theorem 7.2.9 each is positive exactly when 𝜆 is an eigenvalue. The geometric multiplicity is computed by row reduction, through a rank.

Lemma 7.4.13 (Geometric multiplicity through a rank). For every basis ℬ︀ of 𝑉 and every 𝜆∈𝐹,

gm𝑇(𝜆)=𝑛−rank([𝑇]ℬ︀ℬ︀−𝜆𝐼𝑛).

In particular, gm𝑇(𝜆)=gm𝐴(𝜆) for every matrix 𝐴 of 𝑇, and similar matrices have the same geometric multiplicities.

Proof.routine

Uses. Theorem 4.2.11 (Rank–nullity), Proposition 4.2.15 (Rank through any matrix), Definition 2.3.1 (Matrix of a linear map), Theorem 4.5.25 (Similar matrices are one operator in two bases)

Hint 1. gm𝑇(𝜆) is the nullity of 𝑇−𝜆id. What does rank–nullity say?
Hint 2. What is the matrix of 𝑇−𝜆id in the basis ℬ︀?
Sketch. By rank–nullity, dimker(𝑇−𝜆id)=𝑛−rank(𝑇−𝜆id). The matrix of 𝑇−𝜆id in ℬ︀ is [𝑇]ℬ︀ℬ︀−𝜆𝐼𝑛, and a map has the rank of any of its matrices. For a matrix 𝐴 this is the same formula with 𝑇=𝑇𝐴; similar matrices are matrices of one operator.
By Theorem 4.2.11, gm𝑇(𝜆)=dimker(𝑇−𝜆id)=𝑛−rank(𝑇−𝜆id). The matrix of a sum and of a scalar multiple of maps is the sum and the scalar multiple of the matrices, and [id]ℬ︀ℬ︀=𝐼𝑛, so [𝑇−𝜆id]ℬ︀ℬ︀=[𝑇]ℬ︀ℬ︀−𝜆𝐼𝑛. A linear map has the same rank as each of its matrices (Proposition 4.2.15), which gives the formula. For 𝐴∈𝐹𝑛×𝑛, gm𝐴(𝜆)=𝑛−rank(𝐴−𝜆𝐼𝑛) by the same computation with 𝑇=𝑇𝐴 and the standard basis. Finally, two similar matrices are matrices of one operator in two bases (Theorem 4.5.25), so they have the same geometric multiplicities, namely those of that operator. □

Example 7.4.14 (The multiplicities of a diagonal matrix). Let 𝐷=diag(𝑑1,…,𝑑𝑛) and 𝜆∈𝐹, and let 𝑚 be the number of indices 𝑖 with 𝑑𝑖=𝜆. Then 𝑥𝐼𝑛−𝐷=diag(𝑥−𝑑1,…,𝑥−𝑑𝑛) is diagonal, so

𝑝𝐷=(𝑥−𝑑1)(𝑥−𝑑2)⋯(𝑥−𝑑𝑛),

and by Proposition 7.1.34, am𝐷(𝜆)=𝑚. On the other hand 𝐷−𝜆𝐼𝑛=diag(𝑑1−𝜆,…,𝑑𝑛−𝜆) has exactly 𝑛−𝑚 nonzero columns, which are multiples of distinct standard basis vectors, so its rank is 𝑛−𝑚, and gm𝐷(𝜆)=𝑚 by Lemma 7.4.13. So for a diagonal matrix the two multiplicities agree: both count how often 𝜆 occurs on the diagonal. Consequently, if 𝑇 is diagonalizable, the number of times 𝜆 occurs on the diagonal of a diagonal matrix of 𝑇 does not depend on the eigenbasis: it is gm𝑇(𝜆). This is why 𝐷 in a diagonalization 𝐴=𝑆𝐷𝑆−1 is unique up to order.

For the shear 𝐽=(1101), 𝑝𝐽=(𝑥−1)2, so am𝐽(1)=2, while gm𝐽(1)=2−rank(0100)=1.

Intuition (Why the geometric multiplicity cannot exceed the algebraic one). Suppose 𝜆 has 𝑔 independent eigenvectors. Put them first in a basis. Then the first 𝑔 columns of the matrix of 𝑇 are 𝜆𝐞1,…,𝜆𝐞𝑔, so the matrix is block upper triangular with top left block 𝜆𝐼𝑔, and its characteristic polynomial picks up the factor (𝑥−𝜆)𝑔 from that block. The bottom right block may contribute further factors 𝑥−𝜆, which are not visible as eigenvectors. For the shear, the bottom right block is the 1×1 matrix (1), and it contributes the second factor 𝑥−1 without contributing a second eigenvector.

Theorem 7.4.15 (Geometric multiplicity is at most algebraic multiplicity). For every 𝜆∈𝐹,

gm𝑇(𝜆)≤am𝑇(𝜆).

If 𝜆 is an eigenvalue of 𝑇, then 1≤gm𝑇(𝜆)≤am𝑇(𝜆).

Proof.standard

Uses. Theorem 4.1.16 (Extending to a basis), Proposition 4.5.35 (Preserved subspaces give block triangular matrices), Definition 7.2.13 (Eigenspace), Theorem 5.4.26 (Block-triangular matrices), Definition 7.3.4 (Characteristic polynomial), Definition 7.1.24 (Multiplicity of a root), Proposition 4.1.21 (Dimensions of subspaces)

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

Hint 1. Let 𝑔=gm𝑇(𝜆). Choose a basis of 𝐸𝜆 and extend it to a basis of 𝑉. What does the matrix of 𝑇 look like?
Hint 2. 𝐸𝜆 is 𝑇-invariant and 𝑇 acts on it as 𝜆id, so by Proposition 4.5.35 the matrix is (𝜆𝐼𝑔𝐵0𝐶).
Hint 3. Compute det(𝑥𝐼𝑛−[𝑇]ℬ︀ℬ︀) with the block triangular rule. Which power of 𝑥−𝜆 divides it?
Sketch. Extend a basis of 𝐸𝜆, of length 𝑔, to a basis ℬ︀ of 𝑉. Since 𝑇 acts on 𝐸𝜆 as 𝜆id, [𝑇]ℬ︀ℬ︀=(𝜆𝐼𝑔𝐵0𝐶), so 𝑝𝑇=det((𝑥−𝜆)𝐼𝑔)det(𝑥𝐼𝑛−𝑔−𝐶)=(𝑥−𝜆)𝑔𝑝𝐶. Hence (𝑥−𝜆)𝑔 divides 𝑝𝑇, and 𝑔≤am𝑇(𝜆). If 𝜆 is an eigenvalue, 𝐸𝜆≠{𝟎}, so 𝑔≥1.

Let 𝑔=gm𝑇(𝜆)=dim𝐸𝜆. If 𝑔=0 there is nothing to prove. If 𝑔=𝑛, then 𝐸𝜆=𝑉 (Proposition 4.1.21), so 𝑇=𝜆id𝑉, [𝑇]ℬ︀ℬ︀=𝜆𝐼𝑛 in every basis, and 𝑝𝑇=(𝑥−𝜆)𝑛; so am𝑇(𝜆)=𝑛=𝑔.

Let 1≤𝑔<𝑛. Choose a basis 𝐮1,…,𝐮𝑔 of 𝐸𝜆 and extend it to a basis ℬ︀=𝐮1,…,𝐮𝑔,𝐰1,…,𝐰𝑛−𝑔 of 𝑉 (Theorem 4.1.16). For 𝐮∈𝐸𝜆 we have 𝑇𝐮=𝜆𝐮∈𝐸𝜆, so 𝐸𝜆 is mapped into itself, and the restriction of 𝑇 to 𝐸𝜆 is 𝜆id𝐸𝜆, whose matrix in 𝐮1,…,𝐮𝑔 is 𝜆𝐼𝑔. By Proposition 4.5.35,

[𝑇]ℬ︀ℬ︀=(𝜆𝐼𝑔𝐵0𝐶)

for some 𝐵∈𝐹𝑔×(𝑛−𝑔) and 𝐶∈𝐹(𝑛−𝑔)×(𝑛−𝑔). Then

𝑥𝐼𝑛−[𝑇]ℬ︀ℬ︀=((𝑥−𝜆)𝐼𝑔−𝐵0𝑥𝐼𝑛−𝑔−𝐶)

is block upper triangular with entries in 𝐹[𝑥], and by Theorem 5.4.26

𝑝𝑇=det((𝑥−𝜆)𝐼𝑔)⋅det(𝑥𝐼𝑛−𝑔−𝐶)=(𝑥−𝜆)𝑔𝑝𝐶.

So (𝑥−𝜆)𝑔 divides 𝑝𝑇, and 𝑔≤am𝑇(𝜆) by Definition 7.1.24.

If 𝜆 is an eigenvalue, then 𝐸𝜆 contains a nonzero vector, so 𝑔≥1. □

The lecture calls Theorem 7.4.15 the “AM–GM inequality”, as a joke: the inequality between the arithmetic and geometric means has the same letters in the same order. The proof is Technique 4.5.36 applied to the most obvious invariant subspace of 𝑇, an eigenspace.

We can now say exactly when 𝑇 is diagonalizable. Let 𝜆1,…,𝜆𝑚 be the distinct eigenvalues of 𝑇. Every eigenvector lies in some 𝐸𝜆𝑖, the sum of these eigenspaces is direct (Corollary 7.2.23), and 𝑝𝑇 has degree 𝑛. So

gm(𝜆1)+⋯+gm(𝜆𝑚)≤am(𝜆1)+⋯+am(𝜆𝑚)≤deg𝑝𝑇=𝑛,
(7.4.1)

where the first inequality is Theorem 7.4.15 term by term, and the second is Proposition 7.1.27, since the 𝜆𝑖 are distinct roots of 𝑝𝑇. Diagonalizability turns out to mean that the left end equals the right end, which forces both inequalities to be equalities.

Theorem 7.4.16 (Criteria for diagonalizability). Let 𝜆1,…,𝜆𝑚 be the distinct eigenvalues of 𝑇 (𝑚≥0). The following are equivalent.

  1. (I)𝑇 is diagonalizable.
  2. (II)𝑉=𝐸𝜆1⊕𝐸𝜆2⊕⋯⊕𝐸𝜆𝑚.
  3. (III)gm(𝜆1)+gm(𝜆2)+⋯+gm(𝜆𝑚)=𝑛.
  4. (IV)𝑝𝑇 splits over 𝐹, and gm(𝜆𝑖)=am(𝜆𝑖) for 𝑖=1,…,𝑚.

In that case, for every eigenbasis ℬ︀ the matrix [𝑇]ℬ︀ℬ︀ has 𝜆𝑖 on its diagonal exactly gm(𝜆𝑖)=am(𝜆𝑖) times.

Proof, part 1 of 2: (I) ⇒ (II) ⇒ (III) ⇒ (I).standard

Uses. Theorem 7.4.2 (Diagonalizable means an eigenbasis), Corollary 7.2.23 (Eigenspaces are independent), Theorem 4.4.13 (Bases and dimension of a sum), Proposition 4.1.21 (Dimensions of subspaces), Proposition 4.4.2 (Sums are the smallest containing subspaces)

Techniques. Proving TFAE (Technique 2.7.1), Dimension counting (Technique 4.2.24), Glue bases of the summands (Technique 4.4.14)

Hint 1. For (I) ⇒ (II): where does each vector of an eigenbasis live? The sum is direct by Corollary 7.2.23.
Hint 2. For (II) ⇒ (III) and (III) ⇒ (I), use Theorem 4.4.13: the dimension of a direct sum, and bases of a direct sum by concatenation.
Hint 3. In (III) ⇒ (I), the concatenated list is independent of length 𝑛. Why is it a basis of 𝑉?
Sketch. (I) ⇒ (II): each vector of an eigenbasis lies in some 𝐸𝜆𝑖, so the sum of the eigenspaces contains a basis and equals 𝑉; it is direct by Corollary 7.2.23. (II) ⇒ (III): the dimension of a direct sum is the sum of the dimensions. (III) ⇒ (I): the sum of the eigenspaces is direct, so concatenating their bases gives a basis of the sum, of length 𝑛; the sum is then all of 𝑉, and the concatenation is an eigenbasis.

Write 𝑈=𝐸𝜆1+⋯+𝐸𝜆𝑚. By Corollary 7.2.23 the sum is direct, whatever 𝑇 is.

(I) ⇒ (II). Let ℬ︀ be an eigenbasis. Each vector of ℬ︀ is an eigenvector, whose eigenvalue is one of the 𝜆𝑖, so it lies in 𝐸𝜆𝑖⊆𝑈. So 𝑈, which is a subspace (Proposition 4.4.2), contains a basis of 𝑉, hence contains its span 𝑉. Thus 𝑉=𝑈, and the sum is direct.

(II) ⇒ (III). By Theorem 4.4.13(c), 𝑛=dim𝑉=dim𝐸𝜆1+⋯+dim𝐸𝜆𝑚.

(III) ⇒ (I). Choose a basis of each 𝐸𝜆𝑖 and concatenate them. Since 𝑈 is a direct sum, the concatenation is a basis of 𝑈 by Theorem 4.4.13(b), so dim𝑈=∑𝑖gm(𝜆𝑖)=𝑛, and 𝑈=𝑉 by Proposition 4.1.21(b). So the concatenation is a basis of 𝑉 consisting of eigenvectors, and 𝑇 is diagonalizable by Theorem 7.4.2.

So far nothing used the characteristic polynomial: (I), (II), and (III) are statements about eigenvectors only. Condition (IV) brings in the polynomial, through the squeeze (7.4.1).

Proof, part 2 of 2: (III) ⇔ (IV).standard

Uses. Theorem 7.4.15 (Geometric multiplicity is at most algebraic multiplicity), Proposition 7.1.27 (Multiplicities add up to at most the degree), Proposition 7.1.34 (Roots of a split polynomial), Theorem 7.3.6 (Eigenvalues are the roots of the characteristic polynomial), Proposition 7.1.4 (Degree of a product), Example 7.4.14 (The multiplicities of a diagonal matrix)

Hint 1. Look at the chain (7.4.1). Condition (III) says that its two ends are equal.
Hint 2. A chain 𝑎≤𝑏≤𝑐 with 𝑎=𝑐 forces 𝑎=𝑏 and 𝑏=𝑐. The first equality is a sum of termwise inequalities.
Hint 3. For ∑𝑖am(𝜆𝑖)=𝑛 ⇔ “𝑝𝑇 splits”: Proposition 7.1.27 gives a monic divisor of 𝑝𝑇 of degree ∑𝑖am(𝜆𝑖); compare degrees.
Sketch. In ∑gm(𝜆𝑖)≤∑am(𝜆𝑖)≤𝑛, (III) says the ends are equal. That holds iff gm(𝜆𝑖)=am(𝜆𝑖) for each 𝑖 (the terms satisfy gm≤am) and ∑am(𝜆𝑖)=𝑛. The latter means that ∏(𝑥−𝜆𝑖)am(𝜆𝑖), a monic divisor of 𝑝𝑇, has degree 𝑛, so equals 𝑝𝑇, i.e. 𝑝𝑇 splits; conversely, if 𝑝𝑇 splits, its roots are the 𝜆𝑖 and their multiplicities add up to 𝑛.

Write 𝑔𝑖=gm(𝜆𝑖) and 𝑎𝑖=am(𝜆𝑖). By Theorem 7.4.15, 𝑔𝑖≤𝑎𝑖 for each 𝑖, so ∑𝑖𝑔𝑖≤∑𝑖𝑎𝑖, with equality if and only if 𝑔𝑖=𝑎𝑖 for every 𝑖. By Theorem 7.3.6 the 𝜆𝑖 are distinct roots of 𝑝𝑇, so by Proposition 7.1.27,

𝑞=(𝑥−𝜆1)𝑎1⋯(𝑥−𝜆𝑚)𝑎𝑚divides𝑝𝑇,and∑𝑖𝑎𝑖≤𝑛.

We claim that ∑𝑖𝑎𝑖=𝑛 if and only if 𝑝𝑇 splits. If ∑𝑖𝑎𝑖=𝑛, write 𝑝𝑇=𝑞ℎ. Comparing degrees (Proposition 7.1.4), degℎ=0, and comparing leading coefficients, which are 1 for 𝑝𝑇 and 𝑞, gives ℎ=1. So 𝑝𝑇=𝑞 splits. Conversely, if 𝑝𝑇 splits, then by Proposition 7.1.34 the multiplicities of its distinct roots add up to deg𝑝𝑇=𝑛, and its roots are exactly the eigenvalues 𝜆1,…,𝜆𝑚 (Theorem 7.3.6). So ∑𝑖𝑎𝑖=𝑛.

Now (III) says ∑𝑖𝑔𝑖=𝑛. Since ∑𝑖𝑔𝑖≤∑𝑖𝑎𝑖≤𝑛, this holds if and only if ∑𝑖𝑔𝑖=∑𝑖𝑎𝑖 and ∑𝑖𝑎𝑖=𝑛, that is, if and only if 𝑔𝑖=𝑎𝑖 for every 𝑖 and 𝑝𝑇 splits. This is (IV).

For the last claim, if ℬ︀ is an eigenbasis, [𝑇]ℬ︀ℬ︀ is a diagonal matrix of 𝑇, and by Example 7.4.14 and Lemma 7.4.13, 𝜆𝑖 occurs on its diagonal gm(𝜆𝑖) times. □

Remark 7.4.17 (Over ℂ, and over other fields). Over ℂ every polynomial splits (Corollary 7.1.35), so a complex operator is diagonalizable if and only if gm(𝜆)=am(𝜆) for every eigenvalue. Over ℝ, or over ℚ or 𝔽𝑝, splitting is a genuine condition: the real rotation by a right angle has 𝑝=𝑥2+1, both sums in (7.4.1) are 0, and 𝑛=2. Condition (III) handles both failures at once, which makes it the practical test: add up the dimensions of the eigenspaces and compare with 𝑛. Splitting alone is not enough, as the shear shows.

Algorithm 7.4.18 (Diagonalizing a matrix)

Input
𝐴∈𝐹𝑛×𝑛.
Output
An invertible 𝑆 and a diagonal 𝐷 with 𝐴=𝑆𝐷𝑆−1, or the conclusion that 𝐴 is not diagonalizable over 𝐹.

1.

Compute 𝑝𝐴=det(𝑥𝐼𝑛−𝐴) and find its distinct roots 𝜆1,…,𝜆𝑚 in 𝐹.

2.

For each 𝜆𝑖, row reduce 𝐴−𝜆𝑖𝐼𝑛 and read off a basis of 𝐸𝜆𝑖=ker(𝐴−𝜆𝑖𝐼𝑛), of length 𝑔𝑖=𝑛−rank(𝐴−𝜆𝑖𝐼𝑛).

3.

If 𝑔1+⋯+𝑔𝑚<𝑛, stop: 𝐴 is not diagonalizable over 𝐹.

4.

Otherwise let 𝑆 be the matrix whose columns are the basis vectors found in step 2, and 𝐷 the diagonal matrix whose 𝑗-th entry is the eigenvalue of the 𝑗-th column of 𝑆. Check 𝐴𝑆=𝑆𝐷.

The algorithm is correct by Theorem 7.4.16, (I) ⇔ (III), and the proof of (III) ⇒ (I) shows that the columns found in step 2 form a basis. In step 1 one can stop early if 𝑝𝐴 does not split, and in step 2 as soon as some 𝑔𝑖 is less than the multiplicity of 𝜆𝑖 as a root.

Example 7.4.19 (Two matrices that differ in one entry). Let

𝐴=(311020002),𝐴′=(311021002).

Both are upper triangular with diagonal 3,2,2, so 𝑝𝐴=𝑝𝐴′=(𝑥−3)(𝑥−2)2: the eigenvalue 3 is simple, and am(2)=2 for both. For 𝐴,

𝐴−2𝐼3=(111000000)

has rank 1, so gm𝐴(2)=2, with 𝐸2={(𝑥,𝑦,𝑧)|𝑥+𝑦+𝑧=0}=span((1,−1,0),(1,0,−1)). Since 3 is simple, gm𝐴(3)=1, by Theorem 7.4.15, and 𝐸3=span(𝐞1). The dimensions add up to 3, and

𝐴=𝑆(300020002)𝑆−1,𝑆=(1110−1000−1).

For 𝐴′,

𝐴′−2𝐼3=(111001000)

has rank 2, so gm𝐴′(2)=1<2=am𝐴′(2), and 𝐴′ is not diagonalizable. The single entry 𝑎23′=1 destroys one eigenvector for 2. In particular 𝐴 and 𝐴′ are not similar, although they have the same characteristic polynomial (Lemma 7.4.13).

Example 7.4.20 (The Fibonacci matrix over several fields). Let 𝑄=(1110), so that 𝑝𝑄=det(𝑥−1−1−1𝑥)=𝑥2−𝑥−1. Each failure mode of Theorem 7.4.16 occurs over some field.

  • Over 𝔽2, 𝑝𝑄=𝑥2+𝑥+1 has no roots (𝑝𝑄(0)=𝑝𝑄(1)=1). So 𝑄 has no eigenvalues and is not diagonalizable: 𝑝𝑄 does not split.
  • Over 𝔽5, (𝑥+2)2=𝑥2+4𝑥+4=𝑥2−𝑥−1, so 𝑝𝑄 splits with the single root −2=3, of algebraic multiplicity 2. But 𝑄−3𝐼2=(−211−3)=(3112) has rank 1 (its second row is 2 times its first in 𝔽5), so gm(3)=1<2. Not diagonalizable: the multiplicities differ.
  • Over 𝔽11, (𝑥−4)(𝑥+3)=𝑥2−𝑥−12=𝑥2−𝑥−1, so 𝑄 has the distinct eigenvalues 4 and −3=8, and is diagonalizable by Corollary 7.4.5.
  • Over ℝ, the roots 𝜑=(1+5)/2 and 𝜓=(1−5)/2 are distinct, and 𝑄 is diagonalizable; Section 7.4.4 uses this to find the Fibonacci numbers.

The discriminant 5 of 𝑥2−𝑥−1 explains the pattern: the roots are (1±5)/2 whenever 1+1≠0, they coincide when 5=0, and they lie in 𝔽𝑝 when 5 is a square there, as 42=16=5 is in 𝔽11.

Exercise 7.4.21.routine Let

𝐴=(2100020000200003)∈ℚ4×4.

Find am𝐴(𝜆) and gm𝐴(𝜆) for every eigenvalue 𝜆, and decide whether 𝐴 is diagonalizable. Which single entry of 𝐴 above the diagonal can you change to make 𝐴 diagonalizable?

Solution.

Uses. Definition 7.4.11 (Algebraic multiplicity), Lemma 7.4.13 (Geometric multiplicity through a rank), Theorem 7.4.16 (Criteria for diagonalizability), Proposition 5.3.15 (Triangular matrices)

Hint 1. 𝐴 is upper triangular, so 𝑝𝐴 can be read off the diagonal.
Hint 2. For gm(2), compute the rank of 𝐴−2𝐼4.
Sketch. 𝑝𝐴=(𝑥−2)3(𝑥−3), so am(2)=3 and am(3)=1. 𝐴−2𝐼4 has rank 2 (columns 𝐞1 and 𝐞4), so gm(2)=2<3: not diagonalizable. Changing 𝑎12=1 to 0 makes 𝐴 diagonal. Any other change above the diagonal keeps 𝐴 triangular, so am(2)=3, and keeps columns 2 and 4 of 𝐴−2𝐼4 independent.

The matrix 𝑥𝐼4−𝐴 is upper triangular with diagonal 𝑥−2,𝑥−2,𝑥−2,𝑥−3, so 𝑝𝐴=(𝑥−2)3(𝑥−3) by Proposition 5.3.15. Hence the eigenvalues are 2 and 3, with am(2)=3 and am(3)=1. Next,

𝐴−2𝐼4=(0100000000000001)

has rank 2, so gm(2)=4−2=2 by Lemma 7.4.13, and 1≤gm(3)≤am(3)=1 gives gm(3)=1. Since gm(2)=2≠3=am(2), 𝐴 is not diagonalizable (Theorem 7.4.16); equivalently, gm(2)+gm(3)=3<4.

The only nonzero off-diagonal entry is 𝑎12=1, and replacing it by 0 gives a diagonal matrix. Changing any other entry above the diagonal keeps 𝐴 upper triangular with the same diagonal, so still 𝑝𝐴=(𝑥−2)3(𝑥−3) and am(2)=3. It also keeps 𝑎12=1, and column 2 of 𝐴−2𝐼4 is still 𝐞1 (the only entry above the diagonal in column 2 is 𝑎12), while column 4 still has last entry 1. These two columns are independent, so the rank stays at least 2, gm(2)≤2<3, and the matrix stays non-diagonalizable.

Below the diagonal the answer is different, because the characteristic polynomial can change. Setting 𝑎21=1 makes the top left block (2112), with eigenvalues 1 and 3 (Example 7.2.8); the new matrix has the eigenvalues 1,3,2,3 and is diagonalizable, as one checks with Theorem 7.4.16.

7.4.3Computing with a diagonalization

A diagonalization turns questions about 𝐴 into questions about the scalars on the diagonal of 𝐷, because conjugating by 𝑆 respects products.

Technique 7.4.22 (Diagonalize, solve, conjugate back). If 𝐴=𝑆𝐷𝑆−1 with 𝐷 diagonal, then 𝐴𝑘=𝑆𝐷𝑘𝑆−1, and more generally every expression built from 𝐴 by sums, products, and scalar multiples is 𝑆(⋯)𝑆−1 with 𝐷 in place of 𝐴. Diagonal matrices are computed entry by entry, so a problem about 𝐴 splits into 𝑛 problems about scalars. The recipe: diagonalize; solve the problem for 𝐷, one diagonal entry at a time; conjugate back by 𝑆. In operator language, pass to eigenbasis coordinates, where 𝑇 is 𝑛 independent scalings, solve there, and translate back.

Reach for it when a problem involves powers or polynomials of a matrix (closed forms for 𝐴𝑘, square roots, 𝑝(𝐴)), a linear recurrence 𝐱𝑘+1=𝐴𝐱𝑘, or a linear differential equation 𝐱′=𝐴𝐱. If 𝐴 is not diagonalizable, the same recipe works with a triangular or Jordan form in place of 𝐷 (Section 7.5, Section 10.3), at the price of harder computations with that form. The canonical example is Proposition 7.4.23.

Proposition 7.4.23 (Powers and polynomials of a diagonalizable matrix). Let 𝐴=𝑆𝐷𝑆−1 with 𝐷=diag(𝜆1,…,𝜆𝑛). Then for every integer 𝑘≥0 and every 𝑝∈𝐹[𝑥],

𝐴𝑘=𝑆(𝜆1𝑘⋱𝜆𝑛𝑘)𝑆−1,𝑝(𝐴)=𝑆(𝑝(𝜆1)⋱𝑝(𝜆𝑛))𝑆−1.

If every 𝜆𝑖 is nonzero, then 𝐴 is invertible, 𝐴−1=𝑆diag(𝜆1−1,…,𝜆𝑛−1)𝑆−1, and the formula for 𝐴𝑘 holds for all 𝑘∈ℤ.

Proof.routine

Uses. Proposition 2.4.14 (Algebra of matrices), Example 2.4.10 (Diagonal matrices), Definition 7.2.26 (Polynomials in an operator), Proposition 2.5.17 (Algebra of inverses)

Techniques. Diagonalize, solve, conjugate back (Technique 7.4.22)

Hint 1. Multiply out (𝑆𝐷𝑆−1)(𝑆𝐷𝑆−1). What happens in the middle?
Hint 2. Induct on 𝑘, then use linearity for 𝑝(𝐴)=∑𝑘𝑎𝑘𝐴𝑘. Diagonal matrices multiply entry by entry.
Sketch. The inner factors cancel: (𝑆𝐷𝑆−1)(𝑆𝐷𝑘𝑆−1)=𝑆𝐷𝑘+1𝑆−1, so by induction 𝐴𝑘=𝑆𝐷𝑘𝑆−1, and 𝐷𝑘=diag(𝜆𝑖𝑘) since diagonal matrices multiply entrywise. Summing, 𝑝(𝐴)=𝑆𝑝(𝐷)𝑆−1 with 𝑝(𝐷)=diag(𝑝(𝜆𝑖)). For the inverse, check that 𝑆diag(𝜆𝑖−1)𝑆−1 multiplies 𝐴 to 𝐼𝑛 on both sides.

We show 𝐴𝑘=𝑆𝐷𝑘𝑆−1 by induction on 𝑘. For 𝑘=0 both sides are 𝐼𝑛. If it holds for 𝑘, then by associativity (Proposition 2.4.14)

𝐴𝑘+1=𝐴𝐴𝑘=(𝑆𝐷𝑆−1)(𝑆𝐷𝑘𝑆−1)=𝑆𝐷(𝑆−1𝑆)𝐷𝑘𝑆−1=𝑆𝐷𝑘+1𝑆−1.

Diagonal matrices multiply entry by entry (Example 2.4.10), so 𝐷𝑘=diag(𝜆1𝑘,…,𝜆𝑛𝑘). For 𝑝=∑𝑘𝑎𝑘𝑥𝑘, distributivity gives

𝑝(𝐴)=∑𝑘𝑎𝑘𝑆𝐷𝑘𝑆−1=𝑆(∑𝑘𝑎𝑘𝐷𝑘)𝑆−1=𝑆diag(𝑝(𝜆1),…,𝑝(𝜆𝑛))𝑆−1.

If every 𝜆𝑖≠0, let 𝐷′=diag(𝜆1−1,…,𝜆𝑛−1), so that 𝐷𝐷′=𝐷′𝐷=𝐼𝑛. Then 𝐴(𝑆𝐷′𝑆−1)=𝑆𝐷𝐷′𝑆−1=𝐼𝑛 and likewise (𝑆𝐷′𝑆−1)𝐴=𝐼𝑛, so 𝐴 is invertible with 𝐴−1=𝑆𝐷′𝑆−1. Applying the first part to 𝐴−1=𝑆𝐷′𝑆−1 gives 𝐴−𝑘=𝑆𝐷′𝑘𝑆−1=𝑆diag(𝜆𝑖−𝑘)𝑆−1 for 𝑘≥0. □

Example 7.4.24 (A closed form for 𝐴𝑘). For 𝐴=(4−211)=𝑆diag(2,3)𝑆−1 from Example 7.4.4,

𝐴𝑘=(1211)(2𝑘003𝑘)(−121−1)=(2⋅3𝑘−2𝑘2𝑘+1−2⋅3𝑘3𝑘−2𝑘2𝑘+1−3𝑘)

for every 𝑘∈ℤ. At 𝑘=1 this returns 𝐴, and at 𝑘=0 it returns 𝐼2. Every entry of 𝐴𝑘 is a combination 𝛼2𝑘+𝛽3𝑘, and for large 𝑘 the term 3𝑘 dominates.

Often only 𝐴𝑘𝐯 is needed, and then 𝑆−1 is unnecessary. Expand 𝐯 in the eigenbasis, 𝐯=𝑐1𝐬1+𝑐2𝐬2; then 𝐴𝑘𝐯=𝑐12𝑘𝐬1+𝑐23𝑘𝐬2, because 𝐴𝑘𝐬𝑗=𝜆𝑗𝑘𝐬𝑗 (Proposition 7.2.28). For instance (0,1)=2(1,1)−(2,1), so 𝐴𝑘(0,1)=(2𝑘+1−2⋅3𝑘,2𝑘+1−3𝑘), the second column above.

The second paragraph of Example 7.4.24 proved a matrix identity by checking it on eigenvectors. For diagonalizable operators this is a general method.

Technique 7.4.25 (Test an identity on eigenvectors). To show that two operators built from a diagonalizable 𝑇 are equal, such as 𝑝(𝑇)=0 or 𝑆=𝑞(𝑇), it is enough to show that they agree on each vector of an eigenbasis. On an eigenvector 𝐯 with eigenvalue 𝜆, every polynomial in 𝑇 acts as the scalar 𝑝(𝜆) (Proposition 7.2.28), so the check becomes an identity between scalars, one eigenvalue at a time. This is Technique 2.2.7 with the most convenient basis.

Reach for it when the hypothesis says “diagonalizable” and the claim is an identity between polynomials in 𝑇, or between 𝑇 and another operator with the same eigenvectors. The canonical example is Proposition 7.4.26. Without diagonalizability the method fails: the shear 𝐽 satisfies 𝑝(𝐽)𝐯=𝑝(1)𝐯 on its only eigenvector line, but (𝐽−𝐼) kills that line without being 0.

Proposition 7.4.26 (Polynomials that kill a diagonalizable operator). Let 𝑇 be diagonalizable with distinct eigenvalues 𝜆1,…,𝜆𝑚, and let 𝑝∈𝐹[𝑥]. Then 𝑝(𝑇)=0 if and only if 𝑝(𝜆𝑖)=0 for 𝑖=1,…,𝑚. In particular,

(𝑇−𝜆1id)(𝑇−𝜆2id)⋯(𝑇−𝜆𝑚id)=0.

Proof.routine

Uses. Proposition 7.2.28 (Polynomials act on eigenvectors by their values), Theorem 7.4.2 (Diagonalizable means an eigenbasis), Proposition 7.2.27 (Evaluating polynomials at an operator)

Techniques. Test an identity on eigenvectors (Technique 7.4.25), Check it on a basis (Technique 2.2.7)

Hint 1. One direction is the last sentence of Proposition 7.2.28.
Hint 2. For the other, evaluate 𝑝(𝑇) on each vector of an eigenbasis.
Sketch. If 𝑝(𝑇)=0, every eigenvalue is a root of 𝑝 (Proposition 7.2.28). Conversely, on an eigenbasis vector 𝐛 with eigenvalue 𝜆𝑖, 𝑝(𝑇)𝐛=𝑝(𝜆𝑖)𝐛=𝟎; a linear map that vanishes on a basis is 0. The product is 𝑝(𝑇) for 𝑝=∏(𝑥−𝜆𝑖) (Proposition 7.2.27).
If 𝑝(𝑇)=0, then every eigenvalue of 𝑇 is a root of 𝑝, by Proposition 7.2.28. Conversely, suppose 𝑝(𝜆𝑖)=0 for every 𝑖, and let ℬ︀ be an eigenbasis. Each 𝐛∈ℬ︀ has some eigenvalue 𝜆𝑖, so 𝑝(𝑇)𝐛=𝑝(𝜆𝑖)𝐛=𝟎 by Proposition 7.2.28. The linear map 𝑝(𝑇) vanishes on a basis, so 𝑝(𝑇)=0. For the last claim, take 𝑝=(𝑥−𝜆1)⋯(𝑥−𝜆𝑚); by Proposition 7.2.27(a), 𝑝(𝑇) is the product of the operators 𝑇−𝜆𝑖id. □

For instance, a diagonalizable 𝑇 whose eigenvalues lie in {0,1} satisfies 𝑇2−𝑇=𝑇(𝑇−id)=0, so it is a projection; one with eigenvalues in {1,−1} satisfies 𝑇2=id. The converse of the last statement, that an operator killed by a product of distinct linear factors is diagonalizable, is also true. Exercise 7.4.49 proves it, and Section 7.6 builds the minimal polynomial on it.

Exercise 7.4.27.routine Let 𝐴∈𝐹𝑛×𝑛 be diagonalizable. Show that 𝐴𝖳 and 𝑝(𝐴), for every 𝑝∈𝐹[𝑥], are diagonalizable, and that 𝐴−1 is diagonalizable if 𝐴 is invertible. How are the eigenvalues of 𝑝(𝐴) related to those of 𝐴?

Solution.

Uses. Proposition 7.4.23 (Powers and polynomials of a diagonalizable matrix), Proposition 2.4.23 (Transpose of a product), Definition 7.4.1 (Diagonalizable), Theorem 7.4.16 (Criteria for diagonalizability)

Techniques. Diagonalize, solve, conjugate back (Technique 7.4.22)

Hint 1. Write 𝐴=𝑆𝐷𝑆−1 and apply Proposition 7.4.23.
Hint 2. For the transpose, transpose the factorization; the inverse of 𝑆𝖳 is (𝑆−1)𝖳.
Sketch. 𝑝(𝐴)=𝑆𝑝(𝐷)𝑆−1 with 𝑝(𝐷) diagonal; 𝐴−1=𝑆𝐷−1𝑆−1; 𝐴𝖳=(𝑆−1)𝖳𝐷𝑆𝖳=𝑅𝐷𝑅−1 with 𝑅=(𝑆−1)𝖳. The eigenvalues of 𝑝(𝐴) are the diagonal entries 𝑝(𝜆𝑖) of 𝑝(𝐷).

Write 𝐴=𝑆𝐷𝑆−1 with 𝐷=diag(𝜆1,…,𝜆𝑛). By Proposition 7.4.23, 𝑝(𝐴)=𝑆diag(𝑝(𝜆1),…,𝑝(𝜆𝑛))𝑆−1, which is a diagonalization of 𝑝(𝐴). If 𝐴 is invertible, every 𝜆𝑖 is nonzero (otherwise 𝐷, and with it 𝐴=𝑆𝐷𝑆−1, would not be invertible), and Proposition 7.4.23 gives 𝐴−1=𝑆diag(𝜆1−1,…,𝜆𝑛−1)𝑆−1. For the transpose, Proposition 2.4.23 gives

𝐴𝖳=(𝑆−1)𝖳𝐷𝖳𝑆𝖳=𝑅𝐷𝑅−1,𝑅=(𝑆−1)𝖳,

since 𝐷𝖳=𝐷 and 𝑅𝑆𝖳=(𝑆𝑆−1)𝖳=𝐼𝑛, and similarly 𝑆𝖳𝑅=𝐼𝑛.

By the last part of Theorem 7.4.16, the eigenvalues of a diagonalizable matrix are its diagonal entries in any diagonalization. So the eigenvalues of 𝑝(𝐴) are 𝑝(𝜆1),…,𝑝(𝜆𝑛), that is, the values 𝑝(𝜆) for the eigenvalues 𝜆 of 𝐴. Distinct eigenvalues of 𝐴 can have the same image: for 𝑝=𝑥2 and 𝐴=diag(1,−1), 𝑝(𝐴)=𝐼2.

7.4.4Linear recurrences

The Fibonacci numbers 𝑓0=0, 𝑓1=1, 𝑓𝑘+2=𝑓𝑘+1+𝑓𝑘 are generated by a matrix: with 𝐮𝑘=(𝑓𝑘+1,𝑓𝑘),

𝐮𝑘+1=(𝑓𝑘+1+𝑓𝑘𝑓𝑘+1)=(1110)𝐮𝑘=𝑄𝐮𝑘,

so 𝐮𝑘=𝑄𝑘𝐮0 with 𝐮0=(1,0), as in Exercise 2.4.47. Diagonalizing 𝑄 gives a closed formula.

Example 7.4.28 (Binet’s formula). By Example 7.4.20, 𝑄 has the distinct real eigenvalues 𝜑=(1+5)/2 and 𝜓=(1−5)/2, the roots of 𝑥2=𝑥+1. Eigenvectors are (𝜑,1) and (𝜓,1): for instance 𝑄(𝜑,1)=(𝜑+1,𝜑)=(𝜑2,𝜑)=𝜑(𝜑,1). Since 𝜑−𝜓=5,

𝐮0=(10)=15((𝜑1)−(𝜓1)),so𝐮𝑘=𝑄𝑘𝐮0=15(𝜑𝑘(𝜑1)−𝜓𝑘(𝜓1)).

The second coordinate is Binet’s formula

𝑓𝑘=𝜑𝑘−𝜓𝑘5for all𝑘≥0.

Since |𝜓|<1, the term 𝜓𝑘/5 has absolute value at most 1/5<1/2, so 𝑓𝑘 is the integer nearest to 𝜑𝑘/5. For instance 𝜑10/5≈55.004, and 𝑓10=55. The eigenvalue 𝜑≈1.618 of larger absolute value controls the growth.

The same formula appeared in Exercise 1.6.39 without any matrices, by finding a basis of geometric sequences of the space of all solutions. The two arguments are the same argument seen from two sides: the map 𝐬↦(𝑠1,𝑠0) identifies the solution space with ℝ2, and under it, shifting a sequence one step to the left becomes multiplication by 𝑄. The eigenvectors (𝜑,1) and (𝜓,1) correspond to the geometric sequences (𝜑𝑘) and (𝜓𝑘), which the left shift multiplies by 𝜑 and 𝜓 (Example 7.2.20). We now do this for every recurrence of finite order.

Fix 𝑘≥1 and 𝑐0,…,𝑐𝑘−1∈𝐹, and consider sequences 𝐬=(𝑠0,𝑠1,𝑠2,…) in 𝐹 with

𝑠𝑗+𝑘=𝑐𝑘−1𝑠𝑗+𝑘−1+⋯+𝑐1𝑠𝑗+1+𝑐0𝑠𝑗for all𝑗≥0.
(7.4.2)

The characteristic polynomial of the recurrence is 𝑞=𝑥𝑘−𝑐𝑘−1𝑥𝑘−1−⋯−𝑐1𝑥−𝑐0; a geometric sequence (𝑟𝑗)𝑗≥0 satisfies (7.4.2) exactly when 𝑟𝑗𝑞(𝑟)=0 for all 𝑗, that is, when 𝑞(𝑟)=0. Let 𝑊 be the set of solutions of (7.4.2), a subspace of the sequence space 𝐹𝑋, 𝑋={0,1,2,…} (the recurrence is linear in 𝐬). Let 𝐶∈𝐹𝑘×𝑘 be the matrix with

𝐶(𝑥1⋮𝑥𝑘−1𝑥𝑘)=(𝑥2⋮𝑥𝑘𝑐0𝑥1+𝑐1𝑥2+⋯+𝑐𝑘−1𝑥𝑘),that is,𝐶=(01⋱⋱01𝑐0𝑐1⋯𝑐𝑘−1),

a companion matrix of 𝑞 (compare Definition 7.3.23; conventions differ by a transpose and the order of the basis). For 𝑘=2, 𝑐0=𝑐1=1, it is 𝑄 with the basis reversed.

Proposition 7.4.29 (The matrix of a recurrence). With the notation above:

  1. (a)if 𝐬∈𝑊 and 𝐮𝑗=(𝑠𝑗,𝑠𝑗+1,…,𝑠𝑗+𝑘−1), then 𝐮𝑗+1=𝐶𝐮𝑗, so 𝐮𝑗=𝐶𝑗𝐮0;
  2. (b)the eigenvalues of 𝐶 are the roots of 𝑞 in 𝐹, and for each root 𝑟, 𝐸𝑟(𝐶)=span((1,𝑟,𝑟2,…,𝑟𝑘−1));
  3. (c)𝐶 is diagonalizable if and only if 𝑞 has 𝑘 distinct roots in 𝐹.

Proof.standard

Uses. Definition 7.2.6 (Eigenvalue, eigenvector), Definition 7.2.13 (Eigenspace), Theorem 7.4.16 (Criteria for diagonalizability)

Hint 1. (a) is the definition of 𝐶 and the recurrence at index 𝑗.
Hint 2. For (b), write out 𝐶𝐯=𝜆𝐯 coordinate by coordinate. The first 𝑘−1 equations determine 𝐯 from 𝑣1.
Hint 3. For (c), every eigenspace is a line. Use criterion (III) of Theorem 7.4.16.
Sketch. 𝐶𝐯=𝜆𝐯 says 𝑣𝑖+1=𝜆𝑣𝑖 for 𝑖<𝑘, so 𝐯=𝑣1(1,𝜆,…,𝜆𝑘−1), and the last coordinate then says 𝑣1𝑞(𝜆)=0. So eigenvalues are roots of 𝑞, each with a one-dimensional eigenspace. Then ∑gm= the number of distinct roots of 𝑞 in 𝐹, which is 𝑘 exactly when 𝑞 has 𝑘 distinct roots.

(a) The first 𝑘−1 coordinates of 𝐶𝐮𝑗 are 𝑠𝑗+1,…,𝑠𝑗+𝑘−1, and the last is 𝑐0𝑠𝑗+⋯+𝑐𝑘−1𝑠𝑗+𝑘−1=𝑠𝑗+𝑘 by (7.4.2). So 𝐶𝐮𝑗=𝐮𝑗+1, and 𝐮𝑗=𝐶𝑗𝐮0 by induction.

(b) Let 𝜆∈𝐹 and 𝐯=(𝑣1,…,𝑣𝑘)∈𝐹𝑘. Comparing coordinates, 𝐶𝐯=𝜆𝐯 says

𝑣2=𝜆𝑣1,𝑣3=𝜆𝑣2,…,𝑣𝑘=𝜆𝑣𝑘−1,𝑐0𝑣1+𝑐1𝑣2+⋯+𝑐𝑘−1𝑣𝑘=𝜆𝑣𝑘.

The first 𝑘−1 equations hold if and only if 𝑣𝑖=𝜆𝑖−1𝑣1 for all 𝑖, that is, 𝐯=𝑣1𝐰𝜆 with 𝐰𝜆=(1,𝜆,…,𝜆𝑘−1). Substituting into the last equation gives 𝑣1(𝑐0+𝑐1𝜆+⋯+𝑐𝑘−1𝜆𝑘−1)=𝑣1𝜆𝑘, that is, 𝑣1𝑞(𝜆)=0. So if 𝑞(𝜆)≠0, then 𝑣1=0 and 𝐯=𝟎: 𝜆 is not an eigenvalue. If 𝑞(𝜆)=0, the solutions are exactly the multiples of 𝐰𝜆≠𝟎, so 𝜆 is an eigenvalue with 𝐸𝜆(𝐶)=span(𝐰𝜆).

(c) By (b), gm𝐶(𝑟)=1 for each root 𝑟 of 𝑞 in 𝐹, so the sum of the geometric multiplicities is the number of distinct roots of 𝑞 in 𝐹. By Theorem 7.4.16, (I) ⇔ (III), 𝐶 is diagonalizable if and only if this number is 𝑘. □

So 𝐶 never has two independent eigenvectors for the same eigenvalue: a repeated root of 𝑞 always destroys diagonalizability. When the roots are distinct, the solutions are combinations of geometric sequences.

Theorem 7.4.30 (Recurrences with distinct characteristic roots). Suppose 𝑞 has 𝑘 distinct roots 𝑟1,…,𝑟𝑘∈𝐹. Then the geometric sequences (𝑟𝑖𝑗)𝑗≥0, 𝑖=1,…,𝑘, form a basis of 𝑊. That is, every solution of (7.4.2) is

𝑠𝑗=𝑎1𝑟1𝑗+𝑎2𝑟2𝑗+⋯+𝑎𝑘𝑟𝑘𝑗(𝑗≥0)

for unique 𝑎1,…,𝑎𝑘∈𝐹, which are determined by 𝑠0,…,𝑠𝑘−1.

Proof (diagonalize the matrix of the recurrence).standard

Uses. Proposition 7.4.29 (The matrix of a recurrence), Theorem 7.2.21 (Eigenvectors for distinct eigenvalues are independent), Corollary 7.4.5 (Distinct eigenvalues give a diagonalization), Proposition 7.2.28 (Polynomials act on eigenvectors by their values)

Techniques. Diagonalize, solve, conjugate back (Technique 7.4.22)

Hint 1. By Proposition 7.4.29, 𝐮𝑗=𝐶𝑗𝐮0 and 𝐶 has the eigenbasis 𝐰𝑟1,…,𝐰𝑟𝑘. Expand 𝐮0 in it.
Hint 2. Then 𝐶𝑗𝐮0=∑𝑖𝑎𝑖𝑟𝑖𝑗𝐰𝑟𝑖. Read off the first coordinate, and note that each 𝐰𝑟𝑖 has first coordinate 1.
Hint 3. For uniqueness, 𝑠0,…,𝑠𝑘−1 are the coordinates of 𝐮0, and coordinates in a basis are unique.
Sketch. The vectors 𝐰𝑟𝑖=(1,𝑟𝑖,…,𝑟𝑖𝑘−1) are eigenvectors for distinct eigenvalues, so an eigenbasis. Write 𝐮0=∑𝑎𝑖𝐰𝑟𝑖; then 𝐮𝑗=𝐶𝑗𝐮0=∑𝑎𝑖𝑟𝑖𝑗𝐰𝑟𝑖, whose first coordinate is 𝑠𝑗=∑𝑎𝑖𝑟𝑖𝑗. So the geometric sequences span 𝑊; they lie in 𝑊 since 𝑞(𝑟𝑖)=0. If ∑𝑎𝑖𝑟𝑖𝑗=0 for all 𝑗, then ∑𝑎𝑖𝐰𝑟𝑖=𝟎 (the cases 𝑗<𝑘), so all 𝑎𝑖=0.

Write 𝐰𝑖=(1,𝑟𝑖,𝑟𝑖2,…,𝑟𝑖𝑘−1). By Proposition 7.4.29(b), 𝐰𝑖 is an eigenvector of 𝐶 with eigenvalue 𝑟𝑖, and since the 𝑟𝑖 are distinct, 𝐰1,…,𝐰𝑘 is a basis of 𝐹𝑘 (Corollary 7.4.5). Each geometric sequence (𝑟𝑖𝑗)𝑗 lies in 𝑊, because 𝑞(𝑟𝑖)=0.

Let 𝐬∈𝑊, and write 𝐮0=(𝑠0,…,𝑠𝑘−1)=𝑎1𝐰1+⋯+𝑎𝑘𝐰𝑘. By Proposition 7.4.29(a) and Proposition 7.2.28,

𝐮𝑗=𝐶𝑗𝐮0=∑𝑖=1𝑘𝑎𝑖𝐶𝑗𝐰𝑖=∑𝑖=1𝑘𝑎𝑖𝑟𝑖𝑗𝐰𝑖.

The first coordinate of 𝐮𝑗 is 𝑠𝑗, and the first coordinate of each 𝐰𝑖 is 1, so 𝑠𝑗=∑𝑖𝑎𝑖𝑟𝑖𝑗 for all 𝑗. Thus the geometric sequences span 𝑊.

They are independent: if ∑𝑖𝑎𝑖𝑟𝑖𝑗=0 for all 𝑗≥0, then the cases 𝑗=0,…,𝑘−1 say ∑𝑖𝑎𝑖𝐰𝑖=𝟎, so all 𝑎𝑖=0. Hence they form a basis of 𝑊, the 𝑎𝑖 are unique, and they are the coordinates of 𝐮0=(𝑠0,…,𝑠𝑘−1) in the basis 𝐰1,…,𝐰𝑘. □

Proof (count dimensions in the solution space).standard

Uses. Theorem 4.1.10 (Independent lists are short), Example 7.2.24 (Independence for free), Definition 4.5.1 (Isomorphism)

Techniques. Dimension counting (Technique 4.2.24)

Hint 1. Show that 𝐬↦(𝑠0,…,𝑠𝑘−1) is an isomorphism 𝑊→𝐹𝑘. So dim𝑊=𝑘.
Hint 2. The 𝑘 geometric sequences lie in 𝑊. Why are they independent? Think of the left shift.
Hint 3. 𝑘 independent vectors in a 𝑘-dimensional space form a basis.
Sketch. A solution is determined by its first 𝑘 terms, and any first 𝑘 terms extend (define 𝑠𝑗+𝑘 by (7.4.2)), so 𝐬↦(𝑠0,…,𝑠𝑘−1) is a linear bijection and dim𝑊=𝑘. The geometric sequences lie in 𝑊 and are eigenvectors of the left shift with distinct eigenvalues 𝑟𝑖, hence independent; 𝑘 of them form a basis.

Let Φ:𝑊→𝐹𝑘, Φ(𝐬)=(𝑠0,…,𝑠𝑘−1). It is linear. It is injective: if 𝑠0=⋯=𝑠𝑘−1=0, then (7.4.2) and induction on 𝑗 give 𝑠𝑗=0 for all 𝑗. It is surjective: given (𝑡0,…,𝑡𝑘−1), set 𝑠𝑗=𝑡𝑗 for 𝑗<𝑘 and define 𝑠𝑗+𝑘 for 𝑗=0,1,2,… by (7.4.2); the resulting sequence lies in 𝑊. So Φ is an isomorphism (Definition 4.5.1) and dim𝑊=𝑘.

Since 𝑞(𝑟𝑖)=0, each geometric sequence (𝑟𝑖𝑗)𝑗 lies in 𝑊. These 𝑘 sequences are eigenvectors of the left shift on 𝐹𝑋 with the distinct eigenvalues 𝑟1,…,𝑟𝑘, so they are independent (Example 7.2.24(a)). By Theorem 4.1.10(b) they form a basis of the 𝑘-dimensional space 𝑊, and the coefficients 𝑎𝑖 of a solution are its coordinates in this basis, which are unique. Applying Φ, they are determined by 𝑠0,…,𝑠𝑘−1. □

The first proof computes: it finds the 𝑎𝑖 by solving one linear system, 𝐮0=∑𝑎𝑖𝐰𝑖, whose matrix has rows (1,…,1),(𝑟1,…,𝑟𝑘),… (a Vandermonde matrix, Theorem 5.5.21). The second proof needs no matrix at all, only the dimension of 𝑊; it is Technique 4.2.24, and it is the proof that generalizes to linear differential equations (Section 15.2).

Remark 7.4.31 (Repeated roots). If 𝑞 has a repeated root, 𝐶 is not diagonalizable (Proposition 7.4.29(c)), and the geometric sequences no longer span 𝑊. For 𝑠𝑗+2=2𝑠𝑗+1−𝑠𝑗 we have 𝑞=(𝑥−1)2, the only geometric solution is the constant sequence, and 𝑊 consists of the arithmetic sequences 𝑠𝑗=𝑎+𝑏𝑗. In general, a root 𝑟 of multiplicity 𝑚 with 𝑟≠0 contributes the solutions (𝑗𝑖𝑟𝑗)𝑗 for 𝑖<𝑚 (when 𝐹 has characteristic 0); the source of these is the Jordan form of 𝐶 (Section 10.3), and Section 15.2 develops them.

7.4.5Discrete dynamical systems

A discrete linear dynamical system is a sequence of states 𝐱0,𝐱1,𝐱2,… in 𝐹𝑛 with 𝐱𝑘+1=𝐴𝐱𝑘 for a fixed 𝐴∈𝐹𝑛×𝑛, so 𝐱𝑘=𝐴𝑘𝐱0. If 𝐴 has an eigenbasis 𝐬1,…,𝐬𝑛 with eigenvalues 𝜆1,…,𝜆𝑛, and 𝐱0=𝑐1𝐬1+⋯+𝑐𝑛𝐬𝑛, then

𝐱𝑘=𝑐1𝜆1𝑘𝐬1+𝑐2𝜆2𝑘𝐬2+⋯+𝑐𝑛𝜆𝑛𝑘𝐬𝑛.
(7.4.3)

Each eigen-component evolves on its own, by repeated multiplication by its eigenvalue. Over ℝ or ℂ, a component with |𝜆𝑖|<1 dies out, one with |𝜆𝑖|>1 grows, and the components with the largest |𝜆𝑖| eventually dominate.

Example 7.4.32 (Migration between two cities). Each year, 20% of the population of city 1 moves to city 2, and 30% of the population of city 2 moves to city 1; nobody else moves, and the total population is constant. With 𝐱𝑘=(𝑎𝑘,𝑏𝑘) the populations in year 𝑘,

𝐱𝑘+1=𝐴𝐱𝑘,𝐴=(0.80.30.20.7).

Here det(𝐴−𝜆𝐼2)=𝜆2−1.5𝜆+0.5=(𝜆−1)(𝜆−0.5), with eigenvectors (3,2) for 1 and (1,−1) for 0.5. Writing 𝐱0=(𝑎,𝑏)=𝑐1(3,2)+𝑐2(1,−1) gives 𝑐1=(𝑎+𝑏)/5 and 𝑐2=(2𝑎−3𝑏)/5, so by (7.4.3)

𝐱𝑘=𝑎+𝑏5(32)+2𝑎−3𝑏5⋅2−𝑘(1−1)→𝑎+𝑏5(32)(𝑘→∞).

Whatever the initial distribution, the populations approach a 60:40 split of the total 𝑎+𝑏. The limit is the component of 𝐱0 in 𝐸1, and the deviation from it halves every year. The columns of 𝐴 add up to 1, which is why 1 is an eigenvalue (the rows of 𝐴−𝐼2 add up to zero) and why the total is preserved. Such matrices are the subject of Section 11.1.

Proposition 7.4.33 (Long-run behavior). Let 𝐴∈ℂ𝑛×𝑛 be diagonalizable, with an eigenbasis 𝐬1,…,𝐬𝑛 and corresponding eigenvalues 𝜆1,…,𝜆𝑛. Let 𝐱0=𝑐1𝐬1+⋯+𝑐𝑛𝐬𝑛 and 𝐱𝑘=𝐴𝑘𝐱0. Limits of vectors are taken coordinate by coordinate.

  1. (a)If |𝜆𝑖|<1 for every 𝑖, then 𝐱𝑘→𝟎 for every 𝐱0, and 𝐴𝑘→0 entry by entry.
  2. (b)Suppose 𝜇≠0 is an eigenvalue with |𝜆𝑖|<|𝜇| for every 𝜆𝑖≠𝜇. Then

    𝜇−𝑘𝐱𝑘→∑𝑖:𝜆𝑖=𝜇𝑐𝑖𝐬𝑖,

    the component of 𝐱0 in 𝐸𝜇. In particular, if 𝜇=1, then 𝐱𝑘 converges to the component of 𝐱0 in 𝐸1.

Proof.standard

Uses. Proposition 7.2.28 (Polynomials act on eigenvectors by their values), Theorem 7.4.16 (Criteria for diagonalizability)

Techniques. Diagonalize, solve, conjugate back (Technique 7.4.22)

Hint 1. Start from (7.4.3). Each coordinate of 𝐱𝑘 is a fixed linear combination of the 𝜆𝑖𝑘.
Hint 2. For complex 𝑧 with |𝑧|<1, |𝑧𝑘|=|𝑧|𝑘→0.
Hint 3. For (b), divide (7.4.3) by 𝜇𝑘 and split the sum according to whether 𝜆𝑖=𝜇.
Sketch. By (7.4.3), 𝐱𝑘=∑𝑐𝑖𝜆𝑖𝑘𝐬𝑖, a combination with fixed vectors and coefficients 𝜆𝑖𝑘. If all |𝜆𝑖|<1, each 𝜆𝑖𝑘→0. For (b), 𝜇−𝑘𝐱𝑘=∑𝜆𝑖=𝜇𝑐𝑖𝐬𝑖+∑𝜆𝑖≠𝜇𝑐𝑖(𝜆𝑖/𝜇)𝑘𝐬𝑖, and the second sum tends to 𝟎. Apply (a) to the columns 𝐴𝑘𝐞𝑗 to get 𝐴𝑘→0.

Since 𝐴𝑘𝐬𝑖=𝜆𝑖𝑘𝐬𝑖 (Proposition 7.2.28), (7.4.3) holds: 𝐱𝑘=∑𝑖𝑐𝑖𝜆𝑖𝑘𝐬𝑖. Each coordinate of 𝐱𝑘 is therefore ∑𝑖𝑐𝑖𝜆𝑖𝑘𝑠𝑖,𝑙, where 𝑠𝑖,𝑙 is the 𝑙-th coordinate of 𝐬𝑖, a sum of 𝑛 fixed complex numbers times 𝑘-th powers.

(a) If |𝜆𝑖|<1, then |𝜆𝑖𝑘|=|𝜆𝑖|𝑘→0, so 𝜆𝑖𝑘→0, and each coordinate of 𝐱𝑘 tends to 0. The 𝑗-th column of 𝐴𝑘 is 𝐴𝑘𝐞𝑗, which is the case 𝐱0=𝐞𝑗, so every entry of 𝐴𝑘 tends to 0.

(b) Dividing by 𝜇𝑘,

𝜇−𝑘𝐱𝑘=∑𝑖:𝜆𝑖=𝜇𝑐𝑖𝐬𝑖+∑𝑖:𝜆𝑖≠𝜇𝑐𝑖(𝜆𝑖𝜇)𝑘𝐬𝑖.

The first sum does not depend on 𝑘, and in the second |𝜆𝑖/𝜇|<1, so each term tends to 𝟎 as in (a). The first sum lies in 𝐸𝜇, and the second in the sum of the other eigenspaces; since ℂ𝑛 is the direct sum of the eigenspaces (Theorem 7.4.16), the first sum is the component of 𝐱0 in 𝐸𝜇. □

In Example 7.4.32, 𝜇=1 and the other eigenvalue is 0.5. In Example 7.4.28, 𝜇=𝜑 dominates 𝜓, so 𝜑−𝑘𝐮𝑘→(𝜑,1)/5: the ratio 𝑓𝑘+1/𝑓𝑘 of consecutive Fibonacci numbers tends to 𝜑. When 𝐴 is not diagonalizable, (7.4.3) is replaced by terms 𝑘𝑖𝜆𝑘 coming from the Jordan form, and (a) still holds, but the proof needs Section 10.3.

7.4.6Systems of linear differential equations Beyond 21-242

Differential equations are not part of 21-242.

The continuous analog of 𝐱𝑘+1=𝐴𝐱𝑘 is the system 𝐱′=𝐴𝐱, where 𝐱:ℝ→ℝ𝑛 is differentiable and the derivative is taken coordinate by coordinate. For 𝑛=1 this is 𝑥′=𝜆𝑥, whose solutions are 𝑥(𝑡)=𝑥(0)𝑒𝜆𝑡: if 𝑥′=𝜆𝑥, then (𝑒−𝜆𝑡𝑥(𝑡))′=𝑒−𝜆𝑡(𝑥′−𝜆𝑥)=0, so 𝑒−𝜆𝑡𝑥(𝑡) is constant. A diagonalization decouples a system into 𝑛 such equations.

Proposition 7.4.34 (Solving 𝐱′=𝐴𝐱). Let 𝐴=𝑆𝐷𝑆−1∈ℝ𝑛×𝑛 with 𝐷=diag(𝜆1,…,𝜆𝑛) real, and let 𝐬1,…,𝐬𝑛 be the columns of 𝑆. A differentiable 𝐱:ℝ→ℝ𝑛 satisfies 𝐱′=𝐴𝐱 if and only if

𝐱(𝑡)=𝑐1𝑒𝜆1𝑡𝐬1+𝑐2𝑒𝜆2𝑡𝐬2+⋯+𝑐𝑛𝑒𝜆𝑛𝑡𝐬𝑛

for some 𝑐1,…,𝑐𝑛∈ℝ, namely (𝑐1,…,𝑐𝑛)=𝑆−1𝐱(0).

Proof.standard

Uses. Theorem 7.4.2 (Diagonalizable means an eigenbasis)

Techniques. Diagonalize, solve, conjugate back (Technique 7.4.22)

Hint 1. Change variables: let 𝐲(𝑡)=𝑆−1𝐱(𝑡). Differentiation commutes with multiplication by a constant matrix.
Hint 2. Then 𝐲′=𝐷𝐲, which is 𝑛 separate equations 𝑦𝑖′=𝜆𝑖𝑦𝑖.
Sketch. With 𝐲=𝑆−1𝐱, 𝐲′=𝑆−1𝐱′ (each coordinate is a fixed linear combination), so 𝐱′=𝐴𝐱 iff 𝐲′=𝐷𝐲, iff 𝑦𝑖′=𝜆𝑖𝑦𝑖 for all 𝑖, iff 𝑦𝑖(𝑡)=𝑐𝑖𝑒𝜆𝑖𝑡. Then 𝐱=𝑆𝐲=∑𝑐𝑖𝑒𝜆𝑖𝑡𝐬𝑖.

Let 𝐲(𝑡)=𝑆−1𝐱(𝑡). Each coordinate of 𝐲 is a fixed linear combination of the coordinates of 𝐱, so 𝐲 is differentiable with 𝐲′=𝑆−1𝐱′, and similarly 𝐱′=𝑆𝐲′. Hence

𝐱′=𝐴𝐱⟺𝑆𝐲′=𝑆𝐷𝑆−1𝑆𝐲⟺𝐲′=𝐷𝐲,

multiplying by 𝑆−1 or 𝑆. The system 𝐲′=𝐷𝐲 consists of the 𝑛 separate equations 𝑦𝑖′=𝜆𝑖𝑦𝑖, whose solutions are exactly 𝑦𝑖(𝑡)=𝑐𝑖𝑒𝜆𝑖𝑡 with 𝑐𝑖=𝑦𝑖(0), by the case 𝑛=1 above. So 𝐱′=𝐴𝐱 if and only if 𝐱=𝑆𝐲=∑𝑖𝑐𝑖𝑒𝜆𝑖𝑡𝐬𝑖, and then (𝑐1,…,𝑐𝑛)=𝐲(0)=𝑆−1𝐱(0). □

For 𝐴=(4−211) and 𝐱(0)=(1,0), Example 7.4.4 gives 𝑆−1𝐱(0)=(−1,1), so

𝐱(𝑡)=−𝑒2𝑡(11)+𝑒3𝑡(21)=(2𝑒3𝑡−𝑒2𝑡𝑒3𝑡−𝑒2𝑡).

In matrix form, the solution is 𝐱(𝑡)=𝑆diag(𝑒𝜆1𝑡,…,𝑒𝜆𝑛𝑡)𝑆−1𝐱(0), and the matrix in front of 𝐱(0) is the matrix exponential 𝑒𝑡𝐴, which Technique 7.4.22 computes from 𝑒𝑡𝐷. Complex eigenvalues produce oscillations: for the rotation generator 𝐴=(0−110), diagonalizable over ℂ with eigenvalues ±𝑖, the solutions are 𝐱(𝑡)=(𝛼cos𝑡−𝛽sin𝑡,𝛼sin𝑡+𝛽cos𝑡), circles traversed at unit angular speed. A non-diagonalizable 𝐴 produces terms 𝑡𝑗𝑒𝜆𝑡, which come from the Jordan form (Section 10.3).

Exercises

Exercise 7.4.35.routine For each matrix over ℚ, find an invertible 𝑆 and a diagonal 𝐷 with 𝐴=𝑆𝐷𝑆−1, or show that none exist.

  1. (a)𝐴=(1232).
  2. (b)𝐴=(31−11).
  3. (c)𝐴=(001010100).

Solution.

Uses. Algorithm 7.4.18 (Diagonalizing a matrix), Theorem 7.2.9 (Eigenvalues and invertibility), Proposition 7.4.7 (One eigenvalue)

Hint 1. Follow Algorithm 7.4.18: find the eigenvalues from det(𝐴−𝜆𝐼)=0, then a basis of each eigenspace.
Hint 2. In (b) there is only one eigenvalue. In (c), 𝐴 swaps 𝐞1 and 𝐞3 and fixes 𝐞2; guess the eigenvectors.
Sketch. (a) (𝜆−4)(𝜆+1); 𝐸4=span((2,3)), 𝐸−1=span((1,−1)); 𝑆=(213−1), 𝐷=diag(4,−1). (b) (𝜆−2)2, and 𝐴≠2𝐼: not diagonalizable. (c) 𝐸1=span(𝐞2,𝐞1+𝐞3), 𝐸−1=span(𝐞1−𝐞3).

(a) det(𝐴−𝜆𝐼2)=(1−𝜆)(2−𝜆)−6=𝜆2−3𝜆−4=(𝜆−4)(𝜆+1). Then 𝐴−4𝐼2=(−323−2) has kernel span((2,3)), and 𝐴+𝐼2=(2233) has kernel span((1,−1)). The eigenvalues are distinct, so

𝑆=(213−1),𝐷=(400−1).

Check: 𝐴𝑆=(8−1121)=𝑆𝐷.

(b) det(𝐴−𝜆𝐼2)=(3−𝜆)(1−𝜆)+1=𝜆2−4𝜆+4=(𝜆−2)2, so 2 is the only eigenvalue. Since 𝐴≠2𝐼2, 𝐴 is not diagonalizable (Proposition 7.4.7). Indeed 𝐴−2𝐼2=(11−1−1) has rank 1, so 𝐸2=span((1,−1)) is a line.

(c) 𝐴(𝑥,𝑦,𝑧)=(𝑧,𝑦,𝑥). The vectors 𝐞2 and 𝐞1+𝐞3 are fixed, and 𝐴(𝐞1−𝐞3)=𝐞3−𝐞1. The matrix 𝑆 below, whose columns are these three vectors, has determinant 2 (expand along its second row), so they form an eigenbasis and

𝑆=(01110001−1),𝐷=(10001000−1).

To see that nothing was missed, det(𝐴−𝜆𝐼3)=(1−𝜆)(𝜆2−1)=−(𝜆−1)2(𝜆+1) by cofactor expansion along the middle row, and gm(1)=2=am(1).

Exercise 7.4.36.standard Let 𝐴=(7−103−4)∈ℝ2×2.

  1. (a)Find a formula for 𝐴𝑘 valid for every 𝑘∈ℤ.
  2. (b)Find all 𝐵∈ℝ2×2 with 𝐵2=𝐴.

Solution.

Uses. Proposition 7.4.23 (Powers and polynomials of a diagonalizable matrix), Proposition 7.2.30 (Kernels and images of commuting operators are invariant), Proposition 7.2.7 (Invariant lines are spanned by eigenvectors), Theorem 7.4.2 (Diagonalizable means an eigenbasis)

Techniques. Diagonalize, solve, conjugate back (Technique 7.4.22), Look for invariant subspaces among kernels and images of commuting operators (Technique 7.2.31)

Hint 1. Diagonalize: the eigenvalues are 1 and 2.
Hint 2. For (b), any 𝐵 with 𝐵2=𝐴 commutes with 𝐴. What does that say about the eigenlines of 𝐴?
Hint 3. 𝐵 preserves each eigenline of 𝐴, so it is diagonal in the eigenbasis: 𝑆−1𝐵𝑆=diag(𝛽1,𝛽2) with 𝛽12=1, 𝛽22=2.
Sketch. 𝐴=𝑆diag(1,2)𝑆−1 with 𝑆=(5231), 𝑆−1=(−123−5), so 𝐴𝑘=(6⋅2𝑘−510−10⋅2𝑘3⋅2𝑘−36−5⋅2𝑘). A square root 𝐵 commutes with 𝐴, so it preserves the eigenlines 𝐸1,𝐸2 and is diagonal in the eigenbasis with entries ±1,±2: exactly four square roots.

(a) det(𝐴−𝜆𝐼2)=(7−𝜆)(−4−𝜆)+30=𝜆2−3𝜆+2=(𝜆−1)(𝜆−2). Then 𝐴−𝐼2=(6−103−5) has kernel span((5,3)), and 𝐴−2𝐼2=(5−103−6) has kernel span((2,1)). So 𝐴=𝑆𝐷𝑆−1 with

𝑆=(5231),𝐷=(1002),𝑆−1=(−123−5),

and both eigenvalues are nonzero, so by Proposition 7.4.23, for every 𝑘∈ℤ,

𝐴𝑘=(52⋅2𝑘32𝑘)(−123−5)=(6⋅2𝑘−510−10⋅2𝑘3⋅2𝑘−36−5⋅2𝑘).

At 𝑘=1 this is 𝐴, and at 𝑘=−1 it is (−25−3/27/2), which one checks is 𝐴−1.

(b) Let 𝐬1=(5,3) and 𝐬2=(2,1). If 𝐵2=𝐴, then 𝐵𝐴=𝐵3=𝐴𝐵. By Proposition 7.2.30(b), each eigenspace 𝐸1(𝐴)=span(𝐬1) and 𝐸2(𝐴)=span(𝐬2) is 𝐵-invariant, so by Proposition 7.2.7, 𝐵𝐬1=𝛽1𝐬1 and 𝐵𝐬2=𝛽2𝐬2 for some 𝛽1,𝛽2∈ℝ. Then 𝐴𝐬𝑗=𝐵2𝐬𝑗=𝛽𝑗2𝐬𝑗, so 𝛽12=1 and 𝛽22=2. By Theorem 7.4.2, 𝐵=𝑆diag(𝛽1,𝛽2)𝑆−1 with 𝛽1=±1 and 𝛽2=±2. Conversely each of these four matrices squares to 𝑆diag(1,2)𝑆−1=𝐴. So 𝐴 has exactly four real square roots; for 𝛽1=1, 𝛽2=2,

𝐵=(52232)(−123−5)=(62−510−10232−36−52),

and the four square roots are 𝑆diag(𝛽1,𝛽2)𝑆−1 with (𝛽1,𝛽2)∈{±1}×{±2}, namely 𝐵, −𝐵, and the matrices obtained from 𝐵 and −𝐵 by replacing 2 with −2.

Exercise 7.4.37.standard For which 𝑎,𝑏,𝑐∈𝐹 is the matrix

𝐴=(1𝑎𝑏01𝑐002)

diagonalizable?

Solution.

Uses. Theorem 7.4.16 (Criteria for diagonalizability), Lemma 7.4.13 (Geometric multiplicity through a rank), Theorem 7.4.15 (Geometric multiplicity is at most algebraic multiplicity), Proposition 5.3.15 (Triangular matrices)

Hint 1. The characteristic polynomial does not depend on 𝑎,𝑏,𝑐. Which eigenvalue could cause trouble?
Hint 2. Compute gm(1)=3−rank(𝐴−𝐼3) as a function of 𝑎,𝑏,𝑐.
Sketch. 𝑝𝐴=(𝑥−1)2(𝑥−2), so 𝐴 is diagonalizable iff gm(1)=2, iff rank(𝐴−𝐼3)=1. The rows of 𝐴−𝐼3 are (0,𝑎,𝑏), (0,0,𝑐), (0,0,1); the last is nonzero and the first is independent of it iff 𝑎≠0. So 𝐴 is diagonalizable iff 𝑎=0, for all 𝑏,𝑐.

𝐴 is upper triangular, so 𝑝𝐴=(𝑥−1)2(𝑥−2) (Proposition 5.3.15), which splits, with am(1)=2 and am(2)=1. Since 1≤gm(2)≤am(2)=1 (Theorem 7.4.15), gm(2)=1 always. By Theorem 7.4.16, 𝐴 is diagonalizable if and only if gm(1)=2, that is, rank(𝐴−𝐼3)=1 (Lemma 7.4.13). The rows of

𝐴−𝐼3=(0𝑎𝑏00𝑐001)

are (0,𝑎,𝑏), (0,0,𝑐), (0,0,1). The second is a multiple of the third, so the rank is the rank of (0,𝑎,𝑏),(0,0,1), which is 2 if 𝑎≠0 and 1 if 𝑎=0. So 𝐴 is diagonalizable if and only if 𝑎=0, whatever 𝑏 and 𝑐 are. The entries 𝑏 and 𝑐 couple the eigenvalue 1 to the eigenvalue 2, which does no harm because the two eigenvalues are different; the entry 𝑎 couples the eigenvalue 1 to itself, as in a shear.

Exercise 7.4.38.standard Let 𝑛≥2 and 𝐴,𝐵∈𝐹𝑛×𝑛. Prove or disprove:

  1. (a)If 𝐴 is diagonalizable, then 𝐴 is invertible.
  2. (b)If 𝐴 is invertible, then 𝐴 is diagonalizable.
  3. (c)If 𝐴2 is diagonalizable, then 𝐴 is diagonalizable.
  4. (d)If 𝐴 and 𝐵 are diagonalizable, then so is 𝐴+𝐵.
  5. (e)If 𝐴 is diagonalizable and 𝐵 is similar to 𝐴, then 𝐵 is diagonalizable.

Solution.

Uses. Example 7.4.8 (Shears and nilpotent matrices), Proposition 7.4.7 (One eigenvalue), Proposition 4.5.24 (Similarity is an equivalence relation), Corollary 7.4.5 (Distinct eigenvalues give a diagonalization)

Techniques. Finding counterexamples and minimal examples (Technique 1.4.18)

Hint 1. The shear and the nilpotent matrix 𝑁=(0100) are the usual suspects.
Hint 2. For (d), look for two diagonalizable matrices whose sum is 𝑁. For (e), similarity is transitive.
Sketch. (a) False: 0. (b) False: the shear. (c) False: 𝑁2=0 is diagonal, 𝑁 is not diagonalizable. (d) False: (1100)+(−1000)=𝑁. (e) True: 𝐵∼𝐴∼𝐷.

The counterexamples below are 2×2. For 𝑛>2, enlarge each to a block diagonal matrix by adding the block 𝐼𝑛−2 in (b) and a zero block in (c) and (d). A block diagonal matrix with diagonalizable blocks is diagonalizable (combine eigenbases of the blocks), the enlarged shear still has 1 as its only eigenvalue without being 𝐼𝑛, and the enlarged nilpotent matrices are still nilpotent and nonzero.

(a) False: the zero matrix is diagonal, so diagonalizable, but not invertible.

(b) False: the shear (1101) is invertible but not diagonalizable (Example 7.4.8).

(c) False: 𝑁=(0100) has 𝑁2=0, which is diagonal, while 𝑁 is nonzero and nilpotent, so not diagonalizable (Example 7.4.8).

(d) False: (1100) has the distinct eigenvalues 1 and 0 (it is triangular), so it is diagonalizable (Corollary 7.4.5), and (−1000) is diagonal; their sum is 𝑁.

(e) True: if 𝐴 is similar to a diagonal 𝐷 and 𝐵 is similar to 𝐴, then 𝐵 is similar to 𝐷 by transitivity (Proposition 4.5.24).

Exercise 7.4.39.routine Find the error in the following argument.

“Every 𝐴∈ℂ𝑛×𝑛 is diagonalizable. By the Fundamental Theorem of Algebra, 𝑝𝐴=(𝑥−𝜆1)⋯(𝑥−𝜆𝑛) for some 𝜆1,…,𝜆𝑛∈ℂ. Each 𝜆𝑖 is a root of 𝑝𝐴, hence an eigenvalue; choose an eigenvector 𝐯𝑖 for 𝜆𝑖. Eigenvectors for different eigenvalues are independent, so 𝐯1,…,𝐯𝑛 is an independent list of length 𝑛, hence a basis of ℂ𝑛, and it is an eigenbasis.”

Solution.

Uses. Theorem 7.2.21 (Eigenvectors for distinct eigenvalues are independent), Example 7.4.8 (Shears and nilpotent matrices)

Hint 1. The 𝜆𝑖 in the factorization need not be distinct. What does Theorem 7.2.21 assume?
Hint 2. Run the argument on the shear (1101).
Sketch. The roots 𝜆𝑖 are listed with multiplicity, so they can repeat, and Theorem 7.2.21 only applies to distinct eigenvalues. For the shear, 𝜆1=𝜆2=1, and any two eigenvectors are multiples of 𝐞1. The argument proves only that matrices with 𝑛 distinct eigenvalues are diagonalizable.
The factorization lists the roots with multiplicity, so the 𝜆𝑖 need not be distinct, while Theorem 7.2.21 applies only to eigenvectors for pairwise distinct eigenvalues. When 𝜆𝑖=𝜆𝑗 for some 𝑖≠𝑗, nothing prevents 𝐯𝑗 from being a multiple of 𝐯𝑖, and then the list is dependent. For the shear 𝐽=(1101), 𝑝𝐽=(𝑥−1)2, so 𝜆1=𝜆2=1, and every eigenvector is a multiple of 𝐞1; no choice of 𝐯1,𝐯2 is independent, and 𝐽 is not diagonalizable (Example 7.4.8). The argument is correct when the 𝜆𝑖 are distinct, and then it is the proof of Corollary 7.4.5.
Exercise 7.4.40.standard Let 𝐴∈𝐹𝑛×𝑛 have rank 1. Show that 𝐴 is diagonalizable if and only if tr𝐴≠0, and that in this case 𝐴 is similar to diag(tr𝐴,0,…,0).

Solution.

Uses. Proposition 4.3.16 (Rank-one matrices), Example 4.3.18 (Powers of a rank-one matrix), Theorem 4.2.11 (Rank–nullity), Corollary 7.2.23 (Eigenspaces are independent), Theorem 7.4.16 (Criteria for diagonalizability), Example 7.4.8 (Shears and nilpotent matrices)

Techniques. Write a rank-one matrix as 𝐮𝐯𝖳 (Technique 4.3.17)

Hint 1. Write 𝐴=𝐮𝐯𝖳. Compute 𝐴𝐮, and recall 𝐯𝖳𝐮=tr𝐴.
Hint 2. gm(0)=dimker𝐴=𝑛−1. If tr𝐴≠0, you have one more eigenvector. If tr𝐴=0, compute 𝐴2.
Sketch. 𝐴=𝐮𝐯𝖳 with 𝐯𝖳𝐮=tr𝐴, so 𝐴𝐮=(tr𝐴)𝐮, and gm(0)=𝑛−1. If tr𝐴≠0, then gm(0)+gm(tr𝐴)≥𝑛: diagonalizable, with 𝐷=diag(tr𝐴,0,…,0). If tr𝐴=0, 𝐴2=(tr𝐴)𝐴=0 and 𝐴≠0: nilpotent, not diagonalizable.

By Proposition 4.3.16, 𝐴=𝐮𝐯𝖳 with 𝐮,𝐯 nonzero, and 𝐯𝖳𝐮=tr𝐴 (Example 4.3.18). Then

𝐴𝐮=𝐮(𝐯𝖳𝐮)=(tr𝐴)𝐮.

By rank–nullity, gm𝐴(0)=dimker𝐴=𝑛−1.

Suppose tr𝐴≠0. Then 𝐮≠𝟎 is an eigenvector with the nonzero eigenvalue tr𝐴, so gm(0)+gm(tr𝐴)≥(𝑛−1)+1=𝑛. By Corollary 7.2.23 the sum of all geometric multiplicities is at most 𝑛, so it equals 𝑛, and 𝐴 is diagonalizable by Theorem 7.4.16, (III) ⇒ (I). An eigenbasis is 𝐮 followed by a basis of ker𝐴, so 𝐴 is similar to diag(tr𝐴,0,…,0).

Suppose tr𝐴=0. Then 𝐴2=(tr𝐴)𝐴=0 (Example 4.3.18), and 𝐴≠0, so 𝐴 is a nonzero nilpotent matrix, which is not diagonalizable (Example 7.4.8).

Exercise 7.4.41.standard Let 𝑛≥2, and let 𝜏∈L(𝐹𝑛×𝑛) be the transpose map, 𝜏(𝑋)=𝑋𝖳. Show that 𝜏 is diagonalizable if and only if 1+1≠0 in 𝐹. In that case, find gm(1) and gm(−1).

Solution.

Uses. Example 7.2.17 (Reflections, and the field 𝔽2), Theorem 7.4.16 (Criteria for diagonalizability), Proposition 7.4.7 (One eigenvalue), Example 4.5.31 (The trace of transposition)

Techniques. Study a linear map on a space of matrices (Technique 4.2.29)

Hint 1. 𝜏2=id, so every eigenvalue satisfies 𝜆2=1 (Example 7.2.17). What are 𝐸1 and 𝐸−1?
Hint 2. 𝐸1 is the space of symmetric matrices and 𝐸−1 the space of skew-symmetric ones. Count their dimensions when 1+1≠0; when 1+1=0, there is only one eigenvalue.
Sketch. Eigenvalues satisfy 𝜆2=1. If 1+1≠0: 𝐸1 (symmetric, dimension 𝑛(𝑛+1)/2) and 𝐸−1 (skew-symmetric, with zero diagonal, dimension 𝑛(𝑛−1)/2) have dimensions adding up to 𝑛2. If 1+1=0: 1=−1 is the only eigenvalue and 𝜏≠id since 𝑛≥2.

Since 𝜏2=id, every eigenvalue of 𝜏 satisfies 𝜆2=1, so 𝜆=±1 (Example 7.2.17). The eigenspace 𝐸1 is the space of symmetric matrices, 𝑋𝖳=𝑋, and 𝐸−1 is the space of skew-symmetric matrices, 𝑋𝖳=−𝑋.

Suppose 1+1≠0. A symmetric matrix is determined by its entries 𝑥𝑖𝑗 with 𝑖≤𝑗, which are arbitrary, so dim𝐸1=𝑛(𝑛+1)/2; a basis is 𝑈𝑖𝑖 and 𝑈𝑖𝑗+𝑈𝑗𝑖 for 𝑖<𝑗. A skew-symmetric matrix has 𝑥𝑖𝑖=−𝑥𝑖𝑖, so (1+1)𝑥𝑖𝑖=0 and 𝑥𝑖𝑖=0; it is determined by its entries with 𝑖<𝑗, which are arbitrary, so dim𝐸−1=𝑛(𝑛−1)/2, with basis 𝑈𝑖𝑗−𝑈𝑗𝑖 for 𝑖<𝑗. Since 𝑛(𝑛+1)/2+𝑛(𝑛−1)/2=𝑛2=dim𝐹𝑛×𝑛, 𝜏 is diagonalizable by Theorem 7.4.16, with gm(1)=𝑛(𝑛+1)/2 and gm(−1)=𝑛(𝑛−1)/2. (This recovers tr𝜏=𝑛 from Example 4.5.31.)

Suppose 1+1=0. Then −1=1, so 1 is the only eigenvalue of 𝜏. Since 𝑛≥2, 𝜏(𝑈12)=𝑈21≠𝑈12, so 𝜏≠id, and 𝜏 is not diagonalizable by Proposition 7.4.7.

Exercise 7.4.42.standard
  1. (a)Show that if 𝑇∈L(𝑉) is diagonalizable, then 𝑉=ker𝑇⊕im𝑇.
  2. (b)Give an example of a non-diagonalizable 𝑇 with 𝑉=ker𝑇⊕im𝑇, and one with ker𝑇∩im𝑇≠{𝟎}.

Solution.

Uses. Theorem 7.4.16 (Criteria for diagonalizability), Theorem 4.4.7 (Criteria for direct sums), Theorem 4.2.11 (Rank–nullity), Theorem 4.4.19 (Dimension of a sum)

Hint 1. Split an eigenbasis into the vectors with eigenvalue 0 and the others. Where do the others lie?
Hint 2. If 𝑇𝐛=𝜆𝐛 with 𝜆≠0, then 𝐛=𝑇(𝜆−1𝐛)∈im𝑇.
Hint 3. For (b), invertible operators have ker𝑇={𝟎} and im𝑇=𝑉.
Sketch. In an eigenbasis, the vectors with eigenvalue 0 span ker𝑇=𝐸0 and the others lie in im𝑇, since 𝐛=𝑇(𝜆−1𝐛). So ker𝑇+im𝑇=𝑉, and rank–nullity forces the sum to be direct. (b) The invertible shear; and 𝑁=(0100), with ker𝑁=im𝑁=span(𝐞1).

(a) Let ℬ︀ be an eigenbasis, and split it into the vectors with eigenvalue 0 and the vectors 𝐛 with eigenvalue 𝜆≠0. The former lie in ker𝑇. For the latter, 𝐛=𝜆−1𝑇𝐛=𝑇(𝜆−1𝐛)∈im𝑇. So the subspace ker𝑇+im𝑇 contains the basis ℬ︀, and ker𝑇+im𝑇=𝑉. By Theorem 4.4.19 and rank–nullity,

dim(ker𝑇∩im𝑇)=dimker𝑇+dimim𝑇−dim𝑉=0,

so ker𝑇∩im𝑇={𝟎}, and the sum is direct (Theorem 4.4.7).

(b) The shear 𝐽=(1101) is invertible, so ker𝑇𝐽={𝟎} and im𝑇𝐽=𝐹2, and 𝐹2={𝟎}⊕𝐹2; but 𝐽 is not diagonalizable. For the second, 𝑁=(0100) has ker𝑁=span(𝐞1)=im𝑁.

Exercise 7.4.43.standard
  1. (a)Let 𝐴,𝐵∈𝐹𝑛×𝑛 be diagonalizable. Show that 𝐴 and 𝐵 are similar if and only if 𝑝𝐴=𝑝𝐵.
  2. (b)Show that (a) fails if only one of 𝐴,𝐵 is assumed diagonalizable.

Solution.

Uses. Theorem 7.4.16 (Criteria for diagonalizability), Example 7.4.14 (The multiplicities of a diagonal matrix), Proposition 7.1.34 (Roots of a split polynomial), Theorem 7.4.2 (Diagonalizable means an eigenbasis), Proposition 4.5.24 (Similarity is an equivalence relation), Lemma 7.4.13 (Geometric multiplicity through a rank), Example 7.4.8 (Shears and nilpotent matrices)

Techniques. Rule out similarity with an invariant (Technique 4.5.29)

Hint 1. Similar matrices have the same characteristic polynomial (Section 7.3). For the converse, what does 𝑝𝐴 say about the diagonal matrix 𝐷 in 𝐴=𝑆𝐷𝑆−1?
Hint 2. 𝑝𝐴=𝑝𝐷=∏(𝑥−𝑑𝑖), and the 𝑑𝑖 are determined by 𝑝𝐴 up to order (Proposition 7.1.34).
Sketch. If 𝐴=𝑆𝐷𝑆−1 and 𝐵=𝑅𝐸𝑅−1, then 𝑝𝐴=𝑝𝐷 and 𝑝𝐵=𝑝𝐸, and 𝑝𝐷 determines the diagonal of 𝐷 up to order. So 𝑝𝐴=𝑝𝐵 gives 𝐸=𝑃−1𝐷𝑃 for a permutation matrix 𝑃, and 𝐴∼𝐷∼𝐸∼𝐵. (b) 𝐼2 and the shear both have characteristic polynomial (𝑥−1)2.

(a) Similar matrices have the same characteristic polynomial (Section 7.3), which proves “only if”. Conversely, let 𝐴=𝑆𝐷𝑆−1 and 𝐵=𝑅𝐸𝑅−1 with 𝐷,𝐸 diagonal, and suppose 𝑝𝐴=𝑝𝐵. Since similar matrices have equal characteristic polynomials, 𝑝𝐷=𝑝𝐴=𝑝𝐵=𝑝𝐸. By Example 7.4.14, 𝑝𝐷=∏𝑖(𝑥−𝑑𝑖𝑖), and by Proposition 7.1.34 each 𝜆 occurs on the diagonal of 𝐷 as many times as its multiplicity as a root of 𝑝𝐷; likewise for 𝐸. So the diagonal of 𝐸 is a rearrangement of the diagonal of 𝐷. Rearrange the columns of 𝑆 in the same way, obtaining an invertible 𝑆′ whose 𝑗-th column is an eigenvector of 𝐴 with eigenvalue 𝑒𝑗𝑗. By Theorem 7.4.2, 𝐴=𝑆′𝐸𝑆′−1. Hence 𝐴∼𝐸∼𝐵, and 𝐴∼𝐵 by Proposition 4.5.24.

(b) 𝐼2 is diagonal and the shear 𝐽=(1101) is not diagonalizable, and 𝑝𝐼2=𝑝𝐽=(𝑥−1)2. They are not similar, since the only matrix similar to 𝐼2 is 𝐼2. (Alternatively, gm𝐼2(1)=2≠1=gm𝐽(1), by Lemma 7.4.13.)

Exercise 7.4.44.standard Let 𝑓𝑘 be the Fibonacci numbers, and work in 𝔽11.

  1. (a)Using Example 7.4.20, show that 𝑓𝑘=3(8𝑘−4𝑘) in 𝔽11 for all 𝑘≥0.
  2. (b)Deduce that 𝑓𝑘+10=𝑓𝑘 in 𝔽11 for all 𝑘, that is, the Fibonacci numbers modulo 11 repeat with period dividing 10.

Solution.

Uses. Example 7.4.20 (The Fibonacci matrix over several fields), Theorem 7.4.30 (Recurrences with distinct characteristic roots), Exercise 1.1.51

Techniques. Diagonalize, solve, conjugate back (Technique 7.4.22)

Hint 1. Over 𝔽11, 𝑥2−𝑥−1 has the distinct roots 8 and 4. By Theorem 7.4.30, 𝑓𝑘=𝑎8𝑘+𝑏4𝑘; find 𝑎,𝑏 from 𝑓0,𝑓1.
Hint 2. For (b), Fermat’s little theorem (Exercise 1.1.51) gives 810=410=1.
Sketch. The roots of 𝑥2−𝑥−1 in 𝔽11 are 8 and 4. With 𝑓𝑘=𝑎8𝑘+𝑏4𝑘, 𝑓0=0 gives 𝑏=−𝑎, and 𝑓1=1 gives 4𝑎=1, so 𝑎=3. Since 810=410=1, 𝑓𝑘+10=𝑓𝑘.

(a) In 𝔽11, 𝑥2−𝑥−1=(𝑥−4)(𝑥+3)=(𝑥−4)(𝑥−8) (Example 7.4.20), with distinct roots 8 and 4. By Theorem 7.4.30, applied with 𝑘=2 and 𝑐0=𝑐1=1 over 𝐹=𝔽11, there are 𝑎,𝑏∈𝔽11 with 𝑓𝑘=𝑎8𝑘+𝑏4𝑘 for all 𝑘. From 𝑓0=0, 𝑎+𝑏=0, and from 𝑓1=1, 8𝑎+4𝑏=4𝑎=1, so 𝑎=4−1=3 (as 4⋅3=12=1) and 𝑏=−3. So 𝑓𝑘=3(8𝑘−4𝑘). Check: 𝑓2=3(64−16)=3⋅48=144=1 and 𝑓3=3(512−64)=3⋅448=1344=2 in 𝔽11.

(b) By Fermat’s little theorem (Exercise 1.1.51), 810=410=1 in 𝔽11. So 𝑓𝑘+10=3(8𝑘810−4𝑘410)=3(8𝑘−4𝑘)=𝑓𝑘. Indeed, modulo 11 the sequence is 0,1,1,2,3,5,8,2,10,1,0,1,…, with period exactly 10.

Exercise 7.4.45.routine Find a closed formula for the sequence of real numbers with 𝑠0=3, 𝑠1=2, 𝑠2=6, and

𝑠𝑗+3=2𝑠𝑗+2+𝑠𝑗+1−2𝑠𝑗(𝑗≥0).

Solution.

Uses. Theorem 7.4.30 (Recurrences with distinct characteristic roots)

Hint 1. The characteristic polynomial is 𝑥3−2𝑥2−𝑥+2. Look for integer roots.
Hint 2. With roots 2,1,−1, write 𝑠𝑗=𝑎2𝑗+𝑏+𝑐(−1)𝑗 and solve for 𝑎,𝑏,𝑐 from 𝑠0,𝑠1,𝑠2.
Sketch. 𝑥3−2𝑥2−𝑥+2=(𝑥−2)(𝑥−1)(𝑥+1). Solving 𝑎+𝑏+𝑐=3, 2𝑎+𝑏−𝑐=2, 4𝑎+𝑏+𝑐=6 gives 𝑎=𝑏=𝑐=1, so 𝑠𝑗=2𝑗+1+(−1)𝑗.

The characteristic polynomial is 𝑞=𝑥3−2𝑥2−𝑥+2=𝑥2(𝑥−2)−(𝑥−2)=(𝑥−2)(𝑥−1)(𝑥+1), with distinct roots 2,1,−1. By Theorem 7.4.30, 𝑠𝑗=𝑎2𝑗+𝑏+𝑐(−1)𝑗 for unique 𝑎,𝑏,𝑐∈ℝ. The initial values give

𝑎+𝑏+𝑐=3,2𝑎+𝑏−𝑐=2,4𝑎+𝑏+𝑐=6.

Subtracting the first equation from the third gives 3𝑎=3, so 𝑎=1; then 𝑏+𝑐=2 and 𝑏−𝑐=0, so 𝑏=𝑐=1. Hence 𝑠𝑗=2𝑗+1+(−1)𝑗. Check: 𝑠3=8+1−1=8=2⋅6+2−2⋅3.

Exercise 7.4.46.standard Two quantities evolve by

𝑥𝑘+1=0.7𝑥𝑘+0.4𝑦𝑘,𝑦𝑘+1=0.2𝑥𝑘+0.9𝑦𝑘.
  1. (a)Find (𝑥𝑘,𝑦𝑘) in closed form in terms of (𝑥0,𝑦0).
  2. (b)Show that if 𝑥0,𝑦0≥0 are not both 0, then 𝑥𝑘/𝑦𝑘→1 and both quantities eventually grow by about 10% per step.
  3. (c)For which (𝑥0,𝑦0)∈ℝ2 does (𝑥𝑘,𝑦𝑘)→(0,0)?

Solution.

Uses. Proposition 7.4.33 (Long-run behavior), Theorem 7.4.2 (Diagonalizable means an eigenbasis)

Techniques. Diagonalize, solve, conjugate back (Technique 7.4.22)

Hint 1. The matrix is 𝐴=(0.70.40.20.9). Its eigenvalues are 1.1 and 0.5.
Hint 2. Expand (𝑥0,𝑦0) in the eigenbasis (1,1),(−2,1) and use (7.4.3).
Hint 3. The (1,1)-component is (𝑥0+2𝑦0)/3. When is it zero?
Sketch. 𝐴(1,1)=1.1(1,1) and 𝐴(−2,1)=0.5(−2,1). With 𝑐1=(𝑥0+2𝑦0)/3, 𝑐2=(𝑦0−𝑥0)/3, (𝑥𝑘,𝑦𝑘)=𝑐11.1𝑘(1,1)+𝑐20.5𝑘(−2,1). If 𝑥0,𝑦0≥0, not both 0, then 𝑐1>0 and the first term dominates. The state tends to 𝟎 iff 𝑐1=0, iff (𝑥0,𝑦0) is a multiple of (−2,1).

(a) The system is 𝐱𝑘+1=𝐴𝐱𝑘 with 𝐴=(0.70.40.20.9). Here det(𝐴−𝜆𝐼2)=𝜆2−1.6𝜆+0.55=(𝜆−1.1)(𝜆−0.5), and 𝐴(1,1)=(1.1,1.1), 𝐴(−2,1)=(−1,0.5), so (1,1) and (−2,1) form an eigenbasis with eigenvalues 1.1 and 0.5. Solving (𝑥0,𝑦0)=𝑐1(1,1)+𝑐2(−2,1) gives 𝑐2=(𝑦0−𝑥0)/3 and 𝑐1=(𝑥0+2𝑦0)/3. By (7.4.3),

(𝑥𝑘𝑦𝑘)=𝑥0+2𝑦03(1.1)𝑘(11)+𝑦0−𝑥03(0.5)𝑘(−21).

(b) If 𝑥0,𝑦0≥0 are not both 0, then 𝑐1=(𝑥0+2𝑦0)/3>0. By Proposition 7.4.33(b) with 𝜇=1.1, 1.1−𝑘(𝑥𝑘,𝑦𝑘)→𝑐1(1,1), so 1.1−𝑘𝑥𝑘→𝑐1 and 1.1−𝑘𝑦𝑘→𝑐1. Hence 𝑥𝑘/𝑦𝑘→𝑐1/𝑐1=1 (for large 𝑘, 𝑦𝑘>0), and 𝑥𝑘+1/𝑥𝑘=1.1⋅(1.1−𝑘−1𝑥𝑘+1)/(1.1−𝑘𝑥𝑘)→1.1, and likewise for 𝑦𝑘.

(c) If 𝑐1≠0, the 𝑥-coordinate is 𝑐1(1.1)𝑘−2𝑐2(0.5)𝑘, whose absolute value tends to ∞. If 𝑐1=0, then (𝑥𝑘,𝑦𝑘)=𝑐2(0.5)𝑘(−2,1)→(0,0). So the state tends to (0,0) exactly when 𝑥0+2𝑦0=0, that is, when (𝑥0,𝑦0) lies on the eigenline span((−2,1)) of the eigenvalue 0.5.

Exercise 7.4.47.standard Let 𝐴∈ℝ𝑛×𝑛 be diagonalizable over ℂ, and suppose all its complex eigenvalues are real. Show that 𝐴 is diagonalizable over ℝ.

Solution.

Uses. Lemma 7.4.13 (Geometric multiplicity through a rank), Theorem 3.3.17 (Uniqueness of the reduced row echelon form), Theorem 7.4.16 (Criteria for diagonalizability), Example 3.5.3 (A rank that depends on the field)

Hint 1. By Theorem 7.4.16(III), it suffices to compare gm(𝜆) over ℝ and over ℂ for each real eigenvalue 𝜆.
Hint 2. gm(𝜆)=𝑛−rank(𝐴−𝜆𝐼𝑛). Is the rank of a real matrix the same over ℝ and over ℂ?
Sketch. For real 𝜆, 𝐴−𝜆𝐼𝑛 is real, and Gauss–Jordan elimination over ℝ produces an RREF that is also the RREF over ℂ (uniqueness), so the rank, and hence gm(𝜆), is the same over both fields. The complex eigenvalues are exactly the real ones, so ∑gmℝ(𝜆)=∑gmℂ(𝜆)=𝑛.

Let 𝜆1,…,𝜆𝑚∈ℝ be the distinct eigenvalues of 𝐴 over ℂ; by hypothesis they are real, so they are also its eigenvalues over ℝ, and 𝐴 has no other real eigenvalues. For each 𝑖, 𝐵𝑖=𝐴−𝜆𝑖𝐼𝑛 is a real matrix. Row reducing 𝐵𝑖 over ℝ produces a matrix in reduced row echelon form by a sequence of row operations with real scalars. The same operations are complex row operations, so this real matrix is also the RREF of 𝐵𝑖 over ℂ (Theorem 3.3.17). The rank is the number of pivots, so 𝐵𝑖 has the same rank over ℝ and over ℂ, and by Lemma 7.4.13, dimℝkerℝ𝐵𝑖=dimℂkerℂ𝐵𝑖.

Since 𝐴 is diagonalizable over ℂ, ∑𝑖dimℂkerℂ𝐵𝑖=𝑛 by Theorem 7.4.16. Hence ∑𝑖dimℝkerℝ𝐵𝑖=𝑛, and 𝐴 is diagonalizable over ℝ by the same theorem. (Without the hypothesis on the eigenvalues, the rotation of Example 7.4.9 is a counterexample. Contrast Example 3.5.3, where the rank of one integer matrix depends on the field because the fields ℚ and 𝔽𝑝 do different arithmetic; ℝ⊆ℂ do the same arithmetic on real numbers.)

Exercise 7.4.48.standard Find all 𝐴∈ℝ2×2 with 𝐴2=𝐼2. Describe them in terms of their entries, and geometrically.

Solution.

Uses. Proposition 7.2.28 (Polynomials act on eigenvectors by their values), Exercise 7.2.55, Theorem 7.4.16 (Criteria for diagonalizability), Proposition 7.4.7 (One eigenvalue), Proposition 4.5.27 (Similarity invariants), Proposition 5.4.17 (Similar matrices have equal determinants)

Hint 1. Every eigenvalue satisfies 𝜆2=1. Show that 𝐴 is diagonalizable: 𝐯=(𝐯+𝐴𝐯)/2+(𝐯−𝐴𝐯)/2.
Hint 2. If 𝐴≠±𝐼2, then 𝐴 has both eigenvalues 1 and −1. What are its trace and determinant then?
Sketch. Every 𝐯 splits as (𝐯+𝐴𝐯)/2∈𝐸1 plus (𝐯−𝐴𝐯)/2∈𝐸−1, so 𝐴 is diagonalizable with eigenvalues in {1,−1}. So 𝐴=±𝐼2, or 𝐴=𝑆diag(1,−1)𝑆−1, which happens iff tr𝐴=0 and det𝐴=−1: 𝐴=(𝑎𝑏𝑐−𝑎) with 𝑎2+𝑏𝑐=1. Geometrically, the reflections across a line along another line.

If 𝐴𝐯=𝜆𝐯 with 𝐯≠𝟎, then 𝐯=𝐴2𝐯=𝜆2𝐯, so 𝜆=±1. For every 𝐯∈ℝ2,

𝐯=12(𝐯+𝐴𝐯)+12(𝐯−𝐴𝐯),

and 𝐴(𝐯+𝐴𝐯)=𝐴𝐯+𝐯, 𝐴(𝐯−𝐴𝐯)=−(𝐯−𝐴𝐯). So ℝ2=𝐸1+𝐸−1 (this is Exercise 7.2.55). Dropping a summand that is {𝟎} if 1 or −1 is not an eigenvalue, ℝ2 is the sum, hence the direct sum, of the eigenspaces of 𝐴, and 𝐴 is diagonalizable by Theorem 7.4.16, (II) ⇒ (I), with eigenvalues in {1,−1}.

If only one of ±1 is an eigenvalue, then 𝐴=±𝐼2 by Proposition 7.4.7. Otherwise 𝐴=𝑆diag(1,−1)𝑆−1 for an invertible 𝑆, so tr𝐴=trdiag(1,−1)=0 and det𝐴=−1 (Proposition 4.5.27(b), Proposition 5.4.17). Conversely, if 𝐴=(𝑎𝑏𝑐−𝑎) with −𝑎2−𝑏𝑐=det𝐴=−1, then

𝐴2=(𝑎2+𝑏𝑐𝑎𝑏−𝑎𝑏𝑎𝑐−𝑎𝑐𝑏𝑐+𝑎2)=𝐼2.

So the solutions are ±𝐼2 and the matrices (𝑎𝑏𝑐−𝑎) with 𝑎2+𝑏𝑐=1. Geometrically, 𝐴=𝑆diag(1,−1)𝑆−1 fixes the line 𝐸1 pointwise and negates the line 𝐸−1, with ℝ2=𝐸1⊕𝐸−1: it is the reflection across 𝐸1 along 𝐸−1, and every pair of distinct lines through 𝟎 occurs. The orthogonal reflections of Example 2.1.11 are the cases with 𝐸−1⟂𝐸1.

Exercise 7.4.49.hard Let 𝑞=(𝑥−𝑟1)(𝑥−𝑟2)⋯(𝑥−𝑟𝑘)∈𝐹[𝑥] with 𝑟1,…,𝑟𝑘∈𝐹 distinct, and suppose 𝑞(𝑇)=0.

  1. (a)Show that 𝑇 is diagonalizable. (Use the Lagrange polynomials 𝐿1,…,𝐿𝑘 for the nodes 𝑟1,…,𝑟𝑘: show that 𝐿1+⋯+𝐿𝑘=1 and that im𝐿𝑗(𝑇)⊆𝐸𝑟𝑗.)
  2. (b)Deduce that every 𝐴∈ℂ𝑛×𝑛 with 𝐴𝑚=𝐼𝑛 for some 𝑚≥1 is diagonalizable over ℂ.
  3. (c)Let 𝑝 be prime. Deduce that every 𝐴∈𝔽𝑝𝑛×𝑛 with 𝐴𝑝=𝐴 is diagonalizable over 𝔽𝑝.

Solution.

Uses. Definition 3.4.26 (Lagrange polynomials), Proposition 3.4.27 (Lagrange’s formula), Proposition 7.2.27 (Evaluating polynomials at an operator), Theorem 7.1.20 (Factor theorem), Proposition 7.1.27 (Multiplicities add up to at most the degree), Corollary 7.2.23 (Eigenspaces are independent), Theorem 7.4.16 (Criteria for diagonalizability), Exercise 1.1.51, Proposition 7.1.34 (Roots of a split polynomial)

Techniques. Test an identity on eigenvectors (Technique 7.4.25)

Hint 1. 𝐿1+⋯+𝐿𝑘 and 1 both have degree less than 𝑘 and take the value 1 at every node. Evaluate the identity at 𝑇: every 𝐯 is ∑𝑗𝐿𝑗(𝑇)𝐯.
Hint 2. 𝐿𝑗 vanishes at every 𝑟𝑖 with 𝑖≠𝑗, so (𝑥−𝑟𝑗)𝐿𝑗 is a multiple of 𝑞. Apply this at 𝑇 to show (𝑇−𝑟𝑗id)𝐿𝑗(𝑇)=0.
Hint 3. For (b), 𝑥𝑚−1 has 𝑚 distinct roots in ℂ, the 𝑚-th roots of unity. For (c), every 𝑎∈𝔽𝑝 is a root of 𝑥𝑝−𝑥.
Sketch. By uniqueness of interpolation, ∑𝑗𝐿𝑗=1, so 𝐯=∑𝑗𝐿𝑗(𝑇)𝐯 for every 𝐯. Since 𝐿𝑗(𝑟𝑖)=0 for 𝑖≠𝑗, 𝐿𝑗=𝑐𝑗∏𝑖≠𝑗(𝑥−𝑟𝑖), so (𝑥−𝑟𝑗)𝐿𝑗=𝑐𝑗𝑞 and (𝑇−𝑟𝑗id)𝐿𝑗(𝑇)=𝑐𝑗𝑞(𝑇)=0: 𝐿𝑗(𝑇)𝐯∈𝐸𝑟𝑗. So the eigenspaces span 𝑉. (b) 𝑥𝑚−1=∏(𝑥−𝑒2𝜋𝑖𝑙/𝑚). (c) 𝑥𝑝−𝑥=∏𝑎∈𝔽𝑝(𝑥−𝑎) by Fermat.

(a) By Proposition 3.4.27 with all values 𝑑𝑗=1, the unique polynomial of degree at most 𝑘−1 with value 1 at every node is 𝐿1+⋯+𝐿𝑘; the constant 1 is such a polynomial, so 𝐿1+⋯+𝐿𝑘=1. Evaluating at 𝑇 (Proposition 7.2.27(a)),

𝐯=𝐿1(𝑇)𝐯+⋯+𝐿𝑘(𝑇)𝐯for every𝐯∈𝑉.

Next, 𝐿𝑗 vanishes at the 𝑘−1 distinct nodes 𝑟𝑖, 𝑖≠𝑗, so by the factor theorem (applied 𝑘−1 times, as in Proposition 7.1.27) ∏𝑖≠𝑗(𝑥−𝑟𝑖) divides 𝐿𝑗; since deg𝐿𝑗≤𝑘−1, 𝐿𝑗=𝑐𝑗∏𝑖≠𝑗(𝑥−𝑟𝑖) for a scalar 𝑐𝑗. Hence (𝑥−𝑟𝑗)𝐿𝑗=𝑐𝑗𝑞, and evaluating at 𝑇,

(𝑇−𝑟𝑗id)𝐿𝑗(𝑇)=𝑐𝑗𝑞(𝑇)=0.

So 𝐿𝑗(𝑇)𝐯∈ker(𝑇−𝑟𝑗id)=𝐸𝑟𝑗 for every 𝐯. Thus 𝑉=𝐸𝑟1+⋯+𝐸𝑟𝑘. Dropping the 𝑟𝑗 that are not eigenvalues (their eigenspaces are {𝟎}), 𝑉 is the sum of the eigenspaces of 𝑇, so it is their direct sum (Corollary 7.2.23), and 𝑇 is diagonalizable by Theorem 7.4.16, (II) ⇒ (I).

(b) The 𝑚 numbers 𝜔𝑙=𝑒2𝜋𝑖𝑙/𝑚, 𝑙=0,…,𝑚−1, are distinct roots of 𝑥𝑚−1, so by Proposition 7.1.34 (or the factor theorem, repeatedly) 𝑥𝑚−1=∏𝑙(𝑥−𝜔𝑙), comparing degrees and leading coefficients. Since 𝐴𝑚−𝐼𝑛=0, (a) applies with 𝑞=𝑥𝑚−1.

(c) By Fermat’s little theorem (Exercise 1.1.51), 𝑎𝑝−1=1 for 𝑎≠0 in 𝔽𝑝, so 𝑎𝑝=𝑎 for every 𝑎∈𝔽𝑝, including 𝑎=0. So the 𝑝 distinct elements of 𝔽𝑝 are roots of 𝑥𝑝−𝑥, and as in (b), 𝑥𝑝−𝑥=∏𝑎∈𝔽𝑝(𝑥−𝑎). Since 𝐴𝑝−𝐴=0, (a) applies.

Exercise 7.4.50.hard Let 𝑇∈L(𝑉) be diagonalizable with distinct eigenvalues 𝜆1,…,𝜆𝑚, and let 𝑈⊆𝑉 be a 𝑇-invariant subspace. Show that

𝑈=(𝑈∩𝐸𝜆1)⊕⋯⊕(𝑈∩𝐸𝜆𝑚),

and deduce that the restriction of 𝑇 to 𝑈 is diagonalizable.

Solution.

Uses. Definition 3.4.26 (Lagrange polynomials), Proposition 7.2.28 (Polynomials act on eigenvectors by their values), Proposition 7.2.27 (Evaluating polynomials at an operator), Theorem 7.4.16 (Criteria for diagonalizability), Corollary 7.2.23 (Eigenspaces are independent), Theorem 4.4.13 (Bases and dimension of a sum), Proposition 7.4.26 (Polynomials that kill a diagonalizable operator), Exercise 7.4.49

Techniques. Look for invariant subspaces among kernels and images of commuting operators (Technique 7.2.31)

Hint 1. Write 𝐮∈𝑈 as 𝐮=𝐯1+⋯+𝐯𝑚 with 𝐯𝑖∈𝐸𝜆𝑖. You need each 𝐯𝑖 to lie in 𝑈.
Hint 2. Find a polynomial 𝐿𝑖 with 𝐿𝑖(𝑇)𝐮=𝐯𝑖. Polynomials in 𝑇 map 𝑈 into 𝑈.
Hint 3. Alternatively, apply Exercise 7.4.49 to the restriction: it is killed by ∏(𝑥−𝜆𝑖).
Sketch. Let 𝐿𝑖 be the Lagrange polynomials for the nodes 𝜆1,…,𝜆𝑚. For 𝐮=∑𝐯𝑗 with 𝐯𝑗∈𝐸𝜆𝑗, 𝐿𝑖(𝑇)𝐮=∑𝑗𝐿𝑖(𝜆𝑗)𝐯𝑗=𝐯𝑖, and 𝐿𝑖(𝑇)𝐮∈𝑈 since 𝑈 is invariant under polynomials in 𝑇. So 𝐯𝑖∈𝑈∩𝐸𝜆𝑖, and 𝑈 is the (direct) sum of these. Concatenating bases gives an eigenbasis of 𝑈.

First proof. Let 𝐿1,…,𝐿𝑚 be the Lagrange polynomials for the distinct nodes 𝜆1,…,𝜆𝑚 (Definition 3.4.26). Let 𝐮∈𝑈. Since 𝑉=𝐸𝜆1⊕⋯⊕𝐸𝜆𝑚 (Theorem 7.4.16), 𝐮=𝐯1+⋯+𝐯𝑚 with 𝐯𝑗∈𝐸𝜆𝑗, and by Proposition 7.2.28

𝐿𝑖(𝑇)𝐮=∑𝑗=1𝑚𝐿𝑖(𝜆𝑗)𝐯𝑗=𝐯𝑖.

Since 𝑈 is 𝑇-invariant, it is 𝐿𝑖(𝑇)-invariant (Proposition 7.2.27(d)), so 𝐯𝑖=𝐿𝑖(𝑇)𝐮∈𝑈∩𝐸𝜆𝑖. Hence 𝑈⊆∑𝑖(𝑈∩𝐸𝜆𝑖), and the reverse inclusion holds because each 𝑈∩𝐸𝜆𝑖 lies in 𝑈. The sum is direct because 𝐸𝜆1+⋯+𝐸𝜆𝑚 is (Corollary 7.2.23). Concatenating bases of the spaces 𝑈∩𝐸𝜆𝑖 gives a basis of 𝑈 (Theorem 4.4.13) consisting of eigenvectors of 𝑇, which are eigenvectors of the restriction; so the restriction is diagonalizable.

Second proof. By Proposition 7.4.26, 𝑞(𝑇)=0 for 𝑞=(𝑥−𝜆1)⋯(𝑥−𝜆𝑚). For 𝐮∈𝑈, 𝑞(𝑇↾𝑈)𝐮=𝑞(𝑇)𝐮=𝟎, where 𝑇↾𝑈 is regarded as an operator on 𝑈. So the restriction is killed by a product of distinct linear factors, and it is diagonalizable by Exercise 7.4.49(a). Its eigenspaces are ker(𝑇↾𝑈−𝜆𝑖id𝑈)=𝑈∩𝐸𝜆𝑖, and Theorem 7.4.16(II) for the restriction gives the decomposition of 𝑈.

Exercise 7.4.51.hardBeyond 21-242 Let 𝑆,𝑇∈L(𝑉) be diagonalizable with 𝑆𝑇=𝑇𝑆. Show that 𝑆 and 𝑇 are simultaneously diagonalizable: there is a basis of 𝑉 consisting of vectors that are eigenvectors of both 𝑆 and 𝑇. Show by example that the hypothesis 𝑆𝑇=𝑇𝑆 cannot be dropped.

Solution.

Uses. Proposition 7.2.30 (Kernels and images of commuting operators are invariant), Exercise 7.4.50, Theorem 7.4.16 (Criteria for diagonalizability), Theorem 4.4.13 (Bases and dimension of a sum)

Techniques. Look for invariant subspaces among kernels and images of commuting operators (Technique 7.2.31)

Hint 1. Each eigenspace 𝐸𝜆(𝑇) is 𝑆-invariant (Proposition 7.2.30).
Hint 2. The restriction of 𝑆 to 𝐸𝜆(𝑇) is diagonalizable by Exercise 7.4.50. Its eigenvectors are eigenvectors of 𝑇 too.
Hint 3. For the counterexample, a diagonal and a non-diagonal matrix with distinct eigenvalues.
Sketch. 𝑉=⨁𝐸𝜆𝑖(𝑇), and each 𝐸𝜆𝑖(𝑇) is 𝑆-invariant, so 𝑆 restricted to it is diagonalizable (Exercise 7.4.50). An eigenbasis of each restriction consists of common eigenvectors; concatenate. Without commuting: diag(1,2) and (1102) are both diagonalizable but have no common eigenbasis, since common eigenvectors would make them commute.

Let 𝜆1,…,𝜆𝑚 be the distinct eigenvalues of 𝑇, so that 𝑉=𝐸𝜆1(𝑇)⊕⋯⊕𝐸𝜆𝑚(𝑇) by Theorem 7.4.16. Since 𝑆𝑇=𝑇𝑆, each 𝐸𝜆𝑖(𝑇) is 𝑆-invariant (Proposition 7.2.30(b)). By Exercise 7.4.50, the restriction of 𝑆 to 𝐸𝜆𝑖(𝑇) is diagonalizable, so 𝐸𝜆𝑖(𝑇) has a basis ℬ︀𝑖 of eigenvectors of 𝑆. Each vector of ℬ︀𝑖 is nonzero and lies in 𝐸𝜆𝑖(𝑇), so it is also an eigenvector of 𝑇. The concatenation ℬ︀1,…,ℬ︀𝑚 is a basis of 𝑉 (Theorem 4.4.13), consisting of common eigenvectors.

The hypothesis is needed. If ℬ︀ is a common eigenbasis, then [𝑆]ℬ︀ℬ︀ and [𝑇]ℬ︀ℬ︀ are diagonal, hence commute, so 𝑆𝑇=𝑇𝑆. Thus 𝐷=diag(1,2) and 𝐴=(1102), which are both diagonalizable (distinct eigenvalues) but do not commute (𝐷𝐴=(1104), 𝐴𝐷=(1204)), have no common eigenbasis.

Exercise 7.4.52.hard Let 𝐴∈ℂ𝑛×𝑛 satisfy 𝐴3=𝐴. Show that rank𝐴=tr(𝐴2).

Solution.

Uses. Exercise 7.4.49, Proposition 7.4.23 (Powers and polynomials of a diagonalizable matrix), Proposition 4.5.27 (Similarity invariants), Proposition 7.2.28 (Polynomials act on eigenvectors by their values)

Techniques. Diagonalize, solve, conjugate back (Technique 7.4.22), Test an identity on eigenvectors (Technique 7.4.25)

Hint 1. 𝑥3−𝑥=𝑥(𝑥−1)(𝑥+1) has distinct roots, so 𝐴 is diagonalizable (Exercise 7.4.49).
Hint 2. Write 𝐴=𝑆𝐷𝑆−1. Both sides are similarity invariants, so compute them for 𝐷.
Hint 3. The diagonal entries of 𝐷 lie in {0,1,−1}, and 𝑑2=1 exactly when 𝑑≠0.
Sketch. 𝐴 is killed by 𝑥(𝑥−1)(𝑥+1), so 𝐴=𝑆𝐷𝑆−1 with 𝐷 diagonal, entries in {0,±1}. Then rank𝐴=rank𝐷 is the number of nonzero 𝑑𝑖, and tr𝐴2=tr𝐷2=∑𝑑𝑖2 is the same number.

Since 𝐴3−𝐴=0 and 𝑥3−𝑥=𝑥(𝑥−1)(𝑥+1) has the distinct roots 0,1,−1, 𝐴 is diagonalizable by Exercise 7.4.49(a): 𝐴=𝑆𝐷𝑆−1 with 𝐷=diag(𝑑1,…,𝑑𝑛). Each 𝑑𝑖 is an eigenvalue, hence a root of 𝑥3−𝑥 (Proposition 7.2.28), so 𝑑𝑖∈{0,1,−1}. Rank and trace are similarity invariants (Proposition 4.5.27), and 𝐴2=𝑆𝐷2𝑆−1 (Proposition 7.4.23), so

rank𝐴=rank𝐷=#{𝑖|𝑑𝑖≠0},tr(𝐴2)=tr(𝐷2)=∑𝑖=1𝑛𝑑𝑖2.

For 𝑑𝑖∈{0,1,−1}, 𝑑𝑖2=1 if 𝑑𝑖≠0 and 𝑑𝑖2=0 otherwise. So both sides count the nonzero 𝑑𝑖, and rank𝐴=tr(𝐴2).

Exercise 7.4.53.challenge Let 𝑉 be a finite-dimensional complex vector space and 𝑇∈L(𝑉). Show that 𝑇 is diagonalizable if and only if

ker(𝑇−𝜆id)∩im(𝑇−𝜆id)={𝟎}for every𝜆∈ℂ.

Solution.

Uses. Exercise 7.4.42, Proposition 4.4.16 (Complements exist), Definition 4.4.24 (Projection along a complement), Proposition 4.4.26 (Properties of projections), Theorem 7.3.27 (Existence of eigenvalues over ℂ), Corollary 7.2.23 (Eigenspaces are independent), Theorem 7.4.16 (Criteria for diagonalizability)

Techniques. Extremal choice and minimal counterexamples (Technique 1.5.9)

Hint 1. For “only if”, 𝑇−𝜆id is diagonalizable; use Exercise 7.4.42.
Hint 2. For “if”, let 𝑈 be the sum of all eigenspaces and suppose 𝑈≠𝑉. Take a complement 𝑊 and the projection 𝑃 onto 𝑊 along 𝑈. The operator 𝑃𝑇↾𝑊 on 𝑊 has an eigenvector 𝐰 over ℂ: 𝑇𝐰−𝜆𝐰∈𝑈.
Hint 3. Write (𝑇−𝜆id)𝐰=∑𝜇𝐮𝜇 with 𝐮𝜇∈𝐸𝜇. For 𝜇≠𝜆, 𝐮𝜇=(𝑇−𝜆id)(𝐮𝜇/(𝜇−𝜆)). Subtract these from 𝐰 and use the hypothesis.
Sketch. “Only if”: 𝑇−𝜆id is diagonalizable, so its kernel and image meet in 𝟎. “If”: let 𝑈=⨁𝐸𝜇 and suppose 𝑈≠𝑉; choose 𝑉=𝑈⊕𝑊. The operator 𝐰↦𝑃𝑇𝐰 on 𝑊 has an eigenvector 𝐰≠𝟎, so (𝑇−𝜆id)𝐰=∑𝜇𝐮𝜇∈𝑈. Put 𝐰′=𝐰−∑𝜇≠𝜆𝐮𝜇/(𝜇−𝜆); then (𝑇−𝜆id)𝐰′=𝐮𝜆∈ker∩im, so 𝐮𝜆=𝟎 and 𝐰′∈𝐸𝜆⊆𝑈. Then 𝐰∈𝑈∩𝑊={𝟎}, a contradiction.

Only if. If ℬ︀ is an eigenbasis for 𝑇, it is also one for 𝑇−𝜆id, since 𝑇𝐛=𝜇𝐛 gives (𝑇−𝜆id)𝐛=(𝜇−𝜆)𝐛. So 𝑇−𝜆id is diagonalizable, and Exercise 7.4.42(a) gives ker(𝑇−𝜆id)∩im(𝑇−𝜆id)={𝟎}.

If. Let 𝑈 be the sum of the eigenspaces of 𝑇; it is the direct sum of the eigenspaces of the finitely many eigenvalues (Corollary 7.2.23). Suppose, for contradiction, that 𝑈≠𝑉. Choose a complement 𝑊, so 𝑉=𝑈⊕𝑊 and 𝑊≠{𝟎} (Proposition 4.4.16), and let 𝑃=Pr𝑊,𝑈 be the projection onto 𝑊 along 𝑈 (Definition 4.4.24). The map 𝑅:𝑊→𝑊, 𝑅𝐰=𝑃𝑇𝐰, is an operator on the nonzero complex space 𝑊, so it has an eigenvector 𝐰≠𝟎 (Theorem 7.3.27): 𝑃𝑇𝐰=𝜆𝐰 for some 𝜆∈ℂ. Since 𝑃𝐰=𝐰, this says 𝑃(𝑇𝐰−𝜆𝐰)=𝟎, and since ker𝑃=𝑈 and 𝑃 fixes 𝑊 (Proposition 4.4.26(b)), (𝑇−𝜆id)𝐰∈𝑈. Write

(𝑇−𝜆id)𝐰=∑𝜇𝐮𝜇,𝐮𝜇∈𝐸𝜇,

the sum over the eigenvalues 𝜇 of 𝑇. For 𝜇≠𝜆, (𝑇−𝜆id)𝐮𝜇=(𝜇−𝜆)𝐮𝜇, so 𝐮𝜇=(𝑇−𝜆id)(𝐮𝜇/(𝜇−𝜆)). Let

𝐰′=𝐰−∑𝜇≠𝜆𝐮𝜇𝜇−𝜆.

Then (𝑇−𝜆id)𝐰′=∑𝜇𝐮𝜇−∑𝜇≠𝜆𝐮𝜇=𝐮𝜆, where 𝐮𝜆=𝟎 if 𝜆 is not an eigenvalue. In any case (𝑇−𝜆id)𝐰′ lies in ker(𝑇−𝜆id) (it is 𝟎 or in 𝐸𝜆) and in im(𝑇−𝜆id), so it is 𝟎 by hypothesis. Hence 𝐰′∈ker(𝑇−𝜆id)⊆𝑈, and therefore 𝐰=𝐰′+∑𝜇≠𝜆𝐮𝜇/(𝜇−𝜆)∈𝑈. But 𝐰∈𝑊 and 𝑈∩𝑊={𝟎}, so 𝐰=𝟎, a contradiction. Hence 𝑈=𝑉, and 𝑇 is diagonalizable by Theorem 7.4.16, (II) ⇒ (I).

The hypothesis is needed only for the one 𝜆 that the argument produces, but that 𝜆 is not known in advance. Over ℝ the statement fails: the rotation by a right angle has ker(𝑇−𝜆id)={𝟎} for every real 𝜆 and is not diagonalizable over ℝ.

More practice

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