Mathematics · Book 4 · Bachelor Year 2

University Mathematics — Year 2

University Mathematics — Year 2 · Bachelor Year 2

3Reduction of Endomorphisms

To understand an endomorphism, find the directions it merely stretches. This chapter builds the machinery — eigenvalues, characteristic and minimal polynomials, the kernel-decomposition lemma — and its rewards: diagonalization and trigonalization criteria, Cayley–Hamilton, the Dunford decomposition, and the computation of powers and exponentials that Chapter 16 will feed on. Throughout, EE is a finite-dimensional KK-vector space (K=RK = \R or C\C) and uL(E)u \in \mathcal{L}(E), n=dimEn = \dim E.

3.1 Eigenvalues and eigenvectors

Definition 3.1

λK\lambda \in K is an eigenvalue of uu when u(x)=λxu(x) = \lambda x for some x0x \neq 0 (an eigenvector); the eigenspace is Eλ(u)=ker(uλid)E_\lambda(u) = \ker(u - \lambda\,\mathrm{id}). The set of eigenvalues is the spectrum Sp(u)\operatorname{Sp}(u). A subspace FF is stable when u(F)Fu(F) \subseteq F; eigenspaces are stable, and stable subspaces allow induced endomorphisms uFu|_F.

Theorem 3.2 (Independence of eigenspaces)

Eigenvectors associated with pairwise distinct eigenvalues form a free family; equivalently, the sum of the eigenspaces Eλ1++EλrE_{\lambda_1} + \dots + E_{\lambda_r} (distinct λi\lambda_i) is direct. In particular uu has at most nn eigenvalues.

Proof. By induction on rr. Suppose x1++xr=0x_1 + \dots + x_r = 0 with xiEλix_i \in E_{\lambda_i}, the statement known for r1r - 1. Apply uu and subtract λr\lambda_r times the relation:

i=1r1(λiλr)xi=0,\sum_{i=1}^{r-1} (\lambda_i - \lambda_r)\, x_i = 0 ,

so by induction each (λiλr)xi=0(\lambda_i - \lambda_r)x_i = 0, i.e. xi=0x_i = 0 for i<ri < r, then xr=0x_r = 0. Direct sums of nonzero spaces in a space of dimension nn have at most nn summands.

The matrix A = psmallmatrix2 & 1\\ 1 & 2 psmallmatrix acting on the plane: the generic vector e_1 is knocked off its line, but the eigendirections v_1 = (1,1) and v_2 = (1,-1) are merely stretched — by 3 and by 1 (so Av_2 = v_2: the dashed image coincides with v_2). Diagonalization is the change to the basis (v_1, v_2), where A becomes diag(3, 1).
The matrix A=(2112)A = \left(\begin{smallmatrix}2 & 1\\ 1 & 2\end{smallmatrix}\right) acting on the plane: the generic vector e1e_1 is knocked off its line, but the eigendirections v1=(1,1)v_1 = (1,1) and v2=(1,1)v_2 = (1,-1) are merely stretched — by 33 and by 11 (so Av2=v2Av_2 = v_2: the dashed image coincides with v2v_2). Diagonalization is the change to the basis (v1,v2)(v_1, v_2), where AA becomes diag(3,1)\operatorname{diag}(3, 1).

Definition 3.3 (Characteristic polynomial)

χu(X)=det(Xidu)\chi_u(X) = \det(X\,\mathrm{id} - u) — computed in any basis as det(XInA)\det(XI_n - A), a monic polynomial of degree nn, invariant under similarity (Theorem 2.17). Its roots in KK are exactly the eigenvalues (λ\lambda eigenvalue     uλid\iff u - \lambda\,\mathrm{id} not injective     χu(λ)=0\iff \chi_u(\lambda) = 0), and

χu(X)=Xn(tru)Xn1++(1)ndetu.\chi_u(X) = X^n - (\operatorname{tr} u)\, X^{n-1} + \dots + (-1)^n \det u .

The algebraic multiplicity mλm_\lambda of an eigenvalue is its multiplicity as a root of χu\chi_u; the geometric multiplicity is dimEλ\dim E_\lambda, and 1dimEλmλ1 \leq \dim E_\lambda \leq m_\lambda.

Proof of the stated facts. The coefficient claims: expand det(XIA)\det(XI - A) by the permutation formula; the identity permutation contributes i(Xaii)=Xn(aii)Xn1+\prod_i (X - a_{ii}) = X^n - (\sum a_{ii})X^{n-1} + \dots, and every other permutation fixes at most n2n - 2 diagonal positions, contributing degree n2\leq n - 2: the top two coefficients are as stated; X=0X = 0 gives the constant term det(A)=(1)ndetA\det(-A) = (-1)^n\det A.

Geometric \leq algebraic: let d=dimEλd = \dim E_\lambda and complete a basis of EλE_\lambda into a basis of EE; the matrix of uu is block upper-triangular with top-left block λId\lambda I_d, so χu(X)=(Xλ)dχ(bottom block)(X)\chi_u(X) = (X - \lambda)^d\, \chi_{\text{(bottom block)}}(X): the multiplicity of λ\lambda is at least dd.

Example 3.4 (Same χ\chi, different geometry)

The matrices

(2002)and(2102)\begin{pmatrix}2 & 0\\ 0 & 2\end{pmatrix} \qquad\text{and}\qquad \begin{pmatrix}2 & 1\\ 0 & 2\end{pmatrix}

share the characteristic polynomial (X2)2(X - 2)^2, the trace, the determinant, the spectrum — yet are not similar: the first has E2E_2 of dimension 22 (geometric multiplicity 22), the second of dimension 11. The characteristic polynomial sees only algebraic multiplicities; the eigenspace dimensions are the finer invariant, and the minimal polynomial arbitrates (X2X - 2 versus (X2)2(X - 2)^2). Moral for all diagonalizability discussions: χ\chi shortlists the candidates, but kernels cast the votes.

Definition 3.5 (Diagonalizable, trigonalizable)

uu is diagonalizable when EE has a basis of eigenvectors (matrix: similar to a diagonal matrix); trigonalizable when its matrix in some basis is upper triangular.

Theorem 3.6 (Diagonalizability criteria)

The following are equivalent:

  1. uu is diagonalizable;
  2. E=λSpuEλE = \bigoplus_{\lambda \in \operatorname{Sp} u} E_\lambda;
  3. χu\chi_u splits over KK and dimEλ=mλ\dim E_\lambda = m_\lambda for every eigenvalue;
  4. (sufficient, not necessary) χu\chi_u has nn distinct roots in KK.

Proof. (1     \iff 2): a basis of eigenvectors sorts into bases of the EλE_\lambda’s, and conversely concatenating bases of the direct summands gives a basis of EE (Theorem 3.2 makes the sum direct; equality of dimensions makes it everything).

(2     \iff 3): in the diagonal basis, χu=(Xλ)dimEλ\chi_u = \prod (X - \lambda)^{\dim E_\lambda} splits with matching multiplicities. Conversely, suppose χu\chi_u splits with dimEλ=mλ\dim E_\lambda = m_\lambda throughout; then the direct sum of the eigenspaces (direct by Theorem 3.2) has dimension

λdimEλ=λmλ=degχu=n,\sum_{\lambda}\dim E_\lambda = \sum_{\lambda} m_\lambda = \deg\chi_u = n ,

the middle equality because a split polynomial’s degree is the sum of its root multiplicities: the sum is all of EE. Note where each hypothesis worked: splitting filled the degree, equality of multiplicities filled the dimensions.

(4 \Rightarrow 1): nn distinct eigenvalues give nn independent eigenvectors (Theorem 3.2): a basis.

Method 3.7 (Deciding diagonalizability)

In practice, test in this order — each step may finish the job. (1) Does an annihilating polynomial with simple split roots present itself (u2=idu^2 = \mathrm{id}, u2=uu^2 = u, uk=idu^k = \mathrm{id})? If yes: diagonalizable, no computation (Corollary 3.17 below). (2) Compute χu\chi_u; if it has nn distinct roots in KK: diagonalizable (Theorem 3.6 (4)). (3) Otherwise, for each multiple root λ\lambda only, compare dimker(uλid)\dim\ker(u - \lambda\,\mathrm{id}) with the multiplicity mλm_\lambda: any deficit kills diagonalizability; equality everywhere proves it. Never compute eigenspaces of simple roots (their dimension is forced to be 11), and never trigonalize just to decide.

Example 3.8 (Diagonalization put to work)

A=I+J=(211121112)A = I + J = \left(\begin{smallmatrix}2 & 1 & 1\\ 1 & 2 & 1\\ 1 & 1 & 2\end{smallmatrix}\right), with JJ the all-ones matrix: from Sp(J)={3,0}\operatorname{Sp}(J) = \{3, 0\} (Example 2.19), Sp(A)={4,1}\operatorname{Sp}(A) = \{4, 1\}, with eigenspaces R(1,1,1)\R(1,1,1) and the plane {x+y+z=0}\{x + y + z = 0\}: dimensions 1+2=31 + 2 = 3, diagonalizable (Theorem 3.6 (2)). Powers without any change-of-basis matrix: with Π=J/3\Pi = J/3 the projector onto R(1,1,1)\R(1,1,1),

A=4Π+1(IΠ)Ak=4kΠ+(IΠ)=4k13J+I.A = 4\,\Pi + 1\cdot(I - \Pi) \quad\Longrightarrow\quad A^k = 4^k\,\Pi + (I - \Pi) = \frac{4^k - 1}{3}\,J + I .

(Check k=1k = 1: 413J+I=A\frac{4-1}3 J + I = A.) The closing insight: when the eigenspaces are visible, spectral projectors compute powers faster than PDP1PDP^{-1} ever will — and the formula displays the dynamics: AkA^k grows like 4k4^k along (1,1,1)(1,1,1) and stays put on the orthogonal plane.

Theorem 3.9 (Trigonalization)

uu is trigonalizable over KK if and only if χu\chi_u splits over KK. In particular every endomorphism of a C\C-vector space is trigonalizable.

Proof. (\Rightarrow) The characteristic polynomial of a triangular matrix is (Xtii)\prod(X - t_{ii}): split.

(\Leftarrow) Induction on nn. Since χu\chi_u splits, it has a root λ\lambda: pick an eigenvector e1e_1. In a basis starting with e1e_1, the matrix is (λ0B)\begin{pmatrix} \lambda & \ast\\ 0 & B\end{pmatrix}, and χu=(Xλ)χB\chi_u = (X - \lambda)\chi_B: χB\chi_B splits too. By the induction hypothesis applied to the (n1)×(n1)(n-1) \times (n-1) matrix BB, there is an invertible QQ with Q1BQQ^{-1}BQ upper triangular; conjugating the whole matrix by (100Q)\begin{pmatrix}1 & 0\\ 0 & Q\end{pmatrix} triangularizes it.

Example 3.10 (Trigonalizing by hand)

B=(3111)B = \begin{pmatrix}3 & -1\\ 1 & 1\end{pmatrix}: χB=X24X+4=(X2)2\chi_B = X^2 - 4X + 4 = (X - 2)^2, and ker(B2I)=ker(1111)\ker(B - 2I) = \ker\left(\begin{smallmatrix}1 & -1\\ 1 & -1\end{smallmatrix}\right) is the line spanned by e1=(1,1)e_1' = (1, 1): one eigenvalue, a one-dimensional eigenspace — not diagonalizable, but trigonalizable (Theorem 3.9). Complete the basis with e2=(1,0)e_2' = (1, 0) and compute:

u(e1)=(2,2)=2e1,u(e2)=(3,1)=1e1+2e2,u(e_1') = (2, 2) = 2e_1', \qquad u(e_2') = (3, 1) = 1\cdot e_1' + 2\, e_2' ,

so in the basis (e1,e2)(e_1', e_2') the matrix is T=(2102)T = \left(\begin{smallmatrix}2 & 1\\ 0 & 2\end{smallmatrix}\right). The closing insight: the diagonal of TT was forced (both entries must be the double eigenvalue 22); only the corner entry depended on the choice of e2e_2', and rescaling e2e_2' can make it any nonzero value — the resistant “11” is the shadow of the nilpotent part that Dunford will isolate.

3.2 Polynomials of an endomorphism

Definition 3.11

For P=akXkK[X]P = \sum a_k X^k \in K[X], set P(u)=akukL(E)P(u) = \sum a_k u^k \in \mathcal{L}(E). The map PP(u)P \mapsto P(u) is a morphism of algebras K[X]L(E)K[X] \to \mathcal{L}(E) (Definition 1.33); its kernel {P:P(u)=0}\{P : P(u) = 0\} is an ideal of K[X]K[X], nonzero (the family (id,u,,un2)(\mathrm{id}, u, \dots, u^{n^2}) is linked in the n2n^2-dimensional L(E)\mathcal{L}(E)), hence generated by a unique monic polynomial μu\mu_u: the minimal polynomial (Theorem 1.26).

Proposition 3.12

  1. P(u)=0    μuPP(u) = 0 \iff \mu_u \mid P; eigenvalues of uu are roots of every annihilating polynomial, and the roots of μu\mu_u are exactly the eigenvalues.
  2. If FF is stable, μuFμu\mu_{u|_F} \mid \mu_u.

Proof. (1) The divisibility is the definition of a generator. If u(x)=λxu(x) = \lambda x, x0x \neq 0, then 0=P(u)(x)=P(λ)x0 = P(u)(x) = P(\lambda)x, so P(λ)=0P(\lambda) = 0: eigenvalues are roots of annihilators, in particular of μu\mu_u. Conversely, if λ\lambda is a root, μu=(Xλ)Q\mu_u = (X - \lambda)Q with Q(u)0Q(u) \neq 0 (degree of μu\mu_u minimal): pick yy with Q(u)(y)0Q(u)(y) \neq 0; then (uλ)(Q(u)(y))=μu(u)(y)=0(u - \lambda)(Q(u)(y)) = \mu_u(u)(y) = 0 exhibits the eigenvector Q(u)(y)Q(u)(y).

(2) μu(uF)=μu(u)F=0\mu_u(u|_F) = \mu_u(u)|_F = 0, and apply (1) to uFu|_F.

Example 3.13 (Minimal polynomials found by hand)

The minimal polynomial is computed by testing successive degrees. For the all-ones matrix JM3(R)J \in \mathcal{M}_3(\R): JλIJ \neq \lambda I (degree 11 is out), and J2=3JJ^2 = 3J, so

μJ=X23X=X(X3):\mu_J = X^2 - 3X = X(X - 3) :

degree 22, split, simple roots — JJ is diagonalizable with spectrum {0,3}\{0, 3\} (Corollary 3.17 below), confirming Example 2.19 without a single determinant. For the swap matrix AA of Example 3.15: A±IA \neq \pm I and A2=IA^2 = I give μA=X21\mu_A = X^2 - 1. In both cases the pattern is the same: guess a low-degree identity from the structure (rank one forces J2=(trJ)JJ^2 = (\operatorname{tr}J)\,J; an involution forces A2=IA^2 = I), then check that no proper divisor annihilates. Minimal polynomials are usually found, not computed from χ\chi.

Theorem 3.14 (Kernel decomposition lemma)

If P=P1P2PrP = P_1 P_2 \cdots P_r with the PiP_i pairwise coprime, then

kerP(u)=kerP1(u)kerPr(u),\ker P(u) = \ker P_1(u) \oplus \dots \oplus \ker P_r(u),

and the projections onto the summands are polynomials in uu.

Proof. It suffices to treat r=2r = 2 and induct. Bézout in K[X]K[X] (Theorem 1.26): UP1+VP2=1U P_1 + V P_2 = 1, so for every xx,

x=U(u)P1(u)(x)=:x2+V(u)P2(u)(x)=:x1.x = \underbrace{U(u)P_1(u)(x)}_{=:\,x_2} + \underbrace{V(u)P_2(u)(x)}_{=:\,x_1}.

If xkerP(u)x \in \ker P(u): P2(u)(x2)=U(u)P(u)(x)=0P_2(u)(x_2) = U(u)\,P(u)(x) = 0 (polynomials in uu commute), so x2kerP2(u)x_2 \in \ker P_2(u), and symmetrically x1kerP1(u)x_1 \in \ker P_1(u): the sum fills kerP(u)\ker P(u); both summands sit inside kerP(u)\ker P(u) (PiPP_i \mid P). Directness: xkerP1(u)kerP2(u)x \in \ker P_1(u) \cap \ker P_2(u) gives x=U(u)P1(u)x+V(u)P2(u)x=0x = U(u)P_1(u)x + V(u)P_2(u)x = 0. The formulas for x1,x2x_1, x_2 exhibit the projections as V(u)P2(u)V(u)P_2(u) and U(u)P1(u)U(u)P_1(u).

Example 3.15 (The kernel lemma with explicit projectors)

Let A=(010100001)A = \left(\begin{smallmatrix}0 & 1 & 0\\ 1 & 0 & 0\\ 0 & 0 & 1\end{smallmatrix}\right) (swap the first two coordinates). Then A2=IA^2 = I: the polynomial X21=(X1)(X+1)X^2 - 1 = (X - 1)(X + 1) annihilates AA, its factors are coprime, and Bézout is explicit:

12(X+1)12(X1)=1.\frac{1}{2}(X + 1) - \frac12(X - 1) = 1 .

Following the proof of Theorem 3.14, the projections onto ker(AI)\ker(A - I) and ker(A+I)\ker(A + I) are the polynomials in AA

π+=A+I2=12(110110002),π=IA2=12(110110000).\pi_+ = \frac{A + I}{2} = \frac12\begin{pmatrix} 1 & 1 & 0\\ 1 & 1 & 0\\ 0 & 0 & 2\end{pmatrix}, \qquad \pi_- = \frac{I - A}{2} = \frac12\begin{pmatrix} 1 & -1 & 0\\ -1 & 1 & 0\\ 0 & 0 & 0\end{pmatrix}.

Check: π++π=I\pi_+ + \pi_- = I, π+π=0\pi_+\pi_- = 0, π±2=π±\pi_\pm^2 = \pi_\pm, and the images are the plane {x=y}\{x = y\} (symmetric vectors, eigenvalue 11) and the line R(1,1,0)\R(1, -1, 0) (antisymmetric, eigenvalue 1-1). The kernel lemma is not an existence statement: Bézout coefficients are the projector formulas.

Example 3.16 (Projectors compute the exponential too)

The same swap matrix, one dividend further. Since A=π+πA = \pi_+ - \pi_- with orthogonal-in-the-algebraic-sense projectors (π+π=0\pi_+\pi_- = 0), every power obeys Ak=π++(1)kπA^k = \pi_+ + (-1)^k\pi_-, and the exponential series regroups by projector:

etA=ktkk!(π++(1)kπ)=etπ++etπ=(coshtsinht0sinhtcosht000et).\eu^{tA} = \sum_k \frac{t^k}{k!}\bigl(\pi_+ + (-1)^k\pi_-\bigr) = \eu^{t}\,\pi_+ + \eu^{-t}\,\pi_- = \begin{pmatrix} \cosh t & \sinh t & 0\\ \sinh t & \cosh t & 0\\ 0 & 0 & \eu^{t} \end{pmatrix}.

(Check t=0t = 0: the identity; derivative at 00: AA.) The eigendecomposition converts a matrix series into two scalar series — the exact mechanism that Chapter 16 will run on every diagonalizable system, and the reason hyperbolic functions govern symmetric couplings.

Corollary 3.17 (Diagonalizability via the minimal polynomial)

uu is diagonalizable     \iff μu\mu_u splits over KK with simple roots     \iff some annihilating polynomial of uu splits with simple roots.

Proof. If P(u)=0P(u) = 0 with P=i(Xλi)P = \prod_{i}(X - \lambda_i) (distinct λi\lambda_i), the lemma gives E=kerP(u)=iker(uλi)E = \ker P(u) = \bigoplus_i \ker(u - \lambda_i): a direct sum of eigenspaces, so uu is diagonalizable (Theorem 3.6). Conversely, a diagonalizable uu is killed by λSpu(Xλ)\prod_{\lambda \in \operatorname{Sp}u}(X - \lambda) (it kills each eigenspace), which splits with simple roots; and μu\mu_u divides it while having the same roots (Proposition 3.12): μu\mu_u is exactly that product.

Example 3.18

Projections satisfy p2=pp^2 = p: annihilated by X(X1)X(X-1), split simple roots — diagonalizable with spectrum {0,1}\subseteq \{0, 1\}, and E=kerpker(pid)E = \ker p \oplus \ker(p - \mathrm{id}): Year 1’s geometric analysis, re-proved in one line. Symmetries (s2=ids^2 = \mathrm{id}, annihilator X21X^2 - 1): diagonalizable when charK2\operatorname{char} K \neq 2, spectrum {±1}\subseteq \{\pm 1\}. An endomorphism with u3=u2u^3 = u^2 and u2uu^2 \neq u: annihilated by X2(X1)X^2(X - 1), not necessarily diagonalizable — the criterion detects it (double root 00 must be tested: diagonalizable iff moreover keru2=keru\ker u^2 = \ker u).

Example 3.19 (The field decides: a rotation in R3\R^3)

Let RR be the quarter-turn around the zz-axis:

R=(010100001),χR=(X1)(X2+1).R = \begin{pmatrix} 0 & -1 & 0\\ 1 & 0 & 0\\ 0 & 0 & 1 \end{pmatrix}, \qquad \chi_R = (X - 1)(X^2 + 1).

Over R\R: the only eigenvalue is 11, with eigenspace the axis Re3\R e_3 — one line of fixed vectors, and no further reduction: RR is neither diagonalizable nor trigonalizable in M3(R)\mathcal{M}_3(\R) (χR\chi_R does not split). Over C\C: three distinct eigenvalues 1,i,i1, \iu, -\iu, so RR is diagonalizable, with eigenvectors e3e_3 and e1ie2e_1 \mp \iu e_2. The geometry was audible in the algebra: rotations in the plane have no real invariant directions, and the complex eigenvalues ±i\pm\iu of modulus 11 store the angle (±π2\pm\frac\pi2) that the real matrix can only express by mixing coordinates.

Example 3.20 (Minimal versus characteristic)

For D=diag(2,2,3)D = \operatorname{diag}(2, 2, 3): χD=(X2)2(X3)\chi_D = (X - 2)^2(X - 3) but μD=(X2)(X3)\mu_D = (X - 2)(X - 3), since (D2I)(D3I)=0(D - 2I)(D - 3I) = 0 (check on the canonical basis) while neither factor alone kills DD. For the shift block N=(0100)(3)N = \left(\begin{smallmatrix}0 & 1\\ 0 & 0\end{smallmatrix}\right) \oplus (3), i.e. N=(010000003)N' = \left(\begin{smallmatrix}0 & 1 & 0\\ 0 & 0 & 0\\ 0 & 0 & 3\end{smallmatrix}\right): χN=X2(X3)\chi_{N'} = X^2(X - 3) and μN=X2(X3)\mu_{N'} = X^2(X - 3) — the double root is genuinely needed because NN' is not diagonalizable on ker\ker-side (Ne2=e10N'e_2 = e_1 \neq 0). Rule of thumb: μ\mu and χ\chi share their roots (Proposition 3.12); the multiplicity in μ\mu measures the size of the largest nilpotent block, the one in χ\chi the total dimension of the characteristic subspace.

Theorem 3.21 (Cayley–Hamilton)

χu(u)=0\chi_u(u) = 0; consequently μuχu\mu_u \mid \chi_u, and degμun\deg \mu_u \leq n.

Proof. Fix x0x \neq 0 and let dd be maximal with (x,u(x),,ud1(x))(x, u(x), \dots, u^{d-1}(x)) free; write

ud(x)=a0xa1u(x)ad1ud1(x),u^d(x) = -a_0 x - a_1 u(x) - \dots - a_{d-1}u^{d-1}(x),

and set Px=Xd+ad1Xd1++a0P_x = X^d + a_{d-1}X^{d-1} + \dots + a_0, so Px(u)(x)=0P_x(u)(x) = 0. Complete the free family into a basis of EE: on it, uu has block form (C0D)\begin{pmatrix} C & \ast\\ 0 & D\end{pmatrix} where CC is the companion matrix of PxP_x, whose characteristic polynomial is PxP_x (expand det(XIC)\det(XI - C) along the first column, by induction on dd). Hence χu=PxχD\chi_u = P_x \cdot \chi_D, and

χu(u)(x)=χD(u)(Px(u)(x))=0.\chi_u(u)(x) = \chi_D(u)\bigl(P_x(u)(x)\bigr) = 0 .

The argument holds for every xx: χu(u)=0\chi_u(u) = 0.

Example 3.22 (Cayley–Hamilton at work)

A=(1234)A = \begin{pmatrix} 1 & 2\\ 3 & 4\end{pmatrix}: χA=X25X2\chi_A = X^2 - 5X - 2, so A2=5A+2IA^2 = 5A + 2I. Every power of AA collapses to a combination of II and AA:

A4=(5A+2I)2=25A2+20A+4I=145A+54I=(199290435634),A^4 = (5A + 2I)^2 = 25A^2 + 20A + 4I = 145A + 54I = \begin{pmatrix} 199 & 290\\ 435 & 634\end{pmatrix},

and the inverse comes free: A(A5I)=2IA(A - 5I) = 2I gives

A1=12(A5I)=(213/21/2).A^{-1} = \tfrac12(A - 5I) = \begin{pmatrix} -2 & 1\\ 3/2 & -1/2\end{pmatrix}.

The closing insight: Cayley–Hamilton compresses the whole algebra K[A]K[A] into Vect(I,A,,An1)\operatorname{Vect}(I, A, \dots, A^{n-1})dimK[A]=degμAn\dim K[A] = \deg\mu_A \leq n, however large the powers you need.

Remark 3.23 (Common pitfalls)

(i) Eigenvalues do not add: Sp(A+B)\operatorname{Sp}(A + B) is not SpA+SpB\operatorname{Sp}A + \operatorname{Sp}B, and a sum of diagonalizable matrices need not be diagonalizable(1100)+(0001)=(1101)\left(\begin{smallmatrix}1 & 1\\ 0 & 0\end{smallmatrix}\right) + \left(\begin{smallmatrix}0 & 0\\ 0 & 1\end{smallmatrix}\right) = \left(\begin{smallmatrix}1 & 1\\ 0 & 1\end{smallmatrix}\right) is a sum of two diagonalizable matrices (each has distinct eigenvalues) and is not diagonalizable; only commuting families behave (Exercise 3.9). (ii) “χu\chi_u splits” is a hypothesis about the field: a plane rotation has χ=X22cosθX+1\chi = X^2 - 2\cos\theta\,X + 1, split over C\C, not over R\Rdiagonalizable in M2(C)\mathcal{M}_2(\C), not trigonalizable in M2(R)\mathcal{M}_2(\R). (iii) The inequality runs geometric \leq algebraic, never the reverse; testing only dimEλ1\dim E_\lambda \geq 1 proves nothing about diagonalizability. (iv) μu\mu_u is not χu\chi_u: equality holds exactly when each eigenvalue has a single block chain (e.g. companion matrices, this chapter’s weekend problem); using χ\chi where μ\mu is needed inflates every power computation. (v) Dunford’s dd and ν\nu are polynomials in uu — a decomposition u=d+νu = d' + \nu' with the right properties but dννdd'\nu' \neq \nu'd' is not Dunford and is never unique.

Remark 3.24 (Where this chapter is used)

Reduction is the workhorse of the rest of the book: powers and exponentials of matrices drive the linear differential systems of Chapter 16; the spectral theorem of Chapter 12 is diagonalization made orthogonal; generating functions (Chapter 23) re-derive the recurrence asymptotics of this chapter’s weekend problem analytically. In the Year 3 volume the same program runs in infinite dimension: the spectral theory of compact self-adjoint operators, where eigenvalue sequences replace finite spectra, and the Perron–Frobenius theory of positive matrices, which explains why dominant eigenvalues of counting problems are positive and simple.

3.3 Nilpotents and the Dunford decomposition

Proposition 3.25 (Nilpotent endomorphisms)

For uu with χu\chi_u split, the following are equivalent: un=0u^n = 0; uk=0u^k = 0 for some kk; Sp(u)={0}\operatorname{Sp}(u) = \{0\}; χu=Xn\chi_u = X^n; uu is trigonalizable with zero diagonal. A nilpotent endomorphism has μu=X(nilpotence index)\mu_u = X^{\text{(nilpotence index)}}, and index n\leq n.

Proof. uk=0u^k = 0 makes every eigenvalue a root of XkX^k: spectrum {0}\{0\} (nonempty when χ\chi splits — over C\C always). Then χu=Xn\chi_u = X^n (all roots zero) and Cayley–Hamilton gives un=0u^n = 0; trigonalization (Theorem 3.9) puts zeros on the diagonal (the diagonal carries the eigenvalues). Conversely, let AA be strictly upper triangular: aij=0a_{ij} = 0 for jij \leq i. We show by induction that

(Ak)ij=0whenever ji+k1,(A^k)_{ij} = 0 \qquad \text{whenever } j \leq i + k - 1,

i.e. each power pushes the zero region one diagonal higher. For k=1k = 1 this is the hypothesis. For the step,

(Ak+1)ij=(Ak)iaj,(A^{k+1})_{ij} = \sum_{\ell} (A^k)_{i\ell}\,a_{\ell j} ,

and each term vanishes: either i+k1\ell \leq i + k - 1 (the first factor is 00 by induction) or i+k\ell \geq i + k, in which case ji+kj \leq i + k \leq \ell kills the second factor. At k=nk = n the condition ji+n1j \leq i + n - 1 holds for all i,jni, j \leq n: An=0A^n = 0. The minimal polynomial divides XnX^n and annihilation defines the index.

Theorem 3.26 (Dunford decomposition)

Suppose χu\chi_u splits over KK (automatic for K=CK = \C). Then there is a unique pair (d,ν)(d, \nu) with

u=d+ν,d diagonalizable,ν nilpotent,dν=νd,u = d + \nu, \qquad d \text{ diagonalizable}, \quad \nu \text{ nilpotent}, \quad d\nu = \nu d ,

and moreover dd and ν\nu are polynomials in uu.

Proof. Existence. Write χu=i=1r(Xλi)mi\chi_u = \prod_{i=1}^{r} (X - \lambda_i)^{m_i} (distinct λi\lambda_i) and set Ni=ker(uλi)miN_i = \ker(u - \lambda_i)^{m_i}, the characteristic subspaces. By Cayley–Hamilton and the kernel lemma (Theorem 3.14),

E=N1Nr,E = N_1 \oplus \dots \oplus N_r ,

with projections πi\pi_i polynomial in uu; each NiN_i is stable (polynomials in uu commute with uu). Define d=iλiπid = \sum_i \lambda_i \pi_i: a polynomial in uu, diagonalizable (it acts as λi\lambda_i on NiN_i, so EE decomposes into its eigenspaces). Then ν=ud\nu = u - d is a polynomial in uu (hence commutes with dd), and on each NiN_i it acts as uλiu - \lambda_i, with (uλi)mi=0(u - \lambda_i)^{m_i} = 0 there: νmaxmi=0\nu^{\max m_i} = 0 on each summand, so ν\nu is nilpotent.

Uniqueness. Let u=d+νu = d' + \nu' be another such pair. Since dd' and ν\nu' commute with each other, they commute with u=d+νu = d' + \nu', hence with every polynomial in uu — in particular with dd and ν\nu. Then ddd - d' is diagonalizable (two commuting diagonalizable maps are simultaneously diagonalizable: Exercise 3.9) and equals νν\nu' - \nu, which is nilpotent: if νk=0\nu^k = 0 and νk=0\nu'^{k'} = 0, commutation licenses the binomial expansion

(νν)k+k1=j=0k+k1(k+k1j)νj(ν)k+k1j,(\nu' - \nu)^{k + k' - 1} = \sum_{j=0}^{k+k'-1}\binom{k + k' - 1}{j} \,\nu'^{\,j}\,(-\nu)^{k + k' - 1 - j} ,

in which every term dies: either jkj \geq k' (first factor zero) or k+k1jkk + k' - 1 - j \geq k (second factor zero), and one of the two always holds. A diagonalizable nilpotent is zero (its spectrum is {0}\{0\} and it is diagonal in some basis): d=dd = d', ν=ν\nu = \nu'.

Example 3.27 (Powers and exponentials)

A=(3111)A = \begin{pmatrix} 3 & 1\\ -1 & 1\end{pmatrix}: χA=X24X+4=(X2)2\chi_A = X^2 - 4X + 4 = (X-2)^2, single eigenvalue 22, eigenspace of dimension 11: not diagonalizable. Dunford: D=2ID = 2I, N=A2I=(1111)N = A - 2I = \begin{pmatrix} 1 & 1\\ -1 & -1\end{pmatrix}, N2=0N^2 = 0. Then

Ak=(2I+N)k=2kI+k2k1N,etA=e2t(I+tN),A^k = (2I + N)^k = 2^k I + k\,2^{k-1} N , \qquad \eu^{tA} = \eu^{2t}(I + tN),

by the commuting binomial theorem, resp. the exponential series (Chapter 16) split on commuting summands. Reduction turns matrix dynamics into scalar dynamics.

Remark 3.28 (Perspectives within this volume)

Reduction is a hub; here are the four spokes to watch. In Chapter 5, adapted norms turn “all eigenvalues of modulus <1< 1” into “some operator norm <1< 1”, making spectra govern convergence of powers and series. In Chapter 16, the recipe of Example 3.27 becomes the general solution of X=AXX' = AX: Dunford splits etA\eu^{tA} into polynomial-times-exponential blocks, and stability reads off the real parts of eigenvalues. In Chapter 12, a scalar product forces what mere linear algebra cannot: symmetric matrices become orthogonally diagonalizable, with no nilpotent part at all. And in Chapter 23, the dominant-eigenvalue asymptotics of this chapter’s weekend problem reappear analytically, as the smallest singularity of a generating function — two languages for one growth rate.

3.4 Exercises

Exercise 3.1

Diagonalize (eigenvalues, bases of eigenspaces, invertible PP):

A=(1221),B=(011101110).A = \begin{pmatrix} 1 & 2\\ 2 & 1 \end{pmatrix}, \qquad B = \begin{pmatrix} 0 & 1 & 1\\ 1 & 0 & 1\\ 1 & 1 & 0 \end{pmatrix}.
Solution

Solution of Exercise 3.1.

AA: χA=X22X3=(X3)(X+1)\chi_A = X^2 - 2X - 3 = (X - 3)(X + 1). Eigenvectors: for 33: (1,1)(1,1); for 1-1: (1,1)(1,-1). So P=(1111)P = \begin{pmatrix} 1 & 1\\ 1 & -1\end{pmatrix} gives P1AP=diag(3,1)P^{-1}AP = \operatorname{diag}(3, -1).

B=JIB = J - I where JJ is the all-ones matrix. JJ has rank 11 with Jv=3vJv = 3v for v=(1,1,1)v = (1,1,1) and Jw=0Jw = 0 on the plane x+y+z=0x + y + z = 0: spectrum of BB is {2,1}\{2, -1\} with eigenspaces Vect(1,1,1)\operatorname{Vect}(1,1,1) (dimension 11) and {x+y+z=0}\{x + y + z = 0\} (dimension 22, basis (1,1,0),(1,0,1)(1,-1,0), (1,0,-1)). PP with these three columns gives P1BP=diag(2,1,1)P^{-1}BP = \operatorname{diag}(2, -1, -1).

Exercise 3.2

Show that C=(1101)C = \begin{pmatrix} 1 & 1\\ 0 & 1\end{pmatrix} is not diagonalizable, twice: via eigenspaces, and via the minimal polynomial.

Solution

Solution of Exercise 3.2.

Eigenspaces: χC=(X1)2\chi_C = (X-1)^2, single eigenvalue 11; ker(CI)=ker(0100)\ker(C - I) = \ker\begin{pmatrix} 0&1\\ 0&0\end{pmatrix} is the line Vect(e1)\operatorname{Vect}(e_1): dimension 1<2=m11 < 2 = m_1, so not diagonalizable (Theorem 3.6).

Minimal polynomial: μC\mu_C divides (X1)2(X-1)^2 and CIC \neq I, so μC=(X1)2\mu_C = (X-1)^2: a double root, so not diagonalizable (Corollary 3.17).

Exercise 3.3

Let uu satisfy u25u+6id=0u^2 - 5u + 6\,\mathrm{id} = 0. Prove that uu is diagonalizable, determine the possible spectra, and compute uku^k as a combination of id\mathrm{id} and uu.

Solution

Solution of Exercise 3.3.

X25X+6=(X2)(X3)X^2 - 5X + 6 = (X-2)(X-3): split with simple roots, so uu is diagonalizable (Corollary 3.17), with Sp(u){2,3}\operatorname{Sp}(u) \subseteq \{2, 3\}. Possible spectra: {2}\{2\} (u=2idu = 2\,\mathrm{id}), {3}\{3\} (u=3idu = 3\,\mathrm{id}), or {2,3}\{2, 3\}.

Powers: seek uk=akid+bkuu^k = a_k\,\mathrm{id} + b_k\,u. On the eigenspaces, this reads 2k=ak+2bk2^k = a_k + 2b_k and 3k=ak+3bk3^k = a_k + 3b_k: solving, bk=3k2kb_k = 3^k - 2^k, ak=32k23ka_k = 3\cdot2^k - 2\cdot 3^k:

uk=(32k23k)id+(3k2k)u.u^k = (3\cdot 2^k - 2\cdot 3^k)\,\mathrm{id} + (3^k - 2^k)\, u .

(Valid for all three spectra: the identities hold eigenvalue-wise.)

Exercise 3.4 ★★

Let uu be diagonalizable and FF a stable subspace. Prove that uFu|_F is diagonalizable (restrict an annihilating polynomial with simple split roots).

Solution

Solution of Exercise 3.4.

uu diagonalizable: P=λ(Xλ)P = \prod_{\lambda}(X - \lambda) over the spectrum annihilates uu, splits, simple roots. Then P(uF)=P(u)F=0P(u|_F) = P(u)|_F = 0: the restriction is annihilated by a split polynomial with simple roots, hence diagonalizable (Corollary 3.17).

Exercise 3.5 ★★

(Fibonacci) Let A=(1110)A = \begin{pmatrix} 1 & 1\\ 1 & 0\end{pmatrix}. Diagonalize AA over R\R, and deduce Binet’s formula for the Fibonacci sequence (F0=0F_0 = 0, F1=1F_1 = 1, Fn+1=Fn+Fn1F_{n+1} = F_n + F_{n-1}):

Fn=φnψn5,φ=1+52, ψ=152.F_n = \frac{\varphi^n - \psi^n}{\sqrt 5}, \qquad \varphi = \frac{1 + \sqrt5}{2},\ \psi = \frac{1 - \sqrt5}{2}.
Solution

Solution of Exercise 3.5.

χA=X2X1\chi_A = X^2 - X - 1, roots φ\varphi and ψ\psi (distinct): diagonalizable, with eigenvectors (φ,1)(\varphi, 1) and (ψ,1)(\psi, 1). The recurrence gives (Fn+1Fn)=An(10)\begin{pmatrix} F_{n+1}\\ F_n \end{pmatrix} = A^n \begin{pmatrix}1\\ 0\end{pmatrix}. Decompose (1,0)(1, 0) on the eigenvectors: (1,0)=1φψ((φ,1)(ψ,1))(1,0) = \frac{1}{\varphi - \psi}\bigl((\varphi, 1) - (\psi, 1)\bigr) with φψ=5\varphi - \psi = \sqrt5. Applying AnA^n multiplies each eigencomponent by its eigenvalue’s nn-th power; reading the second coordinate:

Fn=φnψn5.F_n = \frac{\varphi^n - \psi^n}{\sqrt 5} .

(Check: n=1n = 1 gives φψ5=1\frac{\varphi - \psi}{\sqrt5} = 1.)

Exercise 3.6 ★★

Let uL(E)u \in \mathcal{L}(E) with u2u^2 diagonalizable and uu invertible (K=CK = \C). Prove that uu is diagonalizable. Give a counterexample when uu is not invertible.

Solution

Solution of Exercise 3.6.

Let P=i(Xμi)P = \prod_i (X - \mu_i) annihilate u2u^2, split with simple roots μi\mu_i (the spectrum of u2u^2). Since uu is invertible, 00 is not an eigenvalue of u2u^2 (detu2=(detu)20\det u^2 = (\det u)^2 \neq 0), so all μi0\mu_i \neq 0. Then

Q(X)=i(X2μi)=i(Xμi)(X+μi)Q(X) = \prod_i (X^2 - \mu_i) = \prod_i (X - \sqrt{\mu_i})(X + \sqrt{\mu_i})

annihilates uu:   Q(u)=i(u2μiid)=P(u2)=0\;Q(u) = \prod_i (u^2 - \mu_i\,\mathrm{id}) = P(u^2) = 0. Its roots ±μi\pm \sqrt{\mu_i} (complex square roots) are pairwise distinct because the μi\mu_i are distinct and nonzero (μi=μj\sqrt{\mu_i} = -\sqrt{\mu_j} would give μi=μj\mu_i = \mu_j). Split + simple roots: uu is diagonalizable.

Counterexample without invertibility: u=(0100)u = \begin{pmatrix} 0 & 1\\ 0 & 0\end{pmatrix}: u2=0u^2 = 0 is diagonalizable, uu is not.

Exercise 3.7 ★★

Compute the Dunford decomposition, AkA^k, and etA\eu^{tA} for

A=(210021002).A = \begin{pmatrix} 2 & 1 & 0\\ 0 & 2 & 1\\ 0 & 0 & 2 \end{pmatrix}.
Solution

Solution of Exercise 3.7.

A=2I+NA = 2I + N with NN the shift (Ne2=e1N e_2 = e_1, Ne3=e2Ne_3 = e_2), N3=0N^3 = 0, N2=E13N^2 = E_{13}: this is the Dunford decomposition (2I2I diagonal, NN nilpotent, they commute; uniqueness makes it the one). Binomial with commuting terms:

Ak=2kI+k2k1N+(k2)2k2N2=(2kk2k1(k2)2k202kk2k1002k),A^k = 2^k I + k 2^{k-1} N + \binom k2 2^{k-2} N^2 = \begin{pmatrix} 2^k & k2^{k-1} & \binom k2 2^{k-2}\\ 0 & 2^k & k2^{k-1}\\ 0 & 0 & 2^k \end{pmatrix},
etA=e2t(I+tN+t22N2)=e2t(1tt2/201t001).\eu^{tA} = \eu^{2t}\Bigl(I + tN + \frac{t^2}{2}N^2\Bigr) = \eu^{2t}\begin{pmatrix} 1 & t & t^2/2\\ 0 & 1 & t\\ 0 & 0 & 1 \end{pmatrix}.

Exercise 3.8 ★★

Let AMn(C)A \in \mathcal{M}_n(\C) with Ak=IA^k = I for some k1k \geq 1. Prove that AA is diagonalizable and its eigenvalues are kk-th roots of unity. Deduce that a finite-order invertible complex matrix similar to a triangular matrix with unit diagonal is the identity.

Solution

Solution of Exercise 3.8.

Xk1X^k - 1 annihilates AA and splits over C\C with the kk distinct roots e2iπj/k\eu^{2\iu\pi j/k}: AA is diagonalizable (Corollary 3.17) and its eigenvalues, roots of Xk1X^k - 1, are kk-th roots of unity.

If moreover AA is similar to a triangular matrix with unit diagonal: all eigenvalues equal 11, and AA, diagonalizable with sole eigenvalue 11, is PIP1=IP\,I\,P^{-1} = I.

Exercise 3.9 ★★★

(Simultaneous diagonalization) Let u,vu, v be diagonalizable and commuting. Prove that they are simultaneously diagonalizable: some basis diagonalizes both. (Each eigenspace of uu is vv-stable; diagonalize the restrictions of vv there, using Exercise 3.4.)

Solution

Solution of Exercise 3.9.

Write E=λEλ(u)E = \bigoplus_\lambda E_\lambda(u) (Theorem 3.6). Each Eλ(u)E_\lambda(u) is vv-stable: for xEλx \in E_\lambda, u(v(x))=v(u(x))=λv(x)u(v(x)) = v(u(x)) = \lambda v(x). The restriction of vv to Eλ(u)E_\lambda(u) is diagonalizable (Exercise 3.4): choose a basis of Eλ(u)E_\lambda(u) made of vv-eigenvectors. Concatenating these bases over all λ\lambda gives a basis of EE whose vectors are eigenvectors of both uu (by membership in Eλ(u)E_\lambda(u)) and vv (by construction).

Exercise 3.10 ★★★

Let uL(Cn)u \in \mathcal{L}(\C^n). Prove that uu is diagonalizable if and only if every uu-stable subspace has a uu-stable supplementary subspace. (For \Leftarrow: apply the property to F=λEλ(u)F = \sum_\lambda E_\lambda(u), the sum of all eigenspaces; if a stable supplement GG were nonzero, trigonalizing uGu|_G would produce an eigenvector of uu inside GG — contradicting GF={0}G \cap F = \{0\}.)

Solution

Solution of Exercise 3.10.

(\Rightarrow) Let uu be diagonalizable and FF stable. Then uFu|_F is diagonalizable (Exercise 3.4): FF has a basis of eigenvectors, which extends, inside each global eigenspace EλE_\lambda, to a basis of EλE_\lambda (incomplete basis theorem within EλE_\lambda, starting from the part of FF’s basis lying there — note F=λ(FEλ)F = \bigoplus_\lambda (F \cap E_\lambda) since uFu|_F is diagonalizable). The added vectors span a stable supplement (each lies in some EλE_\lambda, so their span is uu-stable).

(\Leftarrow) Let F=λEλ(u)F = \sum_\lambda E_\lambda(u) (a stable subspace) and GG a stable supplement. If G{0}G \neq \{0\}: χuG\chi_{u|_G} splits over C\C, so uGu|_G has an eigenvector xGx \in G (Theorem 3.9 or directly the existence of a root); but every eigenvector of uu lies in FF, so xFG={0}x \in F \cap G = \{0\}: contradiction. Hence G={0}G = \{0\} and E=FE = F: the eigenspaces fill EE, i.e. uu is diagonalizable.

Exercise 3.11 ★★★

(Spectral radius by Gelfand-lite, 2×22\times2 taste of analysis to come) Let AM2(C)A \in \mathcal{M}_2(\C) with both eigenvalues of modulus <1< 1. Prove that Ak0A^k \to 0 entrywise as kk \to \infty. (Trigonalize: A=P(T)P1A = P(T)P^{-1} with TT upper triangular; compute TkT^k explicitly — distinguish equal and distinct eigenvalues — and bound.)

Solution

Solution of Exercise 3.11.

Trigonalize: A=PTP1A = PTP^{-1}, T=(λc0μ)T = \begin{pmatrix} \lambda & c\\ 0 & \mu\end{pmatrix}, λ,μ<1\abs\lambda, \abs\mu < 1. Then Ak=PTkP1A^k = PT^kP^{-1}, and it suffices that Tk0T^k \to 0.

Distinct eigenvalues: induction gives

Tk=(λkcλkμkλμ0μk),T^k = \begin{pmatrix} \lambda^k & c\,\dfrac{\lambda^k - \mu^k}{\lambda - \mu}\\[4pt] 0 & \mu^k \end{pmatrix},

and each entry tends to 00 (λk,μk0\abs{\lambda}^k, \abs\mu^k \to 0).

Equal eigenvalues (μ=λ\mu = \lambda): T=λI+cE12T = \lambda I + cE_{12} and Tk=λkI+kλk1cE12T^k = \lambda^k I + k\lambda^{k-1}cE_{12}; the entry kλk10k\lambda^{k-1} \to 0 since λ<1\abs\lambda < 1 (geometric beats polynomial). In both cases Tk0T^k \to 0 entrywise, hence Ak=PTkP10A^k = PT^kP^{-1} \to 0 (matrix multiplication by fixed P,P1P, P^{-1} is continuous in the entries — each entry of the product is a fixed linear combination).

Exercise 3.12 ★★

Let uL(Cn)u \in \mathcal{L}(\C^n) with rku=1\operatorname{rk} u = 1 (n2n \geq 2). Show that χu=Xn1(Xtru)\chi_u = X^{n-1}(X - \operatorname{tr} u), and that uu is diagonalizable if and only if tru0\operatorname{tr} u \neq 0. (Recall from Exercise 2.5 that u2=(tru)uu^2 = (\operatorname{tr} u)\,u.)

Solution

Solution of Exercise 3.12.

keru\ker u has dimension n1n - 1 (rank–nullity), so 00 is an eigenvalue of geometric multiplicity n1n - 1, and χu\chi_u is divisible by Xn1X^{n-1} (Definition 3.3: geometric \leq algebraic). Write χu=Xn1(Xα)\chi_u = X^{n-1}(X - \alpha); the coefficient of Xn1X^{n-1} being tru-\operatorname{tr} u, we get α=tru\alpha = \operatorname{tr} u: χu=Xn1(Xtru)\chi_u = X^{n-1}(X - \operatorname{tr} u).

If tru0\operatorname{tr} u \neq 0: the eigenvalue tru\operatorname{tr} u is a root of χu\chi_u, so it carries an eigenvector; the eigenspaces for 00 and tru\operatorname{tr} u have dimensions n1n - 1 and 1\geq 1, summing to n\geq n: they fill EE, and uu is diagonalizable (Theorem 3.6). If tru=0\operatorname{tr} u = 0: by Exercise 2.5, u2=(tru)u=0u^2 = (\operatorname{tr} u)u = 0 with u0u \neq 0: uu is a nonzero nilpotent, and a diagonalizable nilpotent is zero (Proposition 3.25): not diagonalizable.

3.5 Problem: Linear Recurrences and Companion Matrices

A linear recurrence un+k=ak1un+k1++a0unu_{n+k} = a_{k-1}u_{n+k-1} + \dots + a_0 u_n is a matrix power in disguise, and reduction converts it into closed formulas, growth rates, and error estimates. This weekend problem develops the dictionary — companion matrices on one side, the shift operator on the space of sequences on the other — proves the fundamental theorem of linear recurrences (the general solution is iQi(n)λin\sum_i Q_i(n)\lambda_i^n over the roots of the characteristic polynomial), and spends the dividends on Diophantine approximation of 2\sqrt2, on the counting of walks and words, and on a ring of coupled sequences that only simultaneous diagonalization can untangle.

Problem 3.1

Weekend problem — the fundamental theorem of linear recurrences

Fix k1k \geq 1, scalars a0,,ak1Ca_0, \dots, a_{k-1} \in \C with a00a_0 \neq 0, the monic polynomial P=Xkak1Xk1a1Xa0P = X^k - a_{k-1}X^{k-1} - \dots - a_1 X - a_0, and the recurrence

(R) ⁣:un+k=ak1un+k1++a1un+1+a0un(n0).(\mathcal R)\colon\quad u_{n+k} = a_{k-1}u_{n+k-1} + \dots + a_1 u_{n+1} + a_0 u_n \qquad (n \geq 0).

The companion matrix of PP is

C=(0101a0a1ak1)Mk(C).C = \begin{pmatrix} 0 & 1 & & \\ & \ddots & \ddots & \\ & & 0 & 1\\ a_0 & a_1 & \cdots & a_{k-1} \end{pmatrix} \in \mathcal{M}_k(\C).

Part I — The companion dictionary.

  1. Show that a sequence (un)(u_n) satisfies (R)(\mathcal R) if and only if the vectors vn=(un,un+1,,un+k1)Tv_n = (u_n, u_{n+1}, \dots, u_{n+k-1})^{\mathsf T} satisfy vn+1=Cvnv_{n+1} = Cv_n, hence vn=Cnv0v_n = C^n v_0.
  2. Prove that χC=P\chi_C = P (expand det(XIC)\det(XI - C) along the first column and induct on kk), then that μC=P\mu_C = P as well (pass to CTC^{\mathsf T}, for which e1e_1 is cyclic, and note that a matrix and its transpose have the same minimal polynomial).
  3. Show that for each root λ\lambda of PP, the vector (1,λ,,λk1)T(1, \lambda, \dots, \lambda^{k-1})^{\mathsf T} spans the eigenspace of CC for λ\lambda; deduce that every eigenspace of CC has dimension 11, and that CC is diagonalizable if and only if PP has kk distinct roots.
  4. Suppose PP has distinct roots λ1,,λk\lambda_1, \dots, \lambda_k. Show that the geometric sequences (λin)n(\lambda_i^n)_n form a basis of the space of solutions of (R)(\mathcal R), so every solution is un=iciλinu_n = \sum_i c_i\lambda_i^n for unique constants cic_i.
  5. Solve completely: un+2=un+1+6unu_{n+2} = u_{n+1} + 6u_n, u0=1u_0 = 1, u1=8u_1 = 8.

Part II — The shift operator and the fundamental theorem. Let S\mathcal{S} be the C\C-vector space of all complex sequences and SL(S)S \in \mathcal{L}(\mathcal{S}) the shift, S((un)n)=(un+1)nS\bigl((u_n)_n\bigr) = (u_{n+1})_n.

  1. Show that the solution set of (R)(\mathcal R) is kerP(S)\ker P(S), and that it has dimension exactly kk (map a solution to its initial values).
  2. Explain why the kernel decomposition lemma (Theorem 3.14) applies to SS on the infinite-dimensional S\mathcal{S} without any change, and write the resulting decomposition of kerP(S)\ker P(S) for P=i=1r(Xλi)miP = \prod_{i=1}^{r}(X - \lambda_i)^{m_i} (distinct λi\lambda_i, all nonzero since a00a_0 \neq 0).
  3. For λ0\lambda \neq 0 and m1m \geq 1, show

    ker(Sλid)m={(Q(n)λn)n:QCm1[X]},\ker\,(S - \lambda\,\mathrm{id})^m = \bigl\{\,\bigl(Q(n)\,\lambda^n\bigr)_n : Q \in \C_{m-1}[X]\,\bigr\},

    of dimension mm. (Compute (Sλ)(Q(n)λn)=λn+1(ΔQ)(n)(S - \lambda)\bigl(Q(n)\lambda^n\bigr) = \lambda^{n+1}(\Delta Q)(n) with ΔQ=Q(X+1)Q(X)\Delta Q = Q(X + 1) - Q(X), and use that Δ\Delta drops the degree; for the dimension, bound it by mm via initial values.)

  4. (The fundamental theorem of linear recurrences) Conclude: if P=i=1r(Xλi)miP = \prod_{i=1}^{r}(X - \lambda_i)^{m_i} with the λi\lambda_i distinct and nonzero, the solutions of (R)(\mathcal R) are exactly the sequences

    un=i=1rQi(n)λin,QiCmi1[X],u_n = \sum_{i=1}^{r} Q_i(n)\,\lambda_i^n, \qquad Q_i \in \C_{m_i - 1}[X],

    with uniquely determined polynomials QiQ_i.

  5. Solve completely: un+2=4un+14unu_{n+2} = 4u_{n+1} - 4u_n, u0=1u_0 = 1, u1=0u_1 = 0, and check the answer on u2u_2.

Part III — Dominant roots and Diophantine dividends.

  1. Suppose the roots are simple with λ1>λi\abs{\lambda_1} > \abs{\lambda_i} for i2i \geq 2, and un=iciλinu_n = \sum_i c_i \lambda_i^n with c10c_1 \neq 0. Show unc1λ1nu_n \sim c_1\lambda_1^n and un+1/unλ1u_{n+1}/u_n \to \lambda_1.
  2. (Pell) Define an+1=an+2bna_{n+1} = a_n + 2b_n, bn+1=an+bnb_{n+1} = a_n + b_n, a0=b0=1a_0 = b_0 = 1. Show that q(a,b)=a22b2q(a, b) = a^2 - 2b^2 satisfies q(an+1,bn+1)=q(an,bn)q(a_{n+1}, b_{n+1}) = -q(a_n, b_n), hence an22bn2=(1)n+1a_n^2 - 2b_n^2 = (-1)^{n+1}; relate this to the determinant of M=(1211)M = \left(\begin{smallmatrix}1 & 2\\ 1 & 1\end{smallmatrix}\right).
  3. Deduce the error estimate

    anbn2=1bn(an+2bn)12bn2,\Bigl|\frac{a_n}{b_n} - \sqrt2\Bigr| = \frac{1}{b_n\,(a_n + \sqrt2\,b_n)} \leq \frac{1}{2b_n^2},

    and show it decays geometrically with ratio 3223 - 2\sqrt2 (find the eigenvalues of MM and the growth of bnb_n).

  4. (General growth) From question 9, prove: (a) if every root satisfies λiρ\abs{\lambda_i} \leq \rho, then unCnm1ρn\abs{u_n} \leq C\,n^{m-1}\rho^n with m=maximim = \max_i m_i; (b) if there is a unique root λ1\lambda_1 of maximal modulus and Q10Q_1 \neq 0, then un+1/unλ1u_{n+1}/u_n \to \lambda_1 — check it on question 10’s solution.

Part IV — Counting walks and words. For a finite graph with vertex set {1,,N}\{1, \dots, N\}, the adjacency matrix AA has Aij=1A_{ij} = 1 if ijij is an edge, else 00.

  1. Prove that (An)ij(A^n)_{ij} is the number of walks of length nn from ii to jj (sequences of nn edges, each step along an edge).
  2. (The triangle) For the complete graph on 33 vertices, A=JIA = J - I: using the spectrum of JJ (Example 2.19), show

    (An)ii=2n+2(1)n3,(An)ij=2n(1)n3(ij),(A^n)_{ii} = \frac{2^n + 2(-1)^n}{3}, \qquad (A^n)_{ij} = \frac{2^n - (-1)^n}{3} \quad (i \neq j),

    and check both at n=2n = 2 by listing walks.

  3. (Words without 1111) Let wnw_n be the number of binary words of length nn with no two consecutive 11’s. Encode words by their last letter to get a transfer matrix, show wn+2=wn+1+wnw_{n+2} = w_{n+1} + w_n, deduce wn=Fn+2w_n = F_{n+2} (Fibonacci, Exercise 3.5), and give the growth rate limwn+1/wn\lim w_{n+1}/w_n.
  4. (The path) For the path graph 1231 - 2 - 3, show the eigenvalues of AA are 2,0,2\sqrt2, 0, -\sqrt2 with eigenvectors (1,±2,1)(1, \pm\sqrt2, 1) and (1,0,1)(1, 0, -1), and deduce that the number of walks of length nn from end to end is ((2)n+(2)n)/4\bigl((\sqrt2)^n + (-\sqrt2)^n\bigr)/4: zero for odd nn, and 2n/212^{\,n/2 - 1} for even nn. Check at n=4n = 4.
  5. (Trace formula) Show that the total number of closed walks of length nn (all starting points) is tr(An)=iλin\operatorname{tr}(A^n) = \sum_i \lambda_i^n, and verify it on the triangle.

Part V — A ring of sequences: simultaneous diagonalization. Fix k3k \geq 3, let ω=e2iπ/k\omega = \eu^{2\iu\pi/k}, and let WMk(C)W \in \mathcal{M}_k(\C) be the cyclic shift: Wei=ei+1W e_i = e_{i+1} (indices mod kk, columns indexed 0,,k10, \dots, k-1).

  1. Show that WTW^{\mathsf T} is the companion matrix of Xk1X^k - 1, deduce χW=μW=Xk1\chi_W = \mu_W = X^k - 1, and that WW is diagonalizable with the kk simple eigenvalues ωj\omega^j and eigenvectors fj=(1,ωj,ω2j,,ω(k1)j)Tf_j = (1, \omega^{-j}, \omega^{-2j}, \dots, \omega^{-(k-1)j})^{\mathsf T}.
  2. A circulant matrix is C=c0I+c1W++ck1Wk1C = c_0 I + c_1 W + \dots + c_{k-1}W^{k-1}. Show that all circulants commute, that the basis (f0,,fk1)(f_0, \dots, f_{k-1}) diagonalizes all of them simultaneously, and that the eigenvalues of CC are c^(ωj)=mcmωjm\widehat c(\omega^j) = \sum_m c_m \omega^{jm}, j=0,,k1j = 0, \dots, k-1.
  3. Deduce detC=j=0k1c^(ωj)\det C = \prod_{j=0}^{k-1} \widehat c(\omega^j), and check that k=3k = 3 recovers the factorization of Exercise 2.8.
  4. (The necklace average) Let x(n+1)=Mx(n)x^{(n+1)} = Mx^{(n)} with M=12(W+W1)M = \frac12(W + W^{-1}): each of kk numbers arranged in a ring is replaced by the average of its two neighbours. Show the eigenvalues of MM are cos(2πj/k)\cos(2\pi j/k), and that the coefficient of x(0)x^{(0)} on f0f_0 is the mean 1kmxm(0)\frac1k\sum_m x^{(0)}_m (sum the coordinates of the fjf_j).
  5. Conclude: for odd kk, x(n)x^{(n)} converges to the constant vector whose value is the mean of the initial values; for k=4k = 4, exhibit the eigenvalue responsible for non-convergence and the exact obstruction (an alternating-mean coefficient that must vanish).
  6. (Synthesis) In one sentence each: how the companion matrix converts analysis of (R)(\mathcal R) into reduction; where the kernel decomposition lemma needed no finite dimension; why dominant eigenvalues govern growth rates and Diophantine error; why powers of the adjacency matrix count walks; and what commuting matrices buy. Name the two summits: the fundamental theorem of linear recurrences, and — for the positive matrices of Part IV, in the Year 3 volume — the Perron–Frobenius theorem.
Solution

Solution of Problem 3.1.

1. The first k1k - 1 coordinates of CvnCv_n are un+1,,un+k1u_{n+1}, \dots, u_{n+k-1} (the superdiagonal shifts), and the last is a0un++ak1un+k1a_0 u_n + \dots + a_{k-1}u_{n+k-1}. So vn+1=Cvnv_{n+1} = Cv_n holds for all nn iff the last coordinates match for all nn, i.e. iff (R)(\mathcal R) holds. Iterating, vn=Cnv0v_n = C^nv_0.

2. Expand Dk(X)=det(XIkC)D_k(X) = \det(XI_k - C) along the first column: the two nonzero entries are XX (position (1,1)(1,1)) and a0-a_0 (position (k,1)(k,1)). The first minor is Dk1D_{k-1}-shaped for the coefficients a1,,ak1a_1, \dots, a_{k-1}; the second minor is upper triangular with diagonal 1-1: determinant (1)k1(-1)^{k-1}, with sign (1)k+1(-1)^{k+1} from the position. Induction on kk (base k=1k = 1: Xa0X - a_0) gives

Dk(X)=X(Xk1ak1Xk2a1)a0=P(X).D_k(X) = X\bigl(X^{k-1} - a_{k-1}X^{k-2} - \dots - a_1\bigr) - a_0 = P(X).

For μC\mu_C: since Q(CT)=Q(C)TQ(C^{\mathsf T}) = Q(C)^{\mathsf T} for any polynomial, CC and CTC^{\mathsf T} have the same annihilators, hence the same minimal polynomial. For CTC^{\mathsf T}: the columns read CTe1=e2C^{\mathsf T}e_1 = e_2, …, CTek1=ekC^{\mathsf T}e_{k-1} = e_k, so (e1,CTe1,,(CT)k1e1)(e_1, C^{\mathsf T}e_1, \dots, (C^{\mathsf T})^{k-1}e_1) is the canonical basis: free. A polynomial Q0Q \neq 0 of degree <k< k then has Q(CT)e10Q(C^{\mathsf T})e_1 \neq 0 (it is a nontrivial combination of basis vectors): degμk\deg\mu \geq k. As μχ=P\mu \mid \chi = P with degP=k\deg P = k: μC=P\mu_C = P.

3. For v=(1,λ,,λk1)Tv = (1, \lambda, \dots, \lambda^{k-1})^{\mathsf T}: rows 11 to k1k-1 of CvCv give λ,λ2,,λk1\lambda, \lambda^2, \dots, \lambda^{k-1}, i.e. λ\lambda times the first k1k - 1 entries of vv; the last row gives mamλm=λkP(λ)=λk=λλk1\sum_m a_m\lambda^m = \lambda^k - P(\lambda) = \lambda^k = \lambda\cdot\lambda^{k-1}. So Cv=λvCv = \lambda v. Conversely, the equations (Cx)i=λxi(Cx)_i = \lambda x_i for i<ki < k read xi+1=λxix_{i+1} = \lambda x_i: any eigenvector is proportional to vv — every eigenspace has dimension exactly 11. Diagonalizable iff the eigenspace dimensions sum to kk (Theorem 3.6) iff there are kk distinct eigenvalues iff PP has kk distinct roots (the eigenvalues are the roots of χC=P\chi_C = P).

4. Each (λin)n(\lambda_i^n)_n solves (R)(\mathcal R): λin+k=λinλik=λinmamλim\lambda_i^{n+k} = \lambda_i^n\,\lambda_i^k = \lambda_i^n\sum_m a_m\lambda_i^m. Freeness: a vanishing combination iciλin=0\sum_i c_i\lambda_i^n = 0 for n=0,,k1n = 0, \dots, k-1 is a Vandermonde system (Exercise 2.11) in the cic_i: all ci=0c_i = 0. The solution space has dimension kk (question 6, whose proof is elementary and independent): kk free solutions form a basis, and coordinates are unique.

5. P=X2X6=(X3)(X+2)P = X^2 - X - 6 = (X - 3)(X + 2): general solution un=A3n+B(2)nu_n = A\,3^n + B(-2)^n. Initial conditions: A+B=1A + B = 1, 3A2B=83A - 2B = 8: A=2A = 2, B=1B = -1:

un=23n(2)n.u_n = 2\cdot 3^n - (-2)^n .

(Check: u2=u1+6u0=14u_2 = u_1 + 6u_0 = 14 and 294=142\cdot9 - 4 = 14.)

6. P(S)((un))P(S)\bigl((u_n)\bigr) is the sequence nun+kak1un+k1a0unn \mapsto u_{n+k} - a_{k-1}u_{n+k-1} - \dots - a_0u_n: it vanishes iff (R)(\mathcal R) holds, so the solution set is kerP(S)\ker P(S), a subspace. The map kerP(S)Ck\ker P(S) \to \C^k, u(u0,,uk1)u \mapsto (u_0, \dots, u_{k-1}), is linear, injective (the recurrence determines uk,uk+1,u_{k}, u_{k+1}, \dots from the first kk values, by induction) and surjective (define unu_n recursively from any initial data): dimension kk.

7. The proof of Theorem 3.14 uses only: the Bézout identity in C[X]\C[X], and the fact that polynomials in a fixed endomorphism commute. Neither mentions the dimension of the ambient space: the lemma holds verbatim for SL(S)S \in \mathcal{L}(\mathcal{S}). Hence

kerP(S)=i=1rker(Sλiid)mi.\ker P(S) = \bigoplus_{i=1}^{r} \ker\,(S - \lambda_i\,\mathrm{id})^{m_i}.

8. For QC[X]Q \in \C[X]: (Sλ)(Q(n)λn)n(S - \lambda)\bigl(Q(n)\lambda^n\bigr)_n has nn-th term Q(n+1)λn+1λQ(n)λn=λn+1(ΔQ)(n)Q(n{+}1)\lambda^{n+1} - \lambda Q(n)\lambda^n = \lambda^{n+1}(\Delta Q)(n), with ΔQ=Q(X+1)Q(X)\Delta Q = Q(X{+}1) - Q(X) of degree degQ1\deg Q - 1 (the leading terms cancel). Iterating, (Sλ)m(Q(n)λn)=(λn+m(ΔmQ)(n))n(S - \lambda)^m\bigl(Q(n)\lambda^n\bigr) = \bigl(\lambda^{n+m}(\Delta^m Q)(n)\bigr)_n, and ΔmQ=0\Delta^m Q = 0 when degQm1\deg Q \leq m - 1: the right-hand set is contained in the kernel. It is a subspace of dimension mm: the sequences (njλn)n(n^j\lambda^n)_n, 0j<m0 \leq j < m, are free, since jcjnjλn=0\sum_j c_j n^j\lambda^n = 0 for all nn forces (dividing by λn0\lambda^n \neq 0) the polynomial jcjXj\sum_j c_jX^j to vanish at every nNn \in \N, hence to be zero. Conversely dimker(Sλ)mm\dim\ker(S - \lambda)^m \leq m: expanding (Sλ)m=j(mj)(λ)mjSj(S - \lambda)^m = \sum_j \binom mj(-\lambda)^{m-j}S^j, the equation (Sλ)mu=0(S - \lambda)^m u = 0 is a linear recurrence of order mm (leading coefficient 11), so uu is determined by u0,,um1u_0, \dots, u_{m-1} as in question 6. Equality of dimensions concludes.

9. Combine questions 7 and 8: every solution decomposes uniquely as a sum of elements of the ker(Sλi)mi\ker(S - \lambda_i)^{m_i}, i.e. un=iQi(n)λinu_n = \sum_i Q_i(n)\lambda_i^n with degQimi1\deg Q_i \leq m_i - 1; the QiQ_i are unique because the decomposition is direct and, within each summand, the coefficients of QiQ_i are coordinates in the basis (njλin)j(n^j\lambda_i^n)_j (question 8). Sanity check on dimensions: imi=k\sum_i m_i = k.

10. P=X24X+4=(X2)2P = X^2 - 4X + 4 = (X - 2)^2: solutions (a+bn)2n(a + bn)2^n. Initial data: a=1a = 1, 2(a+b)=02(a + b) = 0, so b=1b = -1:

un=(1n)2n.u_n = (1 - n)\,2^n .

Check: u2=4u14u0=4u_2 = 4u_1 - 4u_0 = -4, and (12)4=4(1 - 2)\cdot4 = -4.

11. Write un=λ1n(c1+i2ci(λi/λ1)n)u_n = \lambda_1^n\bigl(c_1 + \sum_{i\geq2} c_i(\lambda_i/\lambda_1)^n\bigr); each ratio has modulus <1< 1, so the bracket tends to c10c_1 \neq 0: unc1λ1nu_n \sim c_1\lambda_1^n. In particular un0u_n \neq 0 for large nn, and

un+1un=λ1c1+o(1)c1+o(1)λ1.\frac{u_{n+1}}{u_n} = \lambda_1\,\frac{c_1 + o(1)}{c_1 + o(1)} \longrightarrow \lambda_1 .

12. Compute:

q(an+1,bn+1)=(an+2bn)22(an+bn)2=an2+2bn2=q(an,bn).q(a_{n+1}, b_{n+1}) = (a_n + 2b_n)^2 - 2(a_n + b_n)^2 = -a_n^2 + 2b_n^2 = -q(a_n, b_n).

With q(a0,b0)=12=1q(a_0, b_0) = 1 - 2 = -1: an22bn2=(1)n+1a_n^2 - 2b_n^2 = (-1)^{n+1}. Structurally: q(a,b)=(a2b)(a+2b)q(a, b) = (a - \sqrt2\,b)(a + \sqrt2\,b) and the linear map MM multiplies the factor a+2ba + \sqrt2 b by 1+21 + \sqrt2 and the factor a2ba - \sqrt2 b by 121 - \sqrt2 (compute: an+1+2bn+1=(1+2)(an+2bn)a_{n+1} + \sqrt2 b_{n+1} = (1 + \sqrt2)(a_n + \sqrt2 b_n)); the product is multiplied by (1+2)(12)=1=detM(1 + \sqrt2)(1 - \sqrt2) = -1 = \det M each step.

13. Since an22bn2=(an2bn)(an+2bn)=(1)n+1a_n^2 - 2b_n^2 = (a_n - \sqrt2 b_n)(a_n + \sqrt2 b_n) = (-1)^{n+1},

anbn2=an22bn2bn(an+2bn)=1bn(an+2bn)12bn2,\Bigl|\frac{a_n}{b_n} - \sqrt2\Bigr| = \frac{\abs{a_n^2 - 2b_n^2}}{b_n(a_n + \sqrt2 b_n)} = \frac{1}{b_n(a_n + \sqrt2 b_n)} \leq \frac1{2b_n^2},

using anbn1a_n \geq b_n \geq 1 (induction: both increase) so an+2bn(1+2)bn2bna_n + \sqrt2 b_n \geq (1 + \sqrt2)b_n \geq 2b_n. Eigenvalues of MM: χM=X22X1\chi_M = X^2 - 2X - 1, roots 1±21 \pm \sqrt2; since (a0,b0)(a_0, b_0) has a nonzero component on the dominant eigenvector (all entries positive), bnc(1+2)nb_n \sim c(1 + \sqrt2)^n with c>0c > 0 (question 11). Hence the error is (1+2)2n=(3+22)n\asymp (1 + \sqrt2)^{-2n} = (3 + 2\sqrt2)^{-n}: geometric decay with ratio 1/(3+22)=3220.1721/(3 + 2\sqrt2) = 3 - 2\sqrt2 \approx 0.172.

14. (a) From question 9: uniQi(n)λin(iQi(n))ρn\abs{u_n} \leq \sum_i \abs{Q_i(n)}\abs{\lambda_i}^n \leq \bigl(\sum_i \abs{Q_i(n)}\bigr)\rho^n, and each Qi(n)Cinmi1Cinm1\abs{Q_i(n)} \leq C_i n^{m_i - 1} \leq C_i n^{m-1} for n1n \geq 1: sum the constants. (b) Let ρ=maxi2λi<λ1\rho' = \max_{i \geq 2}\abs{\lambda_i} < \abs{\lambda_1} and d=degQ1d = \deg Q_1, leading coefficient c0c \neq 0. Then un=Q1(n)λ1n+Rnu_n = Q_1(n)\lambda_1^n + R_n with RnCnm1ρn\abs{R_n} \leq Cn^{m-1}\rho'^n, and

RnQ1(n)λ1n=O(nm1d(ρ/λ1)n)0\frac{R_n}{Q_1(n)\lambda_1^n} = O\Bigl(n^{m-1-d} \bigl(\rho'/\abs{\lambda_1}\bigr)^n\Bigr) \longrightarrow 0

(geometric beats polynomial). So

unQ1(n)λ1ncndλ1n,un+1unλ1(since Q1(n+1)/Q1(n)1).u_n \sim Q_1(n)\,\lambda_1^n \sim c\,n^d\lambda_1^n, \qquad \frac{u_{n+1}}{u_n} \longrightarrow \lambda_1 \quad\text{(since } Q_1(n{+}1)/Q_1(n) \to 1\text{)}.

Check on question 10: for un=(1n)2nu_n = (1-n)2^n the ratio is

(n)2n+1(1n)2n=2n1n2=λ1.\frac{(-n)2^{n+1}}{(1-n)2^n} = 2\,\frac{-n}{1-n} \longrightarrow 2 = \lambda_1 .

15. Induction on nn. For n=1n = 1, AijA_{ij} counts walks of length 11. Step: a walk of length n+1n + 1 from ii to jj is a walk of length nn from ii to some vertex \ell followed by an edge j\ell j:

#{walks}=(An)iAj=(An+1)ij.\#\{\text{walks}\} = \sum_{\ell} (A^n)_{i\ell}A_{\ell j} = (A^{n+1})_{ij}.

16. J=3ΠJ = 3\Pi where Π=J/3\Pi = J/3 is the projection onto Vect(1,1,1)\operatorname{Vect}(1,1,1) along the plane x+y+z=0x + y + z = 0 (Π2=Π\Pi^2 = \Pi since J2=3JJ^2 = 3J). Then A=JI=2Π(IΠ)A = J - I = 2\Pi - (I - \Pi), and since Π\Pi and IΠI - \Pi are complementary projections,

An=2nΠ+(1)n(IΠ),i.e.(An)ij=2n3+(1)n(δij13),A^n = 2^n\,\Pi + (-1)^n (I - \Pi), \qquad\text{i.e.}\qquad (A^n)_{ij} = \frac{2^n}3 + (-1)^n\Bigl(\delta_{ij} - \frac13\Bigr),

which gives the two displayed formulas. At n=2n = 2: diagonal (4+2)/3=2(4 + 2)/3 = 2 (walks iii \to \ell \to i for the two neighbours \ell); off-diagonal (41)/3=1(4 - 1)/3 = 1 (the single walk iji \to \ell \to j through the third vertex).

17. Let wn(0),wn(1)w_n^{(0)}, w_n^{(1)} count admissible words of length nn ending in 00, resp. 11. Appending a letter: a 00 may follow anything, a 11 only a 00:

(wn+1(0)wn+1(1))=(1110)(wn(0)wn(1)).\begin{pmatrix} w_{n+1}^{(0)}\\ w_{n+1}^{(1)}\end{pmatrix} = \begin{pmatrix} 1 & 1\\ 1 & 0\end{pmatrix} \begin{pmatrix} w_n^{(0)}\\ w_n^{(1)}\end{pmatrix}.

Summing, wn+2=wn+1+wnw_{n+2} = w_{n+1} + w_n (or: condition on the first letter). With w1=2w_1 = 2, w2=3w_2 = 3: wn=Fn+2w_n = F_{n+2} by induction (F3=2F_3 = 2, F4=3F_4 = 3, same recurrence). Growth: the roots of X2X1X^2 - X - 1 are φ>ψ\varphi > \abs\psi (Exercise 3.5), and the φ\varphi-component is nonzero (the wnw_n are positive and ψn0\psi^n \to 0), so question 11 gives wn+1/wnφ=1+52w_{n+1}/w_n \to \varphi = \frac{1 + \sqrt5}2.

18. A=(010101010)A = \left(\begin{smallmatrix} 0&1&0\\ 1&0&1\\ 0&1&0\end{smallmatrix}\right). Check:

A(1,±2,1)T=(±2,2,±2)T=±2(1,±2,1)T,A(1,0,1)T=0:A(1, \pm\sqrt2, 1)^{\mathsf T} = (\pm\sqrt2, 2, \pm\sqrt2)^{\mathsf T} = \pm\sqrt2\,(1, \pm\sqrt2, 1)^{\mathsf T}, \qquad A(1, 0, -1)^{\mathsf T} = 0 :

eigenvalues 2,2,0\sqrt2, -\sqrt2, 0 (=2cosπ4,2cos3π4,2cosπ2= 2\cos\frac\pi4, 2\cos\frac{3\pi}4, 2\cos\frac\pi2). Decompose e1e_1 on the eigenbasis and read the third coordinate, or use symmetry: with v±=(1,±2,1)v_\pm = (1, \pm\sqrt2, 1), v0=(1,0,1)v_0 = (1, 0, -1), one checks e1=14v++14v+12v0e_1 = \frac14 v_+ + \frac14 v_- + \frac12 v_0, so for n1n \geq 1

(An)13=(14(2)nv++14(2)nv+0) ⁣3=(2)n+(2)n4,(A^n)_{13} = \Bigl(\tfrac14(\sqrt2)^n v_+ + \tfrac14(-\sqrt2)^n v_- + 0\Bigr)_{\!3} = \frac{(\sqrt2)^n + (-\sqrt2)^n}{4},

zero for odd nn (bipartite graph: ends are at even distance), and 22n/2/4=2n/212\cdot 2^{n/2}/4 = 2^{n/2 - 1} for even nn. At n=4n = 4: 21=22^{1} = 2, matching the two walks 121231\,2\,1\,2\,3 and 123231\,2\,3\,2\,3.

19. Closed walks of length nn from ii are (An)ii(A^n)_{ii}; summing over ii gives tr(An)\operatorname{tr}(A^n). Trigonalizing AA (over C\C), AnA^n is triangular with diagonal λin\lambda_i^n: tr(An)=iλin\operatorname{tr}(A^n) = \sum_i\lambda_i^n. Triangle: tr(An)=32n+2(1)n3=2n+2(1)n=2n+(1)n+(1)n\operatorname{tr}(A^n) = 3\,\frac{2^n + 2(-1)^n}3 = 2^n + 2(-1)^n = 2^n + (-1)^n + (-1)^n: the spectrum {2,1,1}\{2, -1, -1\}, consistent with question 16.

20. The columns of WTW^{\mathsf T}: WTei=ei1W^{\mathsf T}e_i = e_{i-1} for i1i \geq 1 and WTe0=ek1W^{\mathsf T}e_0 = e_{k-1}; relabelling in the order e0,e1,e_0, e_1, \dots this is exactly the companion matrix of Xk1X^k - 1 (a0=1a_0 = 1, other am=0a_m = 0). Question 2: χW=χWT=Xk1=μW\chi_{W} = \chi_{W^{\mathsf T}} = X^k - 1 = \mu_{W}. The roots ωj\omega^j (j=0,,k1j = 0, \dots, k-1) are the kk distinct kk-th roots of unity: WW is diagonalizable (question 3, or Exercise 3.8: Wk=IW^k = I). Eigenvectors: Wfj=mωjmem+1=mωj(m1)em=ωjfjWf_j = \sum_m \omega^{-jm}e_{m+1} = \sum_{m'}\omega^{-j(m'-1)}e_{m'} = \omega^j f_j.

21. Circulants are polynomials in WW, and polynomials in a fixed matrix commute with one another. Each fjf_j is an eigenvector of every power: Wmfj=ωjmfjW^m f_j = \omega^{jm}f_j, so

Cfj=mcmωjmfj=c^(ωj)fj:Cf_j = \sum_m c_m\omega^{jm} f_j = \widehat c(\omega^j)\,f_j :

the basis (f0,,fk1)(f_0, \dots, f_{k-1}) (free: Vandermonde in the distinct ωj\omega^{-j}, Exercise 2.11) diagonalizes every circulant at once, with the stated eigenvalues.

22. The determinant is the product of the eigenvalues (diagonalize): detC=jc^(ωj)\det C = \prod_{j}\widehat c(\omega^j). For k=3k = 3, c0=ac_0 = a, c1=bc_1 = b, c2=cc_2 = c and ω=j=e2iπ/3\omega = j = \eu^{2\iu\pi/3}:

detC=(a+b+c)(a+bj+cj2)(a+bj2+cj4),\det C = (a + b + c)(a + bj + cj^2)(a + bj^2 + cj^4),

and j4=jj^4 = j: exactly the factorization of Exercise 2.8.

23. M=12(W+W1)M = \frac12(W + W^{-1}) is a circulant (W1=Wk1W^{-1} = W^{k-1}), with eigenvalues 12(ωj+ωj)=cos2πjk\frac12(\omega^j + \omega^{-j}) = \cos\frac{2\pi j}k on the same basis fjf_j. Coordinates: write x(0)=jαjfjx^{(0)} = \sum_j \alpha_j f_j. The coordinates of fjf_j sum to mωjm\sum_m \omega^{-jm}, which is kk for j=0j = 0 and 00 otherwise (geometric sum with ratio ωj1\omega^{-j} \neq 1). Summing the coordinates of x(0)x^{(0)}: mxm(0)=α0k\sum_m x^{(0)}_m = \alpha_0\,k, so α0=1kmxm(0)\alpha_0 = \frac1k\sum_m x^{(0)}_m, the mean.

24. x(n)=Mnx(0)=jαjcosn(2πjk)fjx^{(n)} = M^nx^{(0)} = \sum_j \alpha_j\cos^n\bigl(\tfrac{2\pi j}k\bigr)f_j. For odd kk, cos(2πj/k)<1\abs{\cos(2\pi j/k)} < 1 for every j0j \neq 0 (the angle is never 00 or π\pi), so all terms but j=0j = 0 tend to 00: x(n)α0f0x^{(n)} \to \alpha_0 f_0, the constant vector equal to the mean — averaging on an odd ring equalizes. For k=4k = 4 the eigenvalues are 1,0,1,01, 0, -1, 0: the j=2j = 2 term α2(1)nf2\alpha_2(-1)^nf_2 with f2=(1,1,1,1)Tf_2 = (1, -1, 1, -1)^{\mathsf T} oscillates forever. The obstruction is the alternating mean: multiplying the coordinates of x(0)x^{(0)} by (1)m(-1)^m and summing, the same geometric-sum computation gives m(1)mxm(0)=4α2\sum_m (-1)^mx^{(0)}_m = 4\alpha_2: the process converges iff x0(0)x1(0)+x2(0)x3(0)=0x^{(0)}_0 - x^{(0)}_1 + x^{(0)}_2 - x^{(0)}_3 = 0, and then converges to the mean.

25. The companion matrix converts a scalar recurrence of order kk into a first-order vector recurrence, so that closed formulas become statements about CnC^n — reduction’s home ground (questions 1–5). The kernel decomposition lemma is pure polynomial algebra (Bézout plus commutation), so it splits kerP(S)\ker P(S) even though S\mathcal{S} is infinite-dimensional (questions 7–9). Dominant eigenvalues govern growth because every other contribution is geometrically negligible after normalization — which is also why the Pell error decays at the square of the dominant root (questions 11–14). Powers of the adjacency matrix count walks because matrix multiplication sums over intermediate vertices, so spectra count closed walks (questions 15–19). Commuting matrices share an eigenbasis, and one Fourier basis then diagonalizes the whole circulant algebra in one stroke (questions 20–24). Summits: the fundamental theorem of linear recurrences (question 9); and for nonnegative matrices, the reason dominant roots like φ\varphi or 1+21 + \sqrt2 are automatically real, positive and simple is the Perron–Frobenius theorem, proved in the Year 3 volume.