← mapVector Spaces

Eigenvalues & Eigenvectors

⚗ Dr. Möbius, from the lab

Most directions, when you hit them with a matrix, get knocked off course — rotated, sheared, twisted into something pointing elsewhere. But a precious few directions refuse to turn. The matrix can only stretch them. Those unbreakable directions are the eigenvectors, and they are the transformation's true personality, stripped of all the bullshit. This is the fucking lesson the whole course has been a runway for.

THE BIG IDEA

An eigenvector of A is a nonzero direction v with Av = λv — A merely scales it by the eigenvalue λ — and you find the λ by solving det(A − λI) = 0.

The directions that refuse to turn

A matrix AA is a machine that eats a vector and spits out another. For a typical input vv, the output AvAv points somewhere new. But ask: are there directions vv that AA leaves pointing the same way, only longer or shorter? A direction vv where the output is just a scaled copy of the input:

Av=λv,v0.A v = \lambda v, \qquad v \ne 0.

When this holds, vv is an eigenvector of AA and the scalar λ\lambda is its eigenvalue ("eigen" is German for "own" — these are the matrix's own directions). The condition v0v \ne 0 is non-negotiable: A0=λ0A \cdot 0 = \lambda \cdot 0 holds trivially for every λ\lambda, so the zero vector would tell us nothing. Eigenvectors are nonzero by decree.

Geometrically, eigenvectors are the axes the transformation is built around. Drag a direction and watch where AA sends it — the eigendirections are the ones where the output arrow stays collinear with the input, glowing when you hit them:

eigen lab — hunt the special directions
vAv
A = [2 1 | 1 2]λ = 3, 1

These special directions are the map's skeleton — its actual identity, buried under all that matrix arithmetic. Everything else AA does is just a blend of stretching along them.

Finding the eigenvalues: det(AλI)=0\det(A - \lambda I) = 0

How do we hunt down the λ\lambda? Rearrange the defining equation. Av=λvAv = \lambda v becomes Avλv=0Av - \lambda v = 0, and writing λv=λIv\lambda v = \lambda I v (so we can factor a matrix out):

(AλI)v=0.(A - \lambda I)\,v = 0.

Now think about what this demands. We need a nonzero vv killed by the matrix AλIA - \lambda I. A matrix has a nonzero vector in its kernel exactly when it is not invertible — when it squashes some direction to zero. And from Determinants, a matrix is non-invertible precisely when its determinant is zero. So:

Av=λv for some v0    (AλI) is singular    det(AλI)=0.Av = \lambda v \text{ for some } v \ne 0 \iff (A - \lambda I) \text{ is singular} \iff \det(A - \lambda I) = 0.

That's the whole derivation, and I'll be damned if it isn't beautiful: the determinant is a squash-detector, and an eigenvalue is exactly a λ\lambda that makes AλIA - \lambda I squash. The expression

p(λ)=det(AλI)p(\lambda) = \det(A - \lambda I)

is the characteristic polynomial. For a 2×22\times 2 matrix A=(abcd)A = \begin{pmatrix} a & b \\ c & d \end{pmatrix},

p(λ)=det(aλbcdλ)=(aλ)(dλ)bc=λ2(a+d)λ+(adbc).p(\lambda) = \det\begin{pmatrix} a - \lambda & b \\ c & d - \lambda \end{pmatrix} = (a-\lambda)(d-\lambda) - bc = \lambda^2 - (a+d)\lambda + (ad - bc).

Notice the coefficients are old friends: a+d=trace(A)a + d = \operatorname{trace}(A) and adbc=det(A)ad - bc = \det(A). So p(λ)=λ2(trace)λ+(det)p(\lambda) = \lambda^2 - (\operatorname{trace})\lambda + (\det). Set it to zero, solve the quadratic, and out fall the eigenvalues.

Finding the eigenvectors: the null space of AλIA - \lambda I

Once you have an eigenvalue λ\lambda, plug it back in and find the nonzero vv with (AλI)v=0(A - \lambda I)v = 0. That's just finding the kernel (null space) of the matrix AλIA - \lambda I — a homogeneous system you already know how to solve from Gaussian Elimination.

The full set of solutions, ker(AλI)\ker(A - \lambda I), is called the eigenspace of λ\lambda. It's a subspace (kernels always are — proved back in Linear Maps), so it's closed under scaling and addition: if vv is an eigenvector, so is 5v5v and so is v-v. Eigenvectors come in entire lines (or planes), never as lone points. By convention we report one clean representative — pick any nonzero specimen from the eigenspace and you're done.

When there are no real eigenvectors: rotation

Here's the twist that exposes what eigenvalues really mean, and it's the kind of thing that should shake you. Take the 9090^\circ rotation

R=(0110).R = \begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix}.

Its characteristic polynomial is λ2(0)λ+(1)=λ2+1\lambda^2 - (0)\lambda + (1) = \lambda^2 + 1. Set it to zero: λ2=1\lambda^2 = -1. No real solutions. And of course not — geometrically, a 9090^\circ rotation turns every direction off its line. There is no direction it merely stretches, because it rotates everything. The algebra (no real roots) and the geometry (nothing stays put) are saying the identical thing. This is why "no real eigenvalues" isn't a failure — it's a rotation confessing its nature. (Over the complex numbers λ=±i\lambda = \pm i appear, which is exactly how rotation hides inside the imaginary unit, but real eigenvectors genuinely don't exist here.)

Distinct eigenvalues force independence

One structural fact we'll lean on hard in the next two lessons. Eigenvectors belonging to distinct eigenvalues are linearly independent. Let me prove it for two.

Suppose Av1=λ1v1Av_1 = \lambda_1 v_1 and Av2=λ2v2Av_2 = \lambda_2 v_2 with λ1λ2\lambda_1 \ne \lambda_2, both eigenvectors nonzero. Assume a dependence:

c_1 v_1 + c_2 v_2 = 0. \tag{1}

Apply AA to both sides: c1λ1v1+c2λ2v2=0c_1 \lambda_1 v_1 + c_2 \lambda_2 v_2 = 0. Now also multiply (1) by λ1\lambda_1: c1λ1v1+c2λ1v2=0c_1 \lambda_1 v_1 + c_2 \lambda_1 v_2 = 0. Subtract:

c2(λ2λ1)v2=0.c_2(\lambda_2 - \lambda_1) v_2 = 0.

Since λ2λ10\lambda_2 - \lambda_1 \ne 0 and v20v_2 \ne 0, we're forced to c2=0c_2 = 0. Plug back into (1): c1v1=0c_1 v_1 = 0, and v10v_1 \ne 0 forces c1=0c_1 = 0 too. The only dependence is the trivial one — they're independent. \blacksquare (The general statement: eigenvectors from kk distinct eigenvalues are always independent; same trick, induction.)

This is the seed of the next lesson. Distinct eigenvalues hand you a basis of eigenvectors for free, and a basis of eigenvectors is the magic grid in which the matrix goes diagonal. File it away — it fucking detonates immediately. Go do the gauntlet, you beautiful disaster.

🔬 SPECIMENS (worked examples)

Worked example 1 — a clean symmetric specimen from the lab

Find the eigenvalues and an eigenvector for each, for A=(2112)A = \begin{pmatrix} 2 & 1 \\ 1 & 2 \end{pmatrix}.

Characteristic polynomial. trace=4\operatorname{trace} = 4, det=41=3\det = 4 - 1 = 3, so

p(λ)=λ24λ+3=(λ1)(λ3).p(\lambda) = \lambda^2 - 4\lambda + 3 = (\lambda - 1)(\lambda - 3).

Eigenvalues: λ=1\lambda = 1 and λ=3\lambda = 3.

Eigenvector for λ=3\lambda = 3. Form A3I=(1111)A - 3I = \begin{pmatrix} -1 & 1 \\ 1 & -1 \end{pmatrix}. Solve (A3I)v=0(A-3I)v = 0: the top row gives x+y=0-x + y = 0, i.e. y=xy = x. Pick x=1x = 1:

v1=(11).v_1 = \begin{pmatrix} 1 \\ 1 \end{pmatrix}.

Check: Av1=(2+11+2)=(33)=3v1.A v_1 = \begin{pmatrix} 2+1 \\ 1+2 \end{pmatrix} = \begin{pmatrix} 3 \\ 3 \end{pmatrix} = 3 v_1. \checkmark

Eigenvector for λ=1\lambda = 1. Form AI=(1111)A - I = \begin{pmatrix} 1 & 1 \\ 1 & 1 \end{pmatrix}. The equation is x+y=0x + y = 0, i.e. y=xy = -x. Pick x=1x = 1:

v2=(11).v_2 = \begin{pmatrix} 1 \\ -1 \end{pmatrix}.

Check: Av2=(2112)=(11)=1v2.A v_2 = \begin{pmatrix} 2 - 1 \\ 1 - 2 \end{pmatrix} = \begin{pmatrix} 1 \\ -1 \end{pmatrix} = 1 \cdot v_2. \checkmark Note v1v2v_1 \perp v_2 — no accident, since AA is symmetric. That's a coming attraction.

Worked example 2 — a non-symmetric matrix, full run

Find the eigenvalues and an eigenvector for each, for A=(4123)A = \begin{pmatrix} 4 & 1 \\ 2 & 3 \end{pmatrix}.

Characteristic polynomial. trace=7\operatorname{trace} = 7, det=(4)(3)(1)(2)=122=10\det = (4)(3) - (1)(2) = 12 - 2 = 10, so

p(λ)=λ27λ+10=(λ2)(λ5).p(\lambda) = \lambda^2 - 7\lambda + 10 = (\lambda - 2)(\lambda - 5).

Eigenvalues: λ=2\lambda = 2 and λ=5\lambda = 5.

Eigenvector for λ=5\lambda = 5. A5I=(1122)A - 5I = \begin{pmatrix} -1 & 1 \\ 2 & -2 \end{pmatrix}. Top row: x+y=0y=x-x + y = 0 \Rightarrow y = x. Pick x=1x = 1:

v1=(11),Av1=(55)=5v1.v_1 = \begin{pmatrix} 1 \\ 1 \end{pmatrix}, \qquad Av_1 = \begin{pmatrix} 5 \\ 5 \end{pmatrix} = 5 v_1. \checkmark

Eigenvector for λ=2\lambda = 2. A2I=(2121)A - 2I = \begin{pmatrix} 2 & 1 \\ 2 & 1 \end{pmatrix}. Row: 2x+y=0y=2x2x + y = 0 \Rightarrow y = -2x. Pick x=1x = 1:

v2=(12),Av2=(4226)=(24)=2v2.v_2 = \begin{pmatrix} 1 \\ -2 \end{pmatrix}, \qquad Av_2 = \begin{pmatrix} 4 - 2 \\ 2 - 6 \end{pmatrix} = \begin{pmatrix} 2 \\ -4 \end{pmatrix} = 2 v_2. \checkmark

Two distinct eigenvalues, so v1v_1 and v2v_2 are independent (the theorem we just proved) — they form a basis of R2\mathbb{R}^2. Hold that thought.

Worked example 3 — the trap: a rotation with nothing to grab

Find the real eigenvalues of the 9090^\circ rotation R=(0110)R = \begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix}, and say what the result means geometrically.

Characteristic polynomial. trace=0\operatorname{trace} = 0, det=(0)(0)(1)(1)=1\det = (0)(0) - (-1)(1) = 1, so

p(λ)=λ20λ+1=λ2+1.p(\lambda) = \lambda^2 - 0\cdot\lambda + 1 = \lambda^2 + 1.

Set λ2+1=0λ2=1\lambda^2 + 1 = 0 \Rightarrow \lambda^2 = -1. No real solutions (the discriminant is 024(1)(1)=4<00^2 - 4(1)(1) = -4 < 0).

Meaning. A 9090^\circ rotation turns every vector off its own line — there is no direction it merely stretches, so there's no real eigenvector, and the algebra agrees by refusing to produce a real λ\lambda. This is the trap people fall into: they assume every matrix has real eigenvectors and start fishing for one that isn't there. The absence IS the information — it certifies that RR is a genuine rotation, mixing all directions, never fixing a line. (Complex eigenvalues ±i\pm i do exist, which is how rotations and the imaginary unit secretly shake hands — but over R\mathbb{R}, nothing stays put.)

☠ KNOWN HAZARDS

  • Letting v=0v = 0 count as an eigenvector. It never does — not in my lab, not anywhere. A0=λ0A\cdot 0 = \lambda \cdot 0 for every λ\lambda, so allowing it would make every scalar an eigenvalue. Eigenvectors are nonzero by definition; the eigenvalue, however, is perfectly allowed to be 00.

  • Solving det(AλI)=0\det(A - \lambda I) = 0 but forgetting the eigenvectors. Eigenvalues are only half the answer. For each λ\lambda you must still solve (AλI)v=0(A - \lambda I)v = 0 to get the directions. The λ\lambda tells you how much; the vv tells you which way.

  • Computing det(A)λI\det(A) - \lambda I instead of det(AλI)\det(A - \lambda I). This one makes me want to scream. You must subtract λ\lambda from the diagonal entries first, then take the determinant of the whole modified matrix. λI\lambda I goes inside the determinant.

  • Panicking when a real matrix has no real eigenvalues. That's not an error — it's geometric content (a rotation). Check your characteristic polynomial's discriminant; a negative discriminant is the algebra telling you "this map rotates."

TL;DR

  • Eigenvector: nonzero vv with Av=λvAv = \lambda v; eigenvalue: the scale factor λ\lambda. The matrix merely stretches its eigenvectors.

  • Eigenvalues solve det(AλI)=0\det(A - \lambda I) = 0, because AλIA - \lambda I must be singular (squash some nonzero vv) for an eigenvector to exist. For 2×22\times2: λ2(trace)λ+det=0\lambda^2 - (\operatorname{trace})\lambda + \det = 0.

  • Eigenvectors for λ\lambda are the nonzero vectors in ker(AλI)\ker(A - \lambda I) — the eigenspace, a subspace.

  • A real rotation has no real eigenvalues (λ2+1=0\lambda^2 + 1 = 0): it turns every direction, leaving none merely stretched.

  • Eigenvectors for distinct eigenvalues are linearly independent — the seed of diagonalization.

unlocks