Books
Theme
Text size
Width

[ and ] go to the previous and next page.

In Section 7.2 we found eigenvalues by solving 𝑇𝐯=𝜆𝐯 by hand, with a separate argument for each operator. Theorem 7.2.9 points to a systematic method: 𝜆 is an eigenvalue of 𝐴∈𝐹𝑛×𝑛 exactly when det(𝐴−𝜆𝐼𝑛)=0. If we leave 𝜆 unknown and call it 𝑥, the Leibniz formula turns det(𝑥𝐼𝑛−𝐴) into a polynomial of degree 𝑛 in 𝑥, and the eigenvalues of 𝐴 are its roots. This section makes that precise, reads the trace and the determinant of 𝐴 off the coefficients, and shows that every monic polynomial occurs. It also settles the existence question raised by the rotation in Example 7.2.18. Over ℂ every operator on a nonzero finite-dimensional space has an eigenvalue, and we prove this twice, once with the characteristic polynomial and once without determinants.

7.3.1Determinants with polynomial entries

The matrix 𝑥𝐼𝑛−𝐴 has entries in 𝐹[𝑥], and 𝐹[𝑥] is not a field: 𝑥 has no multiplicative inverse in 𝐹[𝑥] (Corollary 7.1.7). Chapter 5 defined the determinant only for matrices with entries in a field. So before we can write det(𝑥𝐼𝑛−𝐴), we have to say what it means and which rules it obeys.

Intuition (What the Leibniz formula needs). The Leibniz formula (Theorem 5.3.3) builds det𝐴 out of the entries of 𝐴 with additions, subtractions, and multiplications, and never divides. So it makes sense for entries from any commutative ring (Definition 1.1.20), and for 𝐹[𝑥] in particular (Proposition 7.1.3). The real question is which theorems of Chapter 5 survive. Most facts there were proved twice: once by row reduction, which divides by pivots, and once from the Leibniz formula or from the uniqueness of multilinear alternating functions, which do not divide. The second kind of proof works over a commutative ring word for word. The one statement that changes is the invertibility criterion, because in a ring a nonzero element need not be invertible.

Definition 7.3.1 (Determinant over a commutative ring). Let 𝑅 be a commutative ring and 𝑛≥1. We write 𝑅𝑛×𝑛 for the set of 𝑛×𝑛 arrays 𝑀=(𝑚𝑖𝑗) with entries in 𝑅, and add, multiply, and transpose such arrays by the same formulas as matrices over a field. The determinant of 𝑀∈𝑅𝑛×𝑛 is

det𝑀=∑𝜎∈𝑆𝑛sgn(𝜎)𝑚𝜎(1)1𝑚𝜎(2)2⋯𝑚𝜎(𝑛)𝑛∈𝑅,

where sgn(𝜎)∈{1,−1} is read in 𝑅. When 𝑅=𝐹 is a field, this is the determinant of Definition 5.3.6, by Theorem 5.3.3.

For example, with 𝑅=𝐹[𝑥],

det(𝑥−𝑎−𝑏−𝑐𝑥−𝑑)=(𝑥−𝑎)(𝑥−𝑑)−(−𝑏)(−𝑐)=𝑥2−(𝑎+𝑑)𝑥+(𝑎𝑑−𝑏𝑐).

Proposition 7.3.2 (Determinants over a commutative ring). Let 𝑅 and 𝑅′ be commutative rings, 𝑛≥1, and 𝑀,𝑁∈𝑅𝑛×𝑛.

  1. (a)Multiplication in 𝑅𝑛×𝑛 is associative and distributes over addition, (𝑟𝑀)𝑁=𝑟(𝑀𝑁)=𝑀(𝑟𝑁) for 𝑟∈𝑅, and the array 𝐼𝑛 with 1 on the diagonal and 0 elsewhere is a two-sided identity. In particular the scalar arrays 𝑟𝐼𝑛 commute with every array.
  2. (b)det is linear in each column (with coefficients in 𝑅), it vanishes when two columns are equal, and det𝐼𝑛=1. Every function 𝐷:𝑅𝑛×𝑛→𝑅 with these first two properties satisfies 𝐷(𝑀)=𝐷(𝐼𝑛)det𝑀.
  3. (c)det(𝑀𝑁)=det𝑀det𝑁.
  4. (d)det𝑀𝖳=det𝑀. So (b) holds for rows as well as columns.
  5. (e)If 𝑀=(𝑃𝑄0𝑆) with square blocks 𝑃 and 𝑆, then det𝑀=det𝑃det𝑆. In particular, if 𝑀 is upper or lower triangular, then det𝑀=𝑚11𝑚22⋯𝑚𝑛𝑛.
  6. (f)For 𝑛≥2, cofactor expansion along any row or column computes det𝑀, and 𝑀adj(𝑀)=adj(𝑀)𝑀=(det𝑀)𝐼𝑛.
  7. (g)Let 𝜑:𝑅→𝑅′ satisfy 𝜑(𝑎+𝑏)=𝜑(𝑎)+𝜑(𝑏), 𝜑(𝑎𝑏)=𝜑(𝑎)𝜑(𝑏), and 𝜑(1)=1 for all 𝑎,𝑏∈𝑅, and write 𝜑(𝑀)=(𝜑(𝑚𝑖𝑗)). Then det𝜑(𝑀)=𝜑(det𝑀).

Proof.standard

Uses. Definition 7.3.1 (Determinant over a commutative ring), Proposition 2.4.14 (Algebra of matrices), Theorem 5.3.4 (Existence and uniqueness of the determinant), Theorem 5.3.3 (The Leibniz formula), Theorem 5.4.10 (Multiplicativity), Theorem 5.4.14 (The transpose), Theorem 5.4.26 (Block-triangular matrices), Proposition 5.3.15 (Triangular matrices), Theorem 5.5.3 (Cofactor expansion), Theorem 5.5.9 (The adjugate and the inverse), Definition 1.1.20 (Commutative ring)

Techniques. Cancel terms in pairs (Technique 5.3.5), Recognize the determinant by its axioms (Technique 5.4.11)

Hint 1. For (a)–(f), go back to Chapter 5 and, for each fact, find a proof that never divides. Which proofs used row reduction, and which used only the Leibniz formula or the uniqueness of multilinear alternating functions?
Hint 2. Alternation is the delicate point: the shortcut “𝐿(𝐴)=−𝐿(𝐴)” divides by 1+1. The proof of Theorem 5.3.4 cancelled terms in pairs instead.
Hint 3. For (g), apply 𝜑 to the Leibniz sum term by term. You need 𝜑(0)=0 and 𝜑(−𝑎)=−𝜑(𝑎); derive them from additivity.
Sketch. Each of (a)–(f) has a proof in Chapter 5 that uses only the commutative ring axioms: entrywise computation for (a); the Leibniz formula, pairing 𝜎 with 𝜎𝜏, for (b); uniqueness for (c), since 𝐛↦det(𝑀𝐛) is multilinear and alternating; reindexing by 𝜎−1 for (d); the Leibniz formula for (e) and (f). For (g), 𝜑(0)=0 and 𝜑(−𝑎)=−𝜑(𝑎) follow from additivity, so 𝜑 maps each Leibniz term of 𝑀 to the corresponding term of 𝜑(𝑀).

We go through Chapter 5 and point, for each statement, to an argument that uses only the commutative ring axioms, that is, never divides and never uses 1+1≠0.

(a) The entrywise proof of Proposition 2.4.14 compares the entries of both sides, each of which is a sum of products of entries, using associativity, commutativity, and distributivity in the scalars. These hold in 𝑅.

(b) The proof of Theorem 5.3.4 (parts 2 and 3) shows that the Leibniz formula is normalized, linear in each column, and alternating. Normalization and multilinearity are checked term by term, since each term contains exactly one entry from each column. Alternation is proved by pairing each even 𝜎 with 𝜎𝜏, where 𝜏 swaps the two equal columns, and the two terms cancel (Technique 5.3.5). None of this divides. The last statement is Theorem 5.3.3. Its proof expands 𝐷(𝑀) by multilinearity into a sum over functions [𝑛]→[𝑛], discards the non-injective ones because 𝐷 vanishes on arrays with two equal columns, and evaluates 𝐷 on permuted identity columns with Theorem 5.2.31; these steps use multilinearity, alternation, and the ring axioms only.

(c) We use the second proof of Theorem 5.4.10 (Technique 5.4.11). Fix 𝑀 and let 𝐷(𝑁)=det(𝑀𝑁). Column 𝑗 of 𝑀𝑁 is 𝑀 times column 𝑗 of 𝑁, so 𝐷 is linear in each column of 𝑁 by (a) and (b), and two equal columns of 𝑁 give two equal columns of 𝑀𝑁, where det vanishes. By (b), 𝐷(𝑁)=𝐷(𝐼𝑛)det𝑁=det𝑀det𝑁.

(d) The second proof of Theorem 5.4.14 reindexes the Leibniz sum by 𝜎↦𝜎−1 and uses sgn(𝜎−1)=sgn(𝜎); this is ring arithmetic.

(e) In the Leibniz formula for (𝑃𝑄0𝑆), with 𝑃 of size 𝑘×𝑘, a term can be nonzero only if 𝜎 maps {1,…,𝑘} into itself, and then the term is a product of a Leibniz term of 𝑃 and one of 𝑆. This is the proof of Theorem 5.4.26 via the Leibniz formula, and the triangular case is the Leibniz proof of Proposition 5.3.15. Lower block-triangular arrays follow by (d).

(f) The second proof of Theorem 5.5.3 derives the expansion from the Leibniz formula by regrouping its terms, using the signs of cycles and the multiplicativity of sgn, and gets the row version from (d). The identity 𝑀adj(𝑀)=(det𝑀)𝐼𝑛 is proved in Theorem 5.5.9 by expanding along a row and noting that expansion along the “wrong” row computes the determinant of an array with two equal rows, which is 0 by (b) and (d). Neither step divides.

(g) First, 𝜑(0)=𝜑(0+0)=𝜑(0)+𝜑(0), so 𝜑(0)=0; and 𝜑(𝑎)+𝜑(−𝑎)=𝜑(0)=0, so 𝜑(−𝑎)=−𝜑(𝑎). Hence 𝜑(sgn(𝜎))=sgn(𝜎), and since 𝜑 respects sums and products of any finite length,

𝜑(det𝑀)=∑𝜎∈𝑆𝑛𝜑(sgn(𝜎))𝜑(𝑚𝜎(1)1)⋯𝜑(𝑚𝜎(𝑛)𝑛)=det𝜑(𝑀).

□

The most important instance of (g) is evaluation. For 𝑀∈𝐹[𝑥]𝑛×𝑛 and 𝜆∈𝐹, let 𝑀(𝜆)∈𝐹𝑛×𝑛 be the matrix obtained by evaluating every entry at 𝜆. Evaluation respects sums, products, and 1 (Proposition 7.1.10), so

(det𝑀)(𝜆)=det𝑀(𝜆).
(7.3.1)

In words, it does not matter whether we first compute the determinant as a polynomial and then plug in 𝜆, or first plug in 𝜆 and then compute a determinant of scalars. The same holds for substituting a polynomial 𝑟∈𝐹[𝑥] for 𝑥 in every entry (also Proposition 7.1.10), and for complex conjugation of the coefficients of polynomials in ℂ[𝑥] (Lemma 7.1.37).

Remark 7.3.3 (What does not carry over). Over a field, 𝑀 is invertible if and only if det𝑀≠0 (Theorem 5.4.7). Over 𝐹[𝑥] this fails: det(𝑥𝐼𝑛)=𝑥𝑛≠𝟎, but 𝑥𝐼𝑛 has no inverse in 𝐹[𝑥]𝑛×𝑛: if (𝑥𝐼𝑛)𝑁=𝐼𝑛, then 𝑥𝑛11=1, which is impossible by Corollary 7.1.7(c). The proof of Theorem 5.4.7 used row reduction, which divides by pivots. What survives is (f): if det𝑀 has an inverse in 𝑅, then (det𝑀)−1adj(𝑀) is an inverse of 𝑀. The converse also holds, by (c): 𝑀𝑁=𝐼𝑛 gives det𝑀det𝑁=1. So 𝑀 is invertible in 𝑅𝑛×𝑛 exactly when det𝑀 is invertible in 𝑅, which for 𝑅=𝐹[𝑥] means that det𝑀 is a nonzero constant. We will not need this, but it explains why nothing like “det(𝑥𝐼𝑛−𝐴)≠𝟎, so 𝑥𝐼𝑛−𝐴 is invertible” can be used.
Aside (Rational functions). There is a second way to make sense of det(𝑥𝐼𝑛−𝐴). Fractions 𝑝/𝑞 of polynomials with 𝑞≠𝟎, with the usual rules for fractions, form a field 𝐹(𝑥) that contains 𝐹[𝑥], in the same way that ℚ contains ℤ. All of Chapter 5 applies to matrices over 𝐹(𝑥), and the Leibniz formula shows that a matrix with polynomial entries has a polynomial determinant. This route gives everything in Proposition 7.3.2 for 𝑅=𝐹[𝑥] at once, at the cost of constructing 𝐹(𝑥). We stay with commutative rings, which we need anyway for (g).

7.3.2The characteristic polynomial

Definition 7.3.4 (Characteristic polynomial). Let 𝑛≥1 and 𝐴∈𝐹𝑛×𝑛. The characteristic polynomial of 𝐴 is

𝑝𝐴=det(𝑥𝐼𝑛−𝐴)∈𝐹[𝑥],

the determinant (Definition 7.3.1) of the matrix 𝑥𝐼𝑛−𝐴∈𝐹[𝑥]𝑛×𝑛, whose diagonal entries are 𝑥−𝑎𝑖𝑖 and whose off-diagonal entries are the constants −𝑎𝑖𝑗. When the polynomial is evaluated, we write 𝑝𝐴(𝜆) for its value at 𝜆∈𝐹.

Some texts, including the lecture notes, use det(𝐴−𝑥𝐼𝑛) instead. Multiplying each of the 𝑛 columns by −1 shows det(𝐴−𝑥𝐼𝑛)=(−1)𝑛𝑝𝐴 (Proposition 7.3.2(b)), so the two conventions differ by a sign and have the same roots. With our sign, 𝑝𝐴 is monic (Proposition 7.3.18).

Example 7.3.5 (Small cases).
  1. (a)For 𝐴=(𝑎𝑏𝑐𝑑), the computation after Definition 7.3.1 gives

    𝑝𝐴=𝑥2−(𝑎+𝑑)𝑥+(𝑎𝑑−𝑏𝑐)=𝑥2−(tr𝐴)𝑥+det𝐴.
  2. (b)For 𝐴=(2112) of Example 7.2.8, 𝑝𝐴=𝑥2−4𝑥+3=(𝑥−1)(𝑥−3), and its roots 1 and 3 are the eigenvalues found there.
  3. (c)For the rotation 𝐴=(0−110) of Example 7.2.18, 𝑝𝐴=𝑥2+1. It has no root in ℝ or 𝔽3, the roots ±𝑖 in ℂ, and the roots 2,3 in 𝔽5, matching the eigenvalues found there field by field.
  4. (d)If 𝐴 is upper triangular, so is 𝑥𝐼𝑛−𝐴, with diagonal entries 𝑥−𝑎𝑖𝑖. By Proposition 7.3.2(e),

    𝑝𝐴=(𝑥−𝑎11)(𝑥−𝑎22)⋯(𝑥−𝑎𝑛𝑛).

    For instance, 𝑝𝐴=(𝑥−2)2 for 𝐴=(2102), the matrix of Example 7.2.39 in its chain basis.

Theorem 7.3.6 (Eigenvalues are the roots of the characteristic polynomial). Let 𝐴∈𝐹𝑛×𝑛 and 𝜆∈𝐹. Then 𝜆 is an eigenvalue of 𝐴 if and only if 𝑝𝐴(𝜆)=0.

Proof.routine

Uses. Theorem 7.2.9 (Eigenvalues and invertibility), Proposition 7.3.2 (Determinants over a commutative ring), Proposition 7.1.10 (Evaluation respects sums and products), Definition 7.3.4 (Characteristic polynomial)

Hint 1. Theorem 7.2.9 already says that 𝜆 is an eigenvalue if and only if det(𝐴−𝜆𝐼𝑛)=0. How is that number related to the polynomial 𝑝𝐴?
Hint 2. Evaluate 𝑥𝐼𝑛−𝐴 entry by entry at 𝜆 and use (7.3.1).
Sketch. By (7.3.1), 𝑝𝐴(𝜆)=det(𝜆𝐼𝑛−𝐴)=(−1)𝑛det(𝐴−𝜆𝐼𝑛), which vanishes exactly when 𝜆 is an eigenvalue, by Theorem 7.2.9.
Evaluating every entry of 𝑥𝐼𝑛−𝐴 at 𝜆 gives 𝜆𝐼𝑛−𝐴, so by (7.3.1), 𝑝𝐴(𝜆)=det(𝜆𝐼𝑛−𝐴). Multiplying each column by −1 gives det(𝜆𝐼𝑛−𝐴)=(−1)𝑛det(𝐴−𝜆𝐼𝑛). So 𝑝𝐴(𝜆)=0 if and only if det(𝐴−𝜆𝐼𝑛)=0, which by Theorem 7.2.9 holds if and only if 𝜆 is an eigenvalue of 𝐴. □

The theorem turns the search for eigenvalues into a search for roots of one polynomial of degree 𝑛. It reproves Corollary 7.2.22 for matrices: a nonzero polynomial of degree 𝑛 has at most 𝑛 roots (Proposition 7.1.22).

Algorithm 7.3.7 (Eigenvalues by the characteristic polynomial)

Input
A matrix 𝐴∈𝐹𝑛×𝑛.
Output
The eigenvalues of 𝐴 in 𝐹, and a basis of each eigenspace.

1.

Compute 𝑝𝐴=det(𝑥𝐼𝑛−𝐴), by cofactor expansion along a sparse row or column, or by the Leibniz formula for small 𝑛.

2.

Find the roots 𝜆1,…,𝜆𝑘 of 𝑝𝐴 in 𝐹. These are the eigenvalues.

3.

For each 𝜆𝑗, row reduce 𝐴−𝜆𝑗𝐼𝑛 and read off a basis of its null space 𝐸𝜆𝑗 from the free columns.

The algorithm is correct by Theorem 7.3.6 and Definition 7.2.13. In step 3 the matrix 𝐴−𝜆𝑗𝐼𝑛 must have rank less than 𝑛; if row reduction finds 𝑛 pivots, then 𝜆𝑗 was not a root, and the arithmetic in step 1 or step 2 should be checked. The hard step is step 2: for 𝑛≥5 there is no formula in radicals for the roots in general, and in practice eigenvalues are computed by other methods. For hand computations in this book, the roots will be findable.

Example 7.3.8 (A 3×3 computation). Let 𝐴=(401−210−201)∈ℝ3×3. Then

𝑥𝐼3−𝐴=(𝑥−40−12𝑥−1020𝑥−1).

The middle column has a single nonzero entry, so expanding along it (Proposition 7.3.2(f)),

𝑝𝐴=(𝑥−1)det(𝑥−4−12𝑥−1)=(𝑥−1)((𝑥−4)(𝑥−1)+2)=(𝑥−1)(𝑥2−5𝑥+6)=(𝑥−1)(𝑥−2)(𝑥−3).

So the eigenvalues are 1, 2, 3. For 𝜆=2, row reduction gives

𝐴−2𝐼3=(201−2−10−20−1)⇝(101/201−1000),

so 𝐸2=span((1,−2,−2)); check: 𝐴(1,−2,−2)=(2,−4,−4). In the same way, 𝐸1=span((0,1,0)) and 𝐸3=span((1,−1,−1)). Realize that tr𝐴=6=1+2+3 and det𝐴=6=1⋅2⋅3; Remark 7.3.19 explains this.

(201−2−10−20−1)⇝(101/201−1000)
Row reduction of 𝐴−2𝐼3 from Example 7.3.8. The third column is free, so 𝐸2 is a line. Use “Edit” to try 𝐴−𝜆𝐼3 for other values of 𝜆: you get three pivots, and only the zero vector in the null space, unless 𝜆 is 1, 2, or 3.
Exercise 7.3.9.routine Let 𝐴=(1110). Find 𝑝𝐴, and find the eigenvalues of 𝐴 together with an eigenvector for each, over ℝ, over 𝔽2, over 𝔽5, and over 𝔽11.

Solution.

Uses. Example 7.3.5 (Small cases), Theorem 7.3.6 (Eigenvalues are the roots of the characteristic polynomial)

Hint 1. 𝑝𝐴=𝑥2−(tr𝐴)𝑥+det𝐴. Then find the roots field by field.
Hint 2. For a root 𝜆, the first row of 𝐴−𝜆𝐼2 is (1−𝜆,1), so (1,𝜆−1) is in the kernel. Since 𝜆2=𝜆+1, this is (1,1/𝜆), or more simply the multiple (𝜆,1).
Sketch. 𝑝𝐴=𝑥2−𝑥−1. Over ℝ: (1±5)/2. Over 𝔽2: 𝑥2+𝑥+1 has no root. Over 𝔽5: (𝑥−3)2, one eigenvalue 3. Over 𝔽11: roots 4 and 8. Eigenvectors (𝜆,1) in each case.

By Example 7.3.5(a), 𝑝𝐴=𝑥2−𝑥−1. If 𝜆 is a root, then (𝐴−𝜆𝐼2)(𝜆,1)=(𝜆−𝜆2+1,𝜆−𝜆)=(0,0), because 𝜆2=𝜆+1. So (𝜆,1) is an eigenvector for every root 𝜆, and by Theorem 7.3.6 the eigenvalues are exactly the roots.

Over ℝ, the roots are (1±5)/2, the golden ratio and its conjugate.

Over 𝔽2, 𝑝𝐴=𝑥2+𝑥+1 takes the value 1 at both 0 and 1, so 𝐴 has no eigenvalues.

Over 𝔽5, 𝑝𝐴=𝑥2−𝑥−1=𝑥2+4𝑥+4=(𝑥+2)2=(𝑥−3)2, so 3 is the only eigenvalue, with eigenvector (3,1). Indeed 𝐴(3,1)=(4,3)=3(3,1), as 9=4 in 𝔽5.

Over 𝔽11, 𝑝𝐴(4)=16−4−1=11=0 and 𝑝𝐴(8)=64−8−1=55=0, so the eigenvalues are 4 and 8 (there are at most two), with eigenvectors (4,1) and (8,1). These are the two “square roots of 5 formulas” (1±4)/2 read in 𝔽11, since 42=5 there.

Remark 7.3.10 (A polynomial, not a function). Over a finite field, the function 𝜆↦det(𝜆𝐼𝑛−𝐴) on 𝐹 carries less information than the polynomial 𝑝𝐴. Over 𝔽2, take 𝐴=03×3 and 𝐵=(000001011). Then 𝑝𝐴=𝑥3, and by Proposition 7.3.2(e) and Example 7.3.5(a), 𝑝𝐵=𝑥(𝑥2+𝑥+1)=𝑥3+𝑥2+𝑥. Both polynomials take the value 0 at 0 and 1 at 1, so they define the same function on 𝔽2. But they differ, and so do tr𝐴=0 and tr𝐵=1, which Proposition 7.3.18 reads off the coefficients. This is why we defined 𝑝𝐴 through determinants over 𝐹[𝑥], instead of as the function 𝜆↦det(𝜆𝐼𝑛−𝐴).

7.3.3Similarity and operators

The characteristic polynomial of a matrix should depend only on the operator the matrix represents, as the trace and the determinant do (Definition 4.5.30, Definition 5.4.18).

Proposition 7.3.11 (Similar matrices have the same characteristic polynomial). If 𝐴,𝐵∈𝐹𝑛×𝑛 are similar, then 𝑝𝐴=𝑝𝐵.

Proof (by multiplicativity over 𝐹[𝑥]).routine

Uses. Definition 4.5.23 (Similar matrices), Proposition 7.3.2 (Determinants over a commutative ring), Definition 7.3.4 (Characteristic polynomial)

Hint 1. Write 𝐵=𝑆−1𝐴𝑆. Can you write 𝑥𝐼𝑛−𝐵 as a product of three matrices with entries in 𝐹[𝑥]?
Hint 2. 𝑥𝐼𝑛=𝑆−1(𝑥𝐼𝑛)𝑆. Then take determinants in 𝐹[𝑥], using Proposition 7.3.2(c).
Sketch. 𝑥𝐼𝑛−𝑆−1𝐴𝑆=𝑆−1(𝑥𝐼𝑛−𝐴)𝑆 in 𝐹[𝑥]𝑛×𝑛. Taking determinants with Proposition 7.3.2(c) gives 𝑝𝐵=det𝑆−1𝑝𝐴det𝑆=𝑝𝐴, since det𝑆−1det𝑆=det𝐼𝑛=1.

Let 𝐵=𝑆−1𝐴𝑆 with 𝑆∈𝐹𝑛×𝑛 invertible. In 𝐹[𝑥]𝑛×𝑛, by Proposition 7.3.2(a),

𝑆−1(𝑥𝐼𝑛−𝐴)𝑆=𝑥𝑆−1𝑆−𝑆−1𝐴𝑆=𝑥𝐼𝑛−𝐵.

Taking determinants and using Proposition 7.3.2(c) twice,

𝑝𝐵=det(𝑆−1)det(𝑥𝐼𝑛−𝐴)det(𝑆)=det(𝑆−1)det(𝑆)𝑝𝐴=det(𝑆−1𝑆)𝑝𝐴=𝑝𝐴,

where we used that 𝐹[𝑥] is commutative and det𝐼𝑛=1. □

Proof (pointwise, when 𝐹 is infinite).routine

Uses. Theorem 5.4.10 (Multiplicativity), Corollary 7.1.23 (Polynomials and their values), (7.3.1)

Techniques. Polynomial identity arguments (Technique 5.5.24)

Hint 1. Compare the values 𝑝𝐴(𝜆) and 𝑝𝐵(𝜆) for a scalar 𝜆. Only determinants over 𝐹 are involved.
Hint 2. Two polynomials that agree at infinitely many points are equal.
Sketch. For 𝜆∈𝐹, 𝑝𝐵(𝜆)=det(𝑆−1(𝜆𝐼𝑛−𝐴)𝑆)=𝑝𝐴(𝜆) by multiplicativity over 𝐹. If 𝐹 is infinite, polynomials that agree at every point are equal.

Let 𝐵=𝑆−1𝐴𝑆 and 𝜆∈𝐹. By (7.3.1) and multiplicativity over the field 𝐹 (Theorem 5.4.10),

𝑝𝐵(𝜆)=det(𝑆−1(𝜆𝐼𝑛−𝐴)𝑆)=det(𝑆−1)det(𝜆𝐼𝑛−𝐴)det(𝑆)=𝑝𝐴(𝜆).

So 𝑝𝐴 and 𝑝𝐵 agree at every point of 𝐹. If 𝐹 is infinite, Corollary 7.1.23 gives 𝑝𝐴=𝑝𝐵. □

The second proof avoids determinants over 𝐹[𝑥] but only proves equality of functions, which is weaker than equality of polynomials over a finite field (Remark 7.3.10). Over an infinite field this is no loss. The idea of proving a polynomial identity by checking it at many points returns in Proposition 7.3.22.

Definition 7.3.12 (Characteristic polynomial of an operator). Let 𝑉 be finite-dimensional with dim𝑉≥1, and 𝑇∈L(𝑉). The characteristic polynomial of 𝑇 is

𝑝𝑇=𝑝[𝑇]ℬ︀ℬ︀,

where ℬ︀ is any basis of 𝑉. Matrices of 𝑇 in two bases are similar (Theorem 4.5.25), so by Proposition 7.3.11 this does not depend on ℬ︀.

Corollary 7.3.13 (Eigenvalues of an operator). Let 𝑉 be finite-dimensional with dim𝑉≥1, 𝑇∈L(𝑉), and 𝜆∈𝐹. Then 𝜆 is an eigenvalue of 𝑇 if and only if 𝑝𝑇(𝜆)=0.

Proof.routine

Uses. Theorem 7.3.6 (Eigenvalues are the roots of the characteristic polynomial), Definition 7.3.12 (Characteristic polynomial of an operator), Definition 2.3.1 (Matrix of a linear map), Proposition 1.6.20 (Coordinates translate questions about 𝑉 into 𝐹𝑛)

Techniques. Translate between maps and matrices (Technique 2.3.27)

Hint 1. Fix a basis ℬ︀ and let 𝐴=[𝑇]ℬ︀ℬ︀. Compare the eigenvectors of 𝑇 with those of 𝐴.
Hint 2. [𝑇𝐯]ℬ︀=𝐴[𝐯]ℬ︀, and the coordinate map is a linear bijection.
Sketch. With 𝐴=[𝑇]ℬ︀ℬ︀, 𝑇𝐯=𝜆𝐯 if and only if 𝐴[𝐯]ℬ︀=𝜆[𝐯]ℬ︀, and 𝐯≠𝟎 if and only if [𝐯]ℬ︀≠𝟎. So 𝑇 and 𝐴 have the same eigenvalues, which are the roots of 𝑝𝐴=𝑝𝑇.
Let ℬ︀ be a basis of 𝑉 and 𝐴=[𝑇]ℬ︀ℬ︀, so 𝑝𝑇=𝑝𝐴. For 𝐯∈𝑉, [𝑇𝐯]ℬ︀=𝐴[𝐯]ℬ︀ (Definition 2.3.1), and the coordinate map 𝐯↦[𝐯]ℬ︀ is a linear bijection (Proposition 1.6.20). So 𝑇𝐯=𝜆𝐯 with 𝐯≠𝟎 if and only if 𝐴𝐰=𝜆𝐰 with 𝐰=[𝐯]ℬ︀≠𝟎. Hence 𝑇 and 𝐴 have the same eigenvalues, and Theorem 7.3.6 finishes the proof. □
Example 7.3.14 (Characteristic polynomials of operators).
  1. (a)Differentiation 𝐷 on 𝐹[𝑥]≤𝑛 has, in the basis 1,𝑥,…,𝑥𝑛, a strictly upper triangular matrix (it lowers degrees). So 𝑝𝐷=𝑥𝑛+1 by Example 7.3.5(d), and 0 is its only eigenvalue, as found in Example 7.2.19.
  2. (b)On the real vector space 𝑊=span(cos,sin) of functions ℝ→ℝ, differentiation maps cos↦−sin and sin↦cos. Its matrix in the basis cos,sin is (01−10), with characteristic polynomial 𝑥2+1. So it has no real eigenvalue: no nonzero 𝑎cos+𝑏sin is proportional to its own derivative.
  3. (c)The cyclic permutation 𝑍(𝑥,𝑦,𝑧)=(𝑧,𝑥,𝑦) of Example 7.2.36 has matrix (001100010) in the standard basis. Expanding along the first row,

    𝑝𝑍=det(𝑥0−1−1𝑥00−1𝑥)=𝑥⋅𝑥2+(−1)⋅(−1)1+3⋅((−1)(−1)−𝑥⋅0)=𝑥3−1.

The converse of Proposition 7.3.11 is false. The matrices 𝐼2 and (1101) both have characteristic polynomial (𝑥−1)2, but the only matrix similar to 𝐼2 is 𝑆−1𝐼2𝑆=𝐼2 itself. So the characteristic polynomial is one more similarity invariant to try (Technique 4.5.29), stronger than trace and determinant together, but not a complete one. Chapter 10 finds the complete invariant over ℂ (Section 10.3).

Block-triangular matrices and invariant subspaces

Proposition 7.3.15 (Block-triangular matrices). Let 𝑀=(𝐴𝐵0𝐷)∈𝐹𝑛×𝑛, where 𝐴∈𝐹𝑘×𝑘 and 𝐷∈𝐹(𝑛−𝑘)×(𝑛−𝑘) with 1≤𝑘<𝑛. Then 𝑝𝑀=𝑝𝐴𝑝𝐷. The same holds for (𝐴0𝐶𝐷).

Proof.routine

Uses. Proposition 7.3.2 (Determinants over a commutative ring), Definition 7.3.4 (Characteristic polynomial)

Hint 1. Write 𝑥𝐼𝑛−𝑀 in block form. Is it still block triangular?
Hint 2. Apply Proposition 7.3.2(e) over the ring 𝐹[𝑥].
Sketch. 𝑥𝐼𝑛−𝑀=(𝑥𝐼𝑘−𝐴−𝐵0𝑥𝐼𝑛−𝑘−𝐷) is block upper triangular in 𝐹[𝑥]𝑛×𝑛, so its determinant is 𝑝𝐴𝑝𝐷 by Proposition 7.3.2(e). The lower case is the same, or follows by transposing.

In 𝐹[𝑥]𝑛×𝑛,

𝑥𝐼𝑛−𝑀=(𝑥𝐼𝑘−𝐴−𝐵0𝑥𝐼𝑛−𝑘−𝐷),

which is block upper triangular with diagonal blocks 𝑥𝐼𝑘−𝐴 and 𝑥𝐼𝑛−𝑘−𝐷. By Proposition 7.3.2(e), 𝑝𝑀=det(𝑥𝐼𝑘−𝐴)det(𝑥𝐼𝑛−𝑘−𝐷)=𝑝𝐴𝑝𝐷. For (𝐴0𝐶𝐷), the matrix 𝑥𝐼𝑛−𝑀 is block lower triangular, and Proposition 7.3.2(e) applies again. □

Combined with the block form of an operator with an invariant subspace (Remark 7.2.5), this says that an invariant subspace splits off a factor of the characteristic polynomial.

Corollary 7.3.16 (Invariant subspaces give factors). Let 𝑉 be finite-dimensional, 𝑇∈L(𝑉), and 𝑈 a 𝑇-invariant subspace with dim𝑈≥1. Then 𝑝𝑇↾𝑈 divides 𝑝𝑇.

Proof.routine

Uses. Proposition 4.5.35 (Preserved subspaces give block triangular matrices), Proposition 7.3.15 (Block-triangular matrices), Theorem 4.1.16 (Extending to a basis), Definition 7.3.12 (Characteristic polynomial of an operator)

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

Hint 1. Choose a basis of 𝑉 that starts with a basis of 𝑈. What does the matrix of 𝑇 look like?
Hint 2. Proposition 4.5.35 gives (𝐴𝐵0𝐷) with 𝐴=[𝑇↾𝑈]. Now use Proposition 7.3.15.
Sketch. If 𝑈=𝑉 there is nothing to prove. Otherwise extend a basis of 𝑈 to a basis ℬ︀ of 𝑉; then [𝑇]ℬ︀ℬ︀=(𝐴𝐵0𝐷) with 𝐴 the matrix of 𝑇↾𝑈 (Proposition 4.5.35), and 𝑝𝑇=𝑝𝐴𝑝𝐷=𝑝𝑇↾𝑈𝑝𝐷.

If 𝑈=𝑉, then 𝑇↾𝑈=𝑇 and there is nothing to prove. Otherwise let 𝑘=dim𝑈 and 𝑛=dim𝑉, so 1≤𝑘<𝑛. Extend a basis 𝐛1,…,𝐛𝑘 of 𝑈 to a basis ℬ︀=𝐛1,…,𝐛𝑛 of 𝑉 (Theorem 4.1.16). By Proposition 4.5.35,

[𝑇]ℬ︀ℬ︀=(𝐴𝐵0𝐷)with𝐴=[𝑇↾𝑈]ℬ︀𝑈ℬ︀𝑈.

By Proposition 7.3.15, 𝑝𝑇=𝑝𝐴𝑝𝐷=𝑝𝑇↾𝑈𝑝𝐷. □

Exercise 7.3.17.routine Let 𝐴∈𝐹𝑛×𝑛.

  1. (a)Show that 𝑝𝐴𝖳=𝑝𝐴.
  2. (b)For 𝐹=ℂ, let 𝐴 be the entrywise conjugate of 𝐴, and 𝑓 the polynomial whose coefficients are the conjugates of those of 𝑓∈ℂ[𝑥]. Show that 𝑝𝐴=𝑝𝐴, and deduce that the eigenvalues of 𝐴𝖧 are the conjugates of the eigenvalues of 𝐴.

Solution.

Uses. Proposition 7.3.2 (Determinants over a commutative ring), Lemma 7.1.37 (Conjugating polynomials), Theorem 7.3.6 (Eigenvalues are the roots of the characteristic polynomial), Definition 6.3.26 (Conjugate transpose)

Hint 1. For (a), what is (𝑥𝐼𝑛−𝐴)𝖳? Use Proposition 7.3.2(d).
Hint 2. For (b), coefficientwise conjugation ℂ[𝑥]→ℂ[𝑥] respects sums, products, and 1 (Lemma 7.1.37). Apply Proposition 7.3.2(g) to it.
Sketch. (a) (𝑥𝐼𝑛−𝐴)𝖳=𝑥𝐼𝑛−𝐴𝖳, and transposing does not change det. (b) Conjugating every entry of 𝑥𝐼𝑛−𝐴 gives 𝑥𝐼𝑛−𝐴, so 𝑝𝐴=𝑝𝐴 by Proposition 7.3.2(g). With (a), 𝑝𝐴𝖧=𝑝𝐴, whose roots are the conjugates of the roots of 𝑝𝐴 by Lemma 7.1.37(b).

(a) The transpose of 𝑥𝐼𝑛−𝐴 is 𝑥𝐼𝑛−𝐴𝖳, since 𝑥𝐼𝑛 is symmetric. By Proposition 7.3.2(d), 𝑝𝐴𝖳=det((𝑥𝐼𝑛−𝐴)𝖳)=det(𝑥𝐼𝑛−𝐴)=𝑝𝐴.

(b) Let 𝜑:ℂ[𝑥]→ℂ[𝑥], 𝜑(𝑓)=𝑓. By Lemma 7.1.37(a), 𝜑 respects sums and products, and 𝜑(1)=1. Applying 𝜑 to every entry of 𝑥𝐼𝑛−𝐴 gives 𝑥𝐼𝑛−𝐴, because 𝑥=𝑥 and the constant −𝑎𝑖𝑗 becomes −𝑎𝑖𝑗. By Proposition 7.3.2(g), 𝑝𝐴=det𝜑(𝑥𝐼𝑛−𝐴)=𝜑(𝑝𝐴)=𝑝𝐴.

Now 𝐴𝖧=𝐴𝖳, so by (a), 𝑝𝐴𝖧=𝑝𝐴=𝑝𝐴. By Lemma 7.1.37(b), 𝑝𝐴(𝜆)=𝑝𝐴(𝜆), so 𝜆 is a root of 𝑝𝐴 if and only if 𝜆 is a root of 𝑝𝐴. By Theorem 7.3.6, the eigenvalues of 𝐴𝖧 are exactly the conjugates 𝜆 of the eigenvalues 𝜆 of 𝐴.

7.3.4The coefficients: trace and determinant

Proposition 7.3.18 (The first and last coefficients). Let 𝐴∈𝐹𝑛×𝑛. Then 𝑝𝐴 is monic of degree 𝑛, the coefficient of 𝑥𝑛−1 in 𝑝𝐴 is −tr𝐴, and its constant term is (−1)𝑛det𝐴. That is, for 𝑛≥3,

𝑝𝐴=𝑥𝑛−(tr𝐴)𝑥𝑛−1+𝑐𝑛−2𝑥𝑛−2+⋯+𝑐1𝑥+(−1)𝑛det𝐴

for some 𝑐1,…,𝑐𝑛−2∈𝐹, while 𝑝𝐴=𝑥2−(tr𝐴)𝑥+det𝐴 for 𝑛=2 and 𝑝𝐴=𝑥−𝑎11 for 𝑛=1.

Proof.standard

Uses. Definition 7.3.1 (Determinant over a commutative ring), Proposition 7.3.2 (Determinants over a commutative ring), Definition 2.4.26 (Trace), Proposition 7.1.4 (Degree of a product), Proposition 5.3.19 (det(𝐴+𝑡𝐼𝑛) as a function of 𝑡)

Techniques. Read structure off the Leibniz expansion (Technique 5.3.12)

Hint 1. Each entry of 𝑥𝐼𝑛−𝐴 has degree at most 1, and only the diagonal entries involve 𝑥. In the Leibniz term of 𝜎, which factors have degree 1?
Hint 2. Only the factors 𝑥−𝑎𝑗𝑗 with 𝜎(𝑗)=𝑗. A permutation 𝜎≠id moves at least two points. What does that say about the degree of its term?
Hint 3. So 𝑥𝑛 and 𝑥𝑛−1 come only from (𝑥−𝑎11)⋯(𝑥−𝑎𝑛𝑛). For the constant term, evaluate at 0 with (7.3.1).
Sketch. The Leibniz term of 𝜎 is a product of 𝑛 factors of degree at most 1, of which only those with 𝜎(𝑗)=𝑗 have degree 1. A permutation other than id fixes at most 𝑛−2 points, so its term has degree at most 𝑛−2. Thus the coefficients of 𝑥𝑛 and 𝑥𝑛−1 are those of ∏𝑗(𝑥−𝑎𝑗𝑗)=𝑥𝑛−(∑𝑗𝑎𝑗𝑗)𝑥𝑛−1+⋯. The constant term is 𝑝𝐴(0)=det(−𝐴)=(−1)𝑛det𝐴.

Let 𝑀=𝑥𝐼𝑛−𝐴, so 𝑚𝑗𝑗=𝑥−𝑎𝑗𝑗 and 𝑚𝑖𝑗=−𝑎𝑖𝑗 for 𝑖≠𝑗. The Leibniz term of 𝜎∈𝑆𝑛 is

𝑡𝜎=sgn(𝜎)∏𝑗=1𝑛𝑚𝜎(𝑗)𝑗.

Each factor has degree at most 1, with degree 1 exactly when 𝜎(𝑗)=𝑗. So 𝑡𝜎 is 𝟎 or has degree at most the number of fixed points of 𝜎 (Proposition 7.1.4). If 𝜎≠id, pick 𝑗 with 𝜎(𝑗)≠𝑗; then 𝑖=𝜎(𝑗) is also moved, since 𝜎(𝑖)=𝑖=𝜎(𝑗) would contradict injectivity. So 𝜎 fixes at most 𝑛−2 points, and 𝑡𝜎 contributes nothing to the coefficients of 𝑥𝑛 and 𝑥𝑛−1.

These two coefficients therefore come from

𝑡id=(𝑥−𝑎11)(𝑥−𝑎22)⋯(𝑥−𝑎𝑛𝑛).

Expanding by distributivity, a product of 𝑛 factors 𝑥−𝑎𝑗𝑗 is a sum over the ways to choose 𝑥 or −𝑎𝑗𝑗 from each factor. The only way to get 𝑥𝑛 is to choose 𝑥 every time, and the ways to get 𝑥𝑛−1 are to choose −𝑎𝑗𝑗 from exactly one factor 𝑗. So the coefficient of 𝑥𝑛 is 1 and that of 𝑥𝑛−1 is −(𝑎11+⋯+𝑎𝑛𝑛)=−tr𝐴. In particular 𝑝𝐴 is monic of degree 𝑛.

The constant term of a polynomial is its value at 0. By (7.3.1), 𝑝𝐴(0)=det(−𝐴), and multiplying each of the 𝑛 columns of 𝐴 by −1 multiplies det by (−1)𝑛 (Proposition 7.3.2(b)). So the constant term is (−1)𝑛det𝐴. □

This is the computation in the proof of Proposition 5.3.19, done in 𝐹[𝑥] instead of for each scalar 𝑡. The coefficients in between are sums of minors (Exercise 7.3.53).

Remark 7.3.19 (Trace and determinant from the roots). Suppose 𝑝𝐴 splits over 𝐹, say 𝑝𝐴=(𝑥−𝜆1)⋯(𝑥−𝜆𝑛) with 𝜆𝑖∈𝐹. By Proposition 7.1.36, the coefficient of 𝑥𝑛−1 in this product is −(𝜆1+⋯+𝜆𝑛) and its constant term is (−1)𝑛𝜆1⋯𝜆𝑛. Comparing with Proposition 7.3.18 gives

tr𝐴=𝜆1+⋯+𝜆𝑛,det𝐴=𝜆1𝜆2⋯𝜆𝑛.

The 𝜆𝑖 are the eigenvalues of 𝐴 (Theorem 7.3.6), each repeated as often as its multiplicity as a root of 𝑝𝐴, its algebraic multiplicity (Definition 7.4.11). Repetitions matter: (2102) has the single eigenvalue 2 but trace 4=2+2. Over ℂ, 𝑝𝐴 always splits (Corollary 7.1.35), and Theorem 7.5.22 proves the same statement a second way, by triangularizing 𝐴.

For 2×2 matrices the trace and determinant determine 𝑝𝐴=𝑥2−(tr𝐴)𝑥+det𝐴, and with it the eigenvalues. Over ℝ, the quadratic formula says that the number of real eigenvalues is governed by the discriminant

Δ=(tr𝐴)2−4det𝐴:

two real eigenvalues if Δ>0, one repeated real eigenvalue if Δ=0, and a pair of conjugate nonreal eigenvalues if Δ<0 (Theorem B.1.30).

Figure 7.3.1. Real 2×2 matrices placed by (tr𝐴,det𝐴), with Δ=(tr𝐴)2−4det𝐴. Below the parabola det𝐴=(tr𝐴)2/4 there are two real eigenvalues, and above it none. The rotation 𝑅𝜋/2 sits at (0,1), above the parabola. The matrix (2112) of Example 7.2.8 sits at (4,3) and the reflection diag(1,−1) at (0,−1), both below. The shear (1101) sits on the parabola at (2,1), as do 𝐼2 and the matrix of Example 7.2.39 (at (4,4)). A point on the parabola can be a multiple of 𝐼2 or not: the plane cannot tell them apart.
Exercise 7.3.20.routine Let 𝐴∈ℝ3×3 have tr𝐴=2 and det𝐴=−2, and suppose 𝐴−𝐼3 is not invertible and 𝑝𝐴 splits over ℝ. Find 𝑝𝐴 and the eigenvalues of 𝐴.

Solution.

Uses. Theorem 7.3.6 (Eigenvalues are the roots of the characteristic polynomial), Remark 7.3.19 (Trace and determinant from the roots), Theorem 7.2.9 (Eigenvalues and invertibility)

Hint 1. Since 𝐴−𝐼3 is not invertible, 1 is an eigenvalue. Write 𝑝𝐴=(𝑥−1)(𝑥−𝑟)(𝑥−𝑠).
Hint 2. Remark 7.3.19 gives 1+𝑟+𝑠=2 and 𝑟𝑠=−2.
Sketch. 1 is a root, 𝑟+𝑠=1, 𝑟𝑠=−2, so 𝑟,𝑠 are the roots of 𝑡2−𝑡−2=(𝑡−2)(𝑡+1). Hence 𝑝𝐴=(𝑥−1)(𝑥−2)(𝑥+1).
By Theorem 7.2.9, 1 is an eigenvalue of 𝐴, so 𝑝𝐴(1)=0 by Theorem 7.3.6, and since 𝑝𝐴 splits we can write 𝑝𝐴=(𝑥−1)(𝑥−𝑟)(𝑥−𝑠) with 𝑟,𝑠∈ℝ. By Remark 7.3.19, 1+𝑟+𝑠=tr𝐴=2 and 1⋅𝑟𝑠=det𝐴=−2. So 𝑟+𝑠=1 and 𝑟𝑠=−2, which makes 𝑟 and 𝑠 the roots of (𝑡−𝑟)(𝑡−𝑠)=𝑡2−𝑡−2=(𝑡−2)(𝑡+1). Therefore 𝑝𝐴=(𝑥−1)(𝑥−2)(𝑥+1), and the eigenvalues are 1, 2, −1. (The matrix diag(1,2,−1) shows that such an 𝐴 exists.)

Products in either order

For square matrices, 𝐴𝐵 and 𝐵𝐴 are usually different, and they need not be similar. But their characteristic polynomials always agree.

Technique 7.3.21 (Compare 𝐴𝐵 and 𝐵𝐴). The products 𝐴𝐵 and 𝐵𝐴 are close relatives. If 𝐴 is invertible, they are similar: 𝐵𝐴=𝐴−1(𝐴𝐵)𝐴. In general they are not necessarily similar, but they still share the characteristic polynomial (when square) and the nonzero eigenvalues (always). Three ways to transfer information between them are useful. Conjugate by 𝐴 when 𝐴 is invertible, and reduce the general case to it by a polynomial identity in a perturbation parameter. Or embed both products in block matrices of twice the size that are similar. Or transport eigenvectors directly: if 𝐴𝐵𝐯=𝜆𝐯 with 𝜆≠0, then 𝐵𝐯≠𝟎 and 𝐵𝐴(𝐵𝐯)=𝜆𝐵𝐯.

Reach for it whenever a question involves a product in both orders, such as 𝐴𝖳𝐴 and 𝐴𝐴𝖳, or 𝐮𝐯𝖳 and 𝐯𝖳𝐮. The canonical example is Proposition 7.3.22.

Proposition 7.3.22 (𝐴𝐵 and 𝐵𝐴). For 𝐴,𝐵∈𝐹𝑛×𝑛, 𝑝𝐴𝐵=𝑝𝐵𝐴. In particular 𝐴𝐵 and 𝐵𝐴 have the same eigenvalues.
Intuition. If 𝐴 is invertible, then 𝐵𝐴=𝐴−1(𝐴𝐵)𝐴 is similar to 𝐴𝐵, and we are done by Proposition 7.3.11. If 𝐴 is not invertible, 𝐴𝐵 and 𝐵𝐴 need not be similar: for 𝐴=(1000) and 𝐵=(0100), 𝐴𝐵=𝐵≠0 while 𝐵𝐴=0. Two ideas get around this. The first replaces 𝐴 by the invertible matrices 𝐴+𝑡𝐼𝑛, which differ from 𝐴 by a parameter 𝑡 that we then set to 0; this needs the identity to be polynomial in 𝑡, and enough values of 𝑡 to conclude. The second finds a pair of larger matrices, built from 𝐴𝐵 and 𝐵𝐴, that are similar for every 𝐴.

Proof (by block matrices).standard

Uses. Proposition 7.3.11 (Similar matrices have the same characteristic polynomial), Proposition 7.3.15 (Block-triangular matrices), Theorem 2.4.34 (Block multiplication), Corollary 7.1.7 (No zero divisors, cancellation, and units)

Techniques. Compare 𝐴𝐵 and 𝐵𝐴 (Technique 7.3.21), Block decomposition (Technique 2.4.35)

Hint 1. Look for a 2𝑛×2𝑛 matrix 𝑆 with (𝐴𝐵0𝐵0)𝑆=𝑆(00𝐵𝐵𝐴). Try 𝑆 block upper triangular with identity blocks on the diagonal.
Hint 2. 𝑆=(𝐼𝑛𝐴0𝐼𝑛) works. Both block matrices are block triangular; compute their characteristic polynomials.
Hint 3. You get 𝑥𝑛𝑝𝐴𝐵=𝑥𝑛𝑝𝐵𝐴. Cancel 𝑥𝑛 in 𝐹[𝑥].
Sketch. With 𝑋=(𝐴𝐵0𝐵0), 𝑌=(00𝐵𝐵𝐴), and 𝑆=(𝐼𝑛𝐴0𝐼𝑛), block multiplication gives 𝑋𝑆=(𝐴𝐵𝐴𝐵𝐴𝐵𝐵𝐴)=𝑆𝑌, and 𝑆 is invertible, so 𝑋 and 𝑌 are similar. By Proposition 7.3.15, 𝑝𝑋=𝑝𝐴𝐵𝑥𝑛 and 𝑝𝑌=𝑥𝑛𝑝𝐵𝐴. Cancel 𝑥𝑛 (Corollary 7.1.7).

Consider the 2𝑛×2𝑛 matrices

𝑋=(𝐴𝐵0𝐵0),𝑌=(00𝐵𝐵𝐴),𝑆=(𝐼𝑛𝐴0𝐼𝑛).

By block multiplication (Theorem 2.4.34),

𝑋𝑆=(𝐴𝐵𝐴𝐵𝐴𝐵𝐵𝐴)=𝑆𝑌.

The matrix 𝑆 is invertible, with inverse (𝐼𝑛−𝐴0𝐼𝑛) (multiply out). So 𝑌=𝑆−1𝑋𝑆, and 𝑝𝑋=𝑝𝑌 by Proposition 7.3.11. Both matrices are block triangular, so by Proposition 7.3.15,

𝑝𝑋=𝑝𝐴𝐵𝑝0𝑛×𝑛=𝑝𝐴𝐵𝑥𝑛,𝑝𝑌=𝑝0𝑛×𝑛𝑝𝐵𝐴=𝑥𝑛𝑝𝐵𝐴.

So 𝑥𝑛𝑝𝐴𝐵=𝑥𝑛𝑝𝐵𝐴, and since 𝑥𝑛≠𝟎, cancellation in 𝐹[𝑥] (Corollary 7.1.7(b)) gives 𝑝𝐴𝐵=𝑝𝐵𝐴. The statement about eigenvalues follows from Theorem 7.3.6. □

Proof (by perturbation, when 𝐹 is infinite).hard

Uses. Proposition 7.3.11 (Similar matrices have the same characteristic polynomial), Proposition 5.3.19 (det(𝐴+𝑡𝐼𝑛) as a function of 𝑡), Proposition 7.3.2 (Determinants over a commutative ring), Proposition 7.1.22 (Counting roots), Corollary 7.1.23 (Polynomials and their values), Theorem 5.4.7 (Invertibility and the determinant)

Techniques. Compare 𝐴𝐵 and 𝐵𝐴 (Technique 7.3.21), Polynomial identity arguments (Technique 5.5.24), det is a polynomial in the entries (Technique 5.3.18)

Hint 1. Fix 𝜆∈𝐹. For a scalar 𝑠 with 𝐴+𝑠𝐼𝑛 invertible, compare det(𝜆𝐼𝑛−(𝐴+𝑠𝐼𝑛)𝐵) with det(𝜆𝐼𝑛−𝐵(𝐴+𝑠𝐼𝑛)).
Hint 2. They are equal, by similarity. Now regard 𝑠 as an indeterminate 𝑡: both sides are values of polynomials in 𝑡, and their difference has infinitely many roots.
Hint 3. Set 𝑡=0 to get 𝑝𝐴𝐵(𝜆)=𝑝𝐵𝐴(𝜆) for every 𝜆, and use that 𝐹 is infinite once more.
Sketch. Fix 𝜆. In 𝐹[𝑡] let 𝑔=det(𝜆𝐼𝑛−(𝐴+𝑡𝐼𝑛)𝐵)−det(𝜆𝐼𝑛−𝐵(𝐴+𝑡𝐼𝑛)). For each 𝑠∈𝐹 with 𝐴+𝑠𝐼𝑛 invertible, 𝑔(𝑠)=0 by similarity; this excludes at most 𝑛 values of 𝑠. So 𝑔 has infinitely many roots and 𝑔=𝟎. Then 𝑔(0)=𝑝𝐴𝐵(𝜆)−𝑝𝐵𝐴(𝜆)=0 for every 𝜆∈𝐹, and Corollary 7.1.23 finishes.

Assume 𝐹 is infinite, and fix 𝜆∈𝐹. Work in 𝐹[𝑡], polynomials in a new indeterminate 𝑡, and let

𝑔=det(𝜆𝐼𝑛−(𝐴+𝑡𝐼𝑛)𝐵)−det(𝜆𝐼𝑛−𝐵(𝐴+𝑡𝐼𝑛))∈𝐹[𝑡].

For 𝑠∈𝐹, evaluating every entry at 𝑡=𝑠 (Proposition 7.3.2(g)) gives

𝑔(𝑠)=det(𝜆𝐼𝑛−(𝐴+𝑠𝐼𝑛)𝐵)−det(𝜆𝐼𝑛−𝐵(𝐴+𝑠𝐼𝑛)).

If 𝐴+𝑠𝐼𝑛 is invertible, then 𝐵(𝐴+𝑠𝐼𝑛)=(𝐴+𝑠𝐼𝑛)−1((𝐴+𝑠𝐼𝑛)𝐵)(𝐴+𝑠𝐼𝑛), so the two matrices (𝐴+𝑠𝐼𝑛)𝐵 and 𝐵(𝐴+𝑠𝐼𝑛) are similar, and their characteristic polynomials agree at 𝜆 (Proposition 7.3.11). So 𝑔(𝑠)=0.

By Proposition 5.3.19 and Theorem 5.4.7, 𝐴+𝑠𝐼𝑛 fails to be invertible only for roots 𝑠 of a monic polynomial of degree 𝑛, that is, for at most 𝑛 values of 𝑠 (Proposition 7.1.22). Since 𝐹 is infinite, 𝑔 has infinitely many roots, so 𝑔=𝟎, again by Proposition 7.1.22. In particular 𝑔(0)=0, which says 𝑝𝐴𝐵(𝜆)=𝑝𝐵𝐴(𝜆) by (7.3.1). This holds for every 𝜆∈𝐹, so 𝑝𝐴𝐵=𝑝𝐵𝐴 by Corollary 7.1.23. □

The first proof is shorter and works over every field. The second shows a pattern that recurs: prove an identity for invertible matrices, where it is easy, and extend it to all matrices because both sides are polynomials. Over a finite field there are too few values of 𝑡, and this pattern alone does not suffice.

7.3.5Companion matrices

Which polynomials occur as characteristic polynomials? Every characteristic polynomial of an 𝑛×𝑛 matrix is monic of degree 𝑛. Conversely, every such polynomial occurs, and the matrix realizing it already appeared twice: in Proposition 5.5.20 as a determinant computation, and after Proposition 7.2.35 as the matrix of an operator on a Krylov subspace.

Definition 7.3.23 (Companion matrix). Let 𝑝=𝑥𝑛+𝑎𝑛−1𝑥𝑛−1+⋯+𝑎1𝑥+𝑎0∈𝐹[𝑥] be monic of degree 𝑛≥1. The companion matrix of 𝑝 is

𝐶𝑝=(00⋯0−𝑎010⋯0−𝑎101⋯0−𝑎2⋮⋮⋱⋮⋮00⋯1−𝑎𝑛−1)∈𝐹𝑛×𝑛,

whose column 𝑘 is 𝐞𝑘+1 for 𝑘<𝑛 and whose last column is −(𝑎0,𝑎1,…,𝑎𝑛−1).

For example, 𝐶𝑥−𝑐=(𝑐), the companion matrix of 𝑥2+1 is the rotation (0−110), and the companion matrix of 𝑥3−1 is the matrix of the cyclic permutation 𝑍 in Example 7.3.14. In general 𝐶𝑝𝐞𝑘=𝐞𝑘+1 for 𝑘<𝑛: 𝐶𝑝 shifts the standard basis forward, and the last basis vector goes to 𝐶𝑝𝐞𝑛=−(𝑎0𝐞1+⋯+𝑎𝑛−1𝐞𝑛). Some books use the transpose 𝐶𝑝𝖳 instead; it has the same characteristic polynomial (Exercise 7.3.17).

Theorem 7.3.24 (Every monic polynomial is a characteristic polynomial). Let 𝑝∈𝐹[𝑥] be monic of degree 𝑛≥1. Then 𝑝𝐶𝑝=𝑝.

Proof (by rerunning Proposition 5.5.20).standard

Uses. Proposition 5.5.20 (Companion matrices), Proposition 7.3.2 (Determinants over a commutative ring)

Hint 1. Proposition 5.5.20 proves det(𝑡𝐼𝑛−𝐶)=𝑝(𝑡) for every scalar 𝑡∈𝐹. Over a finite field this is weaker than an identity in 𝐹[𝑥]. Which tools did its first proof use?
Hint 2. It used cofactor expansion along the first row, a triangular determinant, and induction on 𝑛. All of these are available over 𝐹[𝑥] by Proposition 7.3.2.
Sketch. The first proof of Proposition 5.5.20 expands det(𝑡𝐼𝑛−𝐶𝑝) along the first row, recognizes one minor as the same kind of matrix for 𝑞=(𝑝−𝑎0)/𝑥 and the other as triangular with −1′s on the diagonal, and inducts. Each step is valid with 𝑡 replaced by the indeterminate 𝑥, by Proposition 7.3.2(e) and (f). So det(𝑥𝐼𝑛−𝐶𝑝)=𝑥𝑞+𝑎0=𝑝 in 𝐹[𝑥].
The first proof of Proposition 5.5.20 computes det(𝑡𝐼𝑛−𝐶𝑝) by induction on 𝑛. It expands along the first row, identifies the minor at position (1,1) as 𝑡𝐼𝑛−1−𝐶𝑞 for 𝑞=𝑎1+𝑎2𝑥+⋯+𝑎𝑛−1𝑥𝑛−2+𝑥𝑛−1, and the minor at position (1,𝑛) as an upper triangular matrix with diagonal entries −1. It then concludes det(𝑡𝐼𝑛−𝐶𝑝)=𝑡𝑞(𝑡)+𝑎0=𝑝(𝑡). Every step uses only cofactor expansion, the determinant of a triangular matrix, and ring arithmetic. By Proposition 7.3.2(e) and (f) these hold in 𝐹[𝑥]𝑛×𝑛, so the same induction, with the indeterminate 𝑥 in place of the scalar 𝑡, proves det(𝑥𝐼𝑛−𝐶𝑝)=𝑥𝑞+𝑎0=𝑝 in 𝐹[𝑥]. □

Proof (by expanding along the last column).standard

Uses. Proposition 7.3.2 (Determinants over a commutative ring), Definition 7.3.23 (Companion matrix)

Techniques. Read structure off the Leibniz expansion (Technique 5.3.12)

Hint 1. The last column of 𝑥𝐼𝑛−𝐶𝑝 is (𝑎0,𝑎1,…,𝑎𝑛−2,𝑥+𝑎𝑛−1). Expand along it.
Hint 2. Deleting row 𝑘 and column 𝑛 leaves a block-diagonal matrix: a lower triangular block of size 𝑘−1 with 𝑥 on the diagonal, and an upper triangular block of size 𝑛−𝑘 with −1 on the diagonal.
Hint 3. So the minor is 𝑥𝑘−1(−1)𝑛−𝑘, and the sign (−1)𝑘+𝑛 cancels the (−1)𝑛−𝑘.
Sketch. Let 𝑀=𝑥𝐼𝑛−𝐶𝑝. The minor 𝑀𝑘𝑛̂ is block diagonal, with a lower triangular (𝑘−1)×(𝑘−1) block with diagonal 𝑥 and an upper triangular (𝑛−𝑘)×(𝑛−𝑘) block with diagonal −1, so det𝑀𝑘𝑛̂=𝑥𝑘−1(−1)𝑛−𝑘. Expanding along column 𝑛: det𝑀=∑𝑘(−1)𝑘+𝑛𝑚𝑘𝑛𝑥𝑘−1(−1)𝑛−𝑘=∑𝑘𝑚𝑘𝑛𝑥𝑘−1=𝑎0+𝑎1𝑥+⋯+𝑎𝑛−2𝑥𝑛−2+(𝑥+𝑎𝑛−1)𝑥𝑛−1=𝑝.

For 𝑛=1, 𝐶𝑝=(−𝑎0) and det(𝑥𝐼1−𝐶𝑝)=𝑥+𝑎0=𝑝. Let 𝑛≥2 and 𝑀=𝑥𝐼𝑛−𝐶𝑝. Row 𝑘 of 𝑀 has −1 in column 𝑘−1 (if 𝑘≥2), 𝑥 in column 𝑘 (if 𝑘≤𝑛−1), and 𝑚𝑘𝑛=𝑎𝑘−1 in column 𝑛, except that 𝑚𝑛𝑛=𝑥+𝑎𝑛−1; all other entries are 0. Fix 𝑘 and delete row 𝑘 and column 𝑛.

  • Rows 1,…,𝑘−1 have their nonzero entries (outside column 𝑛) in columns 1,…,𝑘−1. They form a (𝑘−1)×(𝑘−1) lower triangular block 𝐿 with diagonal entries 𝑥.
  • Rows 𝑘+1,…,𝑛 have their nonzero entries (outside column 𝑛) in columns 𝑘,…,𝑛−1: row 𝑖 has −1 in column 𝑖−1 and 𝑥 in column 𝑖 (if 𝑖≤𝑛−1). They form an (𝑛−𝑘)×(𝑛−𝑘) upper triangular block 𝑈 with diagonal entries −1.

So 𝑀𝑘𝑛̂=(𝐿00𝑈), and by Proposition 7.3.2(e), det𝑀𝑘𝑛̂=𝑥𝑘−1(−1)𝑛−𝑘. (For 𝑘=1 there is no block 𝐿, and for 𝑘=𝑛 no block 𝑈; the formula still holds.) Expanding along column 𝑛 (Proposition 7.3.2(f)),

det𝑀=∑𝑘=1𝑛(−1)𝑘+𝑛𝑚𝑘𝑛𝑥𝑘−1(−1)𝑛−𝑘=∑𝑘=1𝑛𝑚𝑘𝑛𝑥𝑘−1=𝑎0+𝑎1𝑥+⋯+𝑎𝑛−2𝑥𝑛−2+(𝑥+𝑎𝑛−1)𝑥𝑛−1=𝑝,

since (−1)𝑘+𝑛(−1)𝑛−𝑘=(−1)2𝑛=1. □

The first proof builds 𝑝 from the inside out, as in Horner’s scheme. The second reads all coefficients off at once: the cofactor of the entry in row 𝑘 of the last column is exactly 𝑥𝑘−1, so the last column of 𝑥𝐼𝑛−𝐶𝑝 is the list of coefficients of 𝑝 and the expansion reassembles 𝑝.

Example 7.3.25 (Matrices with prescribed eigenvalues).
  1. (a)The polynomial 𝑝=𝑥3−2 has no root in ℚ (a rational root would be an integer dividing 2, by the argument of Exercise 7.3.42, and none of ±1,±2 works). So 𝐶𝑝=(002100010)∈ℚ3×3 has no eigenvalues in ℚ, although over ℝ it has the eigenvalue 23.
  2. (b)To get a real 2×2 matrix with eigenvalues 3 and −1 that is neither diagonal nor triangular, take 𝑝=(𝑥−3)(𝑥+1)=𝑥2−2𝑥−3 and 𝐶𝑝=(0312).

The connection with Krylov subspaces gives the companion matrix its meaning. By Proposition 7.2.35, if 𝐯≠𝟎 and 𝐾𝑇(𝐯) is finite-dimensional, then 𝐯,𝑇𝐯,…,𝑇𝑑−1𝐯 is a basis of 𝐾𝑇(𝐯), and

𝑇𝑑𝐯=−(𝑐0𝐯+𝑐1𝑇𝐯+⋯+𝑐𝑑−1𝑇𝑑−1𝐯).

In this basis, 𝑇↾𝐾𝑇(𝐯) moves each basis vector to the next and sends the last one to the vector with coordinates −(𝑐0,…,𝑐𝑑−1). That is the description of 𝐶𝑝 for 𝑝=𝑥𝑑+𝑐𝑑−1𝑥𝑑−1+⋯+𝑐0.

Proposition 7.3.26 (Krylov subspaces and the characteristic polynomial). Let 𝑉 be finite-dimensional, 𝑇∈L(𝑉), and 𝐯∈𝑉 with 𝐯≠𝟎. Let 𝑑 and the monic polynomial 𝑝 of degree 𝑑 be as in Proposition 7.2.35, so that 𝑝(𝑇)𝐯=𝟎. Then

  1. (a)the matrix of 𝑇↾𝐾𝑇(𝐯) in the basis 𝐯,𝑇𝐯,…,𝑇𝑑−1𝐯 is 𝐶𝑝;
  2. (b)𝑝𝑇↾𝐾𝑇(𝐯)=𝑝;
  3. (c)𝑝 divides 𝑝𝑇.

Proof.routine

Uses. Proposition 7.2.35 (A basis of a Krylov subspace), Proposition 7.2.34 (The smallest invariant subspace containing a vector), Definition 2.3.1 (Matrix of a linear map), Theorem 7.3.24 (Every monic polynomial is a characteristic polynomial), Corollary 7.3.16 (Invariant subspaces give factors)

Hint 1. Compute 𝑇 on each basis vector 𝑇𝑖𝐯 and write the result in the basis.
Hint 2. For (b) use Theorem 7.3.24, and for (c) use that 𝐾𝑇(𝐯) is 𝑇-invariant.
Sketch. 𝑇(𝑇𝑖𝐯)=𝑇𝑖+1𝐯 for 𝑖<𝑑−1, and 𝑇(𝑇𝑑−1𝐯)=−∑𝑖𝑐𝑖𝑇𝑖𝐯, so the columns are 𝐞2,…,𝐞𝑑 and −(𝑐0,…,𝑐𝑑−1), which is 𝐶𝑝. Then the characteristic polynomial is 𝑝𝐶𝑝=𝑝, and 𝑝 divides 𝑝𝑇 because 𝐾𝑇(𝐯) is invariant (Corollary 7.3.16).
Write 𝑝=𝑥𝑑+𝑐𝑑−1𝑥𝑑−1+⋯+𝑐0 and 𝐛𝑖+1=𝑇𝑖𝐯 for 0≤𝑖≤𝑑−1. For 𝑖<𝑑−1, 𝑇𝐛𝑖+1=𝑇𝑖+1𝐯=𝐛𝑖+2, and 𝑇𝐛𝑑=𝑇𝑑𝐯=−(𝑐0𝐛1+⋯+𝑐𝑑−1𝐛𝑑) because 𝑝(𝑇)𝐯=𝟎. So column 𝑘 of the matrix of 𝑇↾𝐾𝑇(𝐯) is 𝐞𝑘+1 for 𝑘<𝑑, and the last column is −(𝑐0,…,𝑐𝑑−1) (Definition 2.3.1). This is 𝐶𝑝, proving (a). Then (b) follows from (a), Definition 7.3.12, and Theorem 7.3.24. Since 𝐾𝑇(𝐯) is 𝑇-invariant (Proposition 7.2.34) and has dimension 𝑑≥1, (c) follows from (b) and Corollary 7.3.16. □

So every nonzero vector produces a factor of 𝑝𝑇, namely the polynomial that records the first dependence among 𝐯,𝑇𝐯,𝑇2𝐯,…. In Example 7.2.36, the vectors (1,1,1), (1,−1,0), and 𝐞1 produced 𝑥−1, 𝑥2+𝑥+1, and 𝑥3−1, and each divides 𝑝𝑍=𝑥3−1. When 𝐾𝑇(𝐯)=𝑉, part (b) says 𝑝𝑇=𝑝 (Exercise 7.3.50).

7.3.6Eigenvalues exist over ℂ

Example 7.2.18 showed that a real operator need not have an eigenvalue, and Example 7.2.20 that an operator on an infinite-dimensional space need not either. Neither obstruction survives over ℂ in finite dimensions.

Theorem 7.3.27 (Existence of eigenvalues over ℂ). Let 𝑉 be a finite-dimensional complex vector space with 𝑉≠{𝟎}, and let 𝑇∈L(𝑉). Then 𝑇 has an eigenvalue. In particular, every 𝐴∈ℂ𝑛×𝑛 with 𝑛≥1 has an eigenvalue.

The two proofs below locate the eigenvalue in two different polynomials. The first uses 𝑝𝑇, which collects all eigenvalues at once and needs determinants. The second uses the polynomial of a single Krylov sequence, which only needs linear dependence.

Proof (via the characteristic polynomial).routine

Uses. Proposition 7.3.18 (The first and last coefficients), Theorem 7.1.32 (Fundamental Theorem of Algebra), Corollary 7.3.13 (Eigenvalues of an operator)

Hint 1. What is the degree of 𝑝𝑇?
Hint 2. A polynomial of degree at least 1 over ℂ has a root (Theorem 7.1.32).
Sketch. 𝑝𝑇 has degree dim𝑉≥1, so it has a root 𝜆∈ℂ by the FTA, and 𝜆 is an eigenvalue by Corollary 7.3.13.
Let 𝑛=dim𝑉≥1. By Proposition 7.3.18, 𝑝𝑇 is monic of degree 𝑛, in particular nonconstant. By Theorem 7.1.32 it has a root 𝜆∈ℂ, and 𝜆 is an eigenvalue of 𝑇 by Corollary 7.3.13. □

Intuition (A proof without determinants). Take any 𝐯≠𝟎 and look at 𝐯,𝑇𝐯,𝑇2𝐯,…. In a space of dimension 𝑛, the 𝑛+1 vectors 𝐯,…,𝑇𝑛𝐯 are dependent, so some nonzero polynomial 𝑞 has 𝑞(𝑇)𝐯=𝟎. Over ℂ we can factor 𝑞=𝑐(𝑥−𝜆1)⋯(𝑥−𝜆𝑚), so

𝑐(𝑇−𝜆1id)(𝑇−𝜆2id)⋯(𝑇−𝜆𝑚id)𝐯=𝟎

with 𝐯≠𝟎. A composite of injective maps is injective, so some factor 𝑇−𝜆𝑗id is not injective, and 𝜆𝑗 is an eigenvalue. If we use the first dependence, the polynomial 𝑝 of Proposition 7.2.35, we can even say which vector is the eigenvector: split off one linear factor, 𝑝=(𝑥−𝜆)𝑞, and apply 𝑞(𝑇) to 𝐯. The result is nonzero because 𝑞 has smaller degree than 𝑝, and 𝑇−𝜆id kills it.

Technique 7.3.28 (Krylov sequences must become dependent). In a finite-dimensional space, the sequence 𝐯,𝑇𝐯,𝑇2𝐯,… eventually becomes dependent, and the first dependence is a monic polynomial 𝑝 with 𝑝(𝑇)𝐯=𝟎 and no nonzero polynomial of smaller degree killing 𝐯 (Proposition 7.2.35). This converts a question about 𝑇 into a question about one polynomial: factor 𝑝, and use the factors to build vectors. If 𝑝=𝑓𝑞 with deg𝑓≥1, then 𝐰=𝑞(𝑇)𝐯 is nonzero and 𝑓(𝑇)𝐰=𝟎, so a factor 𝑓=𝑥−𝜆 produces an eigenvector, and a quadratic factor produces a two-dimensional invariant subspace (Theorem 7.3.31).

Reach for it when you need an eigenvector or a small invariant subspace and want to avoid determinants, or when 𝑉 is infinite-dimensional but some Krylov subspace is not. The canonical example is the second proof of Theorem 7.3.27.

Proof (via a Krylov sequence, without determinants).standard

Uses. Proposition 7.2.35 (A basis of a Krylov subspace), Theorem 7.1.32 (Fundamental Theorem of Algebra), Theorem 7.1.20 (Factor theorem), Proposition 7.2.27 (Evaluating polynomials at an operator), Definition 7.2.6 (Eigenvalue, eigenvector)

Techniques. Krylov sequences must become dependent (Technique 7.3.28), Too many vectors force a dependence (Technique 4.1.30)

Hint 1. Fix 𝐯≠𝟎. The sequence 𝐯,𝑇𝐯,𝑇2𝐯,… cannot stay independent. Use the first dependence, from Proposition 7.2.35.
Hint 2. This gives a monic 𝑝 of degree 𝑑≥1 with 𝑝(𝑇)𝐯=𝟎, and no nonzero polynomial of degree less than 𝑑 kills 𝐯. Use the FTA to split off a linear factor 𝑥−𝜆.
Hint 3. Write 𝑝=(𝑥−𝜆)𝑞. Then 𝑞(𝑇)𝐯≠𝟎 (why?) and (𝑇−𝜆id)𝑞(𝑇)𝐯=𝑝(𝑇)𝐯.
Sketch. Let 𝐯≠𝟎 and let 𝑝, of degree 𝑑≥1, record the first dependence among 𝐯,𝑇𝐯,… (Proposition 7.2.35). By the FTA and the factor theorem, 𝑝=(𝑥−𝜆)𝑞 with 𝑞 monic of degree 𝑑−1. Then 𝐰=𝑞(𝑇)𝐯≠𝟎 by the minimality of 𝑑, and (𝑇−𝜆id)𝐰=𝑝(𝑇)𝐯=𝟎. So 𝐰 is an eigenvector for 𝜆.

Choose 𝐯∈𝑉 with 𝐯≠𝟎. Since 𝑉 is finite-dimensional, so is 𝐾𝑇(𝐯), and Proposition 7.2.35 gives 𝑑≥1 and a monic polynomial 𝑝∈ℂ[𝑥] of degree 𝑑 with 𝑝(𝑇)𝐯=𝟎, such that no nonzero polynomial of degree less than 𝑑 sends 𝐯 to 𝟎.

As deg𝑝=𝑑≥1, the FTA (Theorem 7.1.32) gives a root 𝜆∈ℂ of 𝑝, and the factor theorem (Theorem 7.1.20) gives 𝑝=(𝑥−𝜆)𝑞 for some 𝑞∈ℂ[𝑥]. Comparing leading coefficients and degrees, 𝑞 is monic of degree 𝑑−1; in particular 𝑞≠𝟎. Let 𝐰=𝑞(𝑇)𝐯. Since 𝑞 is nonzero of degree less than 𝑑, 𝐰≠𝟎. By Proposition 7.2.27(a),

(𝑇−𝜆id)𝐰=((𝑥−𝜆)𝑞)(𝑇)𝐯=𝑝(𝑇)𝐯=𝟎.

So 𝑇𝐰=𝜆𝐰 with 𝐰≠𝟎, and 𝜆 is an eigenvalue of 𝑇 (Definition 7.2.6). □

Example 7.3.29 (Finding an eigenvector from a Krylov sequence). Let 𝐴=(0−110)∈ℂ2×2 and 𝐯=𝐞1. Then 𝐴𝐯=𝐞2, which is not a multiple of 𝐞1, and 𝐴2𝐯=−𝐞1=−𝐯. So the first dependence is 𝐴2𝐯+𝐯=𝟎, and 𝑝=𝑥2+1=(𝑥−𝑖)(𝑥+𝑖). Splitting off 𝑥−𝑖 leaves 𝑞=𝑥+𝑖, and

𝐰=(𝐴+𝑖𝐼2)𝐞1=(𝑖,1),𝐴𝐰=(−1,𝑖)=𝑖(𝑖,1).

Splitting off 𝑥+𝑖 instead gives (𝐴−𝑖𝐼2)𝐞1=(−𝑖,1), an eigenvector for −𝑖. For the cyclic permutation 𝑍 and 𝐯=𝐞1, the first dependence is 𝑥3−1=(𝑥−1)(𝑥2+𝑥+1) (Example 7.2.36), and (𝑍2+𝑍+id)𝐞1=(1,1,1) is the eigenvector for 1.

Intuition (Comparing the two proofs). The first proof is shorter, but it rests on determinants over 𝐹[𝑥], and it finds an eigenvalue without an eigenvector. The second proof needs only linear dependence and polynomial arithmetic. It also uses less: it never needs all of 𝑉, only the finite-dimensional subspace 𝐾𝑇(𝐯). So it shows, for instance, that an operator on an infinite-dimensional complex space has an eigenvalue as soon as some nonzero vector has a finite-dimensional Krylov subspace. On the other hand, the first proof locates every eigenvalue at once, with multiplicities, which the second does not.

Both proofs use exactly one property of ℂ: every nonconstant polynomial over ℂ has a root. So both prove the theorem over any field with this property, and over ℝ both fail at the same step, with the rotation as the witness: its characteristic polynomial is 𝑥2+1, and so is the Krylov polynomial of every 𝐯≠𝟎. Over a finite field the theorem always fails (Exercise 7.3.44).

7.3.7Real operators Beyond 21-242

Real operators and complexification are covered in LADR and LADW (5.8), but not in 21-242 lecture.

Over ℝ an operator may have no eigenvalue, but it cannot avoid small invariant subspaces altogether.

Proposition 7.3.30 (Odd dimension). Let 𝑉 be a real vector space of odd dimension. Then every 𝑇∈L(𝑉) has an eigenvalue.

Proof.routine

Uses. Proposition 7.3.18 (The first and last coefficients), Corollary 7.1.43 (Odd degree gives a real root), Corollary 7.3.13 (Eigenvalues of an operator)

Hint. What is the degree of 𝑝𝑇, and what do you know about real polynomials of that degree?
Sketch. 𝑝𝑇∈ℝ[𝑥] has odd degree dim𝑉, so it has a real root (Corollary 7.1.43), which is an eigenvalue.
By Proposition 7.3.18, 𝑝𝑇∈ℝ[𝑥] has degree dim𝑉, which is odd. By Corollary 7.1.43 it has a real root 𝜆, and 𝜆 is an eigenvalue of 𝑇 by Corollary 7.3.13. □

For example, every operator on ℝ3 maps some line to itself; for a rotation this is its axis. In even dimension there may be no eigenvalue: the operator on ℝ4 with matrix (𝑅𝜋/200𝑅𝜋/2) has characteristic polynomial (𝑥2+1)2 by Proposition 7.3.15. It still has invariant planes, such as span(𝐞1,𝐞2).

Theorem 7.3.31 (Small invariant subspaces). Let 𝑉 be a finite-dimensional real vector space with 𝑉≠{𝟎}, and let 𝑇∈L(𝑉). Then 𝑇 has an invariant subspace of dimension 1 or 2.

Proof (via a Krylov sequence).standard

Uses. Proposition 7.2.35 (A basis of a Krylov subspace), Theorem 7.1.39 (Factorization over ℝ), Proposition 7.2.27 (Evaluating polynomials at an operator), Proposition 7.2.4 (Invariance is checked on a spanning list), Proposition 7.2.7 (Invariant lines are spanned by eigenvectors)

Techniques. Krylov sequences must become dependent (Technique 7.3.28)

Hint 1. Imitate the second proof of Theorem 7.3.27. What replaces the FTA over ℝ?
Hint 2. By Theorem 7.1.39, the polynomial 𝑝 of the first dependence has a monic factor 𝑓 of degree 1 or 2. Write 𝑝=𝑓𝑞 and set 𝐰=𝑞(𝑇)𝐯.
Hint 3. If 𝑓=𝑥2+𝑏𝑥+𝑐, then 𝑇2𝐰=−𝑏𝑇𝐰−𝑐𝐰. Which subspace is invariant?
Sketch. Let 𝑝 record the first dependence of a Krylov sequence. Over ℝ, 𝑝 has a monic factor 𝑓 of degree 1 or 2 (Theorem 7.1.39); write 𝑝=𝑓𝑞. Then 𝐰=𝑞(𝑇)𝐯≠𝟎 and 𝑓(𝑇)𝐰=𝟎. If 𝑓=𝑥−𝜆, 𝐰 spans an invariant line. If 𝑓=𝑥2+𝑏𝑥+𝑐, then 𝑇(𝑇𝐰)=−𝑏𝑇𝐰−𝑐𝐰, so span(𝐰,𝑇𝐰) is invariant, of dimension 1 or 2.

Choose 𝐯≠𝟎 and let 𝑝, monic of degree 𝑑≥1, be as in Proposition 7.2.35. By Theorem 7.1.39, 𝑝 is a product of monic factors of degree 1 and 2 (times its leading coefficient 1), so 𝑝=𝑓𝑞 with 𝑓 monic of degree 1 or 2 and 𝑞 monic of degree less than 𝑑. Let 𝐰=𝑞(𝑇)𝐯. Then 𝐰≠𝟎 by the minimality in Proposition 7.2.35, and 𝑓(𝑇)𝐰=𝑝(𝑇)𝐯=𝟎 by Proposition 7.2.27(a).

If 𝑓=𝑥−𝜆, then 𝑇𝐰=𝜆𝐰, so span(𝐰) is an invariant line (Proposition 7.2.7). If 𝑓=𝑥2+𝑏𝑥+𝑐, then 𝑇2𝐰+𝑏𝑇𝐰+𝑐𝐰=𝟎. Let 𝑈=span(𝐰,𝑇𝐰). Then 𝑇𝐰∈𝑈 and 𝑇(𝑇𝐰)=−𝑏𝑇𝐰−𝑐𝐰∈𝑈, so 𝑈 is 𝑇-invariant by Proposition 7.2.4. It has dimension 1 or 2, since 𝐰≠𝟎 and it is spanned by two vectors. □

A second proof, through complex eigenvectors, follows in the next subsection.

7.3.8Complexification Beyond 21-242

To use complex eigenvalues for a real operator, we need a complex vector space for it to act on. For ℝ𝑛 the answer is ℂ𝑛: a real matrix is also a complex matrix. For an abstract real vector space 𝑉, we build the complex space by adjoining “imaginary parts”, just as ℂ is built from ℝ (Definition 1.1.29).

Definition 7.3.32 (Complexification). Let 𝑉 be a real vector space. Its complexification 𝑉ℂ is the set 𝑉×𝑉, whose element (𝐮,𝐯) we write 𝐮+𝑖𝐯, with addition (𝐮+𝑖𝐯)+(𝐮′+𝑖𝐯′)=(𝐮+𝐮′)+𝑖(𝐯+𝐯′) and scalar multiplication by 𝑎+𝑏𝑖∈ℂ (𝑎,𝑏∈ℝ) given by

(𝑎+𝑏𝑖)(𝐮+𝑖𝐯)=(𝑎𝐮−𝑏𝐯)+𝑖(𝑏𝐮+𝑎𝐯).

We identify 𝐮∈𝑉 with 𝐮+𝑖𝟎∈𝑉ℂ. For 𝑇∈L(𝑉), its complexification is 𝑇ℂ:𝑉ℂ→𝑉ℂ, 𝑇ℂ(𝐮+𝑖𝐯)=𝑇𝐮+𝑖𝑇𝐯.

The formula for scalar multiplication is what we get by expanding (𝑎+𝑏𝑖)(𝐮+𝑖𝐯) formally with 𝑖2=−1. For 𝑉=ℝ𝑛, the map 𝐮+𝑖𝐯↦𝐮+𝑖𝐯∈ℂ𝑛 (with the right side computed coordinatewise) identifies 𝑉ℂ with ℂ𝑛, and 𝑇ℂ with the complex matrix of 𝑇.

Proposition 7.3.33 (Basic properties). Let 𝑉 be a real vector space and 𝑇∈L(𝑉).

  1. (a)𝑉ℂ is a complex vector space.
  2. (b)If 𝐛1,…,𝐛𝑛 is a basis of 𝑉 over ℝ, it is a basis of 𝑉ℂ over ℂ. In particular dimℂ𝑉ℂ=dimℝ𝑉.
  3. (c)𝑇ℂ is ℂ-linear, and [𝑇ℂ]ℬ︀ℬ︀=[𝑇]ℬ︀ℬ︀ for every basis ℬ︀ of 𝑉. In particular 𝑝𝑇ℂ=𝑝𝑇 when 𝑉 is finite-dimensional and nonzero.

Proof.standard

Uses. Definition 7.3.32 (Complexification), Definition 1.2.1 (Vector space), Definition 1.6.1 (Basis), Definition 2.3.1 (Matrix of a linear map), Definition 7.3.12 (Characteristic polynomial of an operator)

Hint 1. For (a), most axioms hold coordinatewise. The interesting one is (𝛼𝛽)𝐳=𝛼(𝛽𝐳); expand both sides with 𝛼=𝑎+𝑏𝑖, 𝛽=𝑐+𝑑𝑖.
Hint 2. For (b), write 𝐮+𝑖𝐯 with 𝐮=∑𝛼𝑘𝐛𝑘 and 𝐯=∑𝛽𝑘𝐛𝑘, and compute ∑(𝛼𝑘+𝑖𝛽𝑘)𝐛𝑘.
Hint 3. For (c), check 𝑇ℂ((𝑎+𝑏𝑖)𝐳)=(𝑎+𝑏𝑖)𝑇ℂ(𝐳) directly, and compute 𝑇ℂ𝐛𝑗.
Sketch. (a) Addition is that of 𝑉×𝑉; scalar multiplication satisfies 1𝐳=𝐳, both distributive laws (by linearity of the formula in (𝑎,𝑏) and in (𝐮,𝐯)), and associativity, which reduces to the multiplication rule of ℂ. (b) ∑(𝛼𝑘+𝑖𝛽𝑘)𝐛𝑘=∑𝛼𝑘𝐛𝑘+𝑖∑𝛽𝑘𝐛𝑘, so spanning and independence over ℂ reduce to the same statements over ℝ for the real and imaginary parts. (c) Expanding gives linearity; 𝑇ℂ𝐛𝑗=𝑇𝐛𝑗 has the same real coordinates.

(a) 𝑉×𝑉 is a real vector space under coordinatewise operations, so addition in 𝑉ℂ satisfies the additive axioms. For scalar multiplication, 1(𝐮+𝑖𝐯)=𝐮+𝑖𝐯 is the case 𝑎=1, 𝑏=0. The formula is additive in 𝑎+𝑏𝑖 and in 𝐮+𝑖𝐯 separately, which gives the two distributive laws. For associativity, let 𝛼=𝑎+𝑏𝑖 and 𝛽=𝑐+𝑑𝑖. Then

𝛼(𝛽(𝐮+𝑖𝐯))=𝛼((𝑐𝐮−𝑑𝐯)+𝑖(𝑑𝐮+𝑐𝐯))=(𝑎(𝑐𝐮−𝑑𝐯)−𝑏(𝑑𝐮+𝑐𝐯))+𝑖(𝑏(𝑐𝐮−𝑑𝐯)+𝑎(𝑑𝐮+𝑐𝐯))=((𝑎𝑐−𝑏𝑑)𝐮−(𝑎𝑑+𝑏𝑐)𝐯)+𝑖((𝑎𝑑+𝑏𝑐)𝐮+(𝑎𝑐−𝑏𝑑)𝐯),

which is (𝛼𝛽)(𝐮+𝑖𝐯) since 𝛼𝛽=(𝑎𝑐−𝑏𝑑)+(𝑎𝑑+𝑏𝑐)𝑖.

(b) For 𝛼𝑘,𝛽𝑘∈ℝ, (𝛼𝑘+𝑖𝛽𝑘)𝐛𝑘=(𝛼𝑘+𝑖𝛽𝑘)(𝐛𝑘+𝑖𝟎)=𝛼𝑘𝐛𝑘+𝑖𝛽𝑘𝐛𝑘, so

∑𝑘(𝛼𝑘+𝑖𝛽𝑘)𝐛𝑘=∑𝑘𝛼𝑘𝐛𝑘+𝑖∑𝑘𝛽𝑘𝐛𝑘.
(7.3.2)

Every 𝐮+𝑖𝐯∈𝑉ℂ is of this form, by expanding 𝐮 and 𝐯 in the basis, so the 𝐛𝑘 span 𝑉ℂ. If the left side of (7.3.2) is 𝟎, then ∑𝛼𝑘𝐛𝑘=𝟎 and ∑𝛽𝑘𝐛𝑘=𝟎, so all 𝛼𝑘 and 𝛽𝑘 vanish. Hence the 𝐛𝑘 are independent over ℂ.

(c) 𝑇ℂ is additive because 𝑇 is. For 𝑎+𝑏𝑖∈ℂ,

𝑇ℂ((𝑎+𝑏𝑖)(𝐮+𝑖𝐯))=𝑇(𝑎𝐮−𝑏𝐯)+𝑖𝑇(𝑏𝐮+𝑎𝐯)=(𝑎𝑇𝐮−𝑏𝑇𝐯)+𝑖(𝑏𝑇𝐮+𝑎𝑇𝐯),

which is (𝑎+𝑏𝑖)(𝑇𝐮+𝑖𝑇𝐯). Finally, 𝑇ℂ𝐛𝑗=𝑇𝐛𝑗+𝑖𝟎=∑𝑘𝑡𝑘𝑗𝐛𝑘, where [𝑇]ℬ︀ℬ︀=(𝑡𝑘𝑗), so 𝑇ℂ has the same matrix (Definition 2.3.1) and the same characteristic polynomial. □

Conjugation 𝐮+𝑖𝐯=𝐮−𝑖𝐯 on 𝑉ℂ commutes with 𝑇ℂ, and this pairs up the eigenvectors of 𝑇ℂ.

Proposition 7.3.34 (Eigenvalues of a complexification). Let 𝑉≠{𝟎} be a finite-dimensional real vector space and 𝑇∈L(𝑉).

  1. (a)A real number 𝜆 is an eigenvalue of 𝑇 if and only if it is an eigenvalue of 𝑇ℂ.
  2. (b)If 𝑇ℂ(𝐮+𝑖𝐯)=𝜆(𝐮+𝑖𝐯), then 𝑇ℂ(𝐮−𝑖𝐯)=𝜆(𝐮−𝑖𝐯). So the nonreal eigenvalues of 𝑇ℂ come in conjugate pairs.
  3. (c)If 𝜆=𝑎+𝑏𝑖 with 𝑏≠0 and 𝐮+𝑖𝐯 is an eigenvector of 𝑇ℂ for 𝜆, then 𝐮,𝐯 is linearly independent, 𝑈=span(𝐮,𝐯) is 𝑇-invariant, and the matrix of 𝑇↾𝑈 in the basis 𝐮,𝐯 is (𝑎𝑏−𝑏𝑎).

Proof.standard

Uses. Proposition 7.3.33 (Basic properties), Corollary 7.3.13 (Eigenvalues of an operator), Proposition 7.2.4 (Invariance is checked on a spanning list), Proposition 7.2.7 (Invariant lines are spanned by eigenvectors), Proposition 1.2.35 (Zero products)

Techniques. Complexify, then take real and imaginary parts (Technique 7.3.35)

Hint 1. For (a), compare characteristic polynomials.
Hint 2. For (b) and (c), write out 𝑇ℂ(𝐮+𝑖𝐯)=(𝑎+𝑏𝑖)(𝐮+𝑖𝐯) and compare the parts in front of 𝑖 and outside it.
Hint 3. You should get 𝑇𝐮=𝑎𝐮−𝑏𝐯 and 𝑇𝐯=𝑏𝐮+𝑎𝐯. If 𝐮,𝐯 were dependent, they would lie on a 𝑇-invariant line.
Sketch. (a) 𝑝𝑇ℂ=𝑝𝑇. (b), (c) Comparing real and imaginary parts, 𝑇𝐮=𝑎𝐮−𝑏𝐯 and 𝑇𝐯=𝑏𝐮+𝑎𝐯, which gives both the conjugate eigenvector and the invariant plane with matrix (𝑎𝑏−𝑏𝑎). If 𝐮,𝐯 were dependent, they would lie on an invariant line on which 𝑇 acts by a real 𝜇, forcing 𝜆=𝜇.

(a) By Proposition 7.3.33(c), 𝑝𝑇ℂ=𝑝𝑇. For real 𝜆, both statements say 𝑝𝑇(𝜆)=0, by Corollary 7.3.13 applied to 𝑇 and to 𝑇ℂ.

(b) Let 𝜆=𝑎+𝑏𝑖. The equation 𝑇𝐮+𝑖𝑇𝐯=(𝑎𝐮−𝑏𝐯)+𝑖(𝑏𝐮+𝑎𝐯) is an equation in 𝑉×𝑉, so

𝑇𝐮=𝑎𝐮−𝑏𝐯,𝑇𝐯=𝑏𝐮+𝑎𝐯.
(7.3.3)

Then 𝑇ℂ(𝐮−𝑖𝐯)=(𝑎𝐮−𝑏𝐯)−𝑖(𝑏𝐮+𝑎𝐯), and the definition of scalar multiplication gives the same vector for (𝑎−𝑏𝑖)(𝐮+𝑖(−𝐯)). Since 𝐮+𝑖𝐯≠𝟎, also 𝐮−𝑖𝐯≠𝟎, so 𝜆 is an eigenvalue whenever 𝜆 is.

(c) By (7.3.3) and Proposition 7.2.4, 𝑈 is 𝑇-invariant, and the coordinates of 𝑇𝐮 and 𝑇𝐯 in 𝐮,𝐯 are (𝑎,−𝑏) and (𝑏,𝑎), giving the matrix (𝑎𝑏−𝑏𝑎), once we know that 𝐮,𝐯 is independent. Suppose it is not. Since 𝐮 and 𝐯 are not both 𝟎, 𝑈 is then a line span(𝐳), and it is 𝑇-invariant, so 𝑇𝐳=𝜇𝐳 for some 𝜇∈ℝ (Proposition 7.2.7). Writing 𝐮=𝛼𝐳 and 𝐯=𝛽𝐳 with 𝛼,𝛽∈ℝ, we get 𝑇ℂ(𝐮+𝑖𝐯)=𝜇(𝐮+𝑖𝐯). Comparing with 𝑇ℂ(𝐮+𝑖𝐯)=𝜆(𝐮+𝑖𝐯) gives (𝜆−𝜇)(𝐮+𝑖𝐯)=𝟎, so 𝜆=𝜇∈ℝ by Proposition 1.2.35, contradicting 𝑏≠0. □

Technique 7.3.35 (Complexify, then take real and imaginary parts). To prove something about an operator 𝑇 on a real vector space, pass to 𝑇ℂ on 𝑉ℂ, where the FTA is available, use a complex result (an eigenvector, later a triangular form), and then split the conclusion into real and imaginary parts. A complex eigenvector 𝐮+𝑖𝐯 for 𝑎+𝑏𝑖 becomes the pair of real equations (7.3.3), and conjugation exchanges 𝜆 with 𝜆.

Reach for it when a statement about real operators is easy over ℂ, and the obstruction over ℝ is only that 𝑝𝑇 may have nonreal roots. The canonical example is the second proof of Theorem 7.3.31.

Proof of Theorem 7.3.31 (via complexification).standard

Uses. Theorem 7.3.27 (Existence of eigenvalues over ℂ), Proposition 7.3.33 (Basic properties), Proposition 7.3.34 (Eigenvalues of a complexification), Proposition 7.2.7 (Invariant lines are spanned by eigenvectors)

Techniques. Complexify, then take real and imaginary parts (Technique 7.3.35)

Hint 1. 𝑇ℂ acts on a nonzero finite-dimensional complex space. What does Theorem 7.3.27 give?
Hint 2. Split into the cases 𝜆∈ℝ and 𝜆∉ℝ, and use Proposition 7.3.34.
Sketch. 𝑇ℂ has an eigenvalue 𝜆. If 𝜆 is real, it is an eigenvalue of 𝑇, and an eigenvector spans an invariant line. Otherwise an eigenvector 𝐮+𝑖𝐯 gives the invariant plane span(𝐮,𝐯) of Proposition 7.3.34(c).
By Proposition 7.3.33(b), 𝑉ℂ is a nonzero finite-dimensional complex vector space, so 𝑇ℂ has an eigenvalue 𝜆 by Theorem 7.3.27. If 𝜆∈ℝ, it is an eigenvalue of 𝑇 by Proposition 7.3.34(a), and any eigenvector spans a 𝑇-invariant line (Proposition 7.2.7). If 𝜆∉ℝ, an eigenvector 𝐮+𝑖𝐯 of 𝑇ℂ gives the two-dimensional 𝑇-invariant subspace span(𝐮,𝐯) of Proposition 7.3.34(c). □

The Krylov proof stays inside 𝑉 and uses the real factorization of polynomials. The complexification proof uses only the complex theorem, but has to leave 𝑉 and come back, and it is the version that generalizes, to the block triangular form of real operators (Section 7.5) and to real normal operators in Chapter 8.

Example 7.3.36 (Decomplexification). The opposite construction forgets structure. A complex vector space is also a real vector space, by restricting scalars to ℝ, and ℂ𝑛 becomes ℝ2𝑛 with real basis 𝐞1,…,𝐞𝑛,𝑖𝐞1,…,𝑖𝐞𝑛. A complex matrix 𝐴=𝑃+𝑖𝑄, with 𝑃,𝑄∈ℝ𝑛×𝑛, acts by

(𝑃+𝑖𝑄)(𝐱+𝑖𝐲)=(𝑃𝐱−𝑄𝐲)+𝑖(𝑄𝐱+𝑃𝐲)(𝐱,𝐲∈ℝ𝑛),

so as a real linear map of ℝ2𝑛 its matrix is (𝑃−𝑄𝑄𝑃). For 𝑛=1, multiplication by 𝑎+𝑏𝑖 on ℂ=ℝ2 becomes (𝑎−𝑏𝑏𝑎), as in Proposition B.1.18. Conversely, decomplexifying 𝑉ℂ gives 𝑉×𝑉, on which 𝑇ℂ acts as (𝑇00𝑇).

Exercises

Exercise 7.3.37.routine For each matrix, find 𝑝𝐴, the eigenvalues, and a basis of each eigenspace.

  1. (a)𝐴=(122212221)∈ℝ3×3.
  2. (b)𝐴=(0100012−54)∈ℝ3×3.
  3. (c)𝐴=(1−10110002), over ℝ and over ℂ.

Solution.

Uses. Algorithm 7.3.7 (Eigenvalues by the characteristic polynomial), Proposition 7.3.2 (Determinants over a commutative ring), Proposition 7.3.15 (Block-triangular matrices), Exercise 7.3.17, Theorem 7.3.24 (Every monic polynomial is a characteristic polynomial)

Hint 1. In (a), add the second and third rows of 𝑥𝐼3−𝐴 to the first; a common factor appears. In (b), 𝐴 is the transpose of a companion matrix. In (c), 𝐴 is block diagonal.
Hint 2. The answers are (𝑥−5)(𝑥+1)2, (𝑥−1)2(𝑥−2), and (𝑥−2)(𝑥2−2𝑥+2). Now row reduce 𝐴−𝜆𝐼3 for each root.
Sketch. (a) 𝑝𝐴=(𝑥−5)(𝑥+1)2, 𝐸5=span((1,1,1)), 𝐸−1=span((1,−1,0),(1,0,−1)). (b) 𝑝𝐴=𝑥3−4𝑥2+5𝑥−2=(𝑥−1)2(𝑥−2), 𝐸1=span((1,1,1)), 𝐸2=span((1,2,4)). (c) 𝑝𝐴=(𝑥−2)(𝑥2−2𝑥+2); over ℝ only 𝐸2=span(𝐞3); over ℂ also 𝐸1+𝑖=span((𝑖,1,0)) and 𝐸1−𝑖=span((−𝑖,1,0)).

(a) In 𝑥𝐼3−𝐴, add rows 2 and 3 to row 1; this does not change the determinant (Proposition 7.3.2(b), (d)). Row 1 becomes (𝑥−5,𝑥−5,𝑥−5), and pulling out the factor 𝑥−5 and then adding 2 times row 1 to rows 2 and 3,

𝑝𝐴=(𝑥−5)det(111−2𝑥−1−2−2−2𝑥−1)=(𝑥−5)det(1110𝑥+1000𝑥+1)=(𝑥−5)(𝑥+1)2.

Next, 𝐴−5𝐼3=(−4222−4222−4) has rank 2 and kills (1,1,1), so 𝐸5=span((1,1,1)). Every row of 𝐴+𝐼3 is (2,2,2), so 𝐸−1={(𝑥,𝑦,𝑧)|𝑥+𝑦+𝑧=0}=span((1,−1,0),(1,0,−1)). As a check, 𝐴=2𝐽−𝐼3 with 𝐽 the all-ones matrix, and 𝐽 acts as 3 on (1,1,1) and as 0 on the plane 𝑥+𝑦+𝑧=0.

(b) 𝐴𝖳=(00210−5014) is the companion matrix of 𝑝=𝑥3−4𝑥2+5𝑥−2, so 𝑝𝐴=𝑝𝐴𝖳=𝑝 by Exercise 7.3.17 and Theorem 7.3.24. (Expanding det(𝑥𝐼3−𝐴) along the first row gives the same.) Since 𝑝(1)=0, dividing gives 𝑝=(𝑥−1)(𝑥2−3𝑥+2)=(𝑥−1)2(𝑥−2). Then

𝐴−𝐼3=(−1100−112−53),𝐴−2𝐼3=(−2100−212−52).

The first two rows of each are independent, so both have rank 2. The first two rows of 𝐴−𝐼3 force 𝑦=𝑥 and 𝑧=𝑦, and (1,1,1) satisfies the third row; so 𝐸1=span((1,1,1)). Likewise 𝐸2=span((1,2,4)). Although 1 is a double root of 𝑝𝐴, its eigenspace is only a line.

(c) 𝐴 is block diagonal with blocks (1−111) and (2), so by Proposition 7.3.15 and Example 7.3.5(a), 𝑝𝐴=(𝑥2−2𝑥+2)(𝑥−2). Over ℝ, 𝑥2−2𝑥+2=(𝑥−1)2+1 has no root, so the only eigenvalue is 2, and 𝐴−2𝐼3=(−1−101−10000) gives 𝐸2=span(𝐞3). Over ℂ, the quadratic has roots 1±𝑖. For 𝜆=1+𝑖, the first row of 𝐴−𝜆𝐼3 is (−𝑖,−1,0), so (𝑖,1,0) spans 𝐸1+𝑖; indeed 𝐴(𝑖,1,0)=(𝑖−1,𝑖+1,0)=(1+𝑖)(𝑖,1,0). Since 𝐴 is real, conjugating 𝐴𝐯=𝜆𝐯 entrywise gives 𝐴𝐯=𝜆𝐯, so 𝐸1−𝑖=span((−𝑖,1,0)), and 𝐸2=span(𝐞3) as before.

Exercise 7.3.38.standard Let 𝐴,𝐵∈𝐹𝑛×𝑛. Prove or disprove:

  1. (a)If 𝑝𝐴=𝑝𝐵, then 𝐴 and 𝐵 are similar.
  2. (b)If 𝑝𝐴=𝑝𝐵, then rank𝐴=rank𝐵.
  3. (c)For every 𝑐∈𝐹, 𝑝𝐴+𝑐𝐼𝑛 is obtained from 𝑝𝐴 by substituting 𝑥−𝑐 for 𝑥, that is, 𝑝𝐴+𝑐𝐼𝑛=𝑝𝐴(𝑥−𝑐).
  4. (d)If 𝑝𝐴=𝑝𝐵 and 𝑝𝐴 has 𝑛 distinct roots in 𝐹, then 𝐴 and 𝐵 are similar.

Solution.

Uses. Proposition 7.3.2 (Determinants over a commutative ring), Proposition 7.1.10 (Evaluation respects sums and products), Theorem 7.3.6 (Eigenvalues are the roots of the characteristic polynomial), Theorem 7.2.21 (Eigenvectors for distinct eigenvalues are independent), Corollary 3.5.25 (The invertible matrix theorem), Proposition 4.5.24 (Similarity is an equivalence relation)

Techniques. Finding counterexamples and minimal examples (Technique 1.4.18), Rule out similarity with an invariant (Technique 4.5.29)

Hint 1. For (a) and (b), look at 2×2 matrices with a repeated eigenvalue.
Hint 2. For (c), 𝑥𝐼𝑛−(𝐴+𝑐𝐼𝑛)=(𝑥−𝑐)𝐼𝑛−𝐴. Which map 𝐹[𝑥]→𝐹[𝑥] turns 𝑥𝐼𝑛−𝐴 into this matrix?
Hint 3. For (d), pick an eigenvector of 𝐴 for each root. They form a basis, in which 𝐴 becomes diag(𝜆1,…,𝜆𝑛). Do the same for 𝐵 separately, in its own basis of eigenvectors.
Sketch. (a) False: 𝐼2 and (1101). (b) False: 0 and (0100). (c) True: substitute 𝑥−𝑐 for 𝑥 entrywise and use Proposition 7.3.2(g). (d) True: 𝑛 distinct eigenvalues give a basis of eigenvectors, so 𝐴 and 𝐵 are both similar to diag(𝜆1,…,𝜆𝑛).

(a) False. 𝐼2 and 𝑁=(1101) both have characteristic polynomial (𝑥−1)2 (Example 7.3.5(d)). But 𝑆−1𝐼2𝑆=𝐼2≠𝑁 for every invertible 𝑆, so they are not similar.

(b) False. 02×2 and (0100) both have characteristic polynomial 𝑥2, but their ranks are 0 and 1.

(c) True. Substituting 𝑥−𝑐 for 𝑥 respects sums, products, and 1 (Proposition 7.1.10), and applied to every entry of 𝑥𝐼𝑛−𝐴 it gives (𝑥−𝑐)𝐼𝑛−𝐴=𝑥𝐼𝑛−(𝐴+𝑐𝐼𝑛), since the constant entries are unchanged. By Proposition 7.3.2(g), 𝑝𝐴+𝑐𝐼𝑛=𝑝𝐴(𝑥−𝑐). In particular the eigenvalues of 𝐴+𝑐𝐼𝑛 are the 𝜆+𝑐 for the eigenvalues 𝜆 of 𝐴.

(d) True. Let 𝜆1,…,𝜆𝑛 be the distinct roots. By Theorem 7.3.6 each 𝜆𝑖 is an eigenvalue of 𝐴; choose eigenvectors 𝐯𝑖. By Theorem 7.2.21 they are independent, so they form a basis of 𝐹𝑛, and the matrix 𝑆 with columns 𝐯𝑖 is invertible (Corollary 3.5.25). Column 𝑖 of 𝐴𝑆 is 𝐴𝐯𝑖=𝜆𝑖𝐯𝑖, which is column 𝑖 of 𝑆𝐷 for 𝐷=diag(𝜆1,…,𝜆𝑛). So 𝑆−1𝐴𝑆=𝐷. The same argument for 𝐵 gives 𝑅−1𝐵𝑅=𝐷. Hence 𝐴∼𝐷∼𝐵, and 𝐴∼𝐵 by Proposition 4.5.24.

Exercise 7.3.39.standard
  1. (a)Let 𝜏:𝐹𝑛×𝑛→𝐹𝑛×𝑛 be transposition, 𝜏(𝑋)=𝑋𝖳. Show that

    𝑝𝜏=(𝑥−1)𝑛(𝑛+1)/2(𝑥+1)𝑛(𝑛−1)/2.
  2. (b)Let 𝑊 be the real vector space of functions ℝ→ℝ with basis sin𝑡,cos𝑡,𝑡sin𝑡,𝑡cos𝑡. Find the characteristic polynomial of differentiation 𝐷:𝑊→𝑊, and conclude that no nonzero function in 𝑊 satisfies 𝑓′=𝜆𝑓 with 𝜆∈ℝ.

Solution.

Uses. Lemma 2.4.24 (Matrix units), Definition 2.3.1 (Matrix of a linear map), Proposition 7.3.15 (Block-triangular matrices), Corollary 7.3.13 (Eigenvalues of an operator)

Hint 1. For (a), order the matrix units so that 𝑈𝑖𝑗 and 𝑈𝑗𝑖 are adjacent. What does the matrix of 𝜏 look like?
Hint 2. It is block diagonal, with 𝑛 blocks (1) and 𝑛(𝑛−1)/2 blocks (0110). For (b), compute 𝐷 on each basis function; the matrix is block upper triangular.
Sketch. (a) 𝜏(𝑈𝑖𝑖)=𝑈𝑖𝑖 and 𝜏 swaps 𝑈𝑖𝑗 and 𝑈𝑗𝑖; so 𝑝𝜏=(𝑥−1)𝑛(𝑥2−1)𝑛(𝑛−1)/2. (b) The matrix is (𝑅𝐼20𝑅) with 𝑅=(0−110), so 𝑝𝐷=(𝑥2+1)2, which has no real root.

(a) By Lemma 2.4.24, the matrix units 𝑈𝑖𝑗 form a basis of 𝐹𝑛×𝑛. Order it as 𝑈11,…,𝑈𝑛𝑛, followed by the pairs 𝑈𝑖𝑗,𝑈𝑗𝑖 for 𝑖<𝑗. Since 𝜏(𝑈𝑖𝑖)=𝑈𝑖𝑖, 𝜏(𝑈𝑖𝑗)=𝑈𝑗𝑖, and 𝜏(𝑈𝑗𝑖)=𝑈𝑖𝑗, the matrix of 𝜏 is block diagonal with 𝑛 blocks (1) and 𝑛(𝑛−1)/2 blocks (0110), whose characteristic polynomial is 𝑥2−1=(𝑥−1)(𝑥+1). Applying Proposition 7.3.15 once for each block,

𝑝𝜏=(𝑥−1)𝑛(𝑥2−1)𝑛(𝑛−1)/2=(𝑥−1)𝑛+𝑛(𝑛−1)/2(𝑥+1)𝑛(𝑛−1)/2,

and 𝑛+𝑛(𝑛−1)/2=𝑛(𝑛+1)/2. This fits Example 7.2.17: 𝜏2=id, the symmetric matrices form 𝐸1, of dimension 𝑛(𝑛+1)/2, and the antisymmetric ones form 𝐸−1. When 1+1=0 in 𝐹, the formula reads 𝑝𝜏=(𝑥−1)𝑛2.

(b) Let 𝐛1=sin𝑡, 𝐛2=cos𝑡, 𝐛3=𝑡sin𝑡, 𝐛4=𝑡cos𝑡. Then

𝐷𝐛1=𝐛2,𝐷𝐛2=−𝐛1,𝐷𝐛3=𝐛1+𝐛4,𝐷𝐛4=𝐛2−𝐛3,

so 𝐷 maps 𝑊 to 𝑊, and its matrix (Definition 2.3.1) is

[𝐷]=(0−1101001000−10010).

This is block upper triangular with both diagonal blocks equal to (0−110), so 𝑝𝐷=(𝑥2+1)2 by Proposition 7.3.15. It has no real root, so 𝐷 has no real eigenvalue (Corollary 7.3.13): no nonzero 𝑓∈𝑊 satisfies 𝑓′=𝜆𝑓 with 𝜆∈ℝ. The upper left block reflects the invariant plane span(sin𝑡,cos𝑡).

Exercise 7.3.40.standard Let 𝐴∈𝐹𝑚×𝑛 and 𝐵∈𝐹𝑛×𝑚, where 𝑚,𝑛≥1.

  1. (a)Show that 𝑥𝑛𝑝𝐴𝐵=𝑥𝑚𝑝𝐵𝐴.
  2. (b)Deduce that 𝐴𝐵 and 𝐵𝐴 have the same nonzero eigenvalues.
  3. (c)Prove (b) again without characteristic polynomials: if 𝐴𝐵𝐯=𝜆𝐯 with 𝐯≠𝟎 and 𝜆≠0, show that 𝐵𝐯 is an eigenvector of 𝐵𝐴.
  4. (d)Give an example in which 0 is an eigenvalue of 𝐵𝐴 but not of 𝐴𝐵.

Solution.

Uses. Proposition 7.3.22 (𝐴𝐵 and 𝐵𝐴), Proposition 7.3.11 (Similar matrices have the same characteristic polynomial), Proposition 7.3.15 (Block-triangular matrices), Theorem 2.4.34 (Block multiplication), Theorem 7.3.6 (Eigenvalues are the roots of the characteristic polynomial)

Techniques. Compare 𝐴𝐵 and 𝐵𝐴 (Technique 7.3.21)

Hint 1. Imitate the first proof of Proposition 7.3.22 with 𝑆=(𝐼𝑚𝐴0𝐼𝑛). The zero blocks now have different sizes.
Hint 2. For (b), evaluate the identity of (a) at 𝜆≠0.
Hint 3. For (c), why is 𝐵𝐯≠𝟎? Apply 𝐴 to it.
Sketch. (a) 𝑋=(𝐴𝐵0𝐵0𝑛×𝑛) and 𝑌=(0𝑚×𝑚0𝐵𝐵𝐴) satisfy 𝑋𝑆=𝑆𝑌 with 𝑆=(𝐼𝑚𝐴0𝐼𝑛), and 𝑝𝑋=𝑝𝐴𝐵𝑥𝑛, 𝑝𝑌=𝑥𝑚𝑝𝐵𝐴. (b) For 𝜆≠0, 𝜆𝑛𝑝𝐴𝐵(𝜆)=𝜆𝑚𝑝𝐵𝐴(𝜆). (c) 𝐴(𝐵𝐯)=𝜆𝐯≠𝟎, and 𝐵𝐴(𝐵𝐯)=𝐵(𝜆𝐯). (d) 𝐴=(1,0), 𝐵=(1,0)𝖳.

(a) Let 𝑆=(𝐼𝑚𝐴0𝐼𝑛), which is invertible with inverse (𝐼𝑚−𝐴0𝐼𝑛), and let

𝑋=(𝐴𝐵0𝐵0𝑛×𝑛),𝑌=(0𝑚×𝑚0𝐵𝐵𝐴).

By block multiplication (Theorem 2.4.34), 𝑋𝑆=(𝐴𝐵𝐴𝐵𝐴𝐵𝐵𝐴)=𝑆𝑌. So 𝑋 and 𝑌 are similar and 𝑝𝑋=𝑝𝑌 (Proposition 7.3.11). Both are block lower triangular, so by Proposition 7.3.15, 𝑝𝑋=𝑝𝐴𝐵𝑝0𝑛×𝑛=𝑝𝐴𝐵𝑥𝑛 and 𝑝𝑌=𝑥𝑚𝑝𝐵𝐴.

(b) Let 𝜆≠0. Evaluating (a) at 𝜆 gives 𝜆𝑛𝑝𝐴𝐵(𝜆)=𝜆𝑚𝑝𝐵𝐴(𝜆), and 𝜆𝑛,𝜆𝑚≠0. So 𝑝𝐴𝐵(𝜆)=0 if and only if 𝑝𝐵𝐴(𝜆)=0, and Theorem 7.3.6 finishes.

(c) Suppose 𝐴𝐵𝐯=𝜆𝐯 with 𝐯≠𝟎 and 𝜆≠0. Then 𝐴(𝐵𝐯)=𝜆𝐯≠𝟎, so 𝐵𝐯≠𝟎, and (𝐵𝐴)(𝐵𝐯)=𝐵(𝐴𝐵𝐯)=𝜆𝐵𝐯. So every nonzero eigenvalue of 𝐴𝐵 is one of 𝐵𝐴, and exchanging the roles of 𝐴 and 𝐵 gives the converse.

(d) Let 𝐴=(10)∈𝐹1×2 and 𝐵=(10)∈𝐹2×1. Then 𝐴𝐵=(1) has the single eigenvalue 1, while 𝐵𝐴=(1000) also has the eigenvalue 0. Consistently with (a), 𝑥2(𝑥−1)=𝑥⋅𝑥(𝑥−1).

Exercise 7.3.41.standard Let 𝑛≥2 and 𝐮,𝐯∈𝐹𝑛.

  1. (a)Show that 𝑝𝐮𝐯𝖳=𝑥𝑛−1(𝑥−𝐯𝖳𝐮).
  2. (b)Let 𝐽∈𝐹𝑛×𝑛 be the matrix whose entries are all 1. Find 𝑝𝐽 and the eigenvalues of 𝐽 over ℚ, and over 𝔽𝑝 for a prime 𝑝 dividing 𝑛.

Solution.

Uses. Exercise 7.3.40, Corollary 7.1.7 (No zero divisors, cancellation, and units), Proposition 4.5.35 (Preserved subspaces give block triangular matrices), Proposition 4.5.27 (Similarity invariants), Theorem 4.2.11 (Rank–nullity), Example 7.3.5 (Small cases)

Techniques. Compare 𝐴𝐵 and 𝐵𝐴 (Technique 7.3.21), Choose a basis adapted to a preserved subspace (Technique 4.5.36)

Hint 1. Write 𝐮𝐯𝖳 as a product of an 𝑛×1 and a 1×𝑛 matrix, and use Exercise 7.3.40.
Hint 2. Alternatively, if 𝐯≠𝟎, the hyperplane {𝐰|𝐯𝖳𝐰=0} lies in the kernel. Choose a basis starting with a basis of it; the matrix becomes triangular. Its last diagonal entry is the trace.
Sketch. (a) With 𝐴=𝐮∈𝐹𝑛×1 and 𝐵=𝐯𝖳∈𝐹1×𝑛, Exercise 7.3.40 gives 𝑥𝑝𝐮𝐯𝖳=𝑥𝑛𝑝𝐯𝖳𝐮=𝑥𝑛(𝑥−𝐯𝖳𝐮); cancel 𝑥. (b) 𝐽=𝟏𝟏𝖳 and 𝟏𝖳𝟏=𝑛, so 𝑝𝐽=𝑥𝑛−1(𝑥−𝑛): eigenvalues 0,𝑛 over ℚ, and only 0 over 𝔽𝑝 when 𝑝∣𝑛.

(a) First proof. Let 𝐴=𝐮∈𝐹𝑛×1 and 𝐵=𝐯𝖳∈𝐹1×𝑛, so 𝐴𝐵=𝐮𝐯𝖳 and 𝐵𝐴=(𝐯𝖳𝐮) is 1×1, with 𝑝𝐵𝐴=𝑥−𝐯𝖳𝐮. By Exercise 7.3.40, 𝑥𝑝𝐮𝐯𝖳=𝑥𝑛(𝑥−𝐯𝖳𝐮). Cancelling 𝑥 in 𝐹[𝑥] (Corollary 7.1.7) gives the claim.

Second proof. Let 𝑀=𝐮𝐯𝖳. If 𝐯=𝟎, then 𝑀=0 and 𝑝𝑀=𝑥𝑛, as claimed. Otherwise 𝐻={𝐰∈𝐹𝑛|𝐯𝖳𝐰=0} is the kernel of a nonzero linear map 𝐹𝑛→𝐹, so dim𝐻=𝑛−1 (Theorem 4.2.11), and 𝑀𝐰=𝐮(𝐯𝖳𝐰)=𝟎 for 𝐰∈𝐻. Extend a basis 𝐰1,…,𝐰𝑛−1 of 𝐻 to a basis ℬ︀ of 𝐹𝑛. The first 𝑛−1 columns of [𝑇𝑀]ℬ︀ℬ︀ are 𝟎, so this matrix is upper triangular with diagonal entries 0,…,0,𝑐, and 𝑝𝑀=𝑥𝑛−1(𝑥−𝑐) by Example 7.3.5(d). Since the trace is a similarity invariant (Proposition 4.5.27), 𝑐=tr𝑀=∑𝑖𝑢𝑖𝑣𝑖=𝐯𝖳𝐮.

(b) 𝐽=𝟏𝟏𝖳 with 𝟏=(1,…,1), and 𝟏𝖳𝟏=𝑛. By (a), 𝑝𝐽=𝑥𝑛−1(𝑥−𝑛). Over ℚ the eigenvalues are 0 and 𝑛. Over 𝔽𝑝 with 𝑝∣𝑛, we have 𝑛=0, so 𝑝𝐽=𝑥𝑛 and 0 is the only eigenvalue; indeed 𝐽𝟏=𝑛𝟏=𝟎 there.

Exercise 7.3.42.standard Let 𝐴∈ℚ𝑛×𝑛 have integer entries.

  1. (a)Show that all coefficients of 𝑝𝐴 are integers.
  2. (b)Show that every rational eigenvalue of 𝐴 is an integer that divides det𝐴.
  3. (c)Show that (1234) has no rational eigenvalue.

Solution.

Uses. Definition 7.3.1 (Determinant over a commutative ring), Definition 7.1.1 (Product of polynomials), Proposition 7.3.18 (The first and last coefficients), Theorem 7.3.6 (Eigenvalues are the roots of the characteristic polynomial)

Hint 1. For (a), every entry of 𝑥𝐼𝑛−𝐴 is a polynomial with integer coefficients. What happens to such polynomials under sums and products?
Hint 2. For (b), write the eigenvalue as 𝑎/𝑏 in lowest terms with 𝑏≥1, substitute into 𝑝𝐴, and multiply by 𝑏𝑛. Which terms are divisible by 𝑏?
Sketch. (a) The Leibniz formula is a signed sum of products of polynomials with integer coefficients. (b) If 𝑝𝐴(𝑎/𝑏)=0 with gcd(𝑎,𝑏)=1, then 𝑎𝑛=−𝑏(𝑐𝑛−1𝑎𝑛−1+⋯+𝑐0𝑏𝑛−1), so every prime factor of 𝑏 divides 𝑎; hence 𝑏=1. Then 𝑎 divides the constant term (−1)𝑛det𝐴. (c) 𝑝𝐴=𝑥2−5𝑥−2 has no root among ±1,±2.

(a) Every entry of 𝑥𝐼𝑛−𝐴 is 𝑥−𝑎𝑖𝑖 or −𝑎𝑖𝑗, a polynomial with integer coefficients. By the formula for the product (Definition 7.1.1), sums and products of polynomials with integer coefficients have integer coefficients. So every Leibniz term, and hence 𝑝𝐴, has integer coefficients.

(b) Write 𝑝𝐴=𝑥𝑛+𝑐𝑛−1𝑥𝑛−1+⋯+𝑐0 with 𝑐𝑖∈ℤ, and suppose 𝑝𝐴(𝑎/𝑏)=0 with 𝑎,𝑏∈ℤ, 𝑏≥1, and gcd(𝑎,𝑏)=1. Multiplying by 𝑏𝑛,

𝑎𝑛=−𝑏(𝑐𝑛−1𝑎𝑛−1+𝑐𝑛−2𝑎𝑛−2𝑏+⋯+𝑐0𝑏𝑛−1).

If 𝑏>1, a prime 𝑞 dividing 𝑏 divides 𝑎𝑛, hence divides 𝑎, contradicting gcd(𝑎,𝑏)=1. So 𝑏=1 and the eigenvalue 𝑎 is an integer. Then 0=𝑎𝑛+𝑐𝑛−1𝑎𝑛−1+⋯+𝑐1𝑎+𝑐0 gives 𝑐0=−𝑎(𝑎𝑛−1+𝑐𝑛−1𝑎𝑛−2+⋯+𝑐1), so 𝑎 divides 𝑐0=(−1)𝑛det𝐴 (Proposition 7.3.18), and therefore divides det𝐴.

(c) Here 𝑝𝐴=𝑥2−5𝑥−2 and det𝐴=−2. By (b), a rational eigenvalue would be one of ±1,±2, but 𝑝𝐴(1)=−6, 𝑝𝐴(−1)=4, 𝑝𝐴(2)=−8, and 𝑝𝐴(−2)=12. So there is none, and the eigenvalues (5±33)/2 are irrational.

Exercise 7.3.43.standard Let 𝐴∈ℝ𝑛×𝑛 satisfy 𝐴2=−𝐼𝑛.

  1. (a)Show that 𝑛 is even, using determinants.
  2. (b)Give a second proof using 𝑝𝐴 and the eigenvalues of 𝐴 regarded as a complex matrix.
  3. (c)Show that, conversely, 𝐴2=−𝐼𝑛 has a solution in ℝ𝑛×𝑛 for every even 𝑛, and in ℂ𝑛×𝑛 for every 𝑛.

Solution.

Uses. Theorem 5.4.10 (Multiplicativity), Proposition 7.2.28 (Polynomials act on eigenvectors by their values), Corollary 7.1.35 (Factorization over ℂ), Theorem 7.3.6 (Eigenvalues are the roots of the characteristic polynomial), Proposition 7.1.38 (Conjugate roots)

Hint 1. For (a), take determinants of both sides. What sign does a real square have?
Hint 2. For (b), an eigenvalue 𝜆 of 𝐴 over ℂ satisfies 𝜆2=−1. Then 𝑝𝐴=(𝑥−𝑖)𝑎(𝑥+𝑖)𝑏. Use that 𝑝𝐴 has real coefficients.
Sketch. (a) (det𝐴)2=det(−𝐼𝑛)=(−1)𝑛, and a nonzero real square is positive, so (−1)𝑛=1. (b) All complex eigenvalues are ±𝑖, so 𝑝𝐴=(𝑥−𝑖)𝑎(𝑥+𝑖)𝑏 with 𝑎+𝑏=𝑛; real coefficients force 𝑎=𝑏 (Proposition 7.1.38). (c) Block-diagonal rotations by 𝜋/2 over ℝ; 𝑖𝐼𝑛 over ℂ.

(a) By Theorem 5.4.10, (det𝐴)2=det(𝐴2)=det(−𝐼𝑛)=(−1)𝑛. Since det𝐴∈ℝ and (det𝐴)2=(−1)𝑛≠0, we get (−1)𝑛=(det𝐴)2>0, so 𝑛 is even.

(b) Regard 𝐴 as a complex matrix. If 𝐴𝐯=𝜆𝐯 with 𝐯≠𝟎, then −𝐯=𝐴2𝐯=𝜆2𝐯 (Proposition 7.2.28), so 𝜆2=−1 and 𝜆=±𝑖. By Corollary 7.1.35, 𝑝𝐴 is a product of 𝑛 factors 𝑥−𝜆𝑘, and each 𝜆𝑘 is a root, hence an eigenvalue (Theorem 7.3.6), hence ±𝑖. So 𝑝𝐴=(𝑥−𝑖)𝑎(𝑥+𝑖)𝑏 with 𝑎+𝑏=𝑛. The entries of 𝐴 are real, so 𝑝𝐴∈ℝ[𝑥], and by Proposition 7.1.38 the roots 𝑖 and −𝑖=𝑖 have the same multiplicity: 𝑎=𝑏. Hence 𝑛=2𝑎 is even.

(c) Over ℂ, (𝑖𝐼𝑛)2=−𝐼𝑛. Over ℝ with 𝑛=2𝑚, the block-diagonal matrix with 𝑚 blocks 𝑅=(0−110) works, since 𝑅2=−𝐼2.

Exercise 7.3.44.standard Let 𝐹 be a finite field with elements 𝑐1,…,𝑐𝑞.

  1. (a)Show that 𝑝=(𝑥−𝑐1)(𝑥−𝑐2)⋯(𝑥−𝑐𝑞)+1 has no root in 𝐹.
  2. (b)Deduce that some operator on 𝐹𝑞 has no eigenvalue. So Theorem 7.3.27 fails over every finite field.
  3. (c)Write down this operator for 𝐹=𝔽2 and check directly that it has no eigenvalue.

Solution.

Uses. Theorem 7.3.24 (Every monic polynomial is a characteristic polynomial), Theorem 7.3.6 (Eigenvalues are the roots of the characteristic polynomial), Definition 7.3.23 (Companion matrix)

Hint 1. Evaluate 𝑝 at each 𝑐𝑗.
Hint 2. Every monic polynomial is a characteristic polynomial (Theorem 7.3.24).
Sketch. (a) 𝑝(𝑐𝑗)=0+1=1. (b) 𝐶𝑝∈𝐹𝑞×𝑞 has 𝑝𝐶𝑝=𝑝, which has no root. (c) Over 𝔽2, 𝑝=𝑥(𝑥−1)+1=𝑥2+𝑥+1 and 𝐶𝑝=(0111); neither 𝐶𝑝 nor 𝐶𝑝−𝐼2 is singular.

(a) For each 𝑗, the product (𝑐𝑗−𝑐1)⋯(𝑐𝑗−𝑐𝑞) contains the factor 𝑐𝑗−𝑐𝑗=0, so 𝑝(𝑐𝑗)=0+1=1≠0.

(b) 𝑝 is monic of degree 𝑞, so by Theorem 7.3.24 the companion matrix 𝐶𝑝∈𝐹𝑞×𝑞 has 𝑝𝐶𝑝=𝑝. By (a) and Theorem 7.3.6, 𝐶𝑝 has no eigenvalue in 𝐹.

(c) Over 𝔽2, 𝑝=𝑥(𝑥−1)+1=𝑥2−𝑥+1=𝑥2+𝑥+1, and 𝐶𝑝=(0−11−1)=(0111). Its determinant is −1=1≠0, so 0 is not an eigenvalue, and 𝐶𝑝−𝐼2=(1110) has determinant −1=1≠0, so 1 is not one either.

Exercise 7.3.45.standard Let 𝑝=𝑥𝑛+𝑎𝑛−1𝑥𝑛−1+⋯+𝑎0∈𝐹[𝑥] with 𝑛≥2, and let 𝐶=𝐶𝑝.

  1. (a)Show that rank(𝐶−𝜆𝐼𝑛)≥𝑛−1 for every 𝜆∈𝐹. So every eigenspace of 𝐶 is a line.
  2. (b)Let 𝜆 be a root of 𝑝 and write 𝑝=(𝑥−𝜆)𝑞 with 𝑞=𝑥𝑛−1+𝑏𝑛−2𝑥𝑛−2+⋯+𝑏0. Show that (𝑏0,𝑏1,…,𝑏𝑛−2,1) is an eigenvector of 𝐶 for 𝜆.
  3. (c)Show that 𝐹𝑛 has a basis of eigenvectors of 𝐶 if and only if 𝑝 has 𝑛 distinct roots in 𝐹.

Solution.

Uses. Definition 7.3.23 (Companion matrix), Theorem 4.2.11 (Rank–nullity), Remark 7.1.11 (Polynomials in a matrix), Theorem 7.3.24 (Every monic polynomial is a characteristic polynomial), Theorem 7.3.6 (Eigenvalues are the roots of the characteristic polynomial), Theorem 7.2.21 (Eigenvectors for distinct eigenvalues are independent), Proposition 7.1.22 (Counting roots)

Techniques. Krylov sequences must become dependent (Technique 7.3.28)

Hint 1. For (a), look at the first 𝑛−1 columns of 𝐶−𝜆𝐼𝑛.
Hint 2. For (b), 𝐶𝑗𝐞1=𝐞𝑗+1 for 𝑗≤𝑛−1. Compute 𝑞(𝐶)𝐞1 and 𝑝(𝐶)𝐞1.
Hint 3. For (c), a basis of eigenvectors has at most one vector in each eigenspace, by (a).
Sketch. (a) Columns 𝑘<𝑛 of 𝐶−𝜆𝐼𝑛 are 𝐞𝑘+1−𝜆𝐞𝑘, which are independent (look at the last coordinate first). (b) 𝑞(𝐶)𝐞1 is the given vector, and (𝐶−𝜆𝐼𝑛)𝑞(𝐶)𝐞1=𝑝(𝐶)𝐞1=𝐶𝐞𝑛+∑𝑎𝑗𝐞𝑗+1=𝟎. (c) Distinct roots give independent eigenvectors; conversely a basis of eigenvectors needs 𝑛 different eigenvalues, all roots of 𝑝.

(a) For 𝑘≤𝑛−1, column 𝑘 of 𝐶−𝜆𝐼𝑛 is 𝐞𝑘+1−𝜆𝐞𝑘. Suppose ∑𝑘=1𝑛−1𝑐𝑘(𝐞𝑘+1−𝜆𝐞𝑘)=𝟎. The coordinate at 𝐞𝑛 is 𝑐𝑛−1, so 𝑐𝑛−1=0. If 𝑐𝑘+1=⋯=𝑐𝑛−1=0, the coordinate at 𝐞𝑘+1 is 𝑐𝑘−𝜆𝑐𝑘+1=𝑐𝑘, so 𝑐𝑘=0. Hence these 𝑛−1 columns are independent, rank(𝐶−𝜆𝐼𝑛)≥𝑛−1, and dim𝐸𝜆=𝑛−rank(𝐶−𝜆𝐼𝑛)≤1 by Theorem 4.2.11.

(b) Since 𝐶𝐞𝑘=𝐞𝑘+1 for 𝑘<𝑛, we have 𝐶𝑗𝐞1=𝐞𝑗+1 for 0≤𝑗≤𝑛−1. So

𝐰=𝑞(𝐶)𝐞1=𝑏0𝐞1+𝑏1𝐞2+⋯+𝑏𝑛−2𝐞𝑛−1+𝐞𝑛=(𝑏0,…,𝑏𝑛−2,1),

which is nonzero. Also 𝐶𝑛𝐞1=𝐶𝐞𝑛=−∑𝑗=0𝑛−1𝑎𝑗𝐞𝑗+1, so 𝑝(𝐶)𝐞1=𝐶𝑛𝐞1+∑𝑗𝑎𝑗𝐞𝑗+1=𝟎. By Remark 7.1.11, (𝐶−𝜆𝐼𝑛)𝐰=((𝑥−𝜆)𝑞)(𝐶)𝐞1=𝑝(𝐶)𝐞1=𝟎.

(c) By Theorem 7.3.24 and Theorem 7.3.6, the eigenvalues of 𝐶 are the roots of 𝑝 in 𝐹. If there are 𝑛 distinct ones, eigenvectors for them are independent (Theorem 7.2.21), hence a basis of 𝐹𝑛. Conversely, suppose 𝐯1,…,𝐯𝑛 is a basis of eigenvectors. Two of them with the same eigenvalue would be two independent vectors in one eigenspace, which is a line by (a). So their eigenvalues are 𝑛 distinct roots of 𝑝, and 𝑝 has at least, hence exactly, 𝑛 distinct roots (Proposition 7.1.22). In the language of Section 7.4, 𝐶𝑝 is diagonalizable over 𝐹 exactly when 𝑝 splits with distinct roots.

Exercise 7.3.46.standard Let 𝐴=(1−10110002)∈ℂ3×3.

  1. (a)For 𝐯=𝐞1+𝐞3, compute 𝐯,𝐴𝐯,𝐴2𝐯,𝐴3𝐯, and find the polynomial 𝑝 of Proposition 7.2.35.
  2. (b)Factor 𝑝 over ℂ, and use the second proof of Theorem 7.3.27 to produce an eigenvector for each root.
  3. (c)Repeat (a) for 𝐯=𝐞1. Which eigenvalue does this Krylov sequence miss, and why is that consistent with Proposition 7.3.26?

Solution.

Uses. Proposition 7.2.35 (A basis of a Krylov subspace), Theorem 7.3.27 (Existence of eigenvalues over ℂ), Proposition 7.3.26 (Krylov subspaces and the characteristic polynomial), Remark 7.1.11 (Polynomials in a matrix)

Techniques. Krylov sequences must become dependent (Technique 7.3.28)

Hint 1. 𝐯,𝐴𝐯,𝐴2𝐯 turn out to be independent. Solve 𝐴3𝐯=𝛼𝐯+𝛽𝐴𝐯+𝛾𝐴2𝐯.
Hint 2. 𝑝=𝑥3−4𝑥2+6𝑥−4 has the root 2. For a root 𝜆, apply 𝑞=𝑝/(𝑥−𝜆) at 𝐴 to 𝐯.
Sketch. (a) 𝐯=(1,0,1), 𝐴𝐯=(1,1,2), 𝐴2𝐯=(0,2,4), 𝐴3𝐯=(−2,2,8)=4𝐯−6𝐴𝐯+4𝐴2𝐯, so 𝑝=𝑥3−4𝑥2+6𝑥−4. (b) 𝑝=(𝑥−2)(𝑥−1−𝑖)(𝑥−1+𝑖), with eigenvectors (0,0,2), (−1−𝑖,−1+𝑖,0), and (−1+𝑖,−1−𝑖,0). (c) 𝐴2𝐞1=2𝐴𝐞1−2𝐞1, so 𝑝=𝑥2−2𝑥+2, which misses 2; it is a proper divisor of 𝑝𝐴.

(a) We get 𝐯=(1,0,1), 𝐴𝐯=(1,1,2), 𝐴2𝐯=(0,2,4), and 𝐴3𝐯=(−2,2,8). The first three are independent: the matrix with these columns has determinant 1⋅(4−4)−1⋅(0−2)+0=2≠0. Solving 𝐴3𝐯=𝛼𝐯+𝛽𝐴𝐯+𝛾𝐴2𝐯 coordinatewise gives 𝛼+𝛽=−2, 𝛽+2𝛾=2, 𝛼+2𝛽+4𝛾=8, so 𝛾=4, 𝛽=−6, 𝛼=4. Hence 𝑑=3 and 𝑝=𝑥3−4𝑥2+6𝑥−4.

(b) 𝑝(2)=0, and dividing, 𝑝=(𝑥−2)(𝑥2−2𝑥+2)=(𝑥−2)(𝑥−(1+𝑖))(𝑥−(1−𝑖)). For each root 𝜆 let 𝑞=𝑝/(𝑥−𝜆); then 𝑞(𝐴)𝐯 is an eigenvector for 𝜆.

  • 𝜆=2: 𝑞=𝑥2−2𝑥+2, and 𝑞(𝐴)𝐯=(0,2,4)−2(1,1,2)+2(1,0,1)=(0,0,2).
  • 𝜆=1+𝑖: 𝑞=(𝑥−2)(𝑥−1+𝑖)=𝑥2−(3−𝑖)𝑥+(2−2𝑖), and

    𝑞(𝐴)𝐯=(0,2,4)−(3−𝑖)(1,1,2)+(2−2𝑖)(1,0,1)=(−1−𝑖,−1+𝑖,0)=(−1+𝑖)(𝑖,1,0).
  • 𝜆=1−𝑖: conjugating the previous computation, since 𝐴 and 𝐯 are real, gives (−1+𝑖,−1−𝑖,0)=(−1−𝑖)(−𝑖,1,0).

Check: 𝐴(0,0,1)=2(0,0,1) and 𝐴(𝑖,1,0)=(𝑖−1,𝑖+1,0)=(1+𝑖)(𝑖,1,0).

(c) For 𝐯=𝐞1: 𝐴𝐞1=(1,1,0), which is not a multiple of 𝐞1, and 𝐴2𝐞1=(0,2,0)=2𝐴𝐞1−2𝐞1. So 𝑑=2 and 𝑝=𝑥2−2𝑥+2, which produces the eigenvectors (𝐴−(1−𝑖)𝐼3)𝐞1=(𝑖,1,0) and (𝐴−(1+𝑖)𝐼3)𝐞1=(−𝑖,1,0) but misses the eigenvalue 2. This is consistent with Proposition 7.3.26: 𝑝 divides 𝑝𝐴=(𝑥−2)(𝑥2−2𝑥+2), but need not equal it, because 𝐾𝐴(𝐞1)=span(𝐞1,𝐞2) is a proper invariant subspace.

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

Claim. For every 𝐴∈𝐹𝑛×𝑛, 𝑝𝐴(𝐴)=0.

“Proof.” By definition, 𝑝𝐴(𝑥)=det(𝑥𝐼𝑛−𝐴). Substituting 𝑥=𝐴 gives 𝑝𝐴(𝐴)=det(𝐴𝐼𝑛−𝐴)=det(𝐴−𝐴)=det0=0.

Solution.

Uses. Definition 7.3.4 (Characteristic polynomial), Proposition 7.3.2 (Determinants over a commutative ring), Remark 7.1.11 (Polynomials in a matrix)

Hint 1. What kind of object is each side? 𝑝𝐴(𝐴) is a matrix. What is det(…)?
Hint 2. Proposition 7.3.2(g) lets us substitute 𝐴 for 𝑥 in every entry of 𝑥𝐼𝑛−𝐴. What does that array look like? Is it the matrix 𝐴−𝐴?
Hint 3. Apply the same reasoning to the single entry 𝑥−𝑎11 of 𝑥𝐼𝑛−𝐴, instead of to its determinant.
Sketch. The two sides have different types: 𝑝𝐴(𝐴) is a matrix and a determinant is a scalar. Substituting 𝐴 for 𝑥 entrywise turns 𝑥𝐼𝑛−𝐴 into an 𝑛×𝑛 array of matrices, with 𝐴−𝑎𝑖𝑖𝐼𝑛 on the diagonal and −𝑎𝑖𝑗𝐼𝑛 elsewhere, not into 𝐴−𝐴=0. The same move would show 𝐴=𝑎11𝐼𝑛 for every 𝐴. The claim is true (Cayley–Hamilton) but needs a real proof.

The claim is true; it is the Cayley–Hamilton theorem (Theorem 7.6.20), and Exercise 7.3.51 proves it. The argument is not a proof.

First, the types do not match: 𝑝𝐴(𝐴)=𝐴𝑛+𝑐𝑛−1𝐴𝑛−1+⋯+𝑐0𝐼𝑛 is an 𝑛×𝑛 matrix, while a determinant is a scalar.

Second, the substitution is done wrongly. The identity 𝑝𝐴=det(𝑥𝐼𝑛−𝐴) holds in 𝐹[𝑥], and Proposition 7.3.2(g) does allow us to apply a map like “substitute 𝐴 for 𝑥” to it, but entry by entry. The map 𝑞↦𝑞(𝐴) takes values in the commutative ring of polynomials in 𝐴 (Remark 7.1.11), and applied to the entries of 𝑥𝐼𝑛−𝐴 it gives the 𝑛×𝑛 array whose (𝑖,𝑖) entry is the matrix 𝐴−𝑎𝑖𝑖𝐼𝑛 and whose (𝑖,𝑗) entry is the matrix −𝑎𝑖𝑗𝐼𝑛. The determinant of that array, computed in this ring, is indeed 𝑝𝐴(𝐴). But the array is not the zero matrix “𝐴𝐼𝑛−𝐴”: the scalar entries 𝑎𝑖𝑗 became the matrices 𝑎𝑖𝑗𝐼𝑛, not the scalars 𝑎𝑖𝑗. For 𝐴=(0100), the array is (𝐴−𝐼20𝐴), whose determinant in this ring is 𝐴2=0, although no entry is 0.

Third, the argument proves too much. Apply the same move to the polynomial 𝑓=𝑥−𝑎11, the (1,1) entry of 𝑥𝐼𝑛−𝐴: “substituting 𝑥=𝐴” gives 𝑓(𝐴)=(𝐴𝐼𝑛−𝐴)11=0, that is, 𝐴=𝑎11𝐼𝑛 for every 𝐴, which is false.

Exercise 7.3.48.hard Let 𝐴∈𝐹𝑛×𝑛 with 𝐴𝑘=0 for some 𝑘≥1.

  1. (a)Show that in 𝐹[𝑥]𝑛×𝑛,

    (𝑥𝐼𝑛−𝐴)(𝑥𝑘−1𝐼𝑛+𝑥𝑘−2𝐴+⋯+𝑥𝐴𝑘−2+𝐴𝑘−1)=𝑥𝑘𝐼𝑛.
  2. (b)Deduce that 𝑝𝐴=𝑥𝑛. In particular tr𝐴=0, and 0 is the only eigenvalue of 𝐴.
  3. (c)For 𝐹=ℂ, give a shorter proof of (b) using eigenvalues. Why does it not work over every field?

Solution.

Uses. Proposition 7.3.2 (Determinants over a commutative ring), Theorem 7.1.20 (Factor theorem), Corollary 7.1.7 (No zero divisors, cancellation, and units), Proposition 7.3.18 (The first and last coefficients), Proposition 7.2.28 (Polynomials act on eigenvectors by their values), Corollary 7.1.35 (Factorization over ℂ)

Techniques. Polynomial identity arguments (Technique 5.5.24)

Hint 1. For (a), multiply out; the sum telescopes.
Hint 2. For (b), take determinants: 𝑝𝐴𝑔=𝑥𝑘𝑛 for some 𝑔∈𝐹[𝑥]. So 𝑝𝐴 divides a power of 𝑥.
Hint 3. Write 𝑝𝐴=𝑥𝑚ℎ with ℎ(0)≠0. Show that ℎ divides 1, by evaluating at 0 and cancelling 𝑥 repeatedly.
Sketch. (a) Telescoping: ∑𝑗𝑥𝑘−𝑗𝐴𝑗−∑𝑗𝑥𝑘−1−𝑗𝐴𝑗+1=𝑥𝑘𝐼𝑛−𝐴𝑘. (b) Taking determinants, 𝑝𝐴⋅𝑔=𝑥𝑘𝑛. Write 𝑝𝐴=𝑥𝑚ℎ with ℎ(0)≠0; cancelling 𝑥𝑚 gives ℎ𝑔=𝑥𝑘𝑛−𝑚, and evaluating at 0 repeatedly shows that 𝑔 absorbs every factor 𝑥, so ℎ𝑔′=1 and ℎ is a constant, hence 1. So 𝑝𝐴=𝑥𝑛. (c) Over ℂ, every eigenvalue satisfies 𝜆𝑘=0, and 𝑝𝐴 splits.

(a) Let 𝑁=∑𝑗=0𝑘−1𝑥𝑘−1−𝑗𝐴𝑗. By Proposition 7.3.2(a), and since the scalar matrices 𝑥𝑚𝐼𝑛 commute with every matrix,

(𝑥𝐼𝑛−𝐴)𝑁=∑𝑗=0𝑘−1𝑥𝑘−𝑗𝐴𝑗−∑𝑗=0𝑘−1𝑥𝑘−1−𝑗𝐴𝑗+1=∑𝑗=0𝑘−1𝑥𝑘−𝑗𝐴𝑗−∑𝑗=1𝑘𝑥𝑘−𝑗𝐴𝑗=𝑥𝑘𝐼𝑛−𝐴𝑘,

which is 𝑥𝑘𝐼𝑛 because 𝐴𝑘=0.

(b) Let 𝑔=det𝑁∈𝐹[𝑥]. Taking determinants in (a) with Proposition 7.3.2(c) and (e), 𝑝𝐴𝑔=det(𝑥𝑘𝐼𝑛)=𝑥𝑘𝑛. Since 𝑝𝐴≠𝟎, write 𝑝𝐴=𝑥𝑚ℎ with 𝑚≥0 as large as possible, so that 𝑥 does not divide ℎ, that is, ℎ(0)≠0 (Theorem 7.1.20). Comparing degrees in 𝑥𝑚ℎ𝑔=𝑥𝑘𝑛 gives 𝑚≤𝑘𝑛, and cancelling 𝑥𝑚 (Corollary 7.1.7(b)) gives ℎ𝑔=𝑥𝑘𝑛−𝑚.

We show by induction on 𝑟=𝑘𝑛−𝑚 that an equation ℎ𝑔=𝑥𝑟 forces ℎ to be a nonzero constant. If 𝑟=0, then ℎ𝑔=1, so ℎ is invertible in 𝐹[𝑥], hence a nonzero constant (Corollary 7.1.7(c)). If 𝑟≥1, evaluating at 0 gives ℎ(0)𝑔(0)=0, so 𝑔(0)=0 and 𝑔=𝑥𝑔1 (Theorem 7.1.20). Then 𝑥(ℎ𝑔1)=𝑥⋅𝑥𝑟−1, and cancelling 𝑥 gives ℎ𝑔1=𝑥𝑟−1, to which the induction hypothesis applies.

So ℎ is a constant, and ℎ=1 because 𝑝𝐴 is monic. Thus 𝑝𝐴=𝑥𝑚, and 𝑚=deg𝑝𝐴=𝑛 (Proposition 7.3.18). The coefficient of 𝑥𝑛−1 is 0=−tr𝐴, and 0 is the only root.

(c) Over ℂ: if 𝐴𝐯=𝜆𝐯 with 𝐯≠𝟎, then 𝟎=𝐴𝑘𝐯=𝜆𝑘𝐯 (Proposition 7.2.28), so 𝜆=0. By Corollary 7.1.35, 𝑝𝐴 is a product of factors 𝑥−𝜆𝑗 whose roots 𝜆𝑗 are eigenvalues, hence 0; so 𝑝𝐴=𝑥𝑛. Over a general field this argument only shows that 0 is the only root of 𝑝𝐴 in 𝐹; it cannot rule out an irreducible factor of degree at least 2, because 𝑝𝐴 need not split. The proof of (b) never needs 𝑝𝐴 to split.

Exercise 7.3.49.hard Let 𝐴∈𝐹𝑛×𝑛. For 𝑓∈𝐹[𝑥], write 𝑓(𝑥2) and 𝑓(−𝑥) for the polynomials obtained by substituting 𝑥2 and −𝑥 for 𝑥 (Definition 7.1.9).

  1. (a)Show that 𝑝𝐴2(𝑥2)=(−1)𝑛𝑝𝐴(𝑥)𝑝𝐴(−𝑥).
  2. (b)Suppose 𝑝𝐴=(𝑥−𝜆1)⋯(𝑥−𝜆𝑛). Show that 𝑝𝐴2=(𝑥−𝜆12)⋯(𝑥−𝜆𝑛2).
  3. (c)Deduce that if 𝐴∈ℝ𝑛×𝑛, 𝑝𝐴 splits over ℝ, and tr(𝐴2)=0, then 𝑝𝐴=𝑥𝑛.

Solution.

Uses. Proposition 7.3.2 (Determinants over a commutative ring), Proposition 7.1.10 (Evaluation respects sums and products), Remark 7.3.19 (Trace and determinant from the roots)

Techniques. Polynomial identity arguments (Technique 5.5.24)

Hint 1. Factor 𝑥2𝐼𝑛−𝐴2 in 𝐹[𝑥]𝑛×𝑛.
Hint 2. det(𝑥2𝐼𝑛−𝐴2) is 𝑝𝐴2 with 𝑥2 substituted, by Proposition 7.3.2(g). Similarly, what is det(𝑥𝐼𝑛+𝐴) in terms of 𝑝𝐴?
Hint 3. For (b), (𝑥−𝜆)(𝑥+𝜆)=𝑥2−𝜆2. A polynomial 𝑟 is determined by 𝑟(𝑥2).
Sketch. (a) 𝑥2𝐼𝑛−𝐴2=(𝑥𝐼𝑛−𝐴)(𝑥𝐼𝑛+𝐴). Substituting 𝑥2 into 𝑥𝐼𝑛−𝐴2 gives 𝑥2𝐼𝑛−𝐴2, so its determinant is 𝑝𝐴2(𝑥2); and det(𝑥𝐼𝑛+𝐴)=(−1)𝑛det(−𝑥𝐼𝑛−𝐴)=(−1)𝑛𝑝𝐴(−𝑥). (b) (−1)𝑛𝑝𝐴(−𝑥)=∏(𝑥+𝜆𝑖), so the right side is ∏(𝑥2−𝜆𝑖2)=𝑞(𝑥2) with 𝑞=∏(𝑥−𝜆𝑖2); compare coefficients. (c) tr(𝐴2)=∑𝜆𝑖2=0 with real 𝜆𝑖.

(a) In 𝐹[𝑥]𝑛×𝑛, (𝑥𝐼𝑛−𝐴)(𝑥𝐼𝑛+𝐴)=𝑥2𝐼𝑛+𝑥𝐴−𝑥𝐴−𝐴2=𝑥2𝐼𝑛−𝐴2. Substituting 𝑥2 for 𝑥 respects sums, products, and 1 (Proposition 7.1.10), and applied to every entry of 𝑥𝐼𝑛−𝐴2 it gives 𝑥2𝐼𝑛−𝐴2; so by Proposition 7.3.2(g), det(𝑥2𝐼𝑛−𝐴2)=𝑝𝐴2(𝑥2). In the same way, substituting −𝑥 into 𝑥𝐼𝑛−𝐴 gives −𝑥𝐼𝑛−𝐴, so det(−𝑥𝐼𝑛−𝐴)=𝑝𝐴(−𝑥), and det(𝑥𝐼𝑛+𝐴)=(−1)𝑛det(−𝑥𝐼𝑛−𝐴)=(−1)𝑛𝑝𝐴(−𝑥) by multilinearity. Taking determinants of the factorization (Proposition 7.3.2(c)) gives 𝑝𝐴2(𝑥2)=𝑝𝐴(𝑥)⋅(−1)𝑛𝑝𝐴(−𝑥).

(b) Substituting −𝑥 into the factorization gives 𝑝𝐴(−𝑥)=∏𝑖(−𝑥−𝜆𝑖)=(−1)𝑛∏𝑖(𝑥+𝜆𝑖). So by (a),

𝑝𝐴2(𝑥2)=∏𝑖(𝑥−𝜆𝑖)(𝑥+𝜆𝑖)=∏𝑖(𝑥2−𝜆𝑖2)=𝑞(𝑥2),𝑞=∏𝑖(𝑥−𝜆𝑖2).

For any 𝑟=∑𝑘𝑟𝑘𝑥𝑘, the coefficient of 𝑥2𝑘 in 𝑟(𝑥2)=∑𝑘𝑟𝑘𝑥2𝑘 is 𝑟𝑘. So 𝑟(𝑥2)=𝑠(𝑥2) implies 𝑟=𝑠, and 𝑝𝐴2=𝑞.

(c) Let 𝑝𝐴=∏(𝑥−𝜆𝑖) with 𝜆𝑖∈ℝ. By (b) and Remark 7.3.19 applied to 𝐴2, tr(𝐴2)=𝜆12+⋯+𝜆𝑛2. If this is 0, then every 𝜆𝑖=0, since a sum of squares of reals vanishes only when each term does. So 𝑝𝐴=𝑥𝑛.

Exercise 7.3.50.hard Let dim𝑉=𝑛≥1 and 𝑇∈L(𝑉). A vector 𝐯 is a cyclic vector for 𝑇 if 𝐾𝑇(𝐯)=𝑉.

  1. (a)Show that if 𝑇 has a cyclic vector, then the matrix of 𝑇 in some basis is 𝐶𝑝𝑇.
  2. (b)Suppose 𝐴,𝐵∈𝐹𝑛×𝑛 both have cyclic vectors (for 𝑇𝐴 and 𝑇𝐵). Show that 𝐴 and 𝐵 are similar if and only if 𝑝𝐴=𝑝𝐵.
  3. (c)Show that 𝐼2 has no cyclic vector, and that 𝐶𝑝 and 𝐶𝑝𝖳 are similar for every monic 𝑝.

Solution.

Uses. Proposition 7.2.35 (A basis of a Krylov subspace), Proposition 7.3.26 (Krylov subspaces and the characteristic polynomial), Theorem 4.5.25 (Similar matrices are one operator in two bases), Proposition 4.5.24 (Similarity is an equivalence relation), Proposition 7.3.11 (Similar matrices have the same characteristic polynomial), Exercise 7.3.17

Hint 1. For (a), if 𝐾𝑇(𝐯)=𝑉, what is 𝑑 in Proposition 7.2.35? Then use Proposition 7.3.26.
Hint 2. For (b), both matrices are similar to the same companion matrix.
Hint 3. For (c), 𝐶𝑝 has the cyclic vector 𝐞1. For 𝐶𝑝𝖳, try 𝐞𝑛 and show that (𝐶𝑝𝖳)𝑗𝐞𝑛 is 𝐞𝑛−𝑗 plus a combination of 𝐞𝑛−𝑗+1,…,𝐞𝑛.
Sketch. (a) 𝑑=𝑛, and in the basis 𝐯,…,𝑇𝑛−1𝐯 the matrix is 𝐶𝑝 with 𝑝=𝑝𝑇 (Proposition 7.3.26). (b) Both are similar to 𝐶𝑓 for 𝑓=𝑝𝐴=𝑝𝐵; the converse is Proposition 7.3.11. (c) dim𝐾(𝐯)≤1 for 𝐼2. 𝐶𝑝 is cyclic with 𝐞1, 𝐶𝑝𝖳 with 𝐞𝑛 (triangular argument), and 𝑝𝐶𝑝=𝑝𝐶𝑝𝖳.

(a) Let 𝐯 be cyclic, so 𝐯≠𝟎 (as 𝑉≠{𝟎}). In Proposition 7.2.35, 𝑑=dim𝐾𝑇(𝐯)=𝑛, so 𝐯,𝑇𝐯,…,𝑇𝑛−1𝐯 is a basis of 𝑉, and by Proposition 7.3.26 the matrix of 𝑇 in it is 𝐶𝑝, with 𝑝𝑇=𝑝. So the matrix is 𝐶𝑝𝑇.

(b) If 𝐴∼𝐵, then 𝑝𝐴=𝑝𝐵 by Proposition 7.3.11. Conversely, let 𝑓=𝑝𝐴=𝑝𝐵. By (a), 𝑇𝐴 has matrix 𝐶𝑓 in some basis, so 𝐴∼𝐶𝑓 by Theorem 4.5.25, and likewise 𝐵∼𝐶𝑓. By Proposition 4.5.24, 𝐴∼𝐵.

(c) For every 𝐯∈𝐹2, 𝐾𝐼2(𝐯)=span(𝐯) has dimension at most 1, so 𝐼2 has no cyclic vector. (Indeed 𝐼2 and (1101) have the same characteristic polynomial and are not similar, so (b) needs the hypothesis.)

Now let 𝐶=𝐶𝑝. Since 𝐶𝑗𝐞1=𝐞𝑗+1 for 𝑗<𝑛, 𝐞1 is cyclic for 𝐶. For 𝐶𝖳: column 𝑗 of 𝐶𝖳 is row 𝑗 of 𝐶, which is −𝑎0𝐞𝑛𝖳 for 𝑗=1 and 𝐞𝑗−1𝖳−𝑎𝑗−1𝐞𝑛𝖳 for 𝑗≥2. So 𝐶𝖳𝐞𝑗=𝐞𝑗−1−𝑎𝑗−1𝐞𝑛 for 𝑗≥2. By induction on 𝑗, (𝐶𝖳)𝑗𝐞𝑛=𝐞𝑛−𝑗+𝐫𝑗 with 𝐫𝑗∈span(𝐞𝑛−𝑗+1,…,𝐞𝑛) for 0≤𝑗≤𝑛−1 (with 𝐫0=𝟎): applying 𝐶𝖳 lowers each index by one and adds a multiple of 𝐞𝑛. So (𝐶𝖳)𝑗𝐞𝑛 has coordinate 1 at 𝐞𝑛−𝑗 and coordinate 0 at 𝐞1,…,𝐞𝑛−𝑗−1. Hence the matrix with columns (𝐶𝖳)𝑛−1𝐞𝑛,…,𝐶𝖳𝐞𝑛,𝐞𝑛 is lower triangular with 1′s on the diagonal, so these vectors are independent (Proposition 5.3.15, Theorem 5.4.7), and 𝐞𝑛 is cyclic for 𝐶𝖳. Since 𝑝𝐶𝖳=𝑝𝐶 (Exercise 7.3.17), (b) shows 𝐶∼𝐶𝖳.

Exercise 7.3.51.hardBeyond 21-242 Let 𝑉 be finite-dimensional with dim𝑉≥1, and 𝑇∈L(𝑉).

  1. (a)Use Proposition 7.3.26 to show that 𝑝𝑇(𝑇)=0. This is the Cayley–Hamilton theorem (Theorem 7.6.20).
  2. (b)Deduce that every 𝐴∈𝐹𝑛×𝑛 satisfies a monic polynomial of degree 𝑛, improving the degree 𝑛2 of Proposition 4.1.31.

Solution.

Uses. Proposition 7.3.26 (Krylov subspaces and the characteristic polynomial), Proposition 7.2.27 (Evaluating polynomials at an operator), Proposition 4.1.31 (Every square matrix satisfies a polynomial)

Techniques. Krylov sequences must become dependent (Technique 7.3.28)

Hint 1. It is enough to show 𝑝𝑇(𝑇)𝐯=𝟎 for each 𝐯. Fix 𝐯≠𝟎 and its Krylov polynomial 𝑝.
Hint 2. 𝑝 divides 𝑝𝑇, so 𝑝𝑇=𝑔𝑝. Evaluate at 𝑇 and apply to 𝐯.
Sketch. For 𝐯≠𝟎 with Krylov polynomial 𝑝, 𝑝𝑇=𝑔𝑝 by Proposition 7.3.26(c), so 𝑝𝑇(𝑇)𝐯=𝑔(𝑇)𝑝(𝑇)𝐯=𝟎. For matrices, apply this to 𝑇𝐴 and use [𝑞(𝑇)]=𝑞([𝑇]).

(a) Let 𝐯∈𝑉. If 𝐯=𝟎, then 𝑝𝑇(𝑇)𝐯=𝟎. Otherwise let 𝑝 be the polynomial of Proposition 7.2.35 for 𝐯, so 𝑝(𝑇)𝐯=𝟎. By Proposition 7.3.26(c), 𝑝𝑇=𝑔𝑝 for some 𝑔∈𝐹[𝑥], and by Proposition 7.2.27(a),

𝑝𝑇(𝑇)𝐯=𝑔(𝑇)(𝑝(𝑇)𝐯)=𝑔(𝑇)𝟎=𝟎.

So 𝑝𝑇(𝑇) sends every vector to 𝟎, that is, 𝑝𝑇(𝑇)=0.

(b) Apply (a) to 𝑇𝐴 on 𝐹𝑛. Its matrix in the standard basis is 𝐴, so 𝑝𝑇𝐴=𝑝𝐴, and by Proposition 7.2.27(c), 𝑝𝐴(𝐴) is the matrix of 𝑝𝐴(𝑇𝐴)=0. So 𝑝𝐴(𝐴)=0, and 𝑝𝐴 is monic of degree 𝑛.

Exercise 7.3.52.hard Let 𝑞 be a prime and 𝑓=𝑥2−𝑡𝑥+𝑑∈𝔽𝑞[𝑥].

  1. (a)Show that the number of matrices 𝐴∈𝔽𝑞2×2 with 𝑝𝐴=𝑓 is 𝑞2+𝑞 if 𝑓 has two distinct roots in 𝔽𝑞, 𝑞2 if it has one repeated root, and 𝑞2−𝑞 if it has no root.
  2. (b)Deduce that exactly (𝑞2−𝑞)2/2 matrices in 𝔽𝑞2×2 have no eigenvalue in 𝔽𝑞, and check that the counts in (a) add up to 𝑞4.

Solution.

Uses. Example 7.3.5 (Small cases), Theorem 7.3.6 (Eigenvalues are the roots of the characteristic polynomial), Proposition 7.1.22 (Counting roots)

Techniques. Double counting (Technique 4.1.36)

Hint 1. Write 𝐴=(𝑎𝑏𝑐𝑒). The conditions are 𝑎+𝑒=𝑡 and 𝑎𝑒−𝑏𝑐=𝑑. Choose 𝑎 first.
Hint 2. Then 𝑒=𝑡−𝑎 and 𝑏𝑐=𝑎(𝑡−𝑎)−𝑑=−𝑓(𝑎). How many pairs (𝑏,𝑐) have a given product 𝑚? Distinguish 𝑚=0 and 𝑚≠0.
Hint 3. If 𝑓 has 𝑟 roots in 𝔽𝑞, the total is 𝑟(2𝑞−1)+(𝑞−𝑟)(𝑞−1)=𝑞(𝑞−1)+𝑟𝑞.
Sketch. For each 𝑎, 𝑒=𝑡−𝑎 is forced and 𝑏𝑐=−𝑓(𝑎); there are 2𝑞−1 pairs (𝑏,𝑐) if 𝑓(𝑎)=0 and 𝑞−1 otherwise. So the count is 𝑞(𝑞−1)+𝑟𝑞, where 𝑟∈{0,1,2} is the number of roots. There are (𝑞2−𝑞)/2 monic quadratics without roots, each realized by 𝑞2−𝑞 matrices.

(a) Write 𝐴=(𝑎𝑏𝑐𝑒). By Example 7.3.5(a), 𝑝𝐴=𝑓 if and only if 𝑎+𝑒=𝑡 and 𝑎𝑒−𝑏𝑐=𝑑. Choose 𝑎∈𝔽𝑞 arbitrarily (𝑞 ways). Then 𝑒=𝑡−𝑎 is determined, and the second condition becomes

𝑏𝑐=𝑎(𝑡−𝑎)−𝑑=−(𝑎2−𝑡𝑎+𝑑)=−𝑓(𝑎).

The number of pairs (𝑏,𝑐) with 𝑏𝑐=𝑚 is 𝑞−1 if 𝑚≠0 (any 𝑏≠0, and then 𝑐=𝑚/𝑏), and 2𝑞−1 if 𝑚=0 (𝑏=0 with 𝑐 arbitrary, or 𝑏≠0 with 𝑐=0). So if 𝑓 has 𝑟 roots in 𝔽𝑞, the number of matrices is

𝑟(2𝑞−1)+(𝑞−𝑟)(𝑞−1)=𝑞(𝑞−1)+𝑟𝑞,

which is 𝑞2+𝑞, 𝑞2, 𝑞2−𝑞 for 𝑟=2,1,0.

(b) By Theorem 7.3.6, 𝐴 has no eigenvalue exactly when 𝑝𝐴 has no root. There are 𝑞2 monic quadratics 𝑥2−𝑡𝑥+𝑑. Those with two distinct roots are the products (𝑥−𝛼)(𝑥−𝛽) with 𝛼≠𝛽, and there are (𝑞2)=(𝑞2−𝑞)/2 of them; those with a repeated root are the (𝑥−𝛼)2, and there are 𝑞 of them (a monic quadratic is determined by its roots, by Proposition 7.1.22 and the factor theorem). So 𝑞2−(𝑞2−𝑞)/2−𝑞=(𝑞2−𝑞)/2 have no root. Each is realized by 𝑞2−𝑞 matrices, which gives (𝑞2−𝑞)2/2 matrices without eigenvalues. The total is

𝑞2−𝑞2(𝑞2+𝑞)+𝑞⋅𝑞2+𝑞2−𝑞2(𝑞2−𝑞)=𝑞2−𝑞2⋅2𝑞2+𝑞3=𝑞4−𝑞3+𝑞3=𝑞4,

as it should be. For 𝑞=2 the two matrices without eigenvalues are (0111) and (1110).

Exercise 7.3.53.hardBeyond 21-242 For 𝐴∈𝐹𝑛×𝑛 and 𝑆⊆{1,…,𝑛}, let 𝐴𝑆 be the principal submatrix of 𝐴 with rows and columns in 𝑆, and let 𝑒𝑘(𝐴)=∑|𝑆|=𝑘det𝐴𝑆 for 1≤𝑘≤𝑛, and 𝑒0(𝐴)=1. Show that

𝑝𝐴=∑𝑘=0𝑛(−1)𝑘𝑒𝑘(𝐴)𝑥𝑛−𝑘,

and check that this agrees with Proposition 7.3.18. For 𝑛=3, write out the coefficient of 𝑥.

Solution.

Uses. Proposition 7.3.2 (Determinants over a commutative ring), Proposition 7.3.18 (The first and last coefficients)

Techniques. Expand by multilinearity (Technique 5.1.7)

Hint 1. Column 𝑗 of 𝑥𝐼𝑛−𝐴 is 𝑥𝐞𝑗−𝐚𝑗, a sum of two vectors. Expand det by multilinearity in all columns at once.
Hint 2. You get 2𝑛 terms indexed by the set 𝑆 of columns where −𝐚𝑗 was chosen. In the term for 𝑆, the columns 𝑥𝐞𝑗 with 𝑗∉𝑆 have a single nonzero entry. Expand along them.
Sketch. Multilinearity gives 𝑝𝐴=∑𝑆det𝑀𝑆, where column 𝑗 of 𝑀𝑆 is −𝐚𝑗 for 𝑗∈𝑆 and 𝑥𝐞𝑗 otherwise. Expanding along the columns 𝑥𝐞𝑗, 𝑗∉𝑆, each contributes a factor 𝑥 and deletes row and column 𝑗, leaving −𝐴𝑆. So det𝑀𝑆=𝑥𝑛−|𝑆|(−1)|𝑆|det𝐴𝑆.

Column 𝑗 of 𝑥𝐼𝑛−𝐴 is 𝑥𝐞𝑗−𝐚𝑗, where 𝐚𝑗 is column 𝑗 of 𝐴. Expanding by linearity in each of the 𝑛 columns (Proposition 7.3.2(b), Technique 5.1.7),

𝑝𝐴=∑𝑆⊆{1,…,𝑛}det𝑀𝑆,

where column 𝑗 of 𝑀𝑆 is −𝐚𝑗 for 𝑗∈𝑆 and 𝑥𝐞𝑗 for 𝑗∉𝑆.

We claim that det𝑀𝑆=𝑥𝑛−|𝑆|det(−𝐴𝑆) (where det(−𝐴∅)=1), by induction on 𝑛−|𝑆|. If 𝑆={1,…,𝑛}, then 𝑀𝑆=−𝐴. If 𝑛=1 and 𝑆=∅, then 𝑀𝑆=(𝑥) and the claim holds. Otherwise 𝑛≥2; pick 𝑗∉𝑆. Column 𝑗 of 𝑀𝑆 is 𝑥𝐞𝑗, with the single nonzero entry 𝑥 in position (𝑗,𝑗), so cofactor expansion along it (Proposition 7.3.2(f)) gives det𝑀𝑆=𝑥(−1)𝑗+𝑗det𝑀′, where 𝑀′ is 𝑀𝑆 with row 𝑗 and column 𝑗 deleted. In 𝑀′, each remaining column 𝑥𝐞𝑗′ still has its single entry 𝑥 on the diagonal, and each column −𝐚𝑗′ has lost its 𝑗-th entry. So 𝑀′ is the matrix of the same form for the (𝑛−1)×(𝑛−1) matrix obtained from 𝐴 by deleting row and column 𝑗, with the same set 𝑆, and the induction hypothesis applies. Since det(−𝐴𝑆)=(−1)|𝑆|det𝐴𝑆 by multilinearity,

𝑝𝐴=∑𝑆(−1)|𝑆|𝑥𝑛−|𝑆|det𝐴𝑆=∑𝑘=0𝑛(−1)𝑘𝑒𝑘(𝐴)𝑥𝑛−𝑘.

For 𝑘=1, the principal 1×1 minors are the diagonal entries, so 𝑒1(𝐴)=tr𝐴; for 𝑘=𝑛, 𝑒𝑛(𝐴)=det𝐴. This matches Proposition 7.3.18. For 𝑛=3, the coefficient of 𝑥 is

𝑒2(𝐴)=(𝑎11𝑎22−𝑎12𝑎21)+(𝑎11𝑎33−𝑎13𝑎31)+(𝑎22𝑎33−𝑎23𝑎32).

Exercise 7.3.54.standardBeyond 21-242 Let 𝑉 be a finite-dimensional real vector space, and let 𝑇∈L(𝑉) have no eigenvalue.

  1. (a)Show that every 𝑇-invariant subspace of 𝑉 has even dimension.
  2. (b)Give such a 𝑇 on ℝ4, and find three different 𝑇-invariant planes.

Solution.

Uses. Proposition 7.3.30 (Odd dimension), Definition 7.2.1 (Invariant subspace), Proposition 7.3.15 (Block-triangular matrices), Proposition 7.2.4 (Invariance is checked on a spanning list)

Hint 1. If 𝑈 is invariant of odd dimension, then 𝑇↾𝑈 is an operator on an odd-dimensional real space.
Hint 2. For (b), try two copies of the rotation 𝑅𝜋/2. Look for planes spanned by 𝐰 and 𝑇𝐰.
Sketch. (a) An invariant 𝑈 of odd dimension gives an eigenvector of 𝑇↾𝑈 (Proposition 7.3.30), which is an eigenvector of 𝑇. (b) 𝑇=(𝑅00𝑅) with 𝑅=𝑅𝜋/2 has 𝑝𝑇=(𝑥2+1)2; for any 𝐰≠𝟎, span(𝐰,𝑇𝐰) is invariant since 𝑇2=−id.

(a) Let 𝑈 be 𝑇-invariant. If dim𝑈=0 it is even. If dim𝑈 is odd, then 𝑇↾𝑈∈L(𝑈) has an eigenvalue by Proposition 7.3.30: there are 𝐮∈𝑈, 𝐮≠𝟎, and 𝜆∈ℝ with 𝑇𝐮=𝜆𝐮. Then 𝜆 is an eigenvalue of 𝑇, a contradiction.

(b) Let 𝑅=(0−110) and 𝑇=(𝑅00𝑅) on ℝ4. By Proposition 7.3.15, 𝑝𝑇=(𝑥2+1)2, which has no real root, so 𝑇 has no eigenvalue. Since 𝑅2=−𝐼2, 𝑇2=−𝐼4. So for every 𝐰≠𝟎, 𝑇(𝑇𝐰)=−𝐰, and span(𝐰,𝑇𝐰) is 𝑇-invariant (Proposition 7.2.4); it is a plane by (a). Three different ones: 𝐰=𝐞1 gives span(𝐞1,𝐞2); 𝐰=𝐞3 gives span(𝐞3,𝐞4); and 𝐰=𝐞1+𝐞3 gives span(𝐞1+𝐞3,𝐞2+𝐞4), which contains neither 𝐞1 nor 𝐞3.

Exercise 7.3.55.challenge Let 𝐴,𝐵∈ℂ𝑛×𝑛 satisfy 𝐴𝐵−𝐵𝐴=𝐴.

  1. (a)Show that 𝐴𝑘𝐵−𝐵𝐴𝑘=𝑘𝐴𝑘 for every 𝑘≥1.
  2. (b)Show that 𝐴𝑘=0 for some 𝑘.
  3. (c)Deduce that 𝑝𝐴=𝑥𝑛.

Solution.

Uses. Corollary 7.2.22 (At most dim𝑉 eigenvalues), Exercise 7.3.48, Definition 7.2.6 (Eigenvalue, eigenvector)

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

Hint 1. For (a), induct on 𝑘: write 𝐴𝑘+1𝐵−𝐵𝐴𝑘+1=𝐴(𝐴𝑘𝐵−𝐵𝐴𝑘)+(𝐴𝐵−𝐵𝐴)𝐴𝑘.
Hint 2. For (b), consider the linear map 𝐿(𝑋)=𝑋𝐵−𝐵𝑋 on ℂ𝑛×𝑛. What does (a) say about 𝐿(𝐴𝑘)?
Hint 3. An operator on an 𝑛2-dimensional space has at most 𝑛2 eigenvalues, but 1,2,3,… are infinitely many.
Sketch. (a) 𝐴𝑘+1𝐵−𝐵𝐴𝑘+1=𝐴(𝐴𝑘𝐵−𝐵𝐴𝑘)+(𝐴𝐵−𝐵𝐴)𝐴𝑘=𝑘𝐴𝑘+1+𝐴𝑘+1. (b) With 𝐿(𝑋)=𝑋𝐵−𝐵𝑋, 𝐿(𝐴𝑘)=𝑘𝐴𝑘, so if 𝐴𝑘≠0 then 𝑘 is an eigenvalue of 𝐿. 𝐿 has at most 𝑛2 eigenvalues, so 𝐴𝑘=0 for some 𝑘≤𝑛2+1. (c) Exercise 7.3.48.

(a) For 𝑘=1 this is the hypothesis. If it holds for 𝑘, then

𝐴𝑘+1𝐵−𝐵𝐴𝑘+1=𝐴(𝐴𝑘𝐵−𝐵𝐴𝑘)+(𝐴𝐵−𝐵𝐴)𝐴𝑘=𝐴(𝑘𝐴𝑘)+𝐴⋅𝐴𝑘=(𝑘+1)𝐴𝑘+1,

where the first equality holds because the terms −𝐴𝐵𝐴𝑘 and +𝐴𝐵𝐴𝑘 cancel.

(b) Let 𝐿:ℂ𝑛×𝑛→ℂ𝑛×𝑛, 𝐿(𝑋)=𝑋𝐵−𝐵𝑋. It is linear, and (a) says 𝐿(𝐴𝑘)=𝑘𝐴𝑘. If 𝐴𝑘≠0, then 𝑘 is an eigenvalue of 𝐿 (Definition 7.2.6). By Corollary 7.2.22, 𝐿 has at most 𝑛2 eigenvalues, and 1,2,…,𝑛2+1 are 𝑛2+1 distinct complex numbers. So 𝐴𝑘=0 for some 𝑘≤𝑛2+1.

(c) By (b) and Exercise 7.3.48, 𝑝𝐴=𝑥𝑛. The proof used that the integers 1,2,3,… are distinct in ℂ; over 𝔽𝑝 the scalars 𝑘 repeat, and (b) breaks down.

More practice

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