C10 - Eigenvalues and eigenvectors

C10 - Eigenvalues and eigenvectors

A matrix usually changes both the length and direction of a vector. Eigenvectors identify the special directions in which it acts as a simple scale. You will find eigenvalues from a characteristic equation, find their corresponding vectors in two and three dimensions, and interpret the results geometrically.

Characteristic equation in 2D and 3D

Which directions does A=(3113)\mathbf{A}=\begin{pmatrix}3&1\\1&3\end{pmatrix} stretch without turning them away from their original line? Try two vectors:

A(11)=(44),A(1−1)=(2−2)\mathbf{A}\begin{pmatrix}1\\1\end{pmatrix}=\begin{pmatrix}4\\4\end{pmatrix},\qquad \mathbf{A}\begin{pmatrix}1\\-1\end{pmatrix}=\begin{pmatrix}2\\-2\end{pmatrix}

Each output is a scalar multiple of its input. The scale depends on the direction: 4 along one line and 2 along the other.

Eigenvalue and eigenvector

A nonzero vector v\mathbf{v} is an eigenvector of a square matrix M\mathbf{M} if

Mv=λv\mathbf{M}\mathbf{v}=\lambda\mathbf{v}

The scalar λ\lambda is its corresponding eigenvalue. The zero vector is excluded because M0=λ0\mathbf{M}\mathbf{0}=\lambda\mathbf{0} holds for every scalar and identifies no special direction. An eigenvalue itself may be zero.

How can you find the scales without guessing vectors? Rearrange to (M−λI)v=0(\mathbf{M}-\lambda\mathbf{I})\mathbf{v}=\mathbf{0}. If M−λI\mathbf{M}-\lambda\mathbf{I} were invertible, multiplying by its inverse would force v=0\mathbf{v}=\mathbf{0}. A nonzero solution therefore requires this matrix to be singular. Conversely, a singular square matrix has a nonzero solution to its homogeneous equations.

Characteristic equation

The eigenvalues of a square matrix M\mathbf{M} are the roots of

det⁡(M−λI)=0\det(\mathbf{M}-\lambda\mathbf{I})=0

This is the characteristic equation. Subtract λ\lambda from each diagonal entry only. For a 2×22\times2 matrix it is quadratic; for a 3×33\times3 matrix it is cubic. Using det⁡(λI−M)=0\det(\lambda\mathbf{I}-\mathbf{M})=0 gives the same roots.

Not in the formula booklet: learn it.

Use the determinant methods you already know. A row with zeros can make the three-dimensional calculation much shorter.

Worked example 1: two-dimensional eigenvalues

Question. Find the eigenvalues of A=(3113)\mathbf{A}=\begin{pmatrix}3&1\\1&3\end{pmatrix}.

Plan. Form the determinant, expand it and solve the resulting quadratic.

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

Hence (λ−4)(λ−2)=0(\lambda-4)(\lambda-2)=0.

Answer. The eigenvalues are 44 and 22.

Check. The two direct vector products at the start give scale factors 4 and 2, independently confirming both roots.

Keep that two-dimensional block and add a third coordinate scaled by −1-1. The extra direction adds one factor to the determinant.

Worked example 2: extending the same block to 3D

Question. Find the eigenvalues of B=(31013000−1)\mathbf{B}=\begin{pmatrix}3&1&0\\1&3&0\\0&0&-1\end{pmatrix}.

Plan. Expand along the last row, where two entries are zero.

det⁡(B−λI)=(−1−λ)∣3−λ113−λ∣=(−1−λ)(λ−4)(λ−2)\begin{aligned} \det(\mathbf{B}-\lambda\mathbf{I}) &=(-1-\lambda)\begin{vmatrix}3-\lambda&1\\1&3-\lambda\end{vmatrix}\\ &=(-1-\lambda)(\lambda-4)(\lambda-2) \end{aligned}

Set the complete product equal to zero.

Answer. The eigenvalues are 44, 22 and −1-1.

Check. Direct multiplication gives scales 4, 2 and −1-1 on (1,1,0)T(1,1,0)^{\mathsf T}, (1,−1,0)T(1,-1,0)^{\mathsf T} and (0,0,1)T(0,0,1)^{\mathsf T}, respectively. All three vectors are nonzero.

Real matrices need not have only real eigenvalues. For example, a quarter-turn matrix R=(0−110)\mathbf{R}=\begin{pmatrix}0&-1\\1&0\end{pmatrix} has characteristic equation λ2+1=0\lambda^2+1=0, with roots ±i\pm\mathrm{i}. These are valid complex eigenvalues, but they do not describe real directions in the plane. For a real matrix, non-real roots occur in conjugate pairs because its characteristic polynomial has real coefficients. Geometrical directions below use real eigenvectors and real eigenvalues.

Watch out

Wrong: subtracting λ\lambda from every entry. The expression is M−λI\mathbf{M}-\lambda\mathbf{I}, so off-diagonal entries stay unchanged. If M=cN\mathbf{M}=c\mathbf{N}, use det⁡(cN−λI)=0\det(c\mathbf{N}-\lambda\mathbf{I})=0; dropping the factor cc changes the eigenvalues.

Find the characteristic equation for a fresh matrix.

Finding corresponding eigenspaces

An eigenvalue tells you the scale; solving (M−λI)v=0(\mathbf{M}-\lambda\mathbf{I})\mathbf{v}=\mathbf{0} tells you which vectors receive that scale. These are homogeneous equations, so choose a convenient free variable rather than seeking a unique vector.

Eigenspace

For an eigenvalue λ\lambda, its eigenspace is the set of all solutions of (M−λI)v=0(\mathbf{M}-\lambda\mathbf{I})\mathbf{v}=\mathbf{0}, including 0\mathbf{0}. Its eigenvectors are exactly the nonzero vectors in that set. Every nonzero scalar multiple of an eigenvector is another eigenvector for the same eigenvalue.

The equations are dependent when λ\lambda is an eigenvalue. An equation reducing to 0=00=0 adds no restriction; it does not require a free variable to be zero.

Worked example 3: matching each vector with its scale

Question. Find the eigenspaces of A=(3113)\mathbf{A}=\begin{pmatrix}3&1\\1&3\end{pmatrix}, whose eigenvalues are 4 and 2.

Plan. Substitute one eigenvalue at a time and solve for x,yx,y.

λ=4:(−111−1)(xy)=0 ⟹ y=xλ=2:(1111)(xy)=0 ⟹ y=−x\begin{aligned} \lambda=4:\quad \begin{pmatrix}-1&1\\1&-1\end{pmatrix}\begin{pmatrix}x\\y\end{pmatrix}=\mathbf{0} &\ \Longrightarrow\ y=x\\ \lambda=2:\quad \begin{pmatrix}1&1\\1&1\end{pmatrix}\begin{pmatrix}x\\y\end{pmatrix}=\mathbf{0} &\ \Longrightarrow\ y=-x \end{aligned}

The second row repeats the restriction from the first in each case.

Answer. The eigenspaces are {t(1,1)T:t∈R}\{t(1,1)^{\mathsf T}:t\in\mathbb{R}\} for 4 and {t(1,−1)T:t∈R}\{t(1,-1)^{\mathsf T}:t\in\mathbb{R}\} for 2. Restrict to t≠0t\ne0 when listing eigenvectors.

Check. A(t,t)T=4(t,t)T\mathbf{A}(t,t)^{\mathsf T}=4(t,t)^{\mathsf T} and A(t,−t)T=2(t,−t)T\mathbf{A}(t,-t)^{\mathsf T}=2(t,-t)^{\mathsf T} for every real tt, checking the whole sets rather than just one representative.

Use the same eigenvalues and vectors in the first two coordinates, but now retain all three equations. The extra coordinate cannot simply be left arbitrary.

Worked example 4: the corresponding three-dimensional spaces

Question. Find the eigenspaces of B=(31013000−1)\mathbf{B}=\begin{pmatrix}3&1&0\\1&3&0\\0&0&-1\end{pmatrix}, whose eigenvalues are 4, 2 and −1-1.

Plan. For each eigenvalue, solve (B−λI)(x,y,z)T=0(\mathbf{B}-\lambda\mathbf{I})(x,y,z)^{\mathsf T}=\mathbf{0}.

λ=4:−x+y=0,x−y=0,−5z=0⟹(x,y,z)=(t,t,0)λ=2:x+y=0,x+y=0,−3z=0⟹(x,y,z)=(t,−t,0)λ=−1:4x+y=0,x+4y=0,0=0⟹(x,y,z)=(0,0,t)\begin{aligned} \lambda=4:&\quad -x+y=0,\quad x-y=0,\quad -5z=0 &&\Longrightarrow (x,y,z)=(t,t,0)\\ \lambda=2:&\quad x+y=0,\quad x+y=0,\quad -3z=0 &&\Longrightarrow (x,y,z)=(t,-t,0)\\ \lambda=-1:&\quad 4x+y=0,\quad x+4y=0,\quad 0=0 &&\Longrightarrow (x,y,z)=(0,0,t) \end{aligned}

For the last case, y=−4xy=-4x in x+4y=0x+4y=0 gives −15x=0-15x=0; only zz is free.

Answer. The eigenspaces are all real multiples of (1,1,0)T(1,1,0)^{\mathsf T}, (1,−1,0)T(1,-1,0)^{\mathsf T} and (0,0,1)T(0,0,1)^{\mathsf T}, paired with 4, 2 and −1-1, respectively. Nonzero multiples are eigenvectors.

Check. Substituting each general solution into Bv\mathbf{B}\mathbf{v} gives (4t,4t,0)T(4t,4t,0)^{\mathsf T}, (2t,−2t,0)T(2t,-2t,0)^{\mathsf T} and (0,0,−t)T(0,0,-t)^{\mathsf T}, the required scalar multiples.

An eigenspace can be a plane. For diag⁡(2,2,5)\operatorname{diag}(2,2,5), the repeated root 2 gives only z=0z=0, so every (s,t,0)T(s,t,0)^{\mathsf T} lies in its eigenspace. All except (0,0,0)T(0,0,0)^{\mathsf T} are eigenvectors. The two free variables provide two independent directions. A repeated root does not automatically provide that many independent eigenvectors: you must solve the equations.

If complex vectors are allowed, the quarter-turn matrix from part 1 has eigenvectors too. Its equations give −y=ix-y=\mathrm{i}x for λ=i\lambda=\mathrm{i}, so (1,−i)T(1,-\mathrm{i})^{\mathsf T} works; conjugating gives (1,i)T(1,\mathrm{i})^{\mathsf T} for −i-\mathrm{i}. Directly, R(1,−i)T=(i,1)T=i(1,−i)T\mathbf{R}(1,-\mathrm{i})^{\mathsf T}=(\mathrm{i},1)^{\mathsf T}=\mathrm{i}(1,-\mathrm{i})^{\mathsf T}. These are complex vectors, not arrows representing real points of the plane.

Watch out

From y=3xy=3x, a convenient eigenvector is (1,3)T(1,3)^{\mathsf T}, not (3,1)T(3,1)^{\mathsf T}. Read the coordinate ratio in the right order, keep the eigenvalue paired with its vector, and check Mv=λv\mathbf{M}\mathbf{v}=\lambda\mathbf{v}.

Find a two-dimensional eigenspace inside three-dimensional space.

Geometrical meaning of eigenvalues

For a real eigenvector, the equation Mv=λv\mathbf{M}\mathbf{v}=\lambda\mathbf{v} describes motion along its original line. Length is multiplied by ∣λ∣|\lambda|; the sign tells you what happens to the direction.

Real eigenvalueEffect on a nonzero eigenvector
λ>0\lambda>0Same ray, with length multiplied by λ\lambda
λ<0\lambda<0Opposite ray, with length multiplied by $
λ=0\lambda=0Collapses to the zero vector
λ=1\lambda=1Endpoint stays fixed
λ=−1\lambda=-1Length stays the same and direction reverses

For instance, H=(1221)\mathbf{H}=\begin{pmatrix}1&2\\2&1\end{pmatrix} sends (1,1)T(1,1)^{\mathsf T} to (3,3)T(3,3)^{\mathsf T}, but sends (1,−1)T(1,-1)^{\mathsf T} to (−1,1)T(-1,1)^{\mathsf T}. The same transformation stretches one direction and reverses the other without changing its length. Eigenvalues reveal that mixture far more clearly than the four entries alone.

[DIAGRAM: asset_name: Two directions under one matrix; asset_slug: aqa7367-eigen-directions; caption: One matrix stretches vectors on y equals x by three and reverses vectors on y equals minus x; description: Two equal-scale coordinate panels show the matrix with rows one two and two one. In the left panel a blue arrow from the origin to (1,1) and a clay arrow to (3,3) lie on the dashed line y equals x. In the right panel a blue arrow to (1,minus 1) and a clay arrow to (minus 1,1) lie on the dashed line y equals minus x. Labels identify the original and image and eigenvalues three and minus one.]
One matrix stretches vectors on y equals x by three and reverses vectors on y equals minus x

Linearity extends this description to every point tvt\mathbf{v}: its image is λtv\lambda t\mathbf{v}. If λ≠0\lambda\ne0, the whole line is mapped onto itself. When λ=0\lambda=0, the line collapses to the origin instead, so its image is contained in the original line but is not the whole line. Eigenvalue 1 gives a line of invariant points; another nonzero eigenvalue gives an invariant line on which nonzero points move.

Not every vector is an eigenvector. For the same H\mathbf{H}, (1,0)T↦(1,2)T(1,0)^{\mathsf T}\mapsto(1,2)^{\mathsf T}, which is not a scalar multiple of (1,0)T(1,0)^{\mathsf T}. A quarter-turn has no real eigenvector at all: every nonzero real arrow is turned off its original line. Its complex eigenvalues are consistent with that geometry.

The opening matrix A\mathbf{A} has eigenvectors along y=xy=x and y=−xy=-x. These directions receive scales 4 and 2. Any nonzero multiple of either displayed vector therefore answers the opening question.

Key idea

A real eigenvector identifies a line direction, and its real eigenvalue gives the signed scale along it. Negative eigenvalues reverse the ray; zero collapses the direction; only eigenvalue 1 fixes every point along it.

Distinguish a fixed direction from fixed points.

Quick Recap

  • An eigenvector is nonzero and satisfies Mv=λv\mathbf{M}\mathbf{v}=\lambda\mathbf{v}; an eigenvalue may be zero.
  • Solve det⁡(M−λI)=0\det(\mathbf{M}-\lambda\mathbf{I})=0, then solve the homogeneous equations separately for each root.
  • An eigenspace includes zero; its nonzero members are the corresponding eigenvectors.
  • For real vectors, ∣λ∣|\lambda| gives the length scale and a negative sign reverses direction. Complex eigenvectors do not represent real line directions.

The opening matrix stretches directions (1,1)T(1,1)^{\mathsf T} and (1,−1)T(1,-1)^{\mathsf T} by 4 and 2, respectively, without turning them off their original lines.