Practice / Linear algebra

Eigenvalues and eigenvectors by hand

Ten problems on eigenvalues and eigenvectors worked by hand: 2×2 and 3×3 matrices, the characteristic polynomial, a rotation with complex eigenvalues, a shear with too few eigenvectors, diagonalisation and powers, orthogonal eigenvectors of symmetric matrices, power iteration and singular values, with worked solutions and the classic sign and null-space mistakes.

Before you start

Every eigenvalue problem on this page comes down to two computations: a determinant that gives the eigenvalues, and a small linear system for each one that gives the eigenvectors. The errors come from the same few places: a sign in A−λIA - \lambda I, an eigenvector read off a row instead of solved for, a repeated eigenvalue assumed to bring two eigenvectors. These ten problems start with 2×22\times2 matrices, meet a rotation and a shear that break the easy pattern, and end with the three uses machine learning makes of eigenvectors: powers of a matrix, power iteration and singular values.

  • AA is a real n×nn\times n matrix and II the n×nn\times n identity. A number λ\lambda is an eigenvalue of AA, with eigenvector vv, when Av=λvAv = \lambda v and v≠0v \neq 0.
  • Av=λvAv = \lambda v is the same as (A−λI)v=0(A - \lambda I)v = 0. That has a nonzero solution exactly when A−λIA - \lambda I is singular, that is when det⁡(A−λI)=0\det(A - \lambda I) = 0. det⁡(A−λI)\det(A - \lambda I) is a polynomial of degree nn in λ\lambda, the characteristic polynomial, and its roots are the eigenvalues.
  • An eigenvector is any nonzero solution of (A−λI)v=0(A - \lambda I)v = 0, so the answers below give one representative each; any nonzero multiple is also correct.
  • Vectors are columns; (1,−1)⊤(1, -1)^\top is the column with entries 11 and −1-1, written on one line to save space.
  • For 2×22\times2 matrices: det⁡[abcd]=ad−bc\det\begin{bmatrix}a&b\\c&d\end{bmatrix} = ad - bc, and when that is not 00, [abcd]−1=1ad−bc[d−b−ca]\begin{bmatrix}a&b\\c&d\end{bmatrix}^{-1} = \dfrac{1}{ad-bc}\begin{bmatrix}d&-b\\-c&a\end{bmatrix}. tr⁡A\operatorname{tr}A, the trace, is the sum of the diagonal entries.
  • PP is a matrix whose columns are eigenvectors of AA. DD (Problem 7) and Λ\Lambda (Problem 8) are both diagonal matrices of eigenvalues, listed in the same order as the columns of PP; diag⁡(a,b)\operatorname{diag}(a, b) is the diagonal matrix with aa and bb on the diagonal. Λ\Lambda is used when the eigenvectors are orthonormal, and then the eigenvector matrix is written QQ: an orthogonal matrix, meaning Q⊤Q=IQ^\top Q = I, so Q−1=Q⊤Q^{-1} = Q^\top.
  • ii is the imaginary unit, i2=−1i^2 = -1. σ1,σ2\sigma_1, \sigma_2 are singular values (Problem 10).

Problems

  1. ·

    Find the eigenvalues and eigenvectors of A=[2112]A = \begin{bmatrix}2&1\\1&2\end{bmatrix}.

  2. ·

    Find the characteristic polynomial, eigenvalues and eigenvectors of A=[4123]A = \begin{bmatrix}4&1\\2&3\end{bmatrix}.

  3. ··

    Find the eigenvalues of the rotation R=[0−110]R = \begin{bmatrix}0&-1\\1&0\end{bmatrix}. Does it have a real eigenvector? What does that mean geometrically?

  4. ··

    Find the eigenvalues and eigenvectors of the shear A=[1101]A = \begin{bmatrix}1&1\\0&1\end{bmatrix}. Can it be diagonalised?

  5. ··

    For any 2×22\times2 matrix, show that the trace is the sum of the eigenvalues and the determinant is their product. Verify on Problem 2.

  6. ···

    Find the eigenvalues and eigenvectors of A=[210121012]A = \begin{bmatrix}2&1&0\\1&2&1\\0&1&2\end{bmatrix}.

  7. ··

    Diagonalise the matrix of Problem 2 as A=PDP−1A = PDP^{-1} and use it to compute A5A^5.

  8. ··

    Show that the eigenvectors of Problem 1 are orthogonal, and write A=QΛQ⊤A = Q\Lambda Q^\top with QQ orthogonal. Why does this always happen for symmetric matrices?

  9. ···

    Apply A=[2112]A = \begin{bmatrix}2&1\\1&2\end{bmatrix} repeatedly to x0=(1,0)⊤x_0 = (1, 0)^\top. Find Akx0A^kx_0 in closed form and say what direction it approaches and how fast.

  10. ···

    For A=[3045]A = \begin{bmatrix}3&0\\4&5\end{bmatrix}, find the eigenvalues of A⊤AA^\top A and hence the singular values of AA. Check them against det⁡A\det A.

Worked solutions

Problem 1

Find the eigenvalues and eigenvectors of A=[2112]A = \begin{bmatrix}2&1\\1&2\end{bmatrix}.

  1. A−λI=[2−λ112−λ]A - \lambda I = \begin{bmatrix}2-\lambda&1\\1&2-\lambda\end{bmatrix}, so det⁡(A−λI)=(2−λ)2−1\det(A - \lambda I) = (2-\lambda)^2 - 1.Subtract λ\lambda from each diagonal entry, then ad−bcad - bc.
  2. (2−λ)2=1(2 - \lambda)^2 = 1, so 2−λ=±12 - \lambda = \pm1 and λ=1{\lambda = 1} or λ=3{\lambda = 3}.The eigenvalues are the roots of the characteristic polynomial; this one is a square minus 11, so take square roots rather than expanding.
  3. λ=3\lambda = 3: A−3I=[−111−1]A - 3I = \begin{bmatrix}-1&1\\1&-1\end{bmatrix}, and with v=(x,y)⊤v = (x, y)^\top, (A−3I)v=0(A - 3I)v = 0 says −x+y=0-x + y = 0, so v=(1,1)⊤v = (1, 1)^\top.The eigenvectors for λ=3\lambda = 3 are the null space of A−3IA - 3I. Its two rows are multiples of each other, as they must be when the determinant is 00, so one equation remains; choose x=1x = 1.
  4. λ=1\lambda = 1: A−I=[1111]A - I = \begin{bmatrix}1&1\\1&1\end{bmatrix}, and (A−I)v=0(A - I)v = 0 says x+y=0x + y = 0, so v=(1,−1)⊤v = (1, -1)^\top.Null space of A−IA - I; the equation x+y=0x + y = 0 fixes vv up to scale.
  5. λ=3\lambda = 3 with v=(1,1)⊤v = (1,1)^\top; λ=1\lambda = 1 with v=(1,−1)⊤v = (1,-1)^\topCheck: A(1,1)⊤=(3,3)⊤A(1,1)^\top = (3,3)^\top and A(1,−1)⊤=(1,−1)⊤A(1,-1)^\top = (1,-1)^\top.

Problem 2

Find the characteristic polynomial, eigenvalues and eigenvectors of A=[4123]A = \begin{bmatrix}4&1\\2&3\end{bmatrix}.

  1. det⁡(A−λI)=(4−λ)(3−λ)−1⋅2\det(A - \lambda I) = (4-\lambda)(3-\lambda) - 1\cdot2.ad−bcad - bc for A−λI=[4−λ123−λ]A - \lambda I = \begin{bmatrix}4-\lambda&1\\2&3-\lambda\end{bmatrix}.
  2. (4−λ)(3−λ)−2=λ2−7λ+12−2=λ2−7λ+10(4-\lambda)(3-\lambda) - 2 = \lambda^2 - 7\lambda + 12 - 2 = \lambda^2 - 7\lambda + 10.Expand the product and collect terms.
  3. λ2−7λ+10=(λ−5)(λ−2)\lambda^2 - 7\lambda + 10 = (\lambda - 5)(\lambda - 2).Two numbers with product 1010 and sum 77.
  4. λ=5\lambda = 5: A−5I=[−112−2]A - 5I = \begin{bmatrix}-1&1\\2&-2\end{bmatrix} gives −x+y=0-x + y = 0 for v=(x,y)⊤v = (x, y)^\top, so v=(1,1)⊤v = (1, 1)^\top.Null space of A−5IA - 5I; the second row is −2-2 times the first, so it adds nothing.
  5. λ=2\lambda = 2: A−2I=[2121]A - 2I = \begin{bmatrix}2&1\\2&1\end{bmatrix} gives 2x+y=02x + y = 0, so v=(1,−2)⊤v = (1, -2)^\top.Null space of A−2IA - 2I, with x=1x = 1.
  6. λ2−7λ+10=(λ−5)(λ−2)\lambda^2 - 7\lambda + 10 = (\lambda-5)(\lambda-2); λ=5\lambda = 5, v=(1,1)⊤v = (1,1)^\top; λ=2\lambda = 2, v=(1,−2)⊤v = (1,-2)^\topCheck: A(1,1)⊤=(5,5)⊤A(1,1)^\top = (5,5)^\top and A(1,−2)⊤=(2,−4)⊤A(1,-2)^\top = (2,-4)^\top. The eigenvectors are not perpendicular; AA is not symmetric, and Problem 8 shows why that matters.

Problem 3

Find the eigenvalues of the rotation R=[0−110]R = \begin{bmatrix}0&-1\\1&0\end{bmatrix}. Does it have a real eigenvector? What does that mean geometrically?

  1. det⁡(R−λI)=det⁡[−λ−11−λ]=λ2+1\det(R - \lambda I) = \det\begin{bmatrix}-\lambda&-1\\1&-\lambda\end{bmatrix} = \lambda^2 + 1.ad−bc=(−λ)(−λ)−(−1)(1)ad - bc = (-\lambda)(-\lambda) - (-1)(1).
  2. λ2+1=0\lambda^2 + 1 = 0 has no real root; λ=±i\lambda = \pm i.λ2≥0\lambda^2 \ge 0 for every real λ\lambda, so λ2+1≥1\lambda^2 + 1 \ge 1; the two roots are complex.
  3. No real eigenvector.A real eigenvector vv would make Rv=λvRv = \lambda v with both sides real, so λ\lambda would be real; there is no real eigenvalue. The eigenvectors are complex: λ=i\lambda = i has v=(1,−i)⊤v = (1, -i)^\top, since R(1,−i)⊤=(i,1)⊤=i(1,−i)⊤R(1,-i)^\top = (i, 1)^\top = i(1,-i)^\top.
  4. R(x,y)⊤=(−y,x)⊤R(x, y)^\top = (-y, x)^\top, which is v=(x,y)⊤v = (x, y)^\top turned 90°90° anticlockwise.v⋅Rv=−xy+yx=0v \cdot Rv = -xy + yx = 0 and ∥Rv∥=∥v∥\|Rv\| = \|v\|, so for every nonzero real vv, RvRv is a nonzero vector perpendicular to vv, off the line through vv. An eigenvector is exactly a vector that stays on its own line.
  5. λ2+1=0\lambda^2 + 1 = 0, λ=±i\lambda = \pm i; no real eigenvector: a rotation by 90°90° sends every real direction somewhere elseSteps 2 to 4.

Problem 4

Find the eigenvalues and eigenvectors of the shear A=[1101]A = \begin{bmatrix}1&1\\0&1\end{bmatrix}. Can it be diagonalised?

  1. det⁡(A−λI)=(1−λ)2−1⋅0=(λ−1)2\det(A - \lambda I) = (1-\lambda)^2 - 1\cdot 0 = (\lambda - 1)^2.A−λIA - \lambda I is upper triangular, so its determinant is the product of the diagonal entries.
  2. λ=1\lambda = 1, a double root.The only root, counted twice: its algebraic multiplicity, the number of times it is a root, is 22.
  3. A−I=[0100]A - I = \begin{bmatrix}0&1\\0&0\end{bmatrix} has one nonzero row, so rank⁡(A−I)=1\operatorname{rank}(A - I) = 1.The rank decides how many independent eigenvectors λ=1\lambda = 1 has: the null space of a 2×22\times2 matrix has dimension 2−rank⁡2 - \operatorname{rank}.
  4. For v=(x,y)⊤v = (x, y)^\top, (A−I)v=0(A - I)v = 0 says y=0y = 0, so every eigenvector is (x,0)⊤(x, 0)^\top, a multiple of (1,0)⊤(1, 0)^\top.The null space has dimension 2−1=12 - 1 = 1: the geometric multiplicity, the number of independent eigenvectors, is 11.
  5. Not diagonalisable.A=PDP−1A = PDP^{-1} needs two independent eigenvectors for the columns of an invertible PP, and there is only one direction. Directly: DD would have to be diag⁡(1,1)=I\operatorname{diag}(1,1) = I, and then PIP−1=I≠APIP^{-1} = I \neq A.
  6. (λ−1)2(\lambda - 1)^2, λ=1\lambda = 1 (double); A−I=[0100]A - I = \begin{bmatrix}0&1\\0&0\end{bmatrix} has rank 11, so the only eigenvectors are multiples of (1,0)⊤(1,0)^\top; not diagonalisableSteps 1 to 5.

Problem 5

For any 2×22\times2 matrix, show that the trace is the sum of the eigenvalues and the determinant is their product. Verify on Problem 2.

  1. Let A=[abcd]A = \begin{bmatrix}a&b\\c&d\end{bmatrix}. Then det⁡(A−λI)=(a−λ)(d−λ)−bc\det(A - \lambda I) = (a - \lambda)(d - \lambda) - bc.ad−bcad - bc for A−λIA - \lambda I.
  2. (a−λ)(d−λ)−bc=λ2−(a+d)λ+(ad−bc)(a - \lambda)(d - \lambda) - bc = \lambda^2 - (a + d)\lambda + (ad - bc).Expand and collect powers of λ\lambda.
  3. det⁡(A−λI)=(λ−λ1)(λ−λ2)=λ2−(λ1+λ2)λ+λ1λ2\det(A - \lambda I) = (\lambda - \lambda_1)(\lambda - \lambda_2) = \lambda^2 - (\lambda_1 + \lambda_2)\lambda + \lambda_1\lambda_2.λ1\lambda_1 and λ2\lambda_2 are the two roots, complex or repeated if need be, and a quadratic whose λ2\lambda^2 coefficient is 11 is the product of (λ−root)(\lambda - \text{root}) over its roots.
  4. a+d=λ1+λ2a + d = \lambda_1 + \lambda_2 and ad−bc=λ1λ2ad - bc = \lambda_1\lambda_2.Steps 2 and 3 are the same polynomial, so their coefficients of λ\lambda and their constant terms agree.
  5. tr⁡A=λ1+λ2\operatorname{tr}A = \lambda_1+\lambda_2, det⁡A=λ1λ2\det A = \lambda_1\lambda_2; Problem 2: 7=5+27 = 5+2, 10=5⋅210 = 5\cdot2For Problem 2, tr⁡A=4+3=7\operatorname{tr}A = 4 + 3 = 7 and det⁡A=4⋅3−1⋅2=10\det A = 4\cdot3 - 1\cdot2 = 10. This is a quick check on any eigenvalue computation.

Problem 6

Find the eigenvalues and eigenvectors of A=[210121012]A = \begin{bmatrix}2&1&0\\1&2&1\\0&1&2\end{bmatrix}.

Write t=2−λt = 2 - \lambda, so the diagonal entries of A−λIA - \lambda I are all tt and the off-diagonal entries are those of AA. Each eigenvector is written v=(x,y,z)⊤v = (x, y, z)^\top.

  1. det⁡(A−λI)=tdet⁡[t11t]−1⋅det⁡[110t]+0\det(A - \lambda I) = t\det\begin{bmatrix}t&1\\1&t\end{bmatrix} - 1\cdot\det\begin{bmatrix}1&1\\0&t\end{bmatrix} + 0.Cofactor expansion along the first row, (t,1,0)(t, 1, 0): each entry times the determinant left after deleting its row and column, with signs +,−,++, -, +.
  2. t(t2−1)−t=t3−2t=t(t2−2)t(t^2 - 1) - t = t^3 - 2t = t(t^2 - 2).The two 2×22\times2 determinants are t2−1t^2 - 1 and t−0{t - 0}.
  3. det⁡(A−λI)=(2−λ)((2−λ)2−2)\det(A - \lambda I) = (2 - \lambda)\big((2 - \lambda)^2 - 2\big).Substitute t=2−λt = 2 - \lambda back.
  4. λ=2\lambda = 2, or 2−λ=±22 - \lambda = \pm\sqrt2, that is λ=2∓2\lambda = 2 \mp \sqrt2.Set each factor to 00. Check with the trace: 2+(2−2)+(2+2)=6=tr⁡A2 + (2 - \sqrt2) + (2 + \sqrt2) = 6 = \operatorname{tr}A.
  5. λ=2\lambda = 2 (t=0t = 0): the rows of A−2IA - 2I give y=0y = 0, x+z=0x + z = 0, y=0y = 0, so v=(1,0,−1)⊤v = (1, 0, -1)^\top.The null space of A−2IA - 2I, whose rows are (0,1,0)(0,1,0), (1,0,1)(1,0,1), (0,1,0)(0,1,0).
  6. λ=2+2\lambda = 2 + \sqrt2 (t=−2t = -\sqrt2): the rows give −2 x+y=0-\sqrt2\,x + y = 0, x−2 y+z=0x - \sqrt2\,y + z = 0, y−2 z=0y - \sqrt2\,z = 0.The null space of A−(2+2)IA - (2 + \sqrt2)I, whose diagonal entries are −2-\sqrt2.
  7. y=2 xy = \sqrt2\,x and y=2 zy = \sqrt2\,z, so z=xz = x; with x=1x = 1, v=(1,2,1)⊤v = (1, \sqrt2, 1)^\top.The first and third equations; the middle one then holds, 1−2+1=01 - 2 + 1 = 0, as it must because the determinant is 00.
  8. λ=2−2\lambda = 2 - \sqrt2 (t=2t = \sqrt2): the same equations with 2\sqrt2 in place of −2-\sqrt2 give y=−2 x=−2 zy = -\sqrt2\,x = -\sqrt2\,z, so v=(1,−2,1)⊤v = (1, -\sqrt2, 1)^\top.Only the sign of tt changed, so only the sign of yy changes.
  9. λ=2\lambda = 2 with (1,0,−1)⊤(1,0,-1)^\top; λ=2±2\lambda = 2\pm\sqrt2 with (1,±2,1)⊤(1, \pm\sqrt2, 1)^\topCheck: A(1,2,1)⊤=(2+2, 2+22, 2+2)⊤=(2+2)(1,2,1)⊤A(1,\sqrt2,1)^\top = (2+\sqrt2,\ 2+2\sqrt2,\ 2+\sqrt2)^\top = (2+\sqrt2)(1,\sqrt2,1)^\top. AA is symmetric, and the three eigenvectors are mutually perpendicular (Problem 8).

Problem 7

Diagonalise the matrix of Problem 2 as A=PDP−1A = PDP^{-1} and use it to compute A5A^5.

  1. P=[111−2]P = \begin{bmatrix}1&1\\1&-2\end{bmatrix}, D=diag⁡(5,2)D = \operatorname{diag}(5, 2).The columns of PP are the eigenvectors from Problem 2, and DD lists their eigenvalues in the same order. Then AP=PDAP = PD column by column, which is Av=λvAv = \lambda v twice.
  2. det⁡P=1⋅(−2)−1⋅1=−3\det P = 1\cdot(-2) - 1\cdot1 = -3, so P−1=−13[−2−1−11]=[2/31/31/3−1/3]P^{-1} = -\tfrac13\begin{bmatrix}-2&-1\\-1&1\end{bmatrix} = \begin{bmatrix}2/3&1/3\\1/3&-1/3\end{bmatrix}.The 2×22\times2 inverse: swap the diagonal, negate the off-diagonal, divide by the determinant. det⁡P≠0\det P \neq 0 because the two eigenvectors are independent.
  3. PD=[525−4]PD = \begin{bmatrix}5&2\\5&-4\end{bmatrix} and PDP−1=[4123]=APDP^{-1} = \begin{bmatrix}4&1\\2&3\end{bmatrix} = A.Check the diagonalisation by multiplying it out: the first row of PDP−1PDP^{-1} is (10/3+2/3, 5/3−2/3)=(4,1)(10/3 + 2/3,\ 5/3 - 2/3) = (4, 1).
  4. A5=PDP−1 PDP−1⋯PDP−1=PD5P−1A^5 = PDP^{-1}\,PDP^{-1}\cdots PDP^{-1} = PD^5P^{-1}.Each inner P−1PP^{-1}P is II, so the five copies collapse. Powers of a diagonal matrix are powers of its entries: D5=diag⁡(3125,32)D^5 = \operatorname{diag}(3125, 32).
  5. PD5=[3125323125−64]PD^5 = \begin{bmatrix}3125&32\\3125&-64\end{bmatrix}.Multiplying by a diagonal matrix on the right scales the columns of PP by 555^5 and 252^5.
  6. PD5P−1=13[6250+323125−326250−643125+64]=[2094103120621063]PD^5P^{-1} = \tfrac13\begin{bmatrix}6250+32&3125-32\\6250-64&3125+64\end{bmatrix} = \begin{bmatrix}2094&1031\\2062&1063\end{bmatrix}.Row times column with P−1=13[211−1]P^{-1} = \tfrac13\begin{bmatrix}2&1\\1&-1\end{bmatrix}.
  7. P=[111−2]P = \begin{bmatrix}1&1\\1&-2\end{bmatrix}, D=diag⁡(5,2)D = \operatorname{diag}(5,2), P−1=−13[−2−1−11]P^{-1} = -\tfrac13\begin{bmatrix}-2&-1\\-1&1\end{bmatrix}; A5=PD5P−1=[2094103120621063]A^5 = PD^5P^{-1} = \begin{bmatrix}2094&1031\\2062&1063\end{bmatrix}Check: (1,1)⊤(1,1)^\top is an eigenvector with eigenvalue 55, so A5(1,1)⊤=3125(1,1)⊤A^5(1,1)^\top = 3125(1,1)^\top, and both rows do sum to 31253125.

Problem 8

Show that the eigenvectors of Problem 1 are orthogonal, and write A=QΛQ⊤A = Q\Lambda Q^\top with QQ orthogonal. Why does this always happen for symmetric matrices?

  1. (1,1)⊤⋅(1,−1)⊤=1−1=0(1,1)^\top \cdot (1,-1)^\top = 1 - 1 = 0.Orthogonal means a zero dot product.
  2. Q=12[111−1]Q = \tfrac1{\sqrt2}\begin{bmatrix}1&1\\1&-1\end{bmatrix}, Λ=diag⁡(3,1)\Lambda = \operatorname{diag}(3, 1).Divide each eigenvector by its length 2\sqrt2 so the columns are unit vectors; they are still eigenvectors. Λ\Lambda lists the eigenvalues in column order.
  3. Q⊤Q=12[2002]=IQ^\top Q = \tfrac12\begin{bmatrix}2&0\\0&2\end{bmatrix} = I.The columns are orthogonal unit vectors, so QQ is orthogonal and Q−1=Q⊤Q^{-1} = Q^\top.
  4. QΛQ⊤=12[313−1][111−1]=12[4224]=AQ\Lambda Q^\top = \tfrac12\begin{bmatrix}3&1\\3&-1\end{bmatrix}\begin{bmatrix}1&1\\1&-1\end{bmatrix} = \tfrac12\begin{bmatrix}4&2\\2&4\end{bmatrix} = A.A=QΛQ−1A = Q\Lambda Q^{-1} as in Problem 7, with Q−1Q^{-1} replaced by Q⊤Q^\top; multiplying it out confirms it.
  5. For symmetric AA with Av1=λ1v1Av_1 = \lambda_1v_1, Av2=λ2v2Av_2 = \lambda_2v_2: λ1v1⊤v2=(Av1)⊤v2=v1⊤A⊤v2=v1⊤Av2=λ2v1⊤v2\lambda_1v_1^\top v_2 = (Av_1)^\top v_2 = v_1^\top A^\top v_2 = v_1^\top Av_2 = \lambda_2v_1^\top v_2.Move AA across the dot product: (Av1)⊤=v1⊤A⊤(Av_1)^\top = v_1^\top A^\top, and A⊤=AA^\top = A is where symmetry is used.
  6. (λ1−λ2) v1⊤v2=0(\lambda_1 - \lambda_2)\,v_1^\top v_2 = 0, so v1⊤v2=0v_1^\top v_2 = 0 when λ1≠λ2\lambda_1 \neq \lambda_2.Subtract the two ends of step 5; a product is 00 only if a factor is. A repeated eigenvalue of a symmetric matrix has as many independent eigenvectors as its multiplicity (the spectral theorem), and they can be chosen orthogonal, so A=QΛQ⊤A = Q\Lambda Q^\top holds for every real symmetric matrix.
  7. Q=12[111−1]Q = \tfrac1{\sqrt2}\begin{bmatrix}1&1\\1&-1\end{bmatrix}, Λ=diag⁡(3,1)\Lambda = \operatorname{diag}(3,1), A=QΛQ⊤A = Q\Lambda Q^\top; for symmetric AA, λ1v1⊤v2=(Av1)⊤v2=v1⊤Av2=λ2v1⊤v2\lambda_1 v_1^\top v_2 = (Av_1)^\top v_2 = v_1^\top Av_2 = \lambda_2 v_1^\top v_2, so distinct eigenvalues force v1⊤v2=0v_1^\top v_2 = 0Steps 1 to 6.

Problem 9

Apply A=[2112]A = \begin{bmatrix}2&1\\1&2\end{bmatrix} repeatedly to x0=(1,0)⊤x_0 = (1, 0)^\top. Find Akx0A^kx_0 in closed form and say what direction it approaches and how fast.

  1. x0=12(1,1)⊤+12(1,−1)⊤x_0 = \tfrac12(1,1)^\top + \tfrac12(1,-1)^\top.Write x0x_0 in the eigenvectors from Problem 1, where AA acts by scaling; the coefficients solve c1+c2=1c_1 + c_2 = 1, c1−c2=0c_1 - c_2 = 0.
  2. Ak(1,1)⊤=3k(1,1)⊤A^k(1,1)^\top = 3^k(1,1)^\top and Ak(1,−1)⊤=1k(1,−1)⊤=(1,−1)⊤A^k(1,-1)^\top = 1^k(1,-1)^\top = (1,-1)^\top.Each application of AA multiplies an eigenvector by its eigenvalue, so kk applications multiply it by the kk-th power.
  3. Akx0=3k2(1,1)⊤+12(1,−1)⊤A^kx_0 = \tfrac{3^k}{2}(1,1)^\top + \tfrac12(1,-1)^\top.AkA^k is linear, so apply step 2 to each term of step 1. Check: k=1k = 1 gives (2,1)⊤=Ax0(2, 1)^\top = Ax_0, and k=2k = 2 gives (5,4)⊤(5, 4)^\top.
  4. Akx0=3k2((1,1)⊤+3−k(1,−1)⊤)A^kx_0 = \tfrac{3^k}{2}\big((1,1)^\top + 3^{-k}(1,-1)^\top\big).Factor out the growing term; the positive scale factor does not change the direction.
  5. tan⁡θk=3−k\tan\theta_k = 3^{-k}, where θk\theta_k is the angle between Akx0A^kx_0 and (1,1)⊤(1,1)^\top.In step 4 the two eigenvectors are perpendicular and have the same length 2\sqrt2, so the ratio of their coefficients, 3−k3^{-k}, is the tangent of that angle, and it goes to 00.
  6. Akx0=3k2(1,1)⊤+12(1,−1)⊤A^kx_0 = \tfrac{3^k}{2}(1,1)^\top + \tfrac12(1,-1)^\top; the direction approaches (1,1)/2(1,1)/\sqrt2, the error shrinking like (1/3)k(1/3)^k — the ratio of the two eigenvaluesThis is power iteration: repeated multiplication (with a rescaling each step to keep the length finite, since it grows like 3k3^k) finds the eigenvector of the largest eigenvalue, at a rate set by λ2/λ1=1/3\lambda_2/\lambda_1 = 1/3.

Problem 10

For A=[3045]A = \begin{bmatrix}3&0\\4&5\end{bmatrix}, find the eigenvalues of A⊤AA^\top A and hence the singular values of AA. Check them against det⁡A\det A.

  1. A⊤A=[3405][3045]=[9+16202025]=[25202025]A^\top A = \begin{bmatrix}3&4\\0&5\end{bmatrix}\begin{bmatrix}3&0\\4&5\end{bmatrix} = \begin{bmatrix}9+16&20\\20&25\end{bmatrix} = \begin{bmatrix}25&20\\20&25\end{bmatrix}.Row times column; the result is symmetric, as A⊤AA^\top A always is.
  2. det⁡(A⊤A−λI)=(25−λ)2−400\det(A^\top A - \lambda I) = (25 - \lambda)^2 - 400, so 25−λ=±2025 - \lambda = \pm20 and λ=45\lambda = 45 or 55.The same shape as Problem 1: a square minus a constant. The eigenvectors are (1,1)⊤(1,1)^\top and (1,−1)⊤(1,-1)^\top again.
  3. σ1=45=35\sigma_1 = \sqrt{45} = 3\sqrt5 and σ2=5\sigma_2 = \sqrt5.The singular values are the square roots of the eigenvalues of A⊤AA^\top A: if A=UΣV⊤A = U\Sigma V^\top with UU, VV orthogonal and Σ=diag⁡(σ1,σ2)\Sigma = \operatorname{diag}(\sigma_1, \sigma_2), then A⊤A=VΣ2V⊤A^\top A = V\Sigma^2V^\top, an orthogonal diagonalisation as in Problem 8 with Λ=Σ2\Lambda = \Sigma^2. The eigenvalues are never negative, since v⊤A⊤Av=∥Av∥2≥0v^\top A^\top Av = \|Av\|^2 \ge 0.
  4. σ1σ2=35⋅5=15\sigma_1\sigma_2 = 3\sqrt5\cdot\sqrt5 = 15, and det⁡A=3⋅5−0⋅4=15\det A = 3\cdot5 - 0\cdot4 = 15.det⁡(A⊤A)=(det⁡A)2\det(A^\top A) = (\det A)^2 is the product of its eigenvalues σ12σ22\sigma_1^2\sigma_2^2 (Problem 5), so σ1σ2=∣det⁡A∣\sigma_1\sigma_2 = |\det A|; here 45⋅5=225=15245\cdot5 = 225 = 15^2.
  5. A⊤A=[25202025]A^\top A = \begin{bmatrix}25&20\\20&25\end{bmatrix} with eigenvalues 4545 and 55; singular values 353\sqrt5 and 5\sqrt5; their product 15=∣det⁡A∣15 = |\det A|The eigenvalues of AA itself are 33 and 55, its diagonal entries since it is triangular: same product, but not the singular values. They coincide, up to sign, only for special matrices such as symmetric ones.

Where this goes wrong

1. Wrong sign in the characteristic matrix

Av=λvAv = \lambda v has λ\lambda on the right, and moving it across is where the sign slips.

  1. Av=λvAv = \lambda v for A=[2112]A = \begin{bmatrix}2&1\\1&2\end{bmatrix}Right so far: the eigen-equation of Problem 1.
  2. “Collect everything on one side: AA and λI\lambda I together.”The step that causes the mistake: the two terms are gathered without tracking the sign.
  3. det⁡(A+λI)=(2+λ)2−1=0\det(A + \lambda I) = (2 + \lambda)^2 - 1 = 0, so λ=−1\lambda = -1 or λ=−3\lambda = -3Av=λvAv = \lambda v rearranges to Av−λv=0Av - \lambda v = 0, that is (A−λI)v=0(A - \lambda I)v = 0. With the plus sign every root comes out negated: −1-1 and −3-3 instead of 11 and 33. The trace catches it: the eigenvalues must sum to tr⁡A=4\operatorname{tr}A = 4 (Problem 5), and these sum to −4-4.

2. Reading the eigenvector off a row of A − λI

The rows of A−λIA - \lambda I are the first vectors on the page after the eigenvalue is found, and it is tempting to take one as the answer.

  1. λ=3\lambda = 3 in Problem 1, and A−3I=[−111−1]A - 3I = \begin{bmatrix}-1&1\\1&-1\end{bmatrix}Right so far.
  2. “The eigenvector comes from A−3IA - 3I, so read it off the first row.”The step that causes the mistake: the row is taken as the vector instead of as an equation the vector must satisfy.
  3. v=(−1,1)⊤v = (-1, 1)^\topThe eigenvector is in the null space of A−3IA - 3I, which is orthogonal to its rows: the row (−1,1)(-1, 1) says −x+y=0-x + y = 0, and the solution (1,1)⊤(1, 1)^\top is perpendicular to it. Here the wrong vector differs from the right one only in the sign of one entry, which is why the habit survives, but A(−1,1)⊤=(−1,1)⊤A(-1,1)^\top = (-1,1)^\top: it is the eigenvector for λ=1\lambda = 1, not 33. For a 2×22\times2 row (a,b)(a, b), the correct move is to swap the entries and negate one, giving (b,−a)(b, -a).

3. Assuming a double eigenvalue gives two eigenvectors

Two eigenvalues usually bring two eigenvectors, and a double eigenvalue looks like two eigenvalues.

  1. The shear A=[1101]A = \begin{bmatrix}1&1\\0&1\end{bmatrix} has det⁡(A−λI)=(λ−1)2\det(A - \lambda I) = (\lambda - 1)^2Right so far: Problem 4.
  2. “A 2×22\times2 matrix has two eigenvalues, so it has two eigenvectors.”The count that causes the mistake: it counts roots, not independent solutions.
  3. “λ=1\lambda = 1 twice, so R2\mathbb{R}^2 is spanned by eigenvectors and A=PDP−1A = PDP^{-1}.”The algebraic multiplicity is 22 but the geometric multiplicity is 11: A−IA - I has rank 11, so its null space holds only multiples of (1,0)⊤(1,0)^\top (Problem 4). The rank computation is the test; without it the claim would force A=PIP−1=IA = PIP^{-1} = I.

4. P and P⁻¹ swapped

PP and P−1P^{-1} appear on opposite sides of DD, and which side is which is easy to forget.

  1. P=[111−2]P = \begin{bmatrix}1&1\\1&-2\end{bmatrix}, D=diag⁡(5,2)D = \operatorname{diag}(5, 2) for the matrix of Problem 2Right so far: Problem 7.
  2. “AA is DD with PP on one side and P−1P^{-1} on the other.”The half-memory that causes the mistake: it does not say which side.
  3. A5=P−1D5PA^5 = P^{-1}D^5PA=PDP−1A = PDP^{-1} because AP=PDAP = PD: PP holds the eigenvectors as columns and multiplies on the left. The swapped product is not A5A^5; here it is [2094206210311063]\begin{bmatrix}2094&2062\\1031&1063\end{bmatrix}, the transpose of the right answer — only because this PP happens to be symmetric; in general the swapped product is just a different matrix. The eigenvector check catches it: A5(1,1)⊤A^5(1,1)^\top must be 3125(1,1)⊤3125(1,1)^\top, but the first row of this matrix sums to 41564156.

5. Complex eigenvalues read as no eigenvalues

In school algebra a negative discriminant means “no solution”, and the habit carries over.

  1. The rotation R=[0−110]R = \begin{bmatrix}0&-1\\1&0\end{bmatrix} has det⁡(R−λI)=λ2+1\det(R - \lambda I) = \lambda^2 + 1Right so far: Problem 3.
  2. “λ2+1=0\lambda^2 + 1 = 0 has no solution.”The habit that causes the mistake: no real solution is read as no solution.
  3. “The rotation has no eigenvalues.”The characteristic polynomial of an n×nn\times n matrix always has nn roots over C\mathbb{C}, counted with multiplicity; here they are ±i\pm i. What the rotation lacks is a real eigenvector, because no real direction stays on its own line under a 90°90° turn.

Print this set: eigenvalues-and-eigenvectors.pdf (problems, answers, and worked solutions on separate pages).