MAT1041 · Week 13 · - Eigenvalues and Eigenvectors

Eigenvalues and eigenvectors

You learn to set up the characteristic equation and find eigenvalues, to compute a basis of the eigenspace for each eigenvalue, to tell algebraic and geometric multiplicity apart, and the shortcuts relating eigenvalues to trace, determinant, powers and inverse.

eigenvalueeigenvectorcharacteristic polynomialcharacteristic equationeigenspacealgebraic multiplicitygeometric multiplicitytriangular matrixtracedeterminant

When a matrix multiplies a vector, it generally both rotates and stretches it. But there are some special vectors that the matrix does not rotate, it only multiplies by a number. These are the eigenvectors, and the multiplier is the eigenvalue. This week we learn to find them; on the next page, if there are enough eigenvectors, we will make the matrix diagonal. The computing side is completely familiar: determinants (weeks 4-5) and the null space (week 9).

Golden rules:

  1. First the eigenvalue: det⁡(λI−A)=0\det(\lambda I - A) = 0.
  2. Then for each eigenvalue the eigenspace: solve (A−λI)x=0(A - \lambda I)\mathbf{x} = \mathbf{0}.
  3. A−λIA - \lambda I must come out singular. If there is no zero row, you have found the eigenvalue wrongly.
  4. Check: the sum of the eigenvalues is the trace, their product is the determinant.

Definition

Let AA be n×nn \times n. If, for a nonzero v\mathbf{v} and a number λ\lambda,

Av=λvA\mathbf{v} = \lambda\mathbf{v}

then λ\lambda is an eigenvalue of AA and v\mathbf{v} is an eigenvector corresponding to λ\lambda.

  • v=0\mathbf{v} = \mathbf{0} is never counted as an eigenvector (A0=λ0A\mathbf{0} = \lambda\mathbf{0} for every λ\lambda, which would be meaningless).
  • λ=0\lambda = 0 can be an eigenvalue: it means Av=0A\mathbf{v} = \mathbf{0}.
  • If v\mathbf{v} is an eigenvector, so is every cvc\mathbf{v} (c≠0c \neq 0). Eigenvectors are not single, they come along a line (or subspace).

Geometry. AvA\mathbf{v} lies on the same line as v\mathbf{v}. If λ>1\lambda > 1 it stretches, if 0<λ<10 < \lambda < 1 it shrinks, if λ<0\lambda < 0 it reverses direction. For the projection onto the xx-axis in the plane, e1\mathbf{e}_1 corresponds to eigenvalue 1 and e2\mathbf{e}_2 to eigenvalue 0.

The characteristic equation

Av=λv  ⟺  (λI−A)v=0A\mathbf{v} = \lambda\mathbf{v} \iff (\lambda I - A)\mathbf{v} = \mathbf{0}. For a nonzero solution to exist, λI−A\lambda I - A must be singular:

det⁡(λI−A)=0.\det(\lambda I - A) = 0.

The left side is a polynomial of degree nn in λ\lambda: the characteristic polynomial. Its roots are exactly the eigenvalues. (Using det⁡(A−λI)\det(A - \lambda I) gives the same roots; the sign differs only if nn is odd.)

Example 1. A=[2321]A = \begin{bmatrix} 2 & 3 \\ 2 & 1 \end{bmatrix}.

det⁡(λI−A)=∣λ−2−3−2λ−1∣=(λ−2)(λ−1)−6=λ2−3λ−4=(λ−4)(λ+1).\det(\lambda I - A) = \begin{vmatrix} \lambda - 2 & -3 \\ -2 & \lambda - 1 \end{vmatrix} = (\lambda - 2)(\lambda - 1) - 6 = \lambda^2 - 3\lambda - 4 = (\lambda - 4)(\lambda + 1).

The eigenvalues are 44 and −1-1. The 2×22 \times 2 shortcut: λ2−(tr⁡A)λ+det⁡A\lambda^2 - (\operatorname{tr}A)\lambda + \det A, here trace 3, det 2−6=−42 - 6 = -4.

The eigenspace

All the eigenvectors for a given λ\lambda plus 0\mathbf{0} form the null space of A−λIA - \lambda I. It is called the eigenspace EλE_\lambda of λ\lambda; it is a subspace.

Example 1 continued:

λ=4:A−4I=[−232−3]→R2←R2+R1[−2300]⇒−2x+3y=0,E4=span⁡{(3,2)}.\lambda = 4: \quad A - 4I = \begin{bmatrix} -2 & 3 \\ 2 & -3 \end{bmatrix} \xrightarrow{R_2 \leftarrow R_2 + R_1} \begin{bmatrix} -2 & 3 \\ 0 & 0 \end{bmatrix} \Rightarrow -2x + 3y = 0, \quad E_4 = \operatorname{span}\{(3, 2)\}.

λ=−1:A+I=[3322]→R2←R2−23R1[3300]⇒x+y=0,E−1=span⁡{(1,−1)}.\lambda = -1: \quad A + I = \begin{bmatrix} 3 & 3 \\ 2 & 2 \end{bmatrix} \xrightarrow{R_2 \leftarrow R_2 - \frac23 R_1} \begin{bmatrix} 3 & 3 \\ 0 & 0 \end{bmatrix} \Rightarrow x + y = 0, \quad E_{-1} = \operatorname{span}\{(1, -1)\}.

Check: A(3,2)T=(6+6,6+2)=(12,8)=4(3,2)A(3, 2)^T = (6 + 6, 6 + 2) = (12, 8) = 4(3, 2). A(1,−1)T=(2−3,2−1)=(−1,1)=−1⋅(1,−1)A(1, -1)^T = (2 - 3, 2 - 1) = (-1, 1) = -1 \cdot (1, -1).

PITFALL · getting a unique solution. If solving (A−λI)x=0(A - \lambda I)\mathbf{x} = \mathbf{0} gives only x=0\mathbf{x} = \mathbf{0}, there is a calculation error in the eigenvalue. At a true eigenvalue A−λIA - \lambda I is always singular and at least one zero row appears. This is both a warning and a free check.

A 3×3 example and multiplicity

Example 2. A=[32−223−222−1]A = \begin{bmatrix} 3 & 2 & -2 \\ 2 & 3 & -2 \\ 2 & 2 & -1 \end{bmatrix}.

Characteristic polynomial (the calculation is long; using the common structure in the first two columns or expanding along the first row helps): λ3−5λ2+7λ−3=(λ−1)2(λ−3)\lambda^3 - 5\lambda^2 + 7\lambda - 3 = (\lambda - 1)^2(\lambda - 3). Check: trace 3+3−1=5=1+1+33 + 3 - 1 = 5 = 1 + 1 + 3, det =1⋅1⋅3=3= 1 \cdot 1 \cdot 3 = 3.

λ=1:A−I=[22−222−222−2]→[11−1000000]:x=−y+z,\lambda = 1: \quad A - I = \begin{bmatrix} 2 & 2 & -2 \\ 2 & 2 & -2 \\ 2 & 2 & -2 \end{bmatrix} \to \begin{bmatrix} 1 & 1 & -1 \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{bmatrix}: \quad x = -y + z,

E1=span⁡{(−1,1,0), (1,0,1)}(dimension 2).E_1 = \operatorname{span}\{(-1, 1, 0),\ (1, 0, 1)\} \quad (\text{dimension } 2).

λ=3:A−3I=[02−220−222−4]→R1←12R1, R2←R2−2R1R1↔R3[11−20−2202−2]→R3←R3+R2[11−20−22000]\lambda = 3: \quad A - 3I = \begin{bmatrix} 0 & 2 & -2 \\ 2 & 0 & -2 \\ 2 & 2 & -4 \end{bmatrix} \xrightarrow[R_1 \leftarrow \frac12R_1,\ R_2 \leftarrow R_2 - 2R_1]{R_1 \leftrightarrow R_3} \begin{bmatrix} 1 & 1 & -2 \\ 0 & -2 & 2 \\ 0 & 2 & -2 \end{bmatrix} \xrightarrow{R_3 \leftarrow R_3 + R_2} \begin{bmatrix} 1 & 1 & -2 \\ 0 & -2 & 2 \\ 0 & 0 & 0 \end{bmatrix}

y=zy = z, x=−y+2z=zx = -y + 2z = z: E3=span⁡{(1,1,1)}E_3 = \operatorname{span}\{(1, 1, 1)\}.

The two multiplicities:

Definitionλ=1\lambda = 1 in Example 2λ=2\lambda = 2 in Q7
Algebraicthe multiplicity of the root in the characteristic polynomial22
Geometricdim⁡Eλ\dim E_\lambda21

Always 1≤geometric≤algebraic1 \le \text{geometric} \le \text{algebraic}. Whether the two are equal is the key to diagonalization (next page).

Shortcuts

If the eigenvalues of AA are λ1,…,λn\lambda_1, \dots, \lambda_n (with multiplicities):

InformationResult
λ1+⋯+λn\lambda_1 + \dots + \lambda_ntr⁡A\operatorname{tr}A
λ1λ2⋯λn\lambda_1\lambda_2\cdots\lambda_ndet⁡A\det A
AA triangulareigenvalues on the diagonal
00 is an eigenvalueAA is singular
AkA^keigenvalues λk\lambda^k, same eigenvectors
A−1A^{-1}1/λ1/\lambda, same eigenvectors
A+cIA + cIλ+c\lambda + c
ATA^Tsame eigenvalues (eigenvectors may differ)

One more item for the invertible matrix theorem: AA is invertible   ⟺  \iff 0 is not an eigenvalue of AA.

There may be no real eigenvalue

The characteristic polynomial of the 90∘90^\circ rotation matrix [0−110]\begin{bmatrix} 0 & -1 \\ 1 & 0 \end{bmatrix} is λ2+1\lambda^2 + 1: no real root. In this course we usually work with real eigenvalues; in week 14 we will see that the eigenvalues of symmetric matrices are always real.

Common mistakes

WrongRight
Writing det⁡(A−λ)\det(A - \lambda)det⁡(A−λI)\det(A - \lambda I): λ\lambda only on the diagonal
Counting v=0\mathbf{v} = \mathbf{0} as an eigenvectorAn eigenvector must be nonzero
Expanding the characteristic polynomial and losing the factorsKeep the common factor, take (λ−3)(\lambda - 3) out
Thinking the eigenspace is the null space of AAIt is the null space of A−λIA - \lambda I
Writing a single vector for the eigenspace basis and forgetting the second free variableAs many vectors as free variables
Thinking row operations preserve eigenvaluesRow operations change eigenvalues; don’t reduce AA first
Assuming algebraic and geometric multiplicity are equalCompute; they differ in Q7

Test yourself

Exam-style questions. Solve them on paper first, then reveal the answer.

Question 1easy

Show that v=(1,1)\mathbf{v} = (1, 1) is an eigenvector of A=[3122]A = \begin{bmatrix} 3 & 1 \\ 2 & 2 \end{bmatrix} and find the corresponding eigenvalue.

Show answer
λ=4\lambda = 4

Av=(3+1, 2+2)=(4,4)=4(1,1)A\mathbf{v} = (3 + 1,\ 2 + 2) = (4, 4) = 4(1, 1). So Av=4vA\mathbf{v} = 4\mathbf{v}, v≠0\mathbf{v} \neq \mathbf{0}: an eigenvector, with eigenvalue 4.

Question 2easy

Find the eigenvalues of [25−10−34007]\begin{bmatrix} 2 & 5 & -1 \\ 0 & -3 & 4 \\ 0 & 0 & 7 \end{bmatrix}.

Show answer
2, −3, 72,\ -3,\ 7

A triangular matrix: λI−A\lambda I - A is triangular too, and its determinant is the product of the diagonal (λ−2)(λ+3)(λ−7)(\lambda - 2)(\lambda + 3)(\lambda - 7).
The eigenvalues of a triangular matrix are its diagonal entries.

Question 3medium

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

Show answer
λ2−4λ−5\lambda^2 - 4\lambda - 5; λ=5,−1\lambda = 5, -1

det⁡(λI−A)=(λ−1)(λ−3)−8=λ2−4λ−5=(λ−5)(λ+1)\det(\lambda I - A) = (\lambda - 1)(\lambda - 3) - 8 = \lambda^2 - 4\lambda - 5 = (\lambda - 5)(\lambda + 1).
Shortcut (2×22 \times 2): λ2−(tr⁡A)λ+det⁡A\lambda^2 - (\operatorname{tr}A)\lambda + \det A; trace 4, det 3−8=−53 - 8 = -5.

Question 4medium

For A=[1423]A = \begin{bmatrix} 1 & 4 \\ 2 & 3 \end{bmatrix}, find a basis for each eigenspace.

Show answer
E5E_5: {(1,1)}\{(1, 1)\}; E−1E_{-1}: {(−2,1)}\{(-2, 1)\}

λ=5\lambda = 5: A−5I=[−442−2]A - 5I = \begin{bmatrix} -4 & 4 \\ 2 & -2 \end{bmatrix}, the rows are proportional, x=yx = y: (1,1)(1, 1).λ=−1\lambda = -1: A+I=[2424]A + I = \begin{bmatrix} 2 & 4 \\ 2 & 4 \end{bmatrix}, x+2y=0x + 2y = 0: (−2,1)(-2, 1).
Check: A(−2,1)T=(−2+4,−4+3)=(2,−1)=−1⋅(−2,1)A(-2, 1)^T = (-2 + 4, -4 + 3) = (2, -1) = -1 \cdot (-2, 1).

Question 5medium

Find the eigenvalues of [1224]\begin{bmatrix} 1 & 2 \\ 2 & 4 \end{bmatrix}. What does one of them tell you about the matrix?

Show answer
00 and 55; the matrix is singular

Trace 5, det 4−4=04 - 4 = 0: λ2−5λ=λ(λ−5)\lambda^2 - 5\lambda = \lambda(\lambda - 5).λ=0\lambda = 0 is an eigenvalue   ⟺  \iff Av=0A\mathbf{v} = \mathbf{0} has a nonzero solution   ⟺  \iff AA is singular.

Question 6medium

Find the eigenvalues of A=[102030201]A = \begin{bmatrix} 1 & 0 & 2 \\ 0 & 3 & 0 \\ 2 & 0 & 1 \end{bmatrix} with their algebraic multiplicities.

Show answer
λ=3\lambda = 3 (multiplicity 2), λ=−1\lambda = -1 (multiplicity 1)

Expand det⁡(λI−A)\det(\lambda I - A) along row 2 (its only entry is λ−3\lambda - 3): (λ−3)[(λ−1)2−4](\lambda - 3)\left[(\lambda - 1)^2 - 4\right].(λ−1)2−4=(λ−3)(λ+1)(\lambda - 1)^2 - 4 = (\lambda - 3)(\lambda + 1). Total: (λ−3)2(λ+1)(\lambda - 3)^2(\lambda + 1).
Check: the sum of the eigenvalues 3+3−1=5=tr⁡A3 + 3 - 1 = 5 = \operatorname{tr}A.

Question 7hard

For A=[210020003]A = \begin{bmatrix} 2 & 1 & 0 \\ 0 & 2 & 0 \\ 0 & 0 & 3 \end{bmatrix}, find the algebraic and geometric multiplicity of λ=2\lambda = 2.

Show answer
Algebraic 2, geometric 1

Triangular: the eigenvalues are 2,2,32, 2, 3. Algebraic multiplicity 2.A−2I=[010000001]A - 2I = \begin{bmatrix} 0 & 1 & 0 \\ 0 & 0 & 0 \\ 0 & 0 & 1 \end{bmatrix}: y=0y = 0, z=0z = 0, xx free. E2=span⁡{(1,0,0)}E_2 = \operatorname{span}\{(1, 0, 0)\}, dimension 1.
The geometric multiplicity is always ≤\le the algebraic; here it is smaller.

Question 8medium

λ=3\lambda = 3 is an eigenvalue of an invertible matrix AA with eigenvector v\mathbf{v}. Find the corresponding eigenvalues of A2−2A+IA^2 - 2A + I and of A−1A^{-1}.

Show answer
44 and 1/31/3

A2v=9vA^2\mathbf{v} = 9\mathbf{v}, so (A2−2A+I)v=(9−6+1)v=4v(A^2 - 2A + I)\mathbf{v} = (9 - 6 + 1)\mathbf{v} = 4\mathbf{v}.
Multiply Av=3vA\mathbf{v} = 3\mathbf{v} from the left by A−1A^{-1}: v=3A−1v\mathbf{v} = 3A^{-1}\mathbf{v}, A−1v=13vA^{-1}\mathbf{v} = \frac13\mathbf{v}.
The eigenvector stays the same.

Question 9medium

A 3×33 \times 3 matrix AA has eigenvalues 1,2,−31, 2, -3. Find det⁡A\det A and tr⁡A\operatorname{tr}A.

Show answer
det⁡A=−6\det A = -6, tr⁡A=0\operatorname{tr}A = 0

Determinant = the product of the eigenvalues: 1⋅2⋅(−3)=−61 \cdot 2 \cdot (-3) = -6.
Trace = the sum of the eigenvalues: 1+2−3=01 + 2 - 3 = 0.

Question 10hard

Does A=[0−110]A = \begin{bmatrix} 0 & -1 \\ 1 & 0 \end{bmatrix} (rotation by 90∘90^\circ) have real eigenvalues?

Show answer
No (the eigenvalues are ±i\pm i)

det⁡(λI−A)=λ2+1\det(\lambda I - A) = \lambda^2 + 1. It has no real root.
Geometric meaning: a 90∘90^\circ rotation leaves no nonzero vector along its own direction.

Question 11hard

Find kk such that λ=4\lambda = 4 is an eigenvalue of [2k13]\begin{bmatrix} 2 & k \\ 1 & 3 \end{bmatrix}.

Show answer
k=2k = 2

λ=4\lambda = 4 is an eigenvalue   ⟺  det⁡(A−4I)=0\iff \det(A - 4I) = 0. ∣−2k1−1∣=2−k=0\begin{vmatrix} -2 & k \\ 1 & -1 \end{vmatrix} = 2 - k = 0.
Check k=2k = 2: trace 5, det 6−2=46 - 2 = 4, λ2−5λ+4=(λ−1)(λ−4)\lambda^2 - 5\lambda + 4 = (\lambda - 1)(\lambda - 4).

Verification: The characteristic polynomials were computed exactly in node with the Faddeev-LeVerrier method using BigInt fractions, and their roots were found with their multiplicities by a rational-root scan; each eigenspace basis was obtained as the exact null space of (A − λI), and the equality Av = λv was checked by multiplication for every eigenvector (verify/MAT1041/13.mjs).

Review card

Av = λv, v ≠ 0 → λ is an eigenvalue, v an eigenvector. (v = 0 does NOT count; λ = 0 is allowed.) Characteristic equation: det(λI - A) = 0 (or det(A - λI) = 0, same roots). 2×2 shortcut: λ² - (tr A)λ + det A = 0. Triangular/diagonal → eigenvalues on the diagonal. Eigenspace E_λ = null(A - λI): solve (A - λI)x = 0, write a basis with the free variables. (A - λI) must be SINGULAR; if you get only the trivial solution, you found λ wrongly. Algebraic multiplicity = multiplicity of the root. Geometric = dim E_λ. 1 ≤ geom ≤ algebraic. Σλ = tr A, Πλ = det A. λ = 0 is an eigenvalue ⇔ A is singular. A^k → λ^k; A⁻¹ → 1/λ; A + cI → λ + c; p(A) → p(λ). The eigenvector is the SAME. A and A^T have the same eigenvalues. Similar matrices have the same eigenvalues. There may be no real root (rotation: λ² + 1). For 3×3, expand det(λI - A) along the row with the most zeros, keep the factor.