Skip to main content
© 2026 ePowerAI. All rights reserved.
CoursesContact
eAI
CoursesContact
Week 9: Eigenvalues and Eigenvectors
Linear Algebra
01Week 1: Vectors and Linear Combinations
02Week 2: Linear Transformations and Matrices
03Week 3: Elimination and LU Factorization
04Week 4: Determinants
05Week 5: Vector Spaces, Independence, and Basis
06Week 6: The Four Fundamental Subspaces
07Week 7: Orthogonality and Projections
08Week 8: Least Squares and QR
09Week 9: Eigenvalues and Eigenvectors
10Week 10: Diagonalization and Markov Matrices
11Week 11: Differential Equations
12Week 12: Symmetric and Positive Definite Matrices
13Week 13: The SVD and Complex Matrices
14Week 14: The Fourier Matrix, FFT, and PCA
Week 09· Linear Algebra20 min read

Week 9: Eigenvalues and Eigenvectors

✦Learning Outcomes
  • Define eigenvalues and eigenvectors and set up the characteristic equation
  • Compute eigenvalues via det(A − λI) = 0
  • Find eigenvectors for a given eigenvalue
  • Interpret eigenvectors as directions preserved by a transformation
◆Prerequisites

Background knowledge assumed:

  • Determinants
  • Linear transformations

Recommended: Review Week 4: Determinants and Week 2: Linear Transformations and Matrices before starting.

Eigenvalues and Eigenvectors Defined#

Week 2 treated a square matrix AAA as a linear transformation: every input vector xxx is sent to a new vector AxAxAx. In general AxAxAx points somewhere new — length and direction both change. There is a special class of directions where the transformation is as simple as possible: AAA only scales the vector, without rotating it off its line.

A nonzero vector xxx is an eigenvector of AAA with eigenvalue λ\lambdaλ when

Ax=λx.Ax = \lambda x.Ax=λx.

The left side applies the matrix; the right side multiplies by a scalar. Equality means AxAxAx lies on the same line through the origin as xxx. The number λ\lambdaλ is the stretch factor along that line: λ>1\lambda > 1λ>1 stretches, 0<λ<10 < \lambda < 10<λ<1 shrinks, λ<0\lambda < 0λ<0 flips the direction and scales, and λ=0\lambda = 0λ=0 sends xxx to the zero vector (so xxx is in the null space N(A)N(A)N(A)).

The zero vector is excluded by definition: A0=λ0A0 = \lambda 0A0=λ0 holds for every λ\lambdaλ, so it would not pin down a meaningful direction or eigenvalue.

Geometric picture. Think of AAA as a map of the plane. Most arrows xxx get sent to arrows AxAxAx that are not parallel to xxx. Eigenvectors are the rare arrows that stay parallel: the map stretches or flips them along their own axis. Those invariant axes are the skeleton of the transformation — later weeks (diagonalization, symmetric matrices, SVD) rebuild AAA from them.

Quick check. Take

A=(2003),x=(10).A = \begin{pmatrix} 2 & 0 \\ 0 & 3 \end{pmatrix}, \qquad x = \begin{pmatrix} 1 \\ 0 \end{pmatrix}.A=(20​03​),x=(10​).

Then Ax=(2,0)⊤=2xAx = (2, 0)^\top = 2xAx=(2,0)⊤=2x, so xxx is an eigenvector with λ=2\lambda = 2λ=2. The vector (0,1)⊤(0,1)^\top(0,1)⊤ is an eigenvector with λ=3\lambda = 3λ=3. Every axis-aligned stretch has the standard basis as eigenvectors; the interesting case is when AAA is not diagonal and those special directions are tilted.


The Characteristic Equation#

Rewrite Ax=λxAx = \lambda xAx=λx by moving every term to one side:

Ax−λx=0⟺(A−λI)x=0.Ax - \lambda x = 0 \quad\Longleftrightarrow\quad (A - \lambda I)x = 0.Ax−λx=0⟺(A−λI)x=0.

So xxx is an eigenvector with eigenvalue λ\lambdaλ exactly when xxx is a nonzero vector in the null space of A−λIA - \lambda IA−λI:

x∈N(A−λI),x≠0.x \in N(A - \lambda I),\qquad x \neq 0.x∈N(A−λI),x=0.

Week 4: a square matrix BBB has a nonzero null-space vector if and only if det⁡(B)=0\det(B) = 0det(B)=0 (equivalently, BBB is singular). Apply that test with B=A−λIB = A - \lambda IB=A−λI:

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

This is the characteristic equation of AAA. Its roots are the eigenvalues. Expanding the determinant produces a monic polynomial of degree nnn in λ\lambdaλ — the characteristic polynomial pA(λ)=det⁡(A−λI)p_A(\lambda) = \det(A - \lambda I)pA​(λ)=det(A−λI) — so an n×nn \times nn×n matrix has nnn eigenvalues counting algebraic multiplicity (some may be complex even when AAA is real).

Why determinants reappear. The same singularity test that told you whether Ax=bAx = bAx=b had a unique solution now tells you which stretch factors λ\lambdaλ are compatible with AAA. No new machinery: only “when is A−λIA - \lambda IA−λI singular?”

Trace and determinant shortcuts (2×22 \times 22×2). For a 2×22 \times 22×2 matrix the characteristic polynomial is always

det⁡(A−λI)=λ2−tr⁡(A) λ+det⁡(A).\det(A - \lambda I) = \lambda^2 - \operatorname{tr}(A)\,\lambda + \det(A).det(A−λI)=λ2−tr(A)λ+det(A).

You can verify this by expanding det⁡(a−λbcd−λ)=(a−λ)(d−λ)−bc\det\begin{pmatrix} a-\lambda & b \\ c & d-\lambda \end{pmatrix} = (a-\lambda)(d-\lambda) - bcdet(a−λc​bd−λ​)=(a−λ)(d−λ)−bc. The Vieta relations follow immediately: the sum of the two eigenvalues equals the trace, and their product equals det⁡(A)\det(A)det(A). Those identities hold for every n×nn \times nn×n matrix as well (sum of eigenvalues = trace, product = determinant), but the quadratic formula is special to 2×22 \times 22×2.


Finding Eigenvectors#

Once you have a candidate λ\lambdaλ, finding eigenvectors is a null-space computation from Week 6: row-reduce A−λIA - \lambda IA−λI and solve (A−λI)x=0(A - \lambda I)x = 0(A−λI)x=0. Every nonzero solution is an eigenvector for that λ\lambdaλ. Scaling an eigenvector still yields an eigenvector for the same λ\lambdaλ — you usually report a convenient representative (integer components, or unit length).

Worked example end to end. Let

A=(4123).A = \begin{pmatrix} 4 & 1 \\ 2 & 3 \end{pmatrix}.A=(42​13​).

Step 1 — eigenvalues. Trace tr⁡(A)=4+3=7\operatorname{tr}(A) = 4 + 3 = 7tr(A)=4+3=7 and det⁡(A)=12−2=10\det(A) = 12 - 2 = 10det(A)=12−2=10, so

det⁡(A−λI)=λ2−7λ+10=0.\det(A - \lambda I) = \lambda^2 - 7\lambda + 10 = 0.det(A−λI)=λ2−7λ+10=0.

Factor: (λ−5)(λ−2)=0(\lambda - 5)(\lambda - 2) = 0(λ−5)(λ−2)=0. The eigenvalues are λ1=5\lambda_1 = 5λ1​=5 and λ2=2\lambda_2 = 2λ2​=2.

Step 2 — eigenvector for λ=5\lambda = 5λ=5. Form

A−5I=(4−5123−5)=(−112−2).A - 5I = \begin{pmatrix} 4-5 & 1 \\ 2 & 3-5 \end{pmatrix} = \begin{pmatrix} -1 & 1 \\ 2 & -2 \end{pmatrix}.A−5I=(4−52​13−5​)=(−12​1−2​).

Row 2 is −2-2−2 times row 1, so the rank is 111. The equation (−1)x1+x2=0(-1)x_1 + x_2 = 0(−1)x1​+x2​=0 gives x2=x1x_2 = x_1x2​=x1​. Free variable x1=tx_1 = tx1​=t yields

x=t(11),t≠0.x = t\begin{pmatrix} 1 \\ 1 \end{pmatrix}, \qquad t \neq 0.x=t(11​),t=0.

Take v1=(1,1)⊤v_1 = (1, 1)^\topv1​=(1,1)⊤. Check: Av1=(5,5)⊤=5v1A v_1 = (5, 5)^\top = 5 v_1Av1​=(5,5)⊤=5v1​.

Step 3 — eigenvector for λ=2\lambda = 2λ=2. Form

A−2I=(2121).A - 2I = \begin{pmatrix} 2 & 1 \\ 2 & 1 \end{pmatrix}.A−2I=(22​11​).

The single independent equation is 2x1+x2=02x_1 + x_2 = 02x1​+x2​=0, so x2=−2x1x_2 = -2x_1x2​=−2x1​. Thus

v2=(1−2)v_2 = \begin{pmatrix} 1 \\ -2 \end{pmatrix}v2​=(1−2​)

is an eigenvector. Check: Av2=(4−2, 2−6)⊤=(2,−4)⊤=2v2A v_2 = (4 - 2,\ 2 - 6)^\top = (2, -4)^\top = 2 v_2Av2​=(4−2, 2−6)⊤=(2,−4)⊤=2v2​.

Summary. AAA stretches the line spanned by (1,1)⊤(1,1)^\top(1,1)⊤ by 555 and the line spanned by (1,−2)⊤(1,-2)^\top(1,−2)⊤ by 222. Those two directions are the eigenbasis of this plane map.

Recipe for n×nn \times nn×n.

  1. Form pA(λ)=det⁡(A−λI)p_A(\lambda) = \det(A - \lambda I)pA​(λ)=det(A−λI) and find the roots λ\lambdaλ.
  2. For each root, solve (A−λI)x=0(A - \lambda I)x = 0(A−λI)x=0 by elimination; report a basis of the solution space (nonzero vectors only).
  3. Optionally normalize eigenvectors or choose integer multiples for readability.

For larger nnn, hand expansion of det⁡(A−λI)\det(A - \lambda I)det(A−λI) is painful; production code uses iterative algorithms (QR iteration and variants). The conceptual pipeline — “characteristic roots, then null spaces” — stays the same. The browser lab below uses numpy.linalg.eig for the same 2×22 \times 22×2 matrix and verifies Av=λvAv = \lambda vAv=λv.

Exercise · Fill in the blank

Find the eigenvalues of A = [[4,1],[2,3]]. The characteristic equation is λ² − trace(A)λ + det(A) = 0. Enter the larger eigenvalue.


Eigenspaces and Multiplicity#

For a fixed eigenvalue λ\lambdaλ, the set of all solutions of (A−λI)x=0(A - \lambda I)x = 0(A−λI)x=0 is the null space N(A−λI)N(A - \lambda I)N(A−λI). That set always includes the zero vector, and it is a subspace of Rn\mathbb{R}^nRn. It is called the eigenspace of AAA for λ\lambdaλ:

Eλ=N(A−λI)={x:Ax=λx}.E_\lambda = N(A - \lambda I) = \{ x : Ax = \lambda x \}.Eλ​=N(A−λI)={x:Ax=λx}.

(The zero vector is allowed here because subspaces must contain zero; “eigenvectors” still means the nonzero members of EλE_\lambdaEλ​.)

The geometric multiplicity of λ\lambdaλ is dim⁡Eλ\dim E_\lambdadimEλ​ — how many independent eigenvectors you can find for that eigenvalue. The algebraic multiplicity is the multiplicity of λ\lambdaλ as a root of the characteristic polynomial. Always

1≤geometric multiplicity≤algebraic multiplicity.1 \le \text{geometric multiplicity} \le \text{algebraic multiplicity}.1≤geometric multiplicity≤algebraic multiplicity.

When the two multiplicities match for every eigenvalue and you can assemble a full basis of eigenvectors, AAA is diagonalizable — Week 10's main topic. They need not match.

Repeated root without enough eigenvectors. Consider

J=(1101).J = \begin{pmatrix} 1 & 1 \\ 0 & 1 \end{pmatrix}.J=(10​11​).

Then det⁡(J−λI)=(1−λ)2\det(J - \lambda I) = (1 - \lambda)^2det(J−λI)=(1−λ)2, so λ=1\lambda = 1λ=1 has algebraic multiplicity 222. But

J−I=(0100)J - I = \begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix}J−I=(00​10​)

has rank 111, so dim⁡N(J−I)=1\dim N(J - I) = 1dimN(J−I)=1. Every eigenvector is a multiple of (1,0)⊤(1, 0)^\top(1,0)⊤; there is no second independent eigenvector. Geometric multiplicity is 111, less than algebraic multiplicity 222. Matrices like JJJ are defective: they cannot be diagonalized. Week 10 returns to the exact criterion (a full set of independent eigenvectors) and to the Jordan form that replaces diagonalization when that criterion fails.

What you need this week. Compute eigenvalues from det⁡(A−λI)=0\det(A - \lambda I) = 0det(A−λI)=0, compute eigenspaces as null spaces, and read geometric multiplicity as the dimension of that null space. The diagonalizability question is reserved for next week; the example above only flags that a repeated root does not automatically give you two independent directions.


Knowledge check#

Question 1 of 3

Eigenvalues of A are the roots of which equation?

det(A - λI) = 0
det(A) = λ
trace(A) = λ
Ax = 0

Browser lab#

Compute eigenvalues/eigenvectors with numpy.linalg.eig and verify Av=λvAv = \lambda vAv=λv.

python · runs in browser
import numpy as np

A = np.array([[4., 1.],
              [2., 3.]])

eigvals, eigvecs = np.linalg.eig(A)
print("eigenvalues =", eigvals)
print("eigenvectors (columns) =\n", eigvecs)

for i in range(2):
    lhs = A @ eigvecs[:, i]
    rhs = eigvals[i] * eigvecs[:, i]
    print(f"A v{i+1} ≈ λ{i+1} v{i+1}:", np.allclose(lhs, rhs))
← Previous
Week 8: Least Squares and QR
Next →
Week 10: Diagonalization and Markov Matrices
On this page
  • Eigenvalues and Eigenvectors Defined
  • The Characteristic Equation
  • Finding Eigenvectors
  • Eigenspaces and Multiplicity
  • Knowledge check
  • Browser lab