CalcPath
ProgressProblem setsFeedback
← Question sheet

Linear Algebra · Solutions

Generalized Eigenvectors and Jordan Form

12 questions · about 108 minutes

Work this sheet on paper and show every step — the solutions are written as methods, so you can check your route and not just your final matrix. Name every theorem you use (“by Cayley–Hamilton”, “by the rank formula”, “by the diagonalizability criterion for the minimal polynomial”), and whenever you claim a Jordan form, say which rank or multiplicity data forced each block. No calculator is needed anywhere: every number is a small integer or an exact expression in e. Full worked solutions are on a separate page, and each question has its own reveal so you can check one answer without seeing the rest.

  1. 1.Computation

    Two matrices:

    A=(51​−13​),C=​1−11​131​002​​.

    1. (a)Find the characteristic polynomial of A, both multiplicities of each eigenvalue, and a basis of each eigenspace. Is A diagonalizable?
    2. (b)Do the same for C. (Expanding det(λI−C) along its third column saves work.)
    3. (c)Give a one-line argument, using nothing but the list of eigenvalues, that C cannot be diagonalizable. Name the fact you are using.

    Solution

    (a) Characteristic polynomial and multiplicities. trA=8 and detA=15−(−1)(1)=16, so

    pA​(λ)=λ2−8λ+16=(λ−4)2.

    The only eigenvalue is 4, with ma​(4)=2. For the geometric multiplicity,

    A−4I=(11​−1−1​)

    has rank 1, so mg​(4)=2−1=1. The eigenspace is {(x,y):x−y=0}=span(1,1); check: A(1,1)=(4,4). Since mg​(4)=1<2=ma​(4), the eigenvalue is defective and A is not diagonalizable: there is no basis of C2 made of eigenvectors.

    (b) The 3×3 matrix. The third column of λI−C is (0,0,λ−2), so expanding along it,

    pC​(λ)=(λ−2)det(λ−11​−1λ−3​)=(λ−2)((λ−1)(λ−3)+1)=(λ−2)(λ2−4λ+4)=(λ−2)3.

    So 2 is the only eigenvalue and ma​(2)=3. Now

    C−2I=​−1−11​111​000​​.

    Rows 1 and 2 are equal; rows 1 and 3 are not proportional (their sum is (0,2,0)), so the rank is 2 and mg​(2)=3−2=1. The eigenspace is cut out by −x+y=0 and x+y=0, forcing x=y=0: it is span(0,0,1), and indeed Ce3​=(0,0,2)=2e3​. The deficit is ma​−mg​=2: two chairs are empty. C is not diagonalizable.

    (c) The scalar test. C has the single eigenvalue 2 and C=2I. If C were diagonalizable, P−1CP would be a diagonal matrix whose entries are eigenvalues of C, hence all equal to 2; so P−1CP=2I and C=P(2I)P−1=2I, which is false. This is the proposition that a matrix with only one eigenvalue λ0​ is diagonalizable only if it is λ0​I. The same one-liner disposes of A in (a): its only eigenvalue is 4 and A=4I.

  2. 2.Proof

    Recall that a subspace W⊆Cn is A-invariant if Aw∈W for every w∈W. Let W=span(e1​,e2​)⊆C3 and

    A=​130​240​456​​.

    1. (a)Show that W is A-invariant, write down the matrix A1​ of the restriction A∣W​ in the basis e1​,e2​, and verify directly that pA​(λ)=pA1​​(λ)pA2​​(λ), where A2​ is the remaining diagonal block.
    2. (b)Prove: if W is A-invariant then W is q(A)-invariant for every polynomial q; and if W1​ and W2​ are A-invariant then so are W1​∩W2​ and W1​+W2​.
    3. (c)Prove: if A is invertible and W is A-invariant, then W is A−1-invariant. (Consider the restriction A∣W​:W→W.)

    Solution

    (a) A concrete invariant subspace. Ae1​=(1,3,0) and Ae2​=(2,4,0) both have third coordinate 0, so they lie in W; a general w=xe1​+ye2​ has Aw=xAe1​+yAe2​∈W by linearity. The columns of A1​ are the coordinates of Ae1​,Ae2​ in the basis e1​,e2​:

    A1​=(13​24​),pA1​​(λ)=λ2−5λ+(4−6)=λ2−5λ−2.

    Here P=I already, and A is block upper triangular with A2​=(6), so pA2​​(λ)=λ−6. Directly, expanding det(λI−A) along the third row,

    pA​(λ)=(λ−6)det(λ−1−3​−2λ−4​)=(λ−6)((λ−1)(λ−4)−6)=(λ−6)(λ2−5λ−2),

    which is pA1​​pA2​​, as the block triangular theorem predicts.

    (b) Polynomials, intersections, sums. First, Ajw∈W for every j≥0 and w∈W, by induction: A0w=w∈W, and if Ajw∈W then Aj+1w=A(Ajw)∈W by invariance. For q(λ)=∑j​cj​λj we get q(A)w=∑j​cj​Ajw, a linear combination of vectors of W, hence in W.

    If w∈W1​∩W2​ then Aw∈W1​ (as w∈W1​) and Aw∈W2​ (as w∈W2​), so Aw∈W1​∩W2​. If w=w1​+w2​ with wi​∈Wi​ then Aw=Aw1​+Aw2​ with Awi​∈Wi​, so Aw∈W1​+W2​. ■

    (c) Invariance under the inverse. The restriction A∣W​:W→W is a linear map from W to itself. It is injective, because A itself is injective (A is invertible, so Aw=0 forces w=0). An injective linear map from a finite-dimensional space to itself is surjective, by rank–nullity: dimim=dimW−dimker=dimW.

    Now let w∈W. By surjectivity there is w′∈W with Aw′=w, and then A−1w=w′∈W. So W is A−1-invariant. ■

    (Once Cayley–Hamilton is available there is a second proof: A−1 is a polynomial in A whenever detA=0, so (b) applies directly.)

  3. 3.Computation

    Let A=(24​1−1​), and let B be a 3×3 matrix with characteristic polynomial pB​(λ)=λ3−2λ2−5λ+6.

    1. (a)Write down pA​ and the identity that Cayley–Hamilton gives for A. Use it — without computing A3 or A4 directly — to express A4 in the form αA+βI. Check your answer against each eigenvalue of A.
    2. (b)Use the same identity to write A−1 in the form γA+δI, and verify by multiplying.
    3. (c)Show that B is invertible, find detB, and express B−1 as a polynomial in B of degree 2. Check the expression at each root of pB​.

    Solution

    (a) Reduce the power. trA=1 and detA=−2−4=−6, so pA​(λ)=λ2−λ−6=(λ−3)(λ+2). Cayley–Hamilton says A2−A−6I=0, i.e.

    A2=A+6I.

    Multiply by A: A3=A2+6A=(A+6I)+6A=7A+6I. Multiply again: A4=7A2+6A=7(A+6I)+6A=13A+42I. So α=13, β=42.

    Check at the eigenvalues, which must satisfy every polynomial identity that A does: at λ=3, 34=81 and 13⋅3+42=81; at λ=−2, (−2)4=16 and 13(−2)+42=16. Both agree. (Explicitly A4=(6852​1329​), and squaring A2=(84​15​) gives the same matrix.)

    (b) The inverse. From A2−A=6I we get A(A−I)=6I, so

    A−1=61​(A−I)=61​(14​1−2​),γ=61​, δ=−61​.

    Verify: A(A−I)=(24​1−1​)(14​1−2​)=(2+44−4​2−24+2​)=6I.

    (c) The 3×3 case. The constant term is c0​=pB​(0)=det(0⋅I−B)=det(−B)=(−1)3detB, so detB=−c0​=−6=0 and B is invertible. (Cross-check: pB​(λ)=(λ−1)(λ+2)(λ−3), and the product of the eigenvalues is 1⋅(−2)⋅3=−6.)

    Cayley–Hamilton: B3−2B2−5B+6I=0, so B(B2−2B−5I)=−6I and

    B−1=−61​(B2−2B−5I).

    Check at each eigenvalue, where the identity must read λ−1=−61​(λ2−2λ−5): at λ=1, −61​(1−2−5)=1; at λ=−2, −61​(4+4−5)=−21​; at λ=3, −61​(9−6−5)=31​. All three are the reciprocals of the eigenvalues, as they must be, since applying the identity to an eigenvector v with Bv=λv gives λ−1v on the left and −61​(λ2−2λ−5)v on the right.

  4. 4.Computation

    Throughout, mX​ denotes the minimal polynomial of a matrix X and pX​ its characteristic polynomial.

    1. (a)Let E be the 3×3 matrix every entry of which is 1. Show that E2=3E, find mE​, decide whether E is diagonalizable, and then determine pE​ without computing any determinant.
    2. (b)A matrix B satisfies B3=B. Prove that B is diagonalizable, and list every polynomial that could be mB​, giving a matrix that realises each one.
    3. (c)A 3×3 matrix C has pC​(λ)=(λ−1)2(λ−4). Determine mC​ if rank(C−I)=2, and determine it if instead rank(C−I)=1. Justify each step from the theorems on the minimal polynomial.

    Solution

    (a) The all-ones matrix. Each entry of E2 is a row of ones dotted with a column of ones, namely 3; so E2=3E. Hence q(λ)=λ2−3λ=λ(λ−3) annihilates E, and mE​ divides q. The monic divisors of degree ≥1 are λ, λ−3 and λ(λ−3). mE​=λ would mean E=0 and mE​=λ−3 would mean E=3I; neither holds. So

    mE​(λ)=λ(λ−3),

    a product of distinct linear factors, and by the diagonalizability criterion E is diagonalizable.

    For pE​: the roots of pE​ are the roots of mE​, namely 0 and 3. E has rank 1, so mg​(0)=3−1=2; since E is diagonalizable, ma​(0)=mg​(0)=2, and the remaining eigenvalue is 3 with multiplicity 1. Thus pE​(λ)=λ2(λ−3). Sanity check: trE=3=0+0+3, and E(1,1,1)=(3,3,3).

    (b) B3=B. The polynomial λ3−λ=λ(λ−1)(λ+1) annihilates B, so mB​ divides it. Every monic divisor of degree ≥1 of a product of distinct linear factors is itself a product of distinct linear factors, so mB​ has no repeated root and B is diagonalizable by the criterion. The candidates are the seven nontrivial monic divisors:

    λ,λ−1,λ+1,λ(λ−1),λ(λ+1),(λ−1)(λ+1),λ(λ−1)(λ+1),

    realised respectively by 0, I, −I, diag(0,1), diag(0,−1), diag(1,−1) and diag(0,1,−1). (A diagonal matrix has minimal polynomial ∏(λ−d) over its distinct diagonal entries d, and each of these satisfies B3=B because d3=d for d∈{0,1,−1}.)

    (c) Reading mC​ off a rank. mC​ divides pC​ and has the same roots 1 and 4, so mC​ is either (λ−1)(λ−4) or (λ−1)2(λ−4).

    If rank(C−I)=2 then mg​(1)=3−2=1<2=ma​(1), so C is not diagonalizable. By the criterion mC​ must then have a repeated root, which rules out (λ−1)(λ−4). Hence mC​(λ)=(λ−1)2(λ−4).

    If rank(C−I)=1 then mg​(1)=2=ma​(1), and mg​(4)=1=ma​(4) automatically, so C is diagonalizable and mC​ has distinct roots: mC​(λ)=(λ−1)(λ−4). In both cases the exponent of (λ−1) in mC​ is the size of the largest Jordan block for 1: two blocks of size 1 in the second case, one block of size 2 in the first.

  5. 5.Computation

    Let

    A=​101​110​002​​,N=A−I.

    1. (a)Compute pA​ and both multiplicities of every eigenvalue. Which eigenvalue is defective, and by how much?
    2. (b)Compute N2, describe the generalized eigenspace K1​=kerN2 as a set, and check that its dimension is what the theory predicts.
    3. (c)Find a generalized eigenvector v of height exactly 2 for the eigenvalue 1, write down its Jordan chain Nv,v, verify that Nv is an eigenvector, and compute Av to confirm that v itself is not one.
    4. (d)Find K2​ and exhibit a basis of C3 adapted to the primary decomposition C3=K1​⊕K2​, verifying independence with a determinant.

    Solution

    (a) Eigenvalues and multiplicities. The third column of λI−A is (0,0,λ−2); expanding along it,

    pA​(λ)=(λ−2)det(λ−10​−1λ−1​)=(λ−2)(λ−1)2.

    So ma​(1)=2 and ma​(2)=1. For λ=2, 1≤mg​≤ma​=1 forces mg​(2)=1. For λ=1,

    N=A−I=​001​100​001​​

    has the two independent rows (0,1,0) and (1,0,1), so rankN=2 and mg​(1)=3−2=1. The eigenvalue 1 is defective with deficit 2−1=1; its eigenspace is {y=0, x+z=0}=span(1,0,−1).

    (b) The generalized eigenspace. Row i of N2 is row i of N applied to N: row 1 is (0,1,0)N= row 2 of N =(0,0,0); row 2 is 0; row 3 is (1,0,1)N= row 1 + row 3 =(1,1,1). Thus

    N2=​001​001​001​​,K1​=kerN2={(x,y,z):x+y+z=0}.

    Its dimension is 3−rankN2=3−1=2=ma​(1), exactly as the theorem dimKλ​=ma​(λ) predicts. (The chain of kernels has stabilised: N3=N⋅N2 has rows 0, 0 and (1,0,1)N2=(1,1,1), so N3=N2 and kerN3=kerN2.)

    (c) A chain of length 2. Height exactly 2 means v∈kerN2∖kerN. Take v=(1,−1,0): its coordinates sum to 0, so v∈K1​, while Nv=(−1,0,1)=0. The chain is

    u1​=Nv=(−1,0,1),u2​=v=(1,−1,0).

    Check the bottom: Nu1​=(0,0,−1+1)=(0,0,0), so u1​ is an eigenvector (a multiple of (1,0,−1) from (a)). And Av=(1−1,−1,1)=(0,−1,1), which is not a multiple of v; instead Av=v+Nv=u2​+u1​ — the eigenvalue times v plus something of lower height, which is the whole point of a generalized eigenvector.

    (d) The primary decomposition. ma​(2)=1, so K2​=ker(A−2I). With

    A−2I=​−101​1−10​000​​,

    row 2 gives y=0 and row 3 gives x=0, so K2​=span(0,0,1); indeed Ae3​=(0,0,2).

    Adapted basis: u1​,u2​ for K1​ and e3​ for K2​. Independence: expanding along the third column,

    det​−101​1−10​001​​=1⋅det(−10​1−1​)=1=0.

    So C3=K1​⊕K2​ with dimK1​+dimK2​=2+1=3. In this basis A is block diagonal: on K1​ it acts as I+N∣K1​​, which in the chain basis is J2​(1), and on K2​ as the 1×1 block (2).

  6. 6.Proof

    Let N, M and A be n×n complex matrices. Recall that N is nilpotent if Nk=0 for some k≥1.

    1. (a)Prove that a nilpotent matrix has trace 0 and determinant 0.
    2. (b)Prove that if N is nilpotent and AN=NA, then AN is nilpotent.
    3. (c)Prove that if N and M are nilpotent and NM=MN, then N+M is nilpotent. (If Nk=0 and Ml=0, consider (N+M)k+l−1.)
    4. (d)Show by an explicit 2×2 example that (c) fails without the hypothesis NM=MN, and that the product NM in your example is not nilpotent either.

    Solution

    (a) Trace and determinant. If Nv=λv with v=0 then 0=Nkv=λkv, so λ=0: the only eigenvalue of a nilpotent matrix is 0. Over C the characteristic polynomial splits and all its roots are 0, so pN​(λ)=λn. But for any matrix the coefficient of λn−1 in det(λI−N) is −trN and the constant term is (−1)ndetN. Comparing with λn, both vanish: trN=0 and detN=0. ■

    (For the determinant alone there is a one-liner: (detN)k=detNk=det0=0.)

    (b) A commuting product. Since A and N commute, (AN)j=AjNj for every j≥1, by induction: (AN)j+1=(AN)j(AN)=AjNjAN=AjANjN=Aj+1Nj+1, using NjA=ANj (itself an induction from NA=AN). With Nk=0 this gives (AN)k=AkNk=0. ■

    (c) A commuting sum. Let Nk=0 and Ml=0. Because N and M commute, the binomial theorem applies:

    (N+M)k+l−1=j=0∑k+l−1​(jk+l−1​)NjMk+l−1−j.

    For each term, either j≥k or k+l−1−j≥l: otherwise j≤k−1 and k+l−1−j≤l−1, whose sum is k+l−1≤k+l−2, absurd. In the first case Nj=0, in the second Mk+l−1−j=0; so every term vanishes and (N+M)k+l−1=0. ■

    (d) Without commutativity. Take

    N=(00​10​),M=(01​00​),

    both with square 0. They do not commute: NM=(10​00​) while MN=(00​01​). The sum N+M=(01​10​) satisfies (N+M)2=I, so (N+M)j is I or N+M for every j and never 0 (equivalently, its eigenvalues are ±1, not 0). And NM is a projection, (NM)2=NM=0, with eigenvalue 1; so neither the sum nor the product is nilpotent. The commuting hypothesis in (b) and (c) cannot be dropped.

  7. 7.Sketch

    A nilpotent 9×9 matrix N has

    rankN=5,rankN2=2,rankN3=1,N4=0.

    You will draw its chain diagram: one column per Jordan chain, one dot per basis vector, the eigenvector at the bottom of each column, columns ordered tallest first. Write rk​=rankNk with r0​=9.

    1. (a)Compute b≥k​=rk−1​−rk​ for k=1,…,5 and draw the diagram. Mark with an arrow what N does to each dot.
    2. (b)Read off from the diagram: the block sizes, the Jordan form of N, its index of nilpotency, mN​, and dimkerN.
    3. (c)Explain, using only the picture, why the number of dots in row k (counting from the bottom) is rk−1​−rk​. Begin by saying which dots Nk sends to 0.
    4. (d)Draw the chain diagram of N2 by tracing where N2 sends each dot, read off the Jordan form of N2, and confirm it with the rank formula applied to the sequence rank(N2)k.

    Solution

    (a) The rank differences and the picture. With (r0​,r1​,r2​,r3​,r4​)=(9,5,2,1,0) and r5​=0,

    b≥1​=9−5=4,b≥2​=5−2=3,b≥3​=2−1=1,b≥4​=1−0=1,b≥5​=0.

    b≥k​ is the number of columns of height at least k, i.e. the number of dots in row k. So the diagram has, from the bottom up, rows of 4, 3, 1, 1 dots: four columns of heights 4,2,2,1, drawn tallest first. Every arrow points straight down one square — N sends each dot to the dot beneath it — and each of the four bottom dots is sent to 0 (draw those arrows leaving the diagram).

    (b) Reading the diagram. The column heights are the block sizes: 4,2,2,1 (sum 9, as it must be). So

    N ∼ J4​(0)⊕J2​(0)⊕J2​(0)⊕J1​(0).

    The index of nilpotency is the tallest column, 4 (the dot at the top of that column survives N3 and dies under N4), so mN​(λ)=λ4. The number of columns is the number of eigenvectors in the basis, dimkerN=4=b≥1​. Consistency: b=1​=4−3=1, b=2​=3−1=2, b=3​=1−1=0, b=4​=1−0=1, matching the heights.

    (c) Why the rows count rank drops. Each dot is a basis vector, and Nk moves every dot down k squares. A dot at height h therefore lands on the dot at height h−k if h>k, and is sent to 0 if h≤k. The dots that survive are exactly those in rows k+1 and above, and their images are distinct basis vectors — the bottom h−k dots of each column — hence independent. So rankNk=rk​ is the number of dots in rows k+1,k+2,…; equivalently rk−1​ is the number of dots in rows k,k+1,…. Subtracting,

    rk−1​−rk​=#{dots in row k}=#{columns of height≥k}=b≥k​.

    Check against the data: rows 2 and above hold 3+1+1=5=r1​ dots; rows 3 and above hold 2=r2​; row 4 holds 1=r3​.

    (d) The diagram of N2. N2 moves each dot down two squares, so a column of N splits into the dots of even height and the dots of odd height, each of which is a chain for N2. The column of height 4 splits into two chains of length 2 (heights 4→2 and 3→1); each column of height 2 splits into two chains of length 1; the column of height 1 stays a chain of length 1. So the diagram of N2 has seven columns of heights 2,2,1,1,1,1,1, and

    N2 ∼ J2​(0)⊕J2​(0)⊕J1​(0)⊕5.

    Rank formula check: the powers of N2 have ranks rankN2=2 and rankN4=0, so the sequence is (9,2,0,…) and b≥1​=7, b≥2​=2, b≥3​=0: two blocks of size 2 and five of size 1, total 4+5=9. The picture and the formula agree.

  8. 8.Critique

    A 4×4 matrix A has pA​(λ)=(λ−1)4 and rank(A−I)=2. A student writes:

    Step 1. The number of Jordan blocks for λ=1 is rank(A−I)=2.

    Step 2. The block sizes add up to ma​(1)=4 and there are two blocks, so each block has size 4/2=2.

    Step 3. Hence J=J2​(1)⊕J2​(1) and mA​(λ)=(λ−1)2; since mA​ has a repeated root, A is not diagonalizable.

    The number in Step 1 is correct but the reasoning is not; Step 2 is wrong; the last clause of Step 3 is correct.

    1. (a)State the correct formula for the number of blocks, explain why the student's version happens to give the right number here, and give a 4×4 matrix for which it gives the wrong number.
    2. (b)Explain precisely why Step 2 does not follow, and list every Jordan form consistent with the given data, together with the minimal polynomial of each.
    3. (c)Name one further piece of data that decides between the possibilities in (b), and say what conclusion each of its possible values gives.
    4. (d)Give a correct proof that A is not diagonalizable which uses only the data given, without knowing the Jordan form.

    Solution

    (a) The block count. The number of Jordan blocks for λ is the geometric multiplicity,

    b≥1​=r0​−r1​=n−rank(A−λI),

    not the rank itself. Here n−rank(A−I)=4−2=2, and the student's rank(A−I)=2 agrees only because 4−2=2: the two formulas coincide exactly when the rank is n/2. For A=I4​ the rank of A−I is 0, so the student would count 0 blocks, whereas there are 4 blocks of size 1; for A=J4​(1) the rank is 3 and the student would count 3 blocks, whereas there is 1.

    (b) Equal sizes are not forced. The rank formula with k=1 gives only the number of blocks; it says nothing about their sizes, which are governed by the ranks of the higher powers (A−I)2,(A−I)3,…. Two blocks with sizes summing to 4 can be 3+1 or 2+2 — the two partitions of 4 into two parts — and nothing given so far excludes either. The consistent Jordan forms are

    J3​(1)⊕J1​(1)withmA​(λ)=(λ−1)3,orJ2​(1)⊕J2​(1)withmA​(λ)=(λ−1)2,

    since the exponent in the minimal polynomial is the largest block size.

    (c) What decides it. The rank of (A−I)2. By the rank-of-a-block lemma, rank(Jm​(1)−I)2=max(m−2,0), so r2​=1+0=1 for J3​(1)⊕J1​(1) and r2​=0+0=0 for J2​(1)⊕J2​(1). Thus: if rank(A−I)2=1 then b≥2​=r1​−r2​=1 and b≥3​=r2​−r3​=1, giving J3​(1)⊕J1​(1); if rank(A−I)2=0 then b≥2​=2 and b≥3​=0, giving J2​(1)⊕J2​(1). Equivalently, one may ask whether (A−I)2=0, or which of (λ−1)2 and (λ−1)3 is the minimal polynomial.

    (d) Non-diagonalizability from the data alone. mg​(1)=4−rank(A−I)=2 while ma​(1)=4. The eigenvalue 1 is therefore defective, the eigenvectors of A span only a 2-dimensional subspace of C4, and no basis of eigenvectors exists. (Alternatively: A has the single eigenvalue 1, and A=I because rank(A−I)=2=0, so the scalar test applies.) The student's conclusion stands, but it never needed the Jordan form.

  9. 9.Computation

    Let

    A=​321​020​−1−21​​,

    whose characteristic polynomial is (λ−2)3 (you may take this as given, but check it against trA and detA). Write N=A−2I.

    1. (a)Compute rankN and N2. Deduce the number and sizes of the Jordan blocks, the Jordan form J, and mA​.
    2. (b)Build a Jordan basis, longest chain first: choose a suitable v, form its chain, then complete the basis. Write down P with P−1AP=J (you need not compute P−1).
    3. (c)Verify, for each column uj​ of P, the relation Auj​=2uj​ or Auj​=2uj​+uj−1​ that certifies P−1AP=J.
    4. (d)Show that e2​=(0,1,0) is an eigenvector of A that does not lie in the image of N. Explain what goes wrong for someone who starts from e2​ and tries to solve Nv=e2​ for the top of a chain, and what goes wrong for someone who assigns the eigenvector (1,2,1) to the 1×1 block first.

    Solution

    (a) Block structure. Consistency of the given polynomial: trA=6=3⋅2, and expanding detA along the second column, detA=2det(31​−11​)=2⋅4=8=23. Now

    N=A−2I=​121​000​−1−2−1​​.

    Every row is a multiple of (1,0,−1), so rankN=1. Better, N=uwT with u=(1,2,1)T and w=(1,0,−1)T, and wTu=1+0−1=0, so N2=u(wTu)wT=0. (Row by row: row 1 of N2 is row 1 minus row 3 of N, which is 0, and the other rows are multiples of that.)

    Rank sequence (r0​,r1​,r2​)=(3,1,0): b≥1​=2, b≥2​=1, b≥3​=0. Two blocks, one of size 2 and one of size 1:

    J=J2​(2)⊕J1​(2)=​200​120​002​​,mA​(λ)=(λ−2)2,

    the exponent being the largest block (equivalently, N=0=N2).

    (b) The basis, longest chain first. A chain of length 2 needs a top v∈kerN2∖kerN=C3∖kerN, where kerN={x−z=0}. Take v=e1​=(1,0,0): Nv=(1,2,1)=0. The chain, eigenvector first, is

    u1​=Nv=(1,2,1),u2​=v=(1,0,0).

    The 1×1 block needs an eigenvector independent of u1​. Any vector with x=z is in kerN; u3​=(0,1,0) qualifies and is not a multiple of (1,2,1). Independence of all three: expanding along the third row,

    det​121​100​010​​=1⋅det(10​01​)=1=0.

    So

    P=(u1​∣u2​∣u3​)=​121​100​010​​,P−1AP=J.

    (c) The certifying relations. Au1​=(3−1, 2+4−2, 1+1)=(2,4,2)=2u1​. Au2​=Ae1​=(3,2,1)=2(1,0,0)+(1,2,1)=2u2​+u1​. Au3​=Ae2​=(0,2,0)=2u3​. Column j of P−1AP is the coordinate vector of Auj​ in the basis u1​,u2​,u3​, namely 2e1​, e1​+2e2​ and 2e3​ — precisely the columns of J, with the 1 above the diagonal because the chain is listed eigenvector first.

    (d) The two traps. The second column of N is zero, so Ne2​=0 and e2​ is an eigenvector. Since N=uwT, every vector in the image of N is a multiple of u=(1,2,1); e2​ is not, so e2​∈/imN and the equation Nv=e2​ has no solution (concretely, Nv=(v1​−v3​)(1,2,1) can never equal (0,1,0)). Someone working bottom-up from e2​ would be stuck: an eigenvector is the bottom of a longer chain only if it lies in imN∩kerN, which here is the single line span(1,2,1).

    That is also why the short block must be assigned last. If one first reserved (1,2,1) for the 1×1 block, the length-2 chain would need a different bottom in imN∩kerN, and there is none — the long chain could never be completed. Longest chains first, each started from its top, avoids both problems.

  10. 10.Computation

    Let A be the matrix of Problem 5,

    A=​101​110​002​​,

    and take the Jordan basis u1​=(−1,0,1), u2​=(1,−1,0), u3​=(0,0,1), so that P=(u1​∣u2​∣u3​) satisfies P−1AP=J=J2​(1)⊕J1​(2). (If your chain in Problem 5 was different, either basis works; the solution uses this one.)

    1. (a)Write down Jk for general k≥0, quoting the formula for powers of a Jordan block.
    2. (b)Compute P−1, then Ak=PJkP−1 in closed form. Check your formula at k=1 and k=2 by direct multiplication.
    3. (c)Let xk​=Akx0​. For x0​=(1,−1,0) show that ∥xk​∥ grows linearly in k although the only eigenvalue involved is 1; for x0​=(0,1,0) show that it grows like 2k. Explain both behaviours in terms of the Jordan blocks.

    Solution

    (a) Powers of J. For a block Jm​(λ)=λI+S the entry in position (i,i+j) of Jm​(λ)k is (jk​)λk−j. With m=2, λ=1 this gives J2​(1)k=(10​k1​), and the 1×1 block gives 2k. Powers of a block diagonal matrix are block diagonal, so

    Jk=​100​k10​002k​​.

    (b) Conjugating back. detP=1 (Problem 5), so P−1 is the adjugate. The cofactors of P=​−101​1−10​001​​ are C11​=−1, C12​=0, C13​=1, C21​=−1, C22​=−1, C23​=1, C31​=0, C32​=0, C33​=1, and transposing,

    P−1=​−101​−1−11​001​​.

    Check PP−1: row 1 of P against the columns gives (1+0+0, 1−1+0, 0)=(1,0,0); row 2 gives (0,1,0); row 3 gives (−1+1, −1+1, 1)=(0,0,1). Correct.

    Now PJk has columns Pe1​=u1​, P(ke1​+e2​)=ku1​+u2​ and P(2ke3​)=2ku3​:

    PJk=​−101​1−k−1k​002k​​.

    Multiplying by P−1 row by row: row 1 is −1⋅(−1,−1,0)+(1−k)(0,−1,0)=(1,1−(1−k),0)=(1,k,0); row 2 is −1⋅(0,−1,0)=(0,1,0); row 3 is (−1,−1,0)+k(0,−1,0)+2k(1,1,1)=(2k−1, 2k−1−k, 2k). Hence

    Ak=​102k−1​k12k−1−k​002k​​.

    At k=1 this is ​101​110​002​​=A. At k=2 the formula gives rows (1,2,0), (0,1,0), (3,1,4); directly, row 1 of A2 is row 1 + row 2 of A =(1,2,0), row 2 is (0,1,0), and row 3 is row 1 + 2⋅row 3 =(1,1,0)+(2,0,4)=(3,1,4). Agreed. (Also trAk=2+2k, the sum of the k-th powers of the eigenvalues 1,1,2.)

    (c) Two growth rates. For x0​=(1,−1,0)=u2​, the formula gives

    xk​=Aku2​=(1−k, −1, (2k−1)−(2k−1−k))=(1−k,−1,k),

    so ∥xk​∥2=(1−k)2+1+k2∼2k2 and ∥xk​∥∼2​k: linear growth. In Jordan coordinates x0​=e2​ and Jke2​=(k,1,0), i.e. xk​=ku1​+u2​ — the k is the off-diagonal entry (1k​)1k−1 of J2​(1)k. A diagonalizable matrix with eigenvalue 1 would leave x0​ fixed; the size-2 block at eigenvalue 1 is what produces the factor k.

    For x0​=e2​=(0,1,0), xk​ is the second column of Ak, (k,1,2k−1−k), whose last coordinate dominates: ∥xk​∥∼2k. The reason is that e2​ has a component in K2​: P−1e2​ is the second column of P−1, (−1,−1,1), so e2​=−u1​−u2​+u3​ (check: (1,0,−1)+(−1,1,0)+(0,0,1)=(0,1,0)), and the u3​ component is multiplied by 2k at each step while the K1​ component only grows linearly. Growth is governed by the largest eigenvalue present in x0​'s decomposition, with a polynomial factor km−1 from a block of size m.

  11. 11.Critique

    Let A=(10​12​). A student computes eA as follows.

    Step 1. Split A=D+N with D=(10​02​) and N=(00​10​).

    Step 2. N2=0, so eN=I+N=(10​11​), and eD=(e0​0e2​).

    Step 3. Therefore eA=eDeN=(e0​ee2​).

    Steps 1 and 2 are correct.

    1. (a)Identify the error in Step 3, name the property of the exponential being misused, and verify by a direct computation that its hypothesis fails here.
    2. (b)Compute eAt correctly for all t (diagonalize), state eA, and compare its (1,2) entry with the student's. Check your eAt by confirming that its derivative at t=0 equals A.
    3. (c)Now let B=(20​12​) and split it as B=2I+N. Explain why the same three-step method is valid for B, and write down eBt.

    Solution

    (a) The flaw. Step 3 uses eD+N=eDeN. That identity is the property “if MN=NM then eM+N=eMeN”, and it requires the two summands to commute. Here

    DN=(10​02​)(00​10​)=(00​10​),ND=(00​10​)(10​02​)=(00​20​),

    so DN=ND and the property may not be invoked. A symptom the student could have noticed: the other order gives eNeD=(e0​e2e2​), a different matrix, whereas eA is one definite matrix.

    (b) The correct exponential. A is upper triangular with distinct eigenvalues 1 and 2, so it is diagonalizable. Eigenvectors: for λ=1, A−I=(00​11​) gives (1,0); for λ=2, A−2I=(−10​10​) gives (1,1). So

    P=(10​11​),P−1=(10​−11​),P−1AP=(10​02​).

    Then eAt=PeDtP−1 (the exponential of a diagonal matrix is the diagonal of exponentials, and conjugation passes through the series):

    eAt=(10​11​)(et0​0e2t​)(10​−11​)=(et0​e2te2t​)(10​−11​)=(et0​e2t−ete2t​).

    At t=1, eA=(e0​e2−ee2​). The (1,2) entry is e2−e≈4.67, not the student's e≈2.72. (Note that det does not catch the error: both matrices have determinant e3=etrA.)

    Check: dtd​eAt​t=0​=(et0​2e2t−et2e2t​)t=0​=(10​12​)=A, as dtd​eAt=AeAt requires. The student's matrix, read as eAt=(et0​tete2t​) or in any other way, cannot pass this test and the eigenvector test simultaneously: eA must send (1,1) to e2(1,1), and the student's matrix sends it to (2e,e2).

    (c) When the split is legal. B=2I+N, and 2I commutes with every matrix, in particular with N. So the commuting-sum property applies: eBt=e2tIeNt. Now e2tI=e2tI (the series for a scalar matrix is the scalar series), and eNt=I+tN because N2=0 truncates the series. Hence

    eBt=e2t(10​t1​),

    which is the general formula eJ2​(λ)t=eλt(10​t1​) with λ=2. Check: its derivative at t=0 is (20​12​)=B. The difference between A and B is exactly that A's diagonal part is not a scalar, so it fails to commute with the nilpotent part; the diagonal–plus–nilpotent split only exponentiates termwise inside a single Jordan block.

  12. 12.Proof

    Let A be an n×n complex matrix with eigenvalues λ1​,…,λn​ listed with algebraic multiplicity. You may use that A is similar to an upper triangular matrix T whose diagonal is λ1​,…,λn​, that eP−1MP=P−1eMP, and that the series defining eM converges entry by entry.

    1. (a)Prove that if T is upper triangular with diagonal entries t1​,…,tn​, then Tk is upper triangular with diagonal entries t1k​,…,tnk​, and deduce that eT is upper triangular with diagonal entries et1​,…,etn​.
    2. (b)Prove that the eigenvalues of eA, with multiplicity, are eλ1​,…,eλn​, and that deteA=etrA.
    3. (c)Deduce that eA is invertible for every A, that deteA>0 whenever A is real, and that deteAt=1 for all t whenever trA=0.

    Solution

    (a) Triangular powers and the exponential. Let S and T be upper triangular. For i>j the entry (ST)ij​=∑m​Sim​Tmj​ vanishes: Sim​=0 needs m≥i and Tmj​=0 needs m≤j, and no m satisfies i≤m≤j when i>j. For i=j the only surviving term is m=i, giving (ST)ii​=Sii​Tii​. So a product of upper triangular matrices is upper triangular with the products of the diagonals on its diagonal. By induction on k (with T0=I), Tk is upper triangular with diagonal t1k​,…,tnk​.

    The partial sums EK​=∑k=0K​Tk/k! are therefore upper triangular with diagonal entries ∑k=0K​tik​/k!. Passing to the limit entry by entry: every below-diagonal entry is 0 for every K, hence 0 in the limit, and the i-th diagonal entry tends to eti​. So eT is upper triangular with diagonal et1​,…,etn​. ■

    (b) Eigenvalues and determinant of eA. Write A=PTP−1 with T upper triangular and diagonal λ1​,…,λn​. Then eA=PeTP−1, so eA is similar to eT and has the same characteristic polynomial. By (a), eT is upper triangular with diagonal eλ1​,…,eλn​, and the characteristic polynomial of a triangular matrix is ∏i​(λ−eλi​). So the eigenvalues of eA, with multiplicity, are eλ1​,…,eλn​.

    For the determinant, deteA=detP⋅deteT⋅detP−1=deteT, and the determinant of a triangular matrix is the product of its diagonal:

    deteA=i=1∏n​eλi​=eλ1​+⋯+λn​=etrT=etrA,

    using that tr is a similarity invariant (it is minus the coefficient of λn−1 in the characteristic polynomial, which A and T share). ■

    (c) Consequences. etrA=0 for every complex number trA, so deteA=0 and eA is invertible — consistent with the direct argument that e−A is its inverse. If A is real then trA is real and etrA>0; so, for instance, no real matrix has eA=(−10​01​). Finally tr(At)=ttrA, so if trA=0 then deteAt=e0=1 for every t: the flow of x′=Ax preserves volume exactly when the trace vanishes.

← Back to the question sheet

CalcPath·Free university mathematics

SubjectsProblem setsProgressFeedbackAbout