Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

For an n x n square matrix A, we seek a scalar λ\lambda and vector x such that:

Ax=λx\begin{align} \vv{A} \vv{x} = \lambda \vv{x} \end{align}

We want a nontrivial solution (x0\vv{x} \ne \vv{0}). This pair represents a vector x that, when multiplied into A, does not change its direction. However, it may adopt a new magnitude λ\lambda. Borrowing from a German word for “own”, we call x an eigenvector and λ\lambda an eigenvalue of A.

To find the eigenvalues of A, rearrange:

Axλx=0(AλI)x=0\begin{align} \vv{A} \vv{x} - \lambda \vv{x} &= \vv{0} \\ (\vv{A} - \lambda \vv{I}) \vv{x} &= \vv{0} \end{align}

If AλI\vv{A} - \lambda \vv{I} is invertible, we have only x=0\vv{x} = \vv{0} as a solution. Hence, we require this matrix to be singular! This occurs when its determinant is zero:

AλI=0\begin{align} |\vv{A}- \lambda \vv{I}| = 0 \end{align}

This equation creates a characeristic polynomial of degree n for λ\lambda that can be solved. Then, the eigenvector x that corresponds to each root λ\lambda can be determined. For example, to find the eigenvalues of:

A=[5222]\vv{A} = \begin{bmatrix} -5 & 2\\ 2 & -2 \end{bmatrix}

First compute the determinant:

AλI=5λ222λ=(5λ)(2λ)22=λ2+7λ+6=(λ+1)(λ+6)=0\begin{align} |\vv{A} -\lambda \vv{I}| &= \begin{vmatrix} -5-\lambda & 2 \\ 2 & -2-\lambda \end{vmatrix} \\ &= (-5 -\lambda)(-2-\lambda) - 2 \cdot 2 \\ &= \lambda^2 +7\lambda + 6 \\ &= (\lambda +1)(\lambda +6) = 0 \end{align}

Hence, the eigenvalues of A are λ1=1\lambda_1 = -1 and λ2=6\lambda_2 = -6. Note that the particular ordering of the eigenvalues is not important, and we are only labeling them to make it convenient to refer to a particular eigenvalue later.

Next, we seek the eigenvector x1\vv{x}_1 that corresponds to λ1\lambda_1. This vector must solve the system:

(Aλ1I)x1=0[4221]x1=0\begin{align} (\vv{A} - \lambda_1 \vv{I}) \vv{x}_1 &= \vv{0} \\ \begin{bmatrix} -4 & 2 \\ 2 & -1 \end{bmatrix} \vv{x}_1 &= \vv{0} \end{align}

We can solve for x1\vv{x}_1 using row reduction. Since the last column of the augmented matrix would be only zeros, it is not necessary to include it, and we row reduce only the matrix itself:

[4221]R1+R2/2[4200]÷2R2[11/200]\begin{bmatrix} -4 & 2 \\ 2 & -1 \end{bmatrix} \begin{matrix}\vphantom{R_1} \\ +R_2/2 \end{matrix} \to \begin{bmatrix} -4 & 2 \\ 0 & 0 \end{bmatrix} \begin{matrix} \div -2 \\ \vphantom{R_2} \end{matrix} \to \begin{bmatrix} 1 & -1/2 \\ 0 & 0 \end{bmatrix}

Note that this reduced matrix is equivalent to x1=x2/2x_1 = x_2/2 (using x to represent elements of x1\vv{x}_1 for convenience) with x2x_2 being free. This means that the eigenvector is not unique, which makes sense both from the original equation we are trying to solve and us forcing the matrix to be singular! Equivalently, eigenvectors can be scaled up or down by an arbitrary nonzero constant. We may then choose x2x_2 so that the eigenvector has either nice values (like like integers) or a nice norm (1 is conventional). For example, choosing x2=2x_2 = 2 gives:

x1=[12]\vv{x}_1 = \begin{bmatrix} 1 \\ 2 \end{bmatrix}

Let’s use the same strategy to get the eigenvector for λ2=6\lambda_2 = -6. We’ll jump straight to the matrix we need:

(Aλ2I)=[1224]x2=[21]\begin{align} (\vv{A} - \lambda_2 \vv{I}) = \begin{bmatrix} 1 & 2 \\ 2 & 4 \end{bmatrix} \to \vv{x}_2 = \begin{bmatrix}2 \\ -1 \end{bmatrix} \end{align}

using the swapping trick on the first row.

Multiple eigenvalues

An n x n matrix has n eigenvalues, but they may not be distinct! For example, consider the matrix:

A=[223216120]\vv{A} = \begin{bmatrix} -2 & 2 & -3 \\ 2 & 1 & -6 \\ -1 & -2 & 0 \end{bmatrix}

We need to evaluate the determinant:

AλI=2λ2321λ612λ=1231λ6(2)2λ326λ2λ221λ=[12+3(1λ)]+2[6(2+λ)+6]λ[(λ+2)(λ1)4]=λ3λ2+21λ+45=0\begin{align} |\vv{A} - \lambda \vv{I}| &= \begin{vmatrix} -2-\lambda & 2 & -3 \\ 2 & 1-\lambda & -6 \\ -1 & -2 & -\lambda \end{vmatrix} \\ &= -1 \begin{vmatrix} 2 & -3 \\ 1-\lambda & -6 \end{vmatrix} -(-2) \begin{vmatrix} -2-\lambda & -3 \\ 2 & -6 \end{vmatrix} -\lambda \begin{vmatrix} -2-\lambda & 2 \\ 2 & 1-\lambda \end{vmatrix} \\ &= -[-12+3(1-\lambda)]+2[6(2+\lambda)+6]-\lambda[(\lambda+2)(\lambda-1)-4] \\ &= -\lambda^3-\lambda^2+21\lambda+45 = 0 \end{align}

The roots of this cubic polynomial are λ1=5\lambda_1 = 5 and λ2=λ3=3\lambda_2 = \lambda_3 = -3. The eigenvector for λ1=5\lambda_1 = 5 is obtained from

Aλ1I=[723246123][101012000]\vv{A} - \lambda_1 \vv{I} = \begin{bmatrix} -7 & 2 & -3 \\ 2 & -4 & -6 \\ -1 & -2 & -3 \end{bmatrix} \to\to \begin{bmatrix} 1 & 0 & 1 \\ 0 & 1 & 2 \\ 0 & 0 & 0 \end{bmatrix}

This matrix is equivalent to the system of equations:

x1+x3=0x2+2x3=0\begin{align} x_1 + x_3 &= 0 \\ x_2 + 2 x_3 &= 0 \end{align}

with x3x_3 free. Choosing x3=1x_3 = -1 gives

x1=[121]\vv{x}_1 = \begin{bmatrix} 1 \\ 2 \\ -1 \end{bmatrix}

The eigenvectors for λ2=λ3=3\lambda_2 = \lambda_3 = -3 are obtained from

Aλ2I=[123246123][123000000]\vv{A} - \lambda_2 \vv{I} = \begin{bmatrix} 1 & 2 & -3 \\ 2 & 4 & -6 \\ -1 & -2 & 3 \end{bmatrix} \to\to \begin{bmatrix} 1 & 2 & -3 \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{bmatrix}

This matrix is equivalent to

x1+2x23x3=0x_1 + 2x_2 - 3x_3 = 0

with x2x_2 and x3x_3 free. Choosing x2=1x_2 = 1 and x3=0x_3 = 0 gives one eigenvector

x2=[210]\vv{x}_2 = \begin{bmatrix} -2 \\ 1 \\ 0 \end{bmatrix}

while choosing an independent combination x2=0x_2 = 0 and x3=1x_3 = 1 gives

x3=[301]\vv{x}_3 = \begin{bmatrix} 3 \\ 0 \\ 1 \end{bmatrix}

Complex eigenvalues

Matrices may have complex eigenvaules. They always come in conjugate pairs! For example, for the matrix:

A=[0110]\vv{A} = \begin{bmatrix} 0 & 1 \\ -1 & 0 \end{bmatrix}

We evaluate the determinant:

AλI=λ11λ=λ2+1=0λ1,2=±i\begin{align} |\vv{A} - \lambda \vv{I}| &= \begin{vmatrix} -\lambda & 1 \\ -1 & \lambda \end{vmatrix} \\ &= \lambda^2 + 1 = 0 \\ \lambda_{1,2} = \pm i \end{align}

For λ1=i\lambda_1 = i:

Aλ1I=[i11i][i100]\vv{A} - \lambda_1 \vv{I} = \begin{bmatrix} -i & 1 \\ -1 & i \end{bmatrix} \to \begin{bmatrix} -i & 1 \\ 0 & 0 \end{bmatrix}

This matrix is equivalent to

ix1+x2=0-ix_1 + x_2 = 0

with x2x_2 free. Choosing x2=ix_2 = i gives

x1=[1i]\vv{x}_1 = \begin{bmatrix} 1 \\ i\end{bmatrix}

For λ2=i\lambda_2 = -i:

Aλ1I=[i11i][i100]\vv{A} - \lambda_1 \vv{I} = \begin{bmatrix} i & 1 \\ -1 & i \end{bmatrix} \to \begin{bmatrix} i & 1 \\ 0 & 0 \end{bmatrix}

This matrix is equivalent to

ix1+x2=0i x_1 + x_2 = 0

with x2x_2 free. Choosing x2=ix_2 = -i gives

x2=[1i]\vv{x}_2 =\begin{bmatrix} 1 \\ -i \end{bmatrix}

Useful tricks

Skill builder problems

Find the eigenvalues and eigenvectors for:

Solution to Exercise 1

Since A is upper triangular, the eigenvalues are the diagonal entries, λ1=4\lambda_1 = 4 and λ2=4\lambda_2 = -4. Then, use these eigenvalues to find the eigenvectors that solve (AλI)x=0(\vv{A}-\lambda \vv{I})\vv{x} = \vv{0}.

For λ1=4\lambda_1 = 4,

Aλ1I=[0208]x1=[10]\vv{A}-\lambda_1\vv{I} = \begin{bmatrix} 0 & 2 \\ 0 & -8 \end{bmatrix} \to \vv{x}_1 = \begin{bmatrix} 1 \\ 0 \end{bmatrix}

For λ2=4\lambda_2 = -4,

Aλ2I=[8200]x2=[14]\vv{A}-\lambda_2\vv{I} = \begin{bmatrix} 8 & 2 \\ 0 & 0 \end{bmatrix} \to \vv{x}_2 = \begin{bmatrix} 1 \\ -4 \end{bmatrix}
Solution to Exercise 2

First, find the eigenvalues:

AλI=5λ296λ=(λ5)(λ+6)+18=λ2+λ12=(λ3)(λ+4)=0\begin{align} |\vv{A}-\lambda \vv{I}| &= \begin{vmatrix} 5 - \lambda & -2 \\ 9 & -6 - \lambda \end{vmatrix} \\ &= (\lambda - 5)(\lambda + 6) + 18 \\ &= \lambda^2 + \lambda - 12 \\ &= (\lambda - 3)(\lambda + 4) = 0 \end{align}

So, λ1=3\lambda_1 = 3 and λ2=4\lambda_2 = -4. For λ1=4\lambda_1 = 4,

Aλ1I=[2299]x1=[11]\vv{A}-\lambda_1\vv{I} = \begin{bmatrix} 2 & -2 \\ 9 & -9 \end{bmatrix} \to \vv{x}_1 = \begin{bmatrix} 1 \\ 1 \end{bmatrix}

For λ2=4\lambda_2 = -4,

Aλ2I=[9292]x2=[29]\vv{A}-\lambda_2\vv{I} = \begin{bmatrix} 9 & -2 \\ 9 & -2 \end{bmatrix} \to \vv{x}_2 = \begin{bmatrix} 2 \\ 9 \end{bmatrix}
Solution to Exercise 3

First, find the eigenvalues:

AλI=6λ2225λ0207λ=22207λ+(5λ)6λ227λ=22(7λ)+(5λ)[(λ6)(λ7)4]=λ3+18λ299λ+162=0\begin{align} |\vv{A}-\lambda \vv{I}| &= \begin{vmatrix} 6 - \lambda & 2 & -2 \\ 2 & 5 - \lambda & 0 \\ -2 & 0 & 7 - \lambda \end{vmatrix} \\ &= -2 \begin{vmatrix} 2 & -2 \\ 0 & 7 - \lambda \end{vmatrix} + (5 - \lambda) \begin{vmatrix} 6 - \lambda & -2 \\ -2 & 7 - \lambda \end{vmatrix} \\ &= -2 \cdot 2(7 - \lambda) + (5 - \lambda)[(\lambda - 6)(\lambda - 7) - 4] \\ &= -\lambda^3 + 18\lambda^2 - 99\lambda + 162 = 0 \end{align}

Solving for these roots numerically gives λ1=9\lambda_1 = 9, λ2=6\lambda_2 = 6, and λ3=3\lambda_3 = 3.

For λ1=9\lambda_1 = 9, using row reduction gives

Aλ1I=[322240202][101021000]\vv{A}-\lambda_1\vv{I} = \begin{bmatrix} -3 & 2 & -2 \\ 2 & -4 & 0 \\ -2 & 0 & -2 \end{bmatrix} \to\to \begin{bmatrix} 1 & 0 & 1 \\ 0 & 2 & 1 \\ 0 & 0 & 0 \end{bmatrix}

Choosing x3=2x_3 = 2 as the free variable gives

x1=[212]\vv{x}_1 = \begin{bmatrix} -2 \\ -1 \\ 2 \end{bmatrix}

For λ2=6\lambda_2 = 6, using row reduction gives

Aλ2I=[022210201][201011000]\vv{A}-\lambda_2\vv{I} = \begin{bmatrix} 0 & 2 & -2 \\ 2 & -1 & 0 \\ -2 & 0 & 1 \end{bmatrix} \to\to \begin{bmatrix} -2 & 0 & 1 \\ 0 & 1 & -1 \\ 0 & 0 & 0 \end{bmatrix}

Taking x3=2x_3 = 2 as the free v ariable gives:

x2=[122]\vv{x}_2 = \begin{bmatrix} 1 \\ 2 \\ 2 \end{bmatrix}

For λ3=3\lambda_3 = 3, using row reduction gives

Aλ3I=[322220204][102012000]\vv{A}-\lambda_3\vv{I} = \begin{bmatrix} 3 & 2 & -2 \\ 2 & 2 & 0 \\ -2 & 0 & 4 \end{bmatrix} \to\to \begin{bmatrix} 1 & 0 & -2 \\ 0 & 1 & 2 \\ 0 & 0 & 0 \end{bmatrix}

Choosing x3=1x_3 = 1 as the free variable gives:

x3=[221]\vv{x}_3 = \begin{bmatrix} 2 \\ -2 \\ 1 \end{bmatrix}