Linear Algebra and Vector Geometry 201-SN4-RE • CEGEP, Montreal

Corrected exercises: eigenvalues and diagonalization (201-SN4-RE)

This is the corrected exercise set on diagonalization, eigenvalues and eigenvectors for 201-SN4-RE (formerly 201-NYC), Linear Algebra and Vector Geometry at CEGEP. Part A sets up the mechanics: the definition Av⃗=λv⃗A\vec{v}=\lambda\vec{v} and what it means geometrically, the characteristic polynomial det⁡(A−λI)=0\det(A-\lambda I)=0, eigenspaces found by solving a homogeneous system, and the construction of PP and DD such that A=PDP−1A=PDP^{-1}. Part B rises to exam level: the exact diagonalizability criterion with a matrix that fails the test, AnA^{n} for large powers, a full problem in which two coupled populations decouple in the eigenbasis, the orthogonal diagonalization of symmetric matrices, the Cayley-Hamilton theorem to compute an inverse and powers, and the eigenvalues of the matrices AkA^{k}, A−1A^{-1} and A+cIA+cI. Part C adds five classics: a 3×33\times 3 matrix with a double eigenvalue, eigenvalues read off the geometry, diagonalizability depending on a parameter, the normal modes of two carts on springs and the principal axes of a tilted ellipse.

Exercise 0.1 opens the set with ten eigenvalues to compute directly, as in the textbook drill lists: diagonal and triangular matrices read off their diagonal, the characteristic polynomial λ2−tr⁡(A)λ+det⁡(A)\lambda^{2}-\operatorname{tr}(A)\lambda+\det(A), 2×22\times 2 matrices with integer roots, a symmetric one, one with a double eigenvalue and the dimension of its eigenspace, one with complex eigenvalues such as 2±3i2\pm 3i, a 3×33\times 3 block matrix, and the eigenvalue read off Av⃗A\vec{v}, each with a value to find and checked on the spot.

This set is written for CEGEP students in Science and in the Social Science mathematics profile, whose linear algebra course is 201-105. Depending on the CEGEP, diagonalization is taught at the end of 201-SN4-RE, in 201-105 or only in first-year university: check your teacher's course outline before working through this set.

The guiding idea deserves to be stated up front, because it gives every calculation its meaning: an arbitrary matrix mixes the coordinates, but in the RIGHT basis the same transformation reduces to stretches along independent axes. Diagonalizing means finding that basis. Everything else, the characteristic polynomial included, is only the procedure for getting there.

16 corrected exercises • 160 points • 225 minutes

This chapter is part of 201-SN4-RE, every chapter with its set and sheet Before you start Revision sheet: the traps and the method for this chapter

Self-checking set Type your answers and check them question by question: the page tells you right or wrong before you open the solution. With an account, every answer right on the first try earns points, and your right answers stay filled in when you come back.

Course recap

  • • Definition: λ\lambda is an EIGENVALUE of AA if there is a NON-ZERO vector v⃗\vec{v} such that Av⃗=λv⃗A\vec{v}=\lambda\vec{v}. That v⃗\vec{v} is an EIGENVECTOR associated with λ\lambda. Geometrically, AA only STRETCHES v⃗\vec{v}, without changing its direction.
  • • CHARACTERISTIC POLYNOMIAL: the eigenvalues are the roots of det⁡(A−λI)=0\det(A-\lambda I)=0. For a 2×22\times 2 matrix it reads λ2−tr⁡(A)λ+det⁡(A)=0\lambda^{2}-\operatorname{tr}(A)\lambda+\det(A)=0, which gives an immediate check.
  • • EIGENSPACE associated with λ\lambda: Eλ=ker⁡(A−λI)E_{\lambda}=\ker(A-\lambda I), found by solving the homogeneous system (A−λI)v⃗=0⃗(A-\lambda I)\vec{v}=\vec{0}. It is always a subspace, of dimension at least 1.
  • • ALGEBRAIC multiplicity of λ\lambda: its order as a root of the characteristic polynomial. GEOMETRIC multiplicity: dim⁡Eλ\dim E_{\lambda}. Always, geometric ≤\leq algebraic.
  • • DIAGONALIZABILITY CRITERION: AA (of size nn) is diagonalizable if and only if the dimensions of the eigenspaces add up to nn, that is, if for EVERY eigenvalue the geometric multiplicity equals the algebraic multiplicity.
  • • The easy case: if AA has nn DISTINCT eigenvalues, it is automatically diagonalizable. Every real SYMMETRIC matrix is diagonalizable, with orthogonal eigenvectors.
  • • Diagonalization: A=PDP−1A=PDP^{-1}, where the COLUMNS of PP are the eigenvectors and DD is diagonal with the eigenvalues in the SAME ORDER. Consequence: An=PDnP−1A^{n}=PD^{n}P^{-1}, and DnD^{n} is obtained by raising each diagonal entry to the power nn.

Calculations: eigenvalues to compute

Exercise 0.1: Computing eigenvalues: ten direct calculations

For each line, compute the eigenvalues of the matrix, or the quantity asked for, then give the values indicated. For a complex case, give the real part and the positive imaginary part.

One equation covers every line: det⁡(A−λI)=0\det(A-\lambda I)=0. For a 2×22\times 2 matrix it reads directly λ2−tr⁡(A) λ+det⁡(A)=0\lambda^{2}-\operatorname{tr}(A)\,\lambda+\det(A)=0, and you solve that quadratic. But look at the structure first: a diagonal or triangular matrix carries its eigenvalues on its diagonal, a block matrix is handled block by block, and when an eigenvector is known, Av⃗A\vec{v} gives the eigenvalue without any determinant.

  • a) (−3007)\begin{pmatrix} -3 & 0 \\ 0 & 7 \end{pmatrix}: give the largest and the smallest eigenvalue.
  • b) (6−40−1)\begin{pmatrix} 6 & -4 \\ 0 & -1 \end{pmatrix}: give the largest and the smallest eigenvalue.
  • c) (5−231)\begin{pmatrix} 5 & -2 \\ 3 & 1 \end{pmatrix}: write the characteristic polynomial in the form λ2+bλ+c\lambda^{2}+b\lambda+c, then give bb and cc.
  • d) (1423)\begin{pmatrix} 1 & 4 \\ 2 & 3 \end{pmatrix}: give the largest and the smallest eigenvalue.
  • e) (233−6)\begin{pmatrix} 2 & 3 \\ 3 & -6 \end{pmatrix}: give the largest and the smallest eigenvalue.
  • f) (4−112)\begin{pmatrix} 4 & -1 \\ 1 & 2 \end{pmatrix}: give its eigenvalue, which is a double root, then the dimension of the associated eigenspace.
  • g) (1−523)\begin{pmatrix} 1 & -5 \\ 2 & 3 \end{pmatrix}: its eigenvalues are complex; give their real part and their positive imaginary part.
  • h) (27−10−43005)\begin{pmatrix} 2 & 7 & -1 \\ 0 & -4 & 3 \\ 0 & 0 & 5 \end{pmatrix}: give the largest and the smallest eigenvalue.
  • i) (40002104−1)\begin{pmatrix} 4 & 0 & 0 \\ 0 & 2 & 1 \\ 0 & 4 & -1 \end{pmatrix}: give the largest and the smallest eigenvalue.
  • j) A=(21112−1202)A=\begin{pmatrix} 2 & 1 & 1 \\ 1 & 2 & -1 \\ 2 & 0 & 2 \end{pmatrix} and v⃗=(1−12)\vec{v}=\begin{pmatrix} 1 \\ -1 \\ 2 \end{pmatrix}: compute Av⃗A\vec{v}, then give the eigenvalue associated with v⃗\vec{v}.

Type your answers, the page tells you right or wrong 0/19

a)
b)
c)
d)
e)
f)
g)
h)
i)
j)
Show the solution

Answers

  • a) λ=7\lambda=7 and λ=−3\lambda=-3
  • b) (6−λ)(−1−λ)(6-\lambda)(-1-\lambda); 66 and −1-1
  • c) λ2−6λ+11\lambda^{2}-6\lambda+11; b=−6b=-6, c=11c=11
  • d) (λ−5)(λ+1)(\lambda-5)(\lambda+1); 55 and −1-1
  • e) (λ+7)(λ−3)(\lambda+7)(\lambda-3); 33 and −7-7
  • f) (λ−3)2(\lambda-3)^{2}; λ=3\lambda=3, dim⁡E3=1\dim E_{3}=1
  • g) λ=2±3i\lambda=2\pm 3i; 22 and 33
  • h) 22, −4-4, 55; 55 and −4-4
  • i) (4−λ)(λ−3)(λ+2)(4-\lambda)(\lambda-3)(\lambda+2); 44 and −2-2
  • j) Av⃗=3v⃗A\vec{v}=3\vec{v}; λ=3\lambda=3

The move that decides everything: look at the shape of the matrix BEFORE writing det⁡(A−λI)\det(A-\lambda I). Diagonal, triangular or block, the answer is read off; otherwise the trace and the determinant give the polynomial in one line, and the sum and product of the roots serve as a check.

a)

Diagonal matrix: det⁡(A−λI)=(−3−λ)(7−λ)\det(A-\lambda I)=(-3-\lambda)(7-\lambda), so the eigenvalues are the diagonal entries, 77 and −3-3.

b)

Triangular matrix: det⁡(A−λI)=(6−λ)(−1−λ)\det(A-\lambda I)=(6-\lambda)(-1-\lambda), because the −4-4 does not enter the product of the diagonal. Eigenvalues 66 and −1-1. Pitfall: bringing the −4-4 into the calculation.

c)

tr⁡=5+1=6\operatorname{tr}=5+1=6 and det⁡=5(1)−(−2)(3)=11\det=5(1)-(-2)(3)=11, so λ2−6λ+11\lambda^{2}-6\lambda+11: b=−6b=-6 and c=11c=11. Pitfall: bb carries the MINUS sign of the trace, and det⁡=5+6\det=5+6, not 5−65-6. The discriminant 36−44<036-44<0 announces complex eigenvalues.

d)

tr⁡=4\operatorname{tr}=4 and det⁡=3−8=−5\det=3-8=-5, so λ2−4λ−5=(λ−5)(λ+1)\lambda^{2}-4\lambda-5=(\lambda-5)(\lambda+1): eigenvalues 55 and −1-1. Check: 5+(−1)=45+(-1)=4 and 5×(−1)=−55\times(-1)=-5.

e)

tr⁡=−4\operatorname{tr}=-4 and det⁡=−12−9=−21\det=-12-9=-21, so λ2+4λ−21=(λ+7)(λ−3)\lambda^{2}+4\lambda-21=(\lambda+7)(\lambda-3): eigenvalues 33 and −7-7, real, as for every symmetric matrix.

f)

tr⁡=6\operatorname{tr}=6 and det⁡=8+1=9\det=8+1=9, so λ2−6λ+9=(λ−3)2\lambda^{2}-6\lambda+9=(\lambda-3)^{2}: double eigenvalue 33. Then A−3I=(1−11−1)A-3I=\begin{pmatrix} 1 & -1 \\ 1 & -1 \end{pmatrix}, of rank 11, so dim⁡E3=2−1=1\dim E_{3}=2-1=1. Pitfall: believing that the dimension equals the multiplicity 22; here it is 11, and the matrix is not diagonalizable.

g)

tr⁡=4\operatorname{tr}=4 and det⁡=3+10=13\det=3+10=13, so λ2−4λ+13=0\lambda^{2}-4\lambda+13=0, with discriminant 16−52=−3616-52=-36. Hence λ=4±6i2=2±3i\lambda=\frac{4\pm 6i}{2}=2\pm 3i: real part 22, imaginary part 33. Pitfall: forgetting to divide 6i6i by 22.

h)

Upper triangular matrix: the eigenvalues are 22, −4-4 and 55. The largest is 55, the smallest −4-4.

i)

Block matrix: det⁡(A−λI)=(4−λ)[(2−λ)(−1−λ)−4]=(4−λ)(λ2−λ−6)=(4−λ)(λ−3)(λ+2)\det(A-\lambda I)=(4-\lambda)\left[(2-\lambda)(-1-\lambda)-4\right]=(4-\lambda)(\lambda^{2}-\lambda-6)=(4-\lambda)(\lambda-3)(\lambda+2). Eigenvalues 44, 33 and −2-2: the largest is 44, the smallest −2-2. Pitfall: forgetting the 1×11\times 1 block, which contributes 44.

j)

Av⃗=(2−1+21−2−22+0+4)=(3−36)=3v⃗A\vec{v}=\begin{pmatrix} 2-1+2 \\ 1-2-2 \\ 2+0+4 \end{pmatrix}=\begin{pmatrix} 3 \\ -3 \\ 6 \end{pmatrix}=3\vec{v}: the eigenvalue is 33. No characteristic polynomial is needed when the eigenvector is given.

Check: for a 2×22\times 2 matrix, the sum of the eigenvalues is the trace and their product the determinant. For a 3×33\times 3 matrix the sum is still the trace: in i), 4+3−2=5=4+2−14+3-2=5=4+2-1.

Tick the exercises you have done or want to review: a free account, no password, keeps your ticks and your right answers from one visit to the next, tells you which chapter to tackle next and lets you ask for the exercise you are missing. Create your space, an email is enough.

Part A: the basics (/50)

Exercise 1: Eigenvalues and eigenvectors: the definition

An eigenvector is a vector that the matrix only stretches, without turning it. Everything starts from this geometric picture.

  • a) Check that v⃗=(11)\vec{v}=\begin{pmatrix}1\\1\end{pmatrix} is an eigenvector of A=(3113)A=\begin{pmatrix}3&1\\1&3\end{pmatrix} and give the associated eigenvalue.
  • b) Check the same for w⃗=(1−1)\vec{w}=\begin{pmatrix}1\\-1\end{pmatrix}. What do you notice about the directions of v⃗\vec{v} and w⃗\vec{w}?
  • c) Is the vector (10)\begin{pmatrix}1\\0\end{pmatrix} an eigenvector of AA? Justify by calculation.
  • d) Why does the definition require v⃗≠0⃗\vec{v}\neq\vec{0}? What would happen without this condition?

Type your answers, the page tells you right or wrong 0/5

a)
b)
c)
d)
Show the solution

Answers

  • a) Av⃗=4v⃗A\vec v=4\vec v: λ=4\lambda=4
  • b) Aw⃗=2w⃗A\vec w=2\vec w; v⃗⊥w⃗\vec v\perp\vec w
  • c) No: (3,1)(3, 1) is not parallel to (1,0)(1, 0)
  • d) 0⃗\vec 0 would make every λ\lambda an eigenvalue
a)

Compute Av⃗=(3113)(11)=(3+11+3)=(44)=4(11)=4v⃗A\vec{v}=\begin{pmatrix}3&1\\1&3\end{pmatrix}\begin{pmatrix}1\\1\end{pmatrix}=\begin{pmatrix}3+1\\1+3\end{pmatrix}=\begin{pmatrix}4\\4\end{pmatrix}=4\begin{pmatrix}1\\1\end{pmatrix}=4\vec{v}. The relation Av⃗=λv⃗A\vec{v}=\lambda\vec{v} holds with λ=4\lambda=4: v⃗\vec{v} is indeed an eigenvector, with eigenvalue 4. Geometrically, the matrix multiplies by 4 every vector pointing along the diagonal y=xy=x, without deflecting it.

b)

Aw⃗=(3113)(1−1)=(3−11−3)=(2−2)=2w⃗A\vec{w}=\begin{pmatrix}3&1\\1&3\end{pmatrix}\begin{pmatrix}1\\-1\end{pmatrix}=\begin{pmatrix}3-1\\1-3\end{pmatrix}=\begin{pmatrix}2\\-2\end{pmatrix}=2\vec{w}. So w⃗\vec{w} is an eigenvector with eigenvalue 2. About the directions: v⃗⋅w⃗=(1)(1)+(1)(−1)=0\vec{v}\cdot\vec{w}=(1)(1)+(1)(-1)=0, the two eigenvectors are ORTHOGONAL. This is no accident: AA is symmetric (a12=a21=1a_{12}=a_{21}=1), and the eigenvectors of a symmetric matrix associated with distinct eigenvalues are always orthogonal. Exercise 8 comes back to this.

c)

A(10)=(31)A\begin{pmatrix}1\\0\end{pmatrix}=\begin{pmatrix}3\\1\end{pmatrix}. For this to be an eigenvector we would need (31)=λ(10)=(λ0)\begin{pmatrix}3\\1\end{pmatrix}=\lambda\begin{pmatrix}1\\0\end{pmatrix}=\begin{pmatrix}\lambda\\0\end{pmatrix}, which requires both λ=3\lambda=3 and 1=01=0: impossible. So (10)\begin{pmatrix}1\\0\end{pmatrix} is NOT an eigenvector. Geometrically, AA turns this vector as well as lengthening it, since the image (31)\begin{pmatrix}3\\1\end{pmatrix} is not parallel to the original vector. The practical test is therefore simple: compute Au⃗A\vec{u} and check whether the result is PROPORTIONAL to u⃗\vec{u}.

d)

Without the condition v⃗≠0⃗\vec{v}\neq\vec{0}, the definition would lose all content: A0⃗=0⃗=λ0⃗A\vec{0}=\vec{0}=\lambda\vec{0} holds for EVERY scalar λ\lambda, so any number would be an eigenvalue of any matrix. The notion would no longer distinguish anything. This is also why eigenspaces are defined as NULL SPACES: Eλ=ker⁡(A−λI)E_{\lambda}=\ker(A-\lambda I) necessarily contains 0⃗\vec{0} (every subspace does), and the condition “λ\lambda is an eigenvalue” translates exactly into “this null space contains SOMETHING OTHER than 0⃗\vec{0}”, that is dim⁡Eλ≥1\dim E_{\lambda}\geq 1, or again det⁡(A−λI)=0\det(A-\lambda I)=0. Here you see why the determinant appears: it vanishes precisely when the homogeneous system has non-trivial solutions.

-112345-3-2-112345Av = 4vAw = 2wvwA(1, 0) = (3, 1)

Exercise 2: The characteristic polynomial

Finding the eigenvalues means finding for which λ\lambda the system (A−λI)v⃗=0⃗(A-\lambda I)\vec{v}=\vec{0} has a non-trivial solution. The determinant decides.

  • a) Explain why the condition “λ\lambda is an eigenvalue of AA” is equivalent to det⁡(A−λI)=0\det(A-\lambda I)=0.
  • b) Find the eigenvalues of A=(4213)A=\begin{pmatrix}4&2\\1&3\end{pmatrix}.
  • c) Check your results with the two relations λ1+λ2=tr⁡(A)\lambda_{1}+\lambda_{2}=\operatorname{tr}(A) and λ1λ2=det⁡(A)\lambda_{1}\lambda_{2}=\det(A).
  • d) Find the eigenvalues of B=(20013045−1)B=\begin{pmatrix}2&0&0\\1&3&0\\4&5&-1\end{pmatrix}. What feature makes the calculation immediate?

Type your answers, the page tells you right or wrong 0/6

b)
c)
d)
Show the solution

Answers

  • a) Non-trivial solutions: det⁡(A−λI)=0\det(A-\lambda I)=0
  • b) λ=2\lambda=2 and λ=5\lambda=5
  • c) 2+5=72+5=7, 2×5=102\times 5=10
  • d) Triangular: 22, 33, −1-1
a)

By definition, λ\lambda is an eigenvalue if there is a v⃗≠0⃗\vec{v}\neq\vec{0} with Av⃗=λv⃗A\vec{v}=\lambda\vec{v}, which reads Av⃗−λv⃗=0⃗A\vec{v}-\lambda\vec{v}=\vec{0}, that is (A−λI)v⃗=0⃗(A-\lambda I)\vec{v}=\vec{0} (inserting II is needed for the subtraction to make sense as a matrix operation). The question therefore becomes: does the homogeneous system (A−λI)v⃗=0⃗(A-\lambda I)\vec{v}=\vec{0} have a NON-TRIVIAL solution? A square homogeneous system has non-trivial solutions if and only if its matrix is not invertible, that is, if and only if its determinant is zero. Hence the equivalence λ\lambda eigenvalue   ⟺  det⁡(A−λI)=0\iff\det(A-\lambda I)=0 ✓. The whole chapter rests on this translation, and it reuses directly the invertibility criterion from the determinants chapter.

b)

A−λI=(4−λ213−λ)A-\lambda I=\begin{pmatrix}4-\lambda&2\\1&3-\lambda\end{pmatrix}, with determinant (4−λ)(3−λ)−2=12−4λ−3λ+λ2−2=λ2−7λ+10(4-\lambda)(3-\lambda)-2=12-4\lambda-3\lambda+\lambda^{2}-2=\lambda^{2}-7\lambda+10. Setting it to zero: λ2−7λ+10=0\lambda^{2}-7\lambda+10=0, that is (λ−2)(λ−5)=0(\lambda-2)(\lambda-5)=0. The eigenvalues are therefore λ1=2\lambda_{1}=2 and λ2=5\lambda_{2}=5.

c)

Trace: tr⁡(A)=4+3=7\operatorname{tr}(A)=4+3=7, and λ1+λ2=2+5=7\lambda_{1}+\lambda_{2}=2+5=7 ✓. Determinant: det⁡(A)=4×3−2×1=10\det(A)=4\times 3-2\times 1=10, and λ1λ2=2×5=10\lambda_{1}\lambda_{2}=2\times 5=10 ✓. These two relations are simply the coefficients of the characteristic polynomial read off, since it always reads λ2−tr⁡(A)λ+det⁡(A)\lambda^{2}-\operatorname{tr}(A)\lambda+\det(A) for a 2×22\times 2 matrix. They give a five-second check, and even let you guess the eigenvalues when they are integers: here, looking for two numbers with sum 7 and product 10 gives 2 and 5 at once. This is the check to make every time on the exam.

d)

The matrix BB is LOWER TRIANGULAR (all the entries above the diagonal are zero). The determinant of a triangular matrix is the product of its diagonal entries, so det⁡(B−λI)=(2−λ)(3−λ)(−1−λ)\det(B-\lambda I)=(2-\lambda)(3-\lambda)(-1-\lambda), which vanishes for λ=2\lambda=2, λ=3\lambda=3 and λ=−1\lambda=-1. The eigenvalues of a triangular matrix are therefore exactly its DIAGONAL entries, with no calculation at all. The same holds for upper triangular matrices and, as a limiting case, for diagonal matrices. This fits the goal of the chapter: diagonalizing a matrix means bringing it to a form where its eigenvalues are read directly off the diagonal.

Exercise 3: Eigenspaces

Once the eigenvalues are known, each λ\lambda gives a homogeneous system to solve. The subspace obtained is the eigenspace, and its DIMENSION is what decides everything in the next exercises.

  • a) For A=(4213)A=\begin{pmatrix}4&2\\1&3\end{pmatrix} (eigenvalues 2 and 5), find E2=ker⁡(A−2I)E_{2}=\ker(A-2I) and give a basis of it.
  • b) Likewise, find E5E_{5} and give a basis of it.
  • c) Check directly that your two eigenvectors satisfy Av⃗=λv⃗A\vec{v}=\lambda\vec{v}. Do they form a basis of R2\mathbb{R}^{2}?
  • d) Why is an eigenvector never unique? What can be said about all the eigenvectors associated with the same eigenvalue?

Type your answers, the page tells you right or wrong 0/5

a)
b)
c)
d)
Show the solution

Answers

  • a) E2=span⁡{(−1,1)}E_2=\operatorname{span}\{(-1, 1)\}
  • b) E5=span⁡{(2,1)}E_5=\operatorname{span}\{(2, 1)\}
  • c) Determinant −3-3: a basis of eigenvectors
  • d) The whole direction is an eigendirection
a)

A−2I=(2211)A-2I=\begin{pmatrix}2&2\\1&1\end{pmatrix}. The system (A−2I)v⃗=0⃗(A-2I)\vec{v}=\vec{0} reads 2x+2y=02x+2y=0 and x+y=0x+y=0: the two equations are proportional and reduce to x=−yx=-y. Setting y=ty=t: v⃗=(−t,t)=t(−1,1)\vec{v}=(-t,t)=t(-1,1). So E2=span⁡{(−11)}E_{2}=\operatorname{span}\left\{\begin{pmatrix}-1\\1\end{pmatrix}\right\}, a line, and dim⁡E2=1\dim E_{2}=1. A useful consistency check: the system MUST be degenerate (two proportional equations), otherwise the only solution would be 0⃗\vec{0} and λ=2\lambda=2 would not be an eigenvalue. If the row reduction does not degenerate, there is a mistake in the characteristic polynomial.

b)

A−5I=(−121−2)A-5I=\begin{pmatrix}-1&2\\1&-2\end{pmatrix}. The system gives −x+2y=0-x+2y=0, that is x=2yx=2y. Setting y=ty=t: v⃗=(2t,t)=t(2,1)\vec{v}=(2t,t)=t(2,1). So E5=span⁡{(21)}E_{5}=\operatorname{span}\left\{\begin{pmatrix}2\\1\end{pmatrix}\right\} and dim⁡E5=1\dim E_{5}=1. Here again the two rows are proportional ✓.

c)

Check for λ=2\lambda=2: A(−11)=(−4+2−1+3)=(−22)=2(−11)A\begin{pmatrix}-1\\1\end{pmatrix}=\begin{pmatrix}-4+2\\-1+3\end{pmatrix}=\begin{pmatrix}-2\\2\end{pmatrix}=2\begin{pmatrix}-1\\1\end{pmatrix} ✓. For λ=5\lambda=5: A(21)=(8+22+3)=(105)=5(21)A\begin{pmatrix}2\\1\end{pmatrix}=\begin{pmatrix}8+2\\2+3\end{pmatrix}=\begin{pmatrix}10\\5\end{pmatrix}=5\begin{pmatrix}2\\1\end{pmatrix} ✓. Do they form a basis of R2\mathbb{R}^{2}? There are two of them in a space of dimension 2, so it is enough to check independence: the determinant ∣−1211∣=−1−2=−3≠0\begin{vmatrix}-1&2\\1&1\end{vmatrix}=-1-2=-3\neq 0, they are independent and do form a BASIS. This is the fact that will make AA diagonalizable: there is a basis of R2\mathbb{R}^{2} made entirely of eigenvectors.

d)

An eigenvector is never unique because every non-zero multiple of an eigenvector is still an eigenvector for the same eigenvalue: if Av⃗=λv⃗A\vec{v}=\lambda\vec{v}, then A(cv⃗)=cAv⃗=cλv⃗=λ(cv⃗)A(c\vec{v})=cA\vec{v}=c\lambda\vec{v}=\lambda(c\vec{v}) for every c≠0c\neq 0 ✓. Geometrically, it is the DIRECTION that is an eigendirection, not a particular vector in it. More generally, the set of eigenvectors associated with λ\lambda, together with 0⃗\vec{0}, is exactly the subspace Eλ=ker⁡(A−λI)E_{\lambda}=\ker(A-\lambda I): it is therefore closed under addition and scalar multiplication. The practical consequence on the exam: two students can give (−11)\begin{pmatrix}-1\\1\end{pmatrix} and (3−3)\begin{pmatrix}3\\-3\end{pmatrix} as the eigenvector and both be right. Only the direction counts, which is why the expected answer is a BASIS of the eigenspace rather than one imposed vector.

Exercise 4: Diagonalizing: building P and D

Diagonalizing means rewriting the matrix in the basis of its eigenvectors. The only technical difficulty is keeping the ORDER of the columns.

  • a) For A=(4213)A=\begin{pmatrix}4&2\\1&3\end{pmatrix}, write PP and DD such that A=PDP−1A=PDP^{-1}.
  • b) Compute P−1P^{-1} and check explicitly that PDP−1=APDP^{-1}=A.
  • c) What happens if you swap the two columns of PP without touching DD? Illustrate the problem.
  • d) Interpret the equality A=PDP−1A=PDP^{-1} geometrically: what does each of the three matrices do, read from right to left?

Type your answers, the page tells you right or wrong 0/7

a)
b)
c)
d)
Show the solution

Answers

  • a) P=(−1211)P=\begin{pmatrix}-1&2\\1&1\end{pmatrix}, D=diag⁡(2,5)D=\operatorname{diag}(2, 5)
  • b) P−1=(−13231313)P^{-1}=\begin{pmatrix}-\frac13&\frac23\\\frac13&\frac13\end{pmatrix}; PDP−1=APDP^{-1}=A
  • c) P′DP′−1=(3−2−14)≠AP'DP'^{-1}=\begin{pmatrix}3&-2\\-1&4\end{pmatrix}\neq A
  • d) Change basis, stretch, go back
a)

Put the eigenvectors as the COLUMNS of PP and the corresponding eigenvalues on the diagonal of DD, in the same order. With v⃗2=(−11)\vec{v}_{2}=\begin{pmatrix}-1\\1\end{pmatrix} for λ=2\lambda=2 and v⃗5=(21)\vec{v}_{5}=\begin{pmatrix}2\\1\end{pmatrix} for λ=5\lambda=5: P=(−1211)P=\begin{pmatrix}-1&2\\1&1\end{pmatrix} and D=(2005)D=\begin{pmatrix}2&0\\0&5\end{pmatrix}.

b)

det⁡P=(−1)(1)−(2)(1)=−3\det P=(-1)(1)-(2)(1)=-3, so P−1=1−3(1−2−1−1)=(−13231313)P^{-1}=\frac{1}{-3}\begin{pmatrix}1&-2\\-1&-1\end{pmatrix}=\begin{pmatrix}-\frac{1}{3}&\frac{2}{3}\\\frac{1}{3}&\frac{1}{3}\end{pmatrix}. Check: PD=(−1211)(2005)=(−21025)PD=\begin{pmatrix}-1&2\\1&1\end{pmatrix}\begin{pmatrix}2&0\\0&5\end{pmatrix}=\begin{pmatrix}-2&10\\2&5\end{pmatrix}. Then PDP−1=(−21025)(−13231313)=(23+103−43+103−23+5343+53)=(4213)=APDP^{-1}=\begin{pmatrix}-2&10\\2&5\end{pmatrix}\begin{pmatrix}-\frac{1}{3}&\frac{2}{3}\\\frac{1}{3}&\frac{1}{3}\end{pmatrix}=\begin{pmatrix}\frac{2}{3}+\frac{10}{3}&-\frac{4}{3}+\frac{10}{3}\\-\frac{2}{3}+\frac{5}{3}&\frac{4}{3}+\frac{5}{3}\end{pmatrix}=\begin{pmatrix}4&2\\1&3\end{pmatrix}=A ✓.

c)

If you swap the columns of PP to get P′=(2−111)P'=\begin{pmatrix}2&-1\\1&1\end{pmatrix} while keeping D=(2005)D=\begin{pmatrix}2&0\\0&5\end{pmatrix}, the equality fails: you would be claiming that (21)\begin{pmatrix}2\\1\end{pmatrix} goes with the eigenvalue 2, when it goes with 5. The calculation confirms it, P′DP′−1=(3−2−14)≠AP'DP'^{-1}=\begin{pmatrix}3&-2\\-1&4\end{pmatrix}\neq A. The rule is therefore strict: the ii-th column of PP must be an eigenvector associated with the ii-th diagonal entry of DD. On the other hand, swapping the columns of PP AND the entries of DD together is perfectly legal and gives another diagonalization, just as valid: diagonalization is not unique, only the PAIRING counts.

d)

Read from right to left on a vector x⃗\vec{x}, the equality Ax⃗=PDP−1x⃗A\vec{x}=PDP^{-1}\vec{x} breaks into three steps. First P−1P^{-1} CHANGES THE BASIS: it expresses x⃗\vec{x} in the basis of eigenvectors, that is, it gives the coordinates of x⃗\vec{x} along the two preferred directions. Next DD acts in that basis, and its action is as simple as possible: it multiplies the first coordinate by 2 and the second by 5, with no mixing at all, since DD is diagonal. Finally PP GOES BACK to the standard basis. The complicated transformation AA is therefore a plain double stretch, provided you look at it in the right frame. That is the whole point of the chapter: AA mixes the coordinates only because it is observed in a badly chosen basis, and diagonalizing amounts to choosing the frame in which the transformation becomes transparent.

Exercise 5: Complex eigenvalues: a matrix with no real eigendirection

The characteristic polynomial of a real matrix can have complex roots. The matrix then has no real eigenvector, and that has a very clear geometric reading.

  • a) Find the eigenvalues of the rotation A=(0−110)A=\begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix} (rotation through 90∘90^{\circ}).
  • b) Explain geometrically why this matrix has no real eigenvector.
  • c) Find the eigenvalues of B=(2−112)B=\begin{pmatrix} 2 & -1 \\ 1 & 2 \end{pmatrix}.
  • d) Check on BB that the relations λ1+λ2=tr⁡(B)\lambda_{1}+\lambda_{2}=\operatorname{tr}(B) and λ1λ2=det⁡(B)\lambda_{1}\lambda_{2}=\det(B) still hold with complex conjugate eigenvalues.

Type your answers, the page tells you right or wrong 0/6

a)
b)
c)
d)
Show the solution

Answers

  • a) λ=±i\lambda=\pm i
  • b) Rotation: no direction preserved
  • c) 2±i2\pm i
  • d) Sum 4, product 5
a)

The characteristic polynomial is det⁡(A−λI)=∣−λ−11−λ∣=λ2+1\det(A-\lambda I)=\begin{vmatrix} -\lambda & -1 \\ 1 & -\lambda \end{vmatrix}=\lambda^{2}+1. Setting it to zero: λ2=−1\lambda^{2}=-1, so λ=i\lambda=i and λ=−i\lambda=-i. Both eigenvalues are purely imaginary and conjugate to each other. Note that tr⁡(A)=0=i+(−i)\operatorname{tr}(A)=0=i+(-i) and det⁡(A)=1=i⋅(−i)\det(A)=1=i\cdot(-i).

b)

A real eigenvector would be a non-zero real vector whose direction is preserved by AA. But AA is a rotation through 90∘90^{\circ}: it turns EVERY vector of the plane by a quarter turn, so no non-zero real vector keeps its direction. There can therefore be no real eigenvector, which is consistent with the absence of a real eigenvalue. Eigendirections only appear once you move to the complex plane.

c)

det⁡(B−λI)=∣2−λ−112−λ∣=(2−λ)2+1=λ2−4λ+5\det(B-\lambda I)=\begin{vmatrix} 2-\lambda & -1 \\ 1 & 2-\lambda \end{vmatrix}=(2-\lambda)^{2}+1=\lambda^{2}-4\lambda+5. The discriminant is 16−20=−4<016-20=-4<0, so the roots are complex: λ=4±−42=4±2i2=2±i\lambda=\dfrac{4\pm\sqrt{-4}}{2}=\dfrac{4\pm 2i}{2}=2\pm i. The eigenvalues are 2+i2+i and 2−i2-i.

d)

Sum: (2+i)+(2−i)=4=tr⁡(B)(2+i)+(2-i)=4=\operatorname{tr}(B) ✓. Product: (2+i)(2−i)=4−i2=4+1=5=det⁡(B)(2+i)(2-i)=4-i^{2}=4+1=5=\det(B) ✓. The two relations, proved for real eigenvalues, remain true for complex eigenvalues, because they are only the coefficients of the characteristic polynomial λ2−tr⁡(B)λ+det⁡(B)\lambda^{2}-\operatorname{tr}(B)\lambda+\det(B) read off. Note that the common real part, 22, is half the trace: this is general for a conjugate pair, and it gives a quick check. The matrix BB acts as a rotation combined with a scaling, which the link with complex numbers of modulus 5\sqrt{5} will confirm.

vAvwAwno direction preserved

Part B: exam level (/50)

Exercise 6: The diagonalizability criterion and the matrix that fails it

Not every matrix is diagonalizable, and the exam almost always tests this point. The criterion comes down to comparing two multiplicities.

  • a) State the exact diagonalizability criterion, distinguishing algebraic multiplicity from geometric multiplicity.
  • b) Show that A=(3103)A=\begin{pmatrix}3&1\\0&3\end{pmatrix} is NOT diagonalizable, by computing both multiplicities.
  • c) Show, on the contrary, that B=(3003)B=\begin{pmatrix}3&0\\0&3\end{pmatrix} is diagonalizable, although it also has a double eigenvalue. What does the difference between AA and BB change?
  • d) Prove that if an n×nn\times n matrix has nn DISTINCT eigenvalues, then it is diagonalizable. Use this result to conclude, without calculation, about (570−2)\begin{pmatrix}5&7\\0&-2\end{pmatrix}.

Type your answers, the page tells you right or wrong 0/7

a)
b)
c)
d)
Show the solution

Answers

  • a) Geometric = algebraic for each λ\lambda
  • b) Multiplicities 2 and 1: no
  • c) E3=R2E_3=\mathbb{R}^2: yes
  • d) Distinct eigenvalues: yes
a)

ALGEBRAIC multiplicity of λ\lambda: its multiplicity as a root of the characteristic polynomial. GEOMETRIC multiplicity of λ\lambda: the dimension of the eigenspace Eλ=ker⁡(A−λI)E_{\lambda}=\ker(A-\lambda I), in other words the number of independent eigenvectors associated with λ\lambda. Always, 1≤1\leq geometric ≤\leq algebraic. CRITERION: AA of size nn is diagonalizable if and only if, for EVERY eigenvalue, the geometric multiplicity equals the algebraic multiplicity; equivalently, if and only if the dimensions of the eigenspaces add up to nn, that is, if there is a basis of the space made of eigenvectors.

b)

AA is triangular, so its eigenvalues are its diagonal entries: λ=3\lambda=3, with ALGEBRAIC multiplicity 2 (the characteristic polynomial is (3−λ)2(3-\lambda)^{2}). Geometric multiplicity: A−3I=(0100)A-3I=\begin{pmatrix}0&1\\0&0\end{pmatrix}, and the system (A−3I)v⃗=0⃗(A-3I)\vec{v}=\vec{0} gives the single equation y=0y=0, with xx free. So E3=span⁡{(10)}E_{3}=\operatorname{span}\left\{\begin{pmatrix}1\\0\end{pmatrix}\right\} and dim⁡E3=1\dim E_{3}=1. Comparison: geometric =1<2==1<2= algebraic. The criterion fails, AA is NOT diagonalizable ✓. Concretely, there is only ONE eigendirection in a space of dimension 2: it is impossible to build a basis of eigenvectors, and therefore impossible to form an invertible matrix PP.

c)

B=3IB=3I is also triangular, with eigenvalue 3 of algebraic multiplicity 2. But here B−3I=(0000)B-3I=\begin{pmatrix}0&0\\0&0\end{pmatrix}, the ZERO matrix: the system (B−3I)v⃗=0⃗(B-3I)\vec{v}=\vec{0} imposes NO condition, every vector is a solution, so E3=R2E_{3}=\mathbb{R}^{2} and dim⁡E3=2\dim E_{3}=2. Geometric =2==2= algebraic: the criterion is satisfied and BB is diagonalizable. It is even diagonal already, with P=IP=I. What changes between AA and BB is the off-diagonal entry: B=3IB=3I stretches every direction uniformly, so that EVERY non-zero vector is an eigenvector; AA, on the other hand, adds a SHEAR that leaves only one direction invariant. Two matrices with the same characteristic polynomial can therefore behave in opposite ways: the characteristic polynomial alone never decides diagonalizability, and that is the reason the geometric multiplicity exists.

d)

Proof. Let λ1,…,λn\lambda_{1},\ldots,\lambda_{n} be the nn distinct eigenvalues, and v⃗1,…,v⃗n\vec{v}_{1},\ldots,\vec{v}_{n} associated eigenvectors. The key point is that eigenvectors associated with DISTINCT eigenvalues are linearly independent. Here it is for two: suppose c1v⃗1+c2v⃗2=0⃗c_{1}\vec{v}_{1}+c_{2}\vec{v}_{2}=\vec{0}. Applying AA: c1λ1v⃗1+c2λ2v⃗2=0⃗c_{1}\lambda_{1}\vec{v}_{1}+c_{2}\lambda_{2}\vec{v}_{2}=\vec{0}. Multiplying the first relation by λ1\lambda_{1} and subtracting: c2(λ2−λ1)v⃗2=0⃗c_{2}(\lambda_{2}-\lambda_{1})\vec{v}_{2}=\vec{0}. Since λ2≠λ1\lambda_{2}\neq\lambda_{1} and v⃗2≠0⃗\vec{v}_{2}\neq\vec{0}, we get c2=0c_{2}=0, then c1=0c_{1}=0 ✓ (the general case is handled the same way by induction). We therefore have nn independent vectors in a space of dimension nn: they form a basis of it, made entirely of eigenvectors, and AA is diagonalizable ✓. Application: (570−2)\begin{pmatrix}5&7\\0&-2\end{pmatrix} is triangular, its eigenvalues are 5 and −2-2, which are DISTINCT, two of them in a space of dimension 2. It is therefore diagonalizable, without computing a single eigenvector. Be careful not to reverse the implication: the condition is SUFFICIENT, not necessary, as B=3IB=3I from part c) shows, diagonalizable with a double eigenvalue.

Exercise 7: Problem: powers of a matrix and two coupled populations

Diagonalization is not a formal exercise: it lets you compute AnA^{n} for large nn, and so predict the long-term behaviour of a coupled system. This is the application that justifies the whole chapter.

  • a) For A=(4213)A=\begin{pmatrix}4&2\\1&3\end{pmatrix}, diagonalized with P=(−1211)P=\begin{pmatrix}-1&2\\1&1\end{pmatrix} and D=(2005)D=\begin{pmatrix}2&0\\0&5\end{pmatrix}, prove the formula An=PDnP−1A^{n}=PD^{n}P^{-1}.
  • b) Compute A3A^{3} with this formula, then check by direct multiplication.
  • c) Two populations change each year according to {xk+1=4xk+2ykyk+1=xk+3yk\begin{cases}x_{k+1}=4x_{k}+2y_{k}\\y_{k+1}=x_{k}+3y_{k}\end{cases}, with x0=3x_{0}=3 and y0=1y_{0}=1. Express (xkyk)\begin{pmatrix}x_{k}\\y_{k}\end{pmatrix} in terms of kk.
  • d) What ratio xkyk\frac{x_{k}}{y_{k}} does the system approach as k→∞k\to\infty? Interpret this result in terms of eigenvalues.

Type your answers, the page tells you right or wrong 0/6

b)
c)
d)
Show the solution

Answers

  • a) An=PDnP−1A^n=PD^nP^{-1}
  • b) A3=(86783947)A^3=\begin{pmatrix}86&78\\39&47\end{pmatrix}
  • c) xk=2k+8⋅5k3x_k=\frac{2^k+8\cdot 5^k}{3}, yk=−2k+4⋅5k3y_k=\frac{-2^k+4\cdot 5^k}{3}
  • d) Limiting ratio 2
a)

Proof by induction, or more directly by telescoping. For n=2n=2: A2=(PDP−1)(PDP−1)=PD(P−1P)DP−1=PDIDP−1=PD2P−1A^{2}=\left(PDP^{-1}\right)\left(PDP^{-1}\right)=PD\left(P^{-1}P\right)DP^{-1}=PDIDP^{-1}=PD^{2}P^{-1} ✓, using P−1P=IP^{-1}P=I. The same mechanism repeats: in An=(PDP−1)nA^{n}=\left(PDP^{-1}\right)^{n}, all the inner factors P−1PP^{-1}P cancel in pairs, and only PDnP−1PD^{n}P^{-1} remains. The gain is decisive: DnD^{n} costs nothing, since a diagonal matrix raised to the power nn is just (λ1n00λ2n)\begin{pmatrix}\lambda_{1}^{n}&0\\0&\lambda_{2}^{n}\end{pmatrix}. You therefore replace nn matrix multiplications by two multiplications and two powers of numbers.

b)

D3=(800125)D^{3}=\begin{pmatrix}8&0\\0&125\end{pmatrix}. Then PD3=(−1211)(800125)=(−82508125)PD^{3}=\begin{pmatrix}-1&2\\1&1\end{pmatrix}\begin{pmatrix}8&0\\0&125\end{pmatrix}=\begin{pmatrix}-8&250\\8&125\end{pmatrix}, and A3=PD3P−1=(−82508125)(−13231313)=(83+2503−163+2503−83+1253163+1253)=(86783947)A^{3}=PD^{3}P^{-1}=\begin{pmatrix}-8&250\\8&125\end{pmatrix}\begin{pmatrix}-\frac{1}{3}&\frac{2}{3}\\\frac{1}{3}&\frac{1}{3}\end{pmatrix}=\begin{pmatrix}\frac{8}{3}+\frac{250}{3}&-\frac{16}{3}+\frac{250}{3}\\-\frac{8}{3}+\frac{125}{3}&\frac{16}{3}+\frac{125}{3}\end{pmatrix}=\begin{pmatrix}86&78\\39&47\end{pmatrix}. Direct check: A2=(4213)2=(1814711)A^{2}=\begin{pmatrix}4&2\\1&3\end{pmatrix}^{2}=\begin{pmatrix}18&14\\7&11\end{pmatrix}, then A3=A2A=(1814711)(4213)=(72+1436+4228+1114+33)=(86783947)A^{3}=A^{2}A=\begin{pmatrix}18&14\\7&11\end{pmatrix}\begin{pmatrix}4&2\\1&3\end{pmatrix}=\begin{pmatrix}72+14&36+42\\28+11&14+33\end{pmatrix}=\begin{pmatrix}86&78\\39&47\end{pmatrix} ✓.

c)

The system reads u⃗k+1=Au⃗k\vec{u}_{k+1}=A\vec{u}_{k} with u⃗k=(xkyk)\vec{u}_{k}=\begin{pmatrix}x_{k}\\y_{k}\end{pmatrix}, so u⃗k=Aku⃗0\vec{u}_{k}=A^{k}\vec{u}_{0}. The efficient method is to decompose u⃗0\vec{u}_{0} in the eigenbasis: look for aa and bb such that (31)=a(−11)+b(21)\begin{pmatrix}3\\1\end{pmatrix}=a\begin{pmatrix}-1\\1\end{pmatrix}+b\begin{pmatrix}2\\1\end{pmatrix}. This gives −a+2b=3-a+2b=3 and a+b=1a+b=1; adding, 3b=43b=4, so b=43b=\frac{4}{3} and a=1−43=−13a=1-\frac{4}{3}=-\frac{1}{3}. Since AkA^{k} acts on each eigenvector by multiplying it by λk\lambda^{k}: u⃗k=−13⋅2k(−11)+43⋅5k(21)\vec{u}_{k}=-\frac{1}{3}\cdot 2^{k}\begin{pmatrix}-1\\1\end{pmatrix}+\frac{4}{3}\cdot 5^{k}\begin{pmatrix}2\\1\end{pmatrix}, that is xk=2k3+8⋅5k3x_{k}=\frac{2^{k}}{3}+\frac{8\cdot 5^{k}}{3} and yk=−2k3+4⋅5k3y_{k}=-\frac{2^{k}}{3}+\frac{4\cdot 5^{k}}{3}. Check at k=0k=0: x0=13+83=3x_{0}=\frac{1}{3}+\frac{8}{3}=3 ✓ and y0=−13+43=1y_{0}=-\frac{1}{3}+\frac{4}{3}=1 ✓. Check at k=1k=1: x1=23+403=14x_{1}=\frac{2}{3}+\frac{40}{3}=14 and directly 4(3)+2(1)=144(3)+2(1)=14 ✓.

d)

Compute xkyk=2k3+8⋅5k3−2k3+4⋅5k3=2k+8⋅5k−2k+4⋅5k\frac{x_{k}}{y_{k}}=\frac{\frac{2^{k}}{3}+\frac{8\cdot 5^{k}}{3}}{-\frac{2^{k}}{3}+\frac{4\cdot 5^{k}}{3}}=\frac{2^{k}+8\cdot 5^{k}}{-2^{k}+4\cdot 5^{k}}. Dividing top and bottom by 5k5^{k}: (25)k+8−(25)k+4\frac{\left(\frac{2}{5}\right)^{k}+8}{-\left(\frac{2}{5}\right)^{k}+4}. Since 25<1\frac{2}{5}<1, the term (25)k\left(\frac{2}{5}\right)^{k} tends to 0, and the limit is 84=2\frac{8}{4}=2. The ratio therefore tends to 2, exactly the ratio of the components of the eigenvector (21)\begin{pmatrix}2\\1\end{pmatrix}. Interpretation: the DOMINANT eigenvalue is the one with the largest absolute value, here λ=5\lambda=5, and in the long run its contribution swamps that of λ=2\lambda=2, since 5k5^{k} grows infinitely faster than 2k2^{k}. The system therefore lines up with the dominant eigendirection, whatever the initial condition (provided its component bb along that direction is non-zero). This is the principle behind Markov chains, structured population models and Google's PageRank algorithm: the dominant eigenvalue sets the GROWTH RATE (here, in the long run, the population is multiplied by 5 each period) and its eigenvector sets the equilibrium DISTRIBUTION.

xyX₀X₃ ≈ direction (2, 1)E₅

Exercise 8: Symmetric matrices: orthogonal eigenvectors and orthogonal diagonalization

Real symmetric matrices hold a special place: they are always diagonalizable, and by an ORTHOGONAL change-of-basis matrix. We illustrate it, then prove it.

  • a) For A=(3113)A=\begin{pmatrix} 3 & 1 \\ 1 & 3 \end{pmatrix}, find the eigenvalues and an eigenvector for each, and check that they are orthogonal.
  • b) Normalize the two eigenvectors and form the matrix PP whose columns they are. Check that PP is orthogonal, that is PTP=IP^{T}P=I.
  • c) Check the orthogonal diagonalization A=PDPTA=PDP^{T}, using P−1=PTP^{-1}=P^{T}.
  • d) Prove, for any symmetric matrix, that two eigenvectors associated with distinct eigenvalues are orthogonal.

Type your answers, the page tells you right or wrong 0/3

a)
b)
Show the solution

Answers

  • a) λ=4\lambda=4: (1,1)(1, 1); λ=2\lambda=2: (1,−1)(1, -1)
  • b) P=12(111−1)P=\frac{1}{\sqrt2}\begin{pmatrix}1&1\\1&-1\end{pmatrix}, PTP=IP^TP=I
  • c) PDPT=APDP^T=A
  • d) (λ−μ) u⃗⋅v⃗=0(\lambda-\mu)\,\vec u\cdot\vec v=0
a)

The characteristic polynomial is λ2−6λ+8=(λ−4)(λ−2)\lambda^{2}-6\lambda+8=(\lambda-4)(\lambda-2), so the eigenvalues are λ1=4\lambda_{1}=4 and λ2=2\lambda_{2}=2. For λ=4\lambda=4: A−4I=(−111−1)A-4I=\begin{pmatrix} -1 & 1 \\ 1 & -1 \end{pmatrix} gives x=yx=y, hence v⃗4=(11)\vec{v}_{4}=\begin{pmatrix} 1 \\ 1 \end{pmatrix}. For λ=2\lambda=2: A−2I=(1111)A-2I=\begin{pmatrix} 1 & 1 \\ 1 & 1 \end{pmatrix} gives x=−yx=-y, hence v⃗2=(1−1)\vec{v}_{2}=\begin{pmatrix} 1 \\ -1 \end{pmatrix}. Dot product: v⃗4⋅v⃗2=(1)(1)+(1)(−1)=0\vec{v}_{4}\cdot\vec{v}_{2}=(1)(1)+(1)(-1)=0, the eigenvectors are orthogonal.

b)

Each vector has norm 12+12=2\sqrt{1^{2}+1^{2}}=\sqrt{2}, so the normalized vectors are 12(11)\dfrac{1}{\sqrt{2}}\begin{pmatrix} 1 \\ 1 \end{pmatrix} and 12(1−1)\dfrac{1}{\sqrt{2}}\begin{pmatrix} 1 \\ -1 \end{pmatrix}. Set P=12(111−1)P=\dfrac{1}{\sqrt{2}}\begin{pmatrix} 1 & 1 \\ 1 & -1 \end{pmatrix} and D=(4002)D=\begin{pmatrix} 4 & 0 \\ 0 & 2 \end{pmatrix}. Orthogonality check: PTP=12(111−1)(111−1)=12(2002)=IP^{T}P=\dfrac{1}{2}\begin{pmatrix} 1 & 1 \\ 1 & -1 \end{pmatrix}\begin{pmatrix} 1 & 1 \\ 1 & -1 \end{pmatrix}=\dfrac{1}{2}\begin{pmatrix} 2 & 0 \\ 0 & 2 \end{pmatrix}=I. The columns of PP therefore form an ORTHONORMAL basis, which is exactly what makes a matrix orthogonal.

c)

Since PP is orthogonal, P−1=PTP^{-1}=P^{T}, and the diagonalization A=PDP−1A=PDP^{-1} reads A=PDPTA=PDP^{T}. Check: PD=12(111−1)(4002)=12(424−2)PD=\dfrac{1}{\sqrt{2}}\begin{pmatrix} 1 & 1 \\ 1 & -1 \end{pmatrix}\begin{pmatrix} 4 & 0 \\ 0 & 2 \end{pmatrix}=\dfrac{1}{\sqrt{2}}\begin{pmatrix} 4 & 2 \\ 4 & -2 \end{pmatrix}, then PDPT=12(424−2)(111−1)=12(6226)=(3113)=APDP^{T}=\dfrac{1}{2}\begin{pmatrix} 4 & 2 \\ 4 & -2 \end{pmatrix}\begin{pmatrix} 1 & 1 \\ 1 & -1 \end{pmatrix}=\dfrac{1}{2}\begin{pmatrix} 6 & 2 \\ 2 & 6 \end{pmatrix}=\begin{pmatrix} 3 & 1 \\ 1 & 3 \end{pmatrix}=A ✓. The advantage of orthogonal diagonalization is that it avoids computing an inverse: the transpose, which is immediate, takes its place.

d)

Let AA be symmetric (AT=AA^{T}=A), and u⃗\vec{u}, v⃗\vec{v} eigenvectors for distinct eigenvalues λ≠μ\lambda\neq\mu: Au⃗=λu⃗A\vec{u}=\lambda\vec{u} and Av⃗=μv⃗A\vec{v}=\mu\vec{v}. Compute the dot product (Au⃗)⋅v⃗(A\vec{u})\cdot\vec{v} in two ways. On the one hand, (Au⃗)⋅v⃗=λu⃗⋅v⃗(A\vec{u})\cdot\vec{v}=\lambda\vec{u}\cdot\vec{v}. On the other hand, using symmetry, (Au⃗)⋅v⃗=(Au⃗)Tv⃗=u⃗TATv⃗=u⃗TAv⃗=u⃗⋅(Av⃗)=μ u⃗⋅v⃗(A\vec{u})\cdot\vec{v}=(A\vec{u})^{T}\vec{v}=\vec{u}^{T}A^{T}\vec{v}=\vec{u}^{T}A\vec{v}=\vec{u}\cdot(A\vec{v})=\mu\,\vec{u}\cdot\vec{v}. Equating: λ u⃗⋅v⃗=μ u⃗⋅v⃗\lambda\,\vec{u}\cdot\vec{v}=\mu\,\vec{u}\cdot\vec{v}, that is (λ−μ) u⃗⋅v⃗=0(\lambda-\mu)\,\vec{u}\cdot\vec{v}=0. Since λ≠μ\lambda\neq\mu, we conclude u⃗⋅v⃗=0\vec{u}\cdot\vec{v}=0: the two eigenvectors are orthogonal. This result, together with the fact that a real symmetric matrix always has all its eigenvalues real, is what guarantees orthogonal diagonalization: the spectral theorem.

-3-2-112345-3-2-1123454 (1, 1)2 (1, −1)

Exercise 9: The Cayley-Hamilton theorem: inverse and powers without diagonalizing

Every matrix satisfies its own characteristic polynomial. For a 2×22\times 2 matrix, this gives a relation A2=tr⁡(A) A−det⁡(A) IA^{2}=\operatorname{tr}(A)\,A-\det(A)\,I that lets you compute the inverse and powers without going through the eigenvectors.

  • a) State the Cayley-Hamilton theorem for a 2×22\times 2 matrix, then check it on A=(4213)A=\begin{pmatrix} 4 & 2 \\ 1 & 3 \end{pmatrix}.
  • b) Deduce an expression of A−1A^{-1} in terms of AA and II, and compute it.
  • c) Deduce that A3A^{3} is a linear combination of AA and II, and give its value.
  • d) Explain why every power AnA^{n} and the inverse A−1A^{-1} can be written as a combination of AA and II, and compare this method with diagonalization.

Type your answers, the page tells you right or wrong 0/5

b)
c)
d)
Show the solution

Answers

  • a) A2−7A+10I=0A^2-7A+10I=0
  • b) A−1=110(3−2−14)A^{-1}=\frac{1}{10}\begin{pmatrix}3&-2\\-1&4\end{pmatrix}
  • c) A3=39A−70IA^3=39A-70I
  • d) A recursion that always works against a closed formula
a)

Cayley-Hamilton theorem: a matrix satisfies its characteristic polynomial. For a 2×22\times 2 matrix, that polynomial is λ2−tr⁡(A)λ+det⁡(A)\lambda^{2}-\operatorname{tr}(A)\lambda+\det(A), so A2−tr⁡(A) A+det⁡(A) I=0A^{2}-\operatorname{tr}(A)\,A+\det(A)\,I=0. Here tr⁡(A)=7\operatorname{tr}(A)=7 and det⁡(A)=10\det(A)=10, and the expected relation is A2−7A+10I=0A^{2}-7A+10I=0. Check: A2=(1814711)A^{2}=\begin{pmatrix} 18 & 14 \\ 7 & 11 \end{pmatrix}, 7A=(2814721)7A=\begin{pmatrix} 28 & 14 \\ 7 & 21 \end{pmatrix}, so A2−7A+10I=(18−28+1014−147−711−21+10)=(0000)A^{2}-7A+10I=\begin{pmatrix} 18-28+10 & 14-14 \\ 7-7 & 11-21+10 \end{pmatrix}=\begin{pmatrix} 0 & 0 \\ 0 & 0 \end{pmatrix} ✓.

b)

Isolate the constant term in A2−7A+10I=0A^{2}-7A+10I=0: 10I=7A−A2=A(7I−A)10I=7A-A^{2}=A(7I-A), so A⋅110(7I−A)=IA\cdot\dfrac{1}{10}(7I-A)=I, which identifies the inverse: A−1=110(7I−A)=110((7007)−(4213))=110(3−2−14)A^{-1}=\dfrac{1}{10}(7I-A)=\dfrac{1}{10}\left(\begin{pmatrix} 7 & 0 \\ 0 & 7 \end{pmatrix}-\begin{pmatrix} 4 & 2 \\ 1 & 3 \end{pmatrix}\right)=\dfrac{1}{10}\begin{pmatrix} 3 & -2 \\ -1 & 4 \end{pmatrix}. Check: A A−1=110(4213)(3−2−14)=110(100010)=IA\,A^{-1}=\dfrac{1}{10}\begin{pmatrix} 4 & 2 \\ 1 & 3 \end{pmatrix}\begin{pmatrix} 3 & -2 \\ -1 & 4 \end{pmatrix}=\dfrac{1}{10}\begin{pmatrix} 10 & 0 \\ 0 & 10 \end{pmatrix}=I ✓. This is the usual 2×22\times 2 inverse formula, but derived from the theorem.

c)

From A2=7A−10IA^{2}=7A-10I, multiply by AA: A3=7A2−10A=7(7A−10I)−10A=49A−70I−10A=39A−70IA^{3}=7A^{2}-10A=7(7A-10I)-10A=49A-70I-10A=39A-70I. Numerically: A3=39(4213)−70I=(156−707839117−70)=(86783947)A^{3}=39\begin{pmatrix} 4 & 2 \\ 1 & 3 \end{pmatrix}-70I=\begin{pmatrix} 156-70 & 78 \\ 39 & 117-70 \end{pmatrix}=\begin{pmatrix} 86 & 78 \\ 39 & 47 \end{pmatrix}. This is exactly the A3A^{3} obtained with PD3P−1PD^{3}P^{-1} in the previous problem: two methods, one result.

d)

The relation A2=7A−10IA^{2}=7A-10I lets you replace any A2A^{2} by a combination of AA and II. By induction, each power AnA^{n} reduces step by step to such a combination αnA+βnI\alpha_{n}A+\beta_{n}I: as soon as an A2A^{2} appears, you reduce it. So does the inverse, as seen in b). Comparison with diagonalization: the method An=PDnP−1A^{n}=PD^{n}P^{-1} gives a CLOSED formula (depending explicitly on 2n2^{n} and 5n5^{n}) and works even for symbolic exponents, but it requires diagonalizing, so finding the eigenvectors and the inverse of PP. Cayley-Hamilton only needs the trace and the determinant and applies even to NON-diagonalizable matrices, but it computes the powers step by step rather than by a closed formula. The two approaches complement each other: an explicit formula on one side, a recursion that always works on the other.

Exercise 10: Proofs: eigenvalues of the inverse, of powers and of shifted matrices

Knowing the eigenvalues of AA gives those of many derived matrices for free, with the SAME eigenvectors. These results are proved in one line and save repeated calculations.

  • a) Prove that if λ\lambda is an eigenvalue of AA with eigenvector v⃗\vec{v}, then λk\lambda^{k} is an eigenvalue of AkA^{k}, with the same eigenvector.
  • b) Prove that if AA is invertible and λ\lambda is an eigenvalue, then λ≠0\lambda\neq 0 and 1λ\dfrac{1}{\lambda} is an eigenvalue of A−1A^{-1}, with the same eigenvector.
  • c) Prove that λ+c\lambda+c is an eigenvalue of A+cIA+cI, with the same eigenvector, for every scalar cc.
  • d) Application to A=(4213)A=\begin{pmatrix} 4 & 2 \\ 1 & 3 \end{pmatrix}, with eigenvalues 22 and 55: give, without any new calculation, the eigenvalues of A3A^{3}, of A−1A^{-1} and of A+2IA+2I, and check the trace relation on A+2IA+2I.

Type your answers, the page tells you right or wrong 0/4

d)
Show the solution

Answers

  • a) Akv⃗=λkv⃗A^k\vec v=\lambda^k\vec v
  • b) A−1v⃗=1λv⃗A^{-1}\vec v=\frac1\lambda\vec v
  • c) (A+cI)v⃗=(λ+c)v⃗(A+cI)\vec v=(\lambda+c)\vec v
  • d) 88, 125125; 12\frac12, 15\frac15; 44, 77
a)

Start from Av⃗=λv⃗A\vec{v}=\lambda\vec{v} and apply AA repeatedly. A2v⃗=A(Av⃗)=A(λv⃗)=λAv⃗=λ2v⃗A^{2}\vec{v}=A(A\vec{v})=A(\lambda\vec{v})=\lambda A\vec{v}=\lambda^{2}\vec{v}. Iterating, or by induction, Akv⃗=λkv⃗A^{k}\vec{v}=\lambda^{k}\vec{v}. So λk\lambda^{k} is an eigenvalue of AkA^{k}, with the same eigenvector v⃗\vec{v}: raising the matrix to a power raises the eigenvalues to that power, without touching the eigendirections.

b)

If AA is invertible, no eigenvalue is zero: λ=0\lambda=0 would mean Av⃗=0⃗A\vec{v}=\vec{0} for some v⃗≠0⃗\vec{v}\neq\vec{0}, hence a non-trivial null space, which contradicts invertibility. So 1λ\dfrac{1}{\lambda} makes sense. Applying A−1A^{-1} to Av⃗=λv⃗A\vec{v}=\lambda\vec{v}: v⃗=A−1(λv⃗)=λA−1v⃗\vec{v}=A^{-1}(\lambda\vec{v})=\lambda A^{-1}\vec{v}, hence A−1v⃗=1λv⃗A^{-1}\vec{v}=\dfrac{1}{\lambda}\vec{v}. So 1λ\dfrac{1}{\lambda} is an eigenvalue of A−1A^{-1}, with the same eigenvector.

c)

(A+cI)v⃗=Av⃗+cIv⃗=λv⃗+cv⃗=(λ+c)v⃗(A+cI)\vec{v}=A\vec{v}+cI\vec{v}=\lambda\vec{v}+c\vec{v}=(\lambda+c)\vec{v}. So λ+c\lambda+c is an eigenvalue of A+cIA+cI, with the same eigenvector. Adding cIcI shifts every eigenvalue by cc without changing the eigendirections, since cIcI acts the same way on every vector.

d)

The eigenvalues of AA are 22 and 55. By a), those of A3A^{3} are 23=82^{3}=8 and 53=1255^{3}=125. By b), those of A−1A^{-1} are 12\dfrac{1}{2} and 15\dfrac{1}{5}. By c), those of A+2IA+2I are 2+2=42+2=4 and 5+2=75+2=7. Trace check on A+2I=(6215)A+2I=\begin{pmatrix} 6 & 2 \\ 1 & 5 \end{pmatrix}: its trace is 6+5=116+5=11, and the sum of the eigenvalues 4+7=114+7=11 ✓. All these results come without redoing a single characteristic polynomial: that is the practical value of these three properties, which come up constantly in problems on sequences and dynamical systems. Finally, note that AA and ATA^{T} have exactly the same eigenvalues, since det⁡(AT−λI)=det⁡((A−λI)T)=det⁡(A−λI)\det(A^{T}-\lambda I)=\det((A-\lambda I)^{T})=\det(A-\lambda I), even though their eigenvectors generally differ.

Part C: the classics (/50)

Exercise 11: A 3×3 matrix with a double eigenvalue

Consider A=(211121112)A=\begin{pmatrix} 2 & 1 & 1 \\ 1 & 2 & 1 \\ 1 & 1 & 2 \end{pmatrix}, which can be written A=I+JA=I+J, where JJ is the matrix whose entries are all 11.

  • a) Show that (1,1,1)(1, 1, 1) is an eigenvector of AA and give its eigenvalue.
  • b) Show that Av⃗=v⃗+(x+y+z)(1,1,1)A\vec{v}=\vec{v}+(x+y+z)(1, 1, 1) for v⃗=(x,y,z)\vec{v}=(x, y, z). Deduce a second eigenvalue and the dimension of its eigenspace.
  • c) Check with the trace and the determinant of AA.
  • d) Write PP and DD such that A=PDP−1A=PDP^{-1}, and justify that PP is invertible.
  • e) Given that J2=3JJ^{2}=3J, show that An=I+4n−13JA^{n}=I+\frac{4^{n}-1}{3}J and give the entries of A3A^{3}.

Type your answers, the page tells you right or wrong 0/8

a)
b)
c)
d)
e)
Show the solution

Answers

  • a) λ=4\lambda=4
  • b) λ=1\lambda=1 on the plane x+y+z=0x+y+z=0
  • c) Trace 6, determinant 4
  • d) det⁡P=3\det P=3: diagonalizable
  • e) A3A^3: 22 and 21
a)

A(1,1,1)=(4,4,4)=4(1,1,1)A(1, 1, 1)=(4, 4, 4)=4(1, 1, 1): eigenvalue 44.

b)

Jv⃗=(x+y+z)(1,1,1)J\vec{v}=(x+y+z)(1, 1, 1), so Av⃗=v⃗+Jv⃗A\vec{v}=\vec{v}+J\vec{v} gives the formula. If x+y+z=0x+y+z=0, then Av⃗=v⃗A\vec{v}=\vec{v}: every non-zero vector of the plane x+y+z=0x+y+z=0 is an eigenvector for the eigenvalue 11. This plane has dimension 22, with basis (1,−1,0)(1, -1, 0) and (1,0,−1)(1, 0, -1).

c)

The sum of the eigenvalues, counted with multiplicity, is 4+1+1=6=tr⁡A4+1+1=6=\operatorname{tr}A. Their product is 4×1×1=44\times 1\times 1=4, and det⁡A=2(4−1)−1(2−1)+1(1−2)=6−1−1=4\det A=2(4-1)-1(2-1)+1(1-2)=6-1-1=4. Both checks agree.

d)

P=(1111−1010−1)P=\begin{pmatrix} 1 & 1 & 1 \\ 1 & -1 & 0 \\ 1 & 0 & -1 \end{pmatrix}, whose columns are (1,1,1)(1, 1, 1), (1,−1,0)(1, -1, 0) and (1,0,−1)(1, 0, -1), and D=diag⁡(4,1,1)D=\operatorname{diag}(4, 1, 1). det⁡P=1(1−0)−1(−1−0)+1(0+1)=3≠0\det P=1(1-0)-1(-1-0)+1(0+1)=3\neq 0: the three eigenvectors are independent. The double eigenvalue does have an eigenspace of dimension 22: AA is diagonalizable.

e)

By induction: true for n=1n=1, since 4−13=1\frac{4-1}{3}=1. If An=I+cnJA^{n}=I+c_{n}J, then An+1=(I+cnJ)(I+J)=I+(1+cn+3cn)JA^{n+1}=(I+c_{n}J)(I+J)=I+(1+c_{n}+3c_{n})J, and 1+4cn=1+4n+1−43=4n+1−131+4c_{n}=1+\frac{4^{n+1}-4}{3}=\frac{4^{n+1}-1}{3}. For n=3n=3: 633=21\frac{63}{3}=21, so A3A^{3} has 2222 on the diagonal and 2121 everywhere else.

The thread: a double eigenvalue blocks nothing if its eigenspace has dimension 2; writing A=I+JA=I+J makes the eigenvalues almost obvious.

Exercise 12: Eigenvalues read off the geometry

Consider the orthogonal projection P=12(1111)P=\frac{1}{2}\begin{pmatrix} 1 & 1 \\ 1 & 1 \end{pmatrix} onto the line y=xy=x, the reflection S=(0110)S=\begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix} across the line y=xy=x, the dilation H=3IH=3I and the shear C=(1201)C=\begin{pmatrix} 1 & 2 \\ 0 & 1 \end{pmatrix}.

  • a) Without the characteristic polynomial, give the eigenvalues of PP and associated eigenvectors. Check by calculation.
  • b) Same question for SS.
  • c) Same question for HH.
  • d) Find the eigenvalues of CC and the dimension of the eigenspace. Is CC diagonalizable?
  • e) Express PP in terms of SS and II, and use this to recover the eigenvalues of PP from those of SS.

Type your answers, the page tells you right or wrong 0/8

a)
b)
c)
d)
e)
Show the solution

Answers

  • a) 1 and 0
  • b) 1 and −1-1
  • c) 3, double, every vector
  • d) 1, double, dimension 1: no
  • e) P=S+I2P=\frac{S+I}{2}
a)

A vector on the line y=xy=x is not changed by the projection: (1,1)(1, 1) is an eigenvector for 11. A vector perpendicular to the line is crushed onto the origin: (1,−1)(1, -1) is an eigenvector for 00. Calculation: P(1,1)=(1,1)P(1, 1)=(1, 1) and P(1,−1)=(0,0)P(1, -1)=(0, 0).

b)

The reflection leaves the axis fixed, so (1,1)(1, 1) is an eigenvector for 11, and reverses the perpendicular direction, so (1,−1)(1, -1) is an eigenvector for −1-1. Calculation: S(1,−1)=(−1,1)=−(1,−1)S(1, -1)=(-1, 1)=-(1, -1).

c)

HH multiplies every vector by 33: 33 is the only eigenvalue, a double one, and EVERY non-zero vector is an eigenvector. The eigenspace is the whole plane.

d)

CC is triangular: eigenvalue 11, double. C−I=(0200)C-I=\begin{pmatrix} 0 & 2 \\ 0 & 0 \end{pmatrix} forces y=0y=0: the eigenspace is the xx-axis, of dimension 11. Geometrically, only the horizontal direction is not tilted by the shear. Dimension 1<21<2: CC is not diagonalizable.

e)

S+I2=12(1111)=P\frac{S+I}{2}=\frac{1}{2}\begin{pmatrix} 1 & 1 \\ 1 & 1 \end{pmatrix}=P. The eigenvalues of SS, 11 and −1-1, become those of PP through λ↦λ+12\lambda\mapsto\frac{\lambda+1}{2}: 11 and 00, with the same eigenvectors, as exercise 10 predicts.

The thread: a direction that the transformation only stretches is an eigendirection; the geometry often gives the eigenvalues before any calculation.

Exercise 13: Diagonalizable or not, depending on a parameter

For every real number aa, consider Aa=(1a02)A_{a}=\begin{pmatrix} 1 & a \\ 0 & 2 \end{pmatrix}, Ba=(2a02)B_{a}=\begin{pmatrix} 2 & a \\ 0 & 2 \end{pmatrix} and Ca=(11a1)C_{a}=\begin{pmatrix} 1 & 1 \\ a & 1 \end{pmatrix}.

  • a) Give the eigenvalues of AaA_{a}. For which values of aa is the matrix diagonalizable?
  • b) Same question for BaB_{a}.
  • c) Compute the eigenvalues of CaC_{a} in terms of aa, for a≥0a\geq 0.
  • d) Discuss the diagonalizability of CaC_{a} over R\mathbb{R} according to the sign of aa.
  • e) For a=4a=4, find eigenvectors of C4C_{4}.

Type your answers, the page tells you right or wrong 0/7

a)
b)
c)
d)
e)
Show the solution

Answers

  • a) Always diagonalizable
  • b) Only for a=0a=0
  • c) 1±a1\pm\sqrt a
  • d) Yes if a>0a>0, no otherwise
  • e) (1,2)(1, 2) and (1,−2)(1, -2)
a)

AaA_{a} is triangular: eigenvalues 11 and 22, distinct whatever aa is. Two distinct eigenvalues in dimension 2 are enough: AaA_{a} is diagonalizable for every aa. An eigenvector for 22 is (a,1)(a, 1).

b)

BaB_{a} has the single eigenvalue 22, double. Ba−2I=(0a00)B_{a}-2I=\begin{pmatrix} 0 & a \\ 0 & 0 \end{pmatrix}: if a≠0a\neq 0, the eigenspace is the line y=0y=0, of dimension 11, and BaB_{a} is not diagonalizable; if a=0a=0, B0=2IB_{0}=2I is already diagonal. Diagonalizable if and only if a=0a=0.

c)

det⁡(Ca−λI)=(1−λ)2−a=0\det(C_{a}-\lambda I)=(1-\lambda)^{2}-a=0 gives (1−λ)2=a(1-\lambda)^{2}=a, so λ=1±a\lambda=1\pm\sqrt{a}.

d)

If a>0a>0: two distinct real eigenvalues, CaC_{a} is diagonalizable. If a=0a=0: eigenvalue 11, double, and C0−I=(0100)C_{0}-I=\begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix} has a null space of dimension 11, so C0C_{0} is not diagonalizable. If a<0a<0: complex eigenvalues 1±i−a1\pm i\sqrt{-a}, no real eigendirection, CaC_{a} is not diagonalizable over R\mathbb{R}.

e)

Eigenvalues 33 and −1-1. C4−3I=(−214−2)C_{4}-3I=\begin{pmatrix} -2 & 1 \\ 4 & -2 \end{pmatrix} gives y=2xy=2x: eigenvector (1,2)(1, 2). C4+I=(2142)C_{4}+I=\begin{pmatrix} 2 & 1 \\ 4 & 2 \end{pmatrix} gives y=−2xy=-2x: eigenvector (1,−2)(1, -2). Check: C4(1,2)=(3,6)C_{4}(1, 2)=(3, 6).

The thread: distinct eigenvalues guarantee diagonalization; a repeated eigenvalue requires checking the dimension of the eigenspace.

Exercise 14: Problem: the normal modes of two carts on springs

Two carts of the same mass mm slide without friction between two walls. Three identical springs of stiffness kk connect the left wall to the first cart, the two carts to each other, and the second cart to the right wall. Let x1x_{1} and x2x_{2} be the displacements of the carts from equilibrium, and take km=4 s−2\frac{k}{m}=4\ \mathrm{s^{-2}}.

Newton's second law gives x1′′=4(−2x1+x2)x_{1}''=4(-2x_{1}+x_{2}) and x2′′=4(x1−2x2)x_{2}''=4(x_{1}-2x_{2}), that is X′′=−KXX''=-KX. You may use the fact that if Kv⃗=ω2v⃗K\vec{v}=\omega^{2}\vec{v}, then X(t)=v⃗cos⁡(ωt)X(t)=\vec{v}\cos(\omega t) is a solution: a NORMAL MODE, in which both carts oscillate at the same angular frequency ω\omega.

mmx₁x₂kkk
  • a) Write the matrix KK. What property does it have?
  • b) Compute the eigenvalues of KK.
  • c) Find an eigenvector for each, and describe the motion of the carts in each mode.
  • d) Compute the two angular frequencies ω\omega and the corresponding periods, to the nearest hundredth.
  • e) The first cart is moved 2 cm2\ \mathrm{cm} away from equilibrium, the second staying in place, and both are released from rest. Decompose X(0)=(2,0)X(0)=(2, 0) on the two modes and write x1(t)x_{1}(t) and x2(t)x_{2}(t).

Type your answers, the page tells you right or wrong 0/10

a)
b)
c)
d)
e)
Show the solution

Answers

  • a) K=(8−4−48)K=\begin{pmatrix}8&-4\\-4&8\end{pmatrix}, symmetric
  • b) 4 and 12
  • c) (1,1)(1, 1) in phase, (1,−1)(1, -1) out of phase
  • d) 22 and 3.463.46 rad/s; 3.143.14 s and 1.811.81 s
  • e) x1=cos⁡2t+cos⁡12tx_1=\cos 2t+\cos\sqrt{12}t, x2=cos⁡2t−cos⁡12tx_2=\cos 2t-\cos\sqrt{12}t
a)

K=(8−4−48)K=\begin{pmatrix} 8 & -4 \\ -4 & 8 \end{pmatrix}: it is SYMMETRIC, so diagonalizable with orthogonal eigenvectors.

b)

λ2−16λ+(64−16)=λ2−16λ+48=(λ−4)(λ−12)\lambda^{2}-16\lambda+(64-16)=\lambda^{2}-16\lambda+48=(\lambda-4)(\lambda-12): eigenvalues 44 and 1212.

c)

For 44: K−4I=(4−4−44)K-4I=\begin{pmatrix} 4 & -4 \\ -4 & 4 \end{pmatrix}, so x1=x2x_{1}=x_{2}, vector (1,1)(1, 1): the carts move together, IN PHASE, and the middle spring is never stretched. For 1212: K−12I=(−4−4−4−4)K-12I=\begin{pmatrix} -4 & -4 \\ -4 & -4 \end{pmatrix}, so x1=−x2x_{1}=-x_{2}, vector (1,−1)(1, -1): the carts move OUT OF PHASE, as mirror images.

d)

ω1=4=2 rad/s\omega_{1}=\sqrt{4}=2\ \mathrm{rad/s} and ω2=12≈3.46 rad/s\omega_{2}=\sqrt{12}\approx 3.46\ \mathrm{rad/s}. Periods T=2πωT=\frac{2\pi}{\omega}: π≈3.14 s\pi\approx 3.14\ \mathrm{s} and 2π12≈1.81 s\frac{2\pi}{\sqrt{12}}\approx 1.81\ \mathrm{s}. The out-of-phase mode is faster, because the middle spring works in it too.

e)

(2,0)=1⋅(1,1)+1⋅(1,−1)(2, 0)=1\cdot(1, 1)+1\cdot(1, -1). Each mode evolves at its own angular frequency: x1(t)=cos⁡(2t)+cos⁡(12 t)x_{1}(t)=\cos(2t)+\cos(\sqrt{12}\,t) and x2(t)=cos⁡(2t)−cos⁡(12 t)x_{2}(t)=\cos(2t)-\cos(\sqrt{12}\,t), in centimetres. Since the two angular frequencies are not in a simple whole-number ratio, the motion is not periodic: the energy passes from one cart to the other.

The thread: diagonalizing a coupled system splits it into independent modes; the eigenvalues give the frequencies, the eigenvectors the shapes of the motion.

Exercise 15: Problem: the principal axes of a tilted ellipse

We study the curve with equation 5x2+4xy+2y2=65x^{2}+4xy+2y^{2}=6. Let X=(xy)X=\begin{pmatrix} x \\ y \end{pmatrix} and A=(5222)A=\begin{pmatrix} 5 & 2 \\ 2 & 2 \end{pmatrix}.

-3-2-1123-3-2-1123
  • a) Check that XTAX=5x2+4xy+2y2X^{T}AX=5x^{2}+4xy+2y^{2}.
  • b) Compute the eigenvalues of AA.
  • c) Find an eigenvector for each and check that they are orthogonal.
  • d) In the orthonormal frame carried by the normalized eigenvectors, with coordinate x′x' along the vector of the largest eigenvalue and y′y' along the other, the equation becomes 6x′2+y′2=66x'^{2}+y'^{2}=6. Justify this, then give the semi-axes of the ellipse, to the nearest hundredth.
  • e) Give the angle between the major axis and the xx-axis, to the nearest hundredth of a degree, and a point of the curve at maximum distance from the origin.

Type your answers, the page tells you right or wrong 0/6

b)
c)
d)
e)
Show the solution

Answers

  • a) XTAX=5x2+4xy+2y2X^TAX=5x^2+4xy+2y^2
  • b) 6 and 1
  • c) (2,1)(2, 1) and (1,−2)(1, -2)
  • d) Semi-axes 1 and 6≈2.45\sqrt6\approx 2.45
  • e) −63.43∘-63.43^\circ
a)

AX=(5x+2y,2x+2y)AX=(5x+2y, 2x+2y) and XTAX=x(5x+2y)+y(2x+2y)=5x2+4xy+2y2X^{T}AX=x(5x+2y)+y(2x+2y)=5x^{2}+4xy+2y^{2}. The coefficient 44 of the cross term is split into two 22s off the diagonal, which makes AA symmetric.

b)

λ2−7λ+6=(λ−1)(λ−6)\lambda^{2}-7\lambda+6=(\lambda-1)(\lambda-6): eigenvalues 66 and 11.

c)

For 66: A−6I=(−122−4)A-6I=\begin{pmatrix} -1 & 2 \\ 2 & -4 \end{pmatrix}, so x=2yx=2y: (2,1)(2, 1). For 11: A−I=(4221)A-I=\begin{pmatrix} 4 & 2 \\ 2 & 1 \end{pmatrix}, so y=−2xy=-2x: (1,−2)(1, -2). Dot product 2−2=02-2=0: orthogonal, as for every symmetric matrix.

d)

With PP orthogonal, with columns 15(2,1)\frac{1}{\sqrt{5}}(2, 1) and 15(1,−2)\frac{1}{\sqrt{5}}(1, -2), set X=PX′X=PX': XTAX=X′T(PTAP)X′=X′TDX′=6x′2+y′2X^{T}AX=X'^{T}(P^{T}AP)X'=X'^{T}DX'=6x'^{2}+y'^{2}. The equation x′2+y′26=1x'^{2}+\frac{y'^{2}}{6}=1 is that of an ellipse with semi-axes 11, along (2,1)(2, 1), and 6≈2.45\sqrt{6}\approx 2.45, along (1,−2)(1, -2). The cross term has disappeared: that is the point of orthogonal diagonalization.

e)

The major axis is directed by (1,−2)(1, -2), at the angle arctan⁡(−2)≈−63.43∘\arctan(-2)\approx -63.43^{\circ}. The point 6⋅15(1,−2)\sqrt{6}\cdot\frac{1}{\sqrt{5}}(1, -2) has squared coordinates 65\frac{6}{5} and 245\frac{24}{5} and product of coordinates −125-\frac{12}{5}: 5⋅65+4⋅(−125)+2⋅245=6−485+485=65\cdot\frac{6}{5}+4\cdot\left(-\frac{12}{5}\right)+2\cdot\frac{24}{5}=6-\frac{48}{5}+\frac{48}{5}=6. It lies on the curve, at distance 6\sqrt{6} from the origin.

The thread: a quadratic form is read off its symmetric matrix; its eigenvectors give the axes, its eigenvalues the lengths.

-3-2-1123-3-2-11231√6
Something missing from this set? A type of exercise your teacher gives, a problem you are stuck on, a mistake you spotted: tell me, I read every request.

Sending goes through your free space: that is how I can answer you and tell you when the exercise is added. An email is enough, no password. Create your space. Rather not create an account? Write to me directly.

What is it about?
I read every request myself, nobody else sees it. No need to write your name or your school.
When your request is dealt with, you hear about it by email and in your space.

By sending, you agree that I read your request to answer it. It is kept with your account, erased with it, and never published: privacy policy.

Previous chapter Linear transformations of the plane, 201-SN4-RE Next chapter Matrix models: Markov, Leontief, Leslie, 201-SN4-RE

© Ahmed Squalli Houssaini. Exercise set published at www.letuteurscientifique.ca/en/exercices/cegep-sn4-eigenvalues-diagonalization. Free for personal and classroom use; republishing it elsewhere requires written permission (legal notice).

See also

Looking for a linear algebra tutor in Montreal?

Get in touch for a first session. We work on linear algebra at the real level of CEGEP evaluations, diagonalization included, the chapter that finally gives a geometric meaning to all the matrix computation of the semester.

Site by Studio Squalli