Skip to content

Eigenvalues and Eigenvectors

Attempt these from memory before reading.

  • Define an eigenvalue and an eigenvector.
  • How do you find them? Write the equation whose roots they are.
  • What does it mean for a matrix to be diagonalisable, and what is the condition?
  • State the Cayley-Hamilton theorem.
  • What does the spectral theorem guarantee for a real symmetric matrix, and why is that matrix class special?

Check your answers against the text. Prompt 3 is the one that separates computing eigenvalues from understanding the map: the condition is about eigenvectors, not eigenvalues.

Motivation: eigenvectors are the directions the map does not rotate

Section titled “Motivation: eigenvectors are the directions the map does not rotate”

A linear map mixes coordinates. Almost every vector you feed in comes out pointing somewhere unrelated. The eigenvectors are the exceptions: they come out pointing along themselves, changed only in length.

That is the whole idea, and everything else follows from wanting to exploit it. If a basis of eigenvectors exists, then in that basis the map does nothing but scale each coordinate independently — which is what diagonalisation is. A hard coupled problem becomes nn easy uncoupled ones.

The characteristic polynomial is where the eigenvalues hide. An eigenvector satisfies (A−λI)v=0(A - \lambda I)v = 0 with v≠0v \neq 0, which is possible only when A−λIA - \lambda I is singular, so

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

That determinant is a polynomial of degree nn, so an n×nn \times n matrix has nn eigenvalues counted with multiplicity — over C\mathbb{C}, always. Over R\mathbb{R} it may have none at all, and a rotation is the standard example.

Not every matrix is diagonalisable, and Jordan form is what you get instead. Diagonalisability is a statement about eigenvectors: you need nn of them, independent. A repeated eigenvalue may come with fewer independent eigenvectors than its algebraic multiplicity, and the missing directions are supplied by generalised eigenvectors — which is what a Jordan block records. The distinction between algebraic and geometric multiplicity is the entire content of section 5.5.

Why symmetric matrices are the special case. A real symmetric matrix is always diagonalisable, by an orthogonal matrix, and always has real eigenvalues. That is the spectral theorem, and it is why symmetric matrices are the ones that appear as observables in quantum mechanics and as quadratic forms in geometry: their eigenbases are orthonormal, so the change of basis costs nothing geometrically.

Let A∈Mn×n(F)A \in \mathcal{M}_{n \times n}(F). A scalar λ∈F\lambda \in F is an eigenvalue of AA if there Exists a non-zero vector v∈Fn\mathbf{v} \in F^n such that

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

The vector v\mathbf{v} is called an eigenvector corresponding to λ\lambda.

The eigenspace corresponding to λ\lambda is Eλ=ker⁡(A−λI)E_\lambda = \ker(A - \lambda I). Its dimension is the geometric multiplicity of λ\lambda.

Theorem 5.1. λ\lambda is an eigenvalue of AA if and only if det⁡(A−λI)=0\det(A - \lambda I) = 0.

The polynomial p(λ)=det⁡(A−λI)p(\lambda) = \det(A - \lambda I) is called the characteristic polynomial of AA. Its roots (in the algebraic closure of FF) are the eigenvalues of AA.

If p(λ)=(lambda−lambda1)m1(lambda−lambda2)m2⋯(lambda−lambdak)mkp(\lambda) = (\\lambda - \\lambda_1)^{m_1}(\\lambda - \\lambda_2)^{m_2}\cdots(\\lambda - \\lambda_k)^{m_k} with lambda1,…,lambdak\\lambda_1, \ldots, \\lambda_k distinct, then mim_i is the algebraic multiplicity of lambdai\\lambda_i.

Proposition 5.2. For each eigenvalue λ\lambda, 1≤dim(Eλ)≤mλ1 \leq \mathrm{dim}(E_\lambda) \leq m_\lambda (geometric multiplicity does not exceed algebraic multiplicity).

Definition. AA is diagonalisable if there exists an invertible matrix PP and a diagonal Matrix DD such that

A=PDP−1A = PDP^{-1}

Theorem 5.3. A∈Mn×n(F)A \in \mathcal{M}_{n \times n}(F) is diagonalisable (over FF) if and only if AA has nn linearly independent Eigenvectors (over FF). Equivalently, the sum of the geometric multiplicities equals nn.

Corollary 5.4. If AA has nn distinct eigenvalues, then AA is diagonalisable.

Proof. Eigenvectors corresponding to distinct eigenvalues are linearly independent. With nn distinct Eigenvalues, we obtain nn linearly independent eigenvectors, which form a basis of FnF^n. ■\blacksquare

Theorem 5.5 (Cayley—Hamilton). Every square matrix satisfies its own characteristic polynomial: If p(λ)=det⁡(λI−A)p(\lambda) = \det(\lambda I - A) Then p(A)=0p(A) = 0 (the zero matrix).

Proof sketch. Let p(λ)=λn+cn−1λn−1+⋯+c1λ+c0p(\lambda) = \lambda^n + c_{n-1}\lambda^{n-1} + \cdots + c_1\lambda + c_0. By the adjugate formula (Theorem 3.5), (λI−A)⋅adj(λI−A)=p(λ)⋅I(\lambda I - A) \cdot \mathrm{adj}(\lambda I - A) = p(\lambda) \cdot I. Each entry of adj(λI−A)\mathrm{adj}(\lambda I - A) is a polynomial in λ\lambda of degree at most n−1n - 1 So we can write adj(λI−A)=Bn−1λn−1+⋯+B1λ+B0\mathrm{adj}(\lambda I - A) = B_{n-1}\lambda^{n-1} + \cdots + B_1\lambda + B_0 for Matrices BiB_i. Multiplying out and comparing coefficients of λk\lambda^k:

Bn−1=I,Bn−2−ABn−1=cn−1I,…,−AB0=c0IB_{n-1} = I, \quad B_{n-2} - AB_{n-1} = c_{n-1}I, \quad \ldots, \quad -AB_0 = c_0 I

Multiplying the kk-th equation on the left by AkA^k and summing over kk:

AnBn−1+An−1(Bn−2−ABn−1)+⋯+A0(−AB0)=An+cn−1An−1+⋯+c0I=p(A)A^n B_{n-1} + A^{n-1}(B_{n-2} - AB_{n-1}) + \cdots + A^0(-AB_0) = A^n + c_{n-1}A^{n-1} + \cdots + c_0 I = p(A)

But the left side telescopes to zero, so p(A)=0p(A) = 0. ■\blacksquare

When a matrix is not diagonalisable, the Jordan normal form provides the next-best canonical Representation.

Theorem 5.6. Let A∈Mn×n(C)A \in \mathcal{M}_{n \times n}(\mathbb{C}). Then AA is similar to a block-diagonal Matrix

J=(J1⋱Jk)J = \begin{pmatrix} J_1 & & \\ & \ddots & \\ & & J_k \end{pmatrix}

Where each Jordan block has the form

Ji=(λi1λi⋱⋱1λi)J_i = \begin{pmatrix} \lambda_i & 1 & & \\ & \lambda_i & \ddots & \\ & & \ddots & 1 \\ & & & \lambda_i \end{pmatrix}

The Jordan form is unique up to permutation of the blocks.

Intuition. Each Jordan block corresponds to one eigenvalue. The size of the block equals the Number of steps in the chain v,(A−λI)v,(A−λI)2v,…\mathbf{v}, (A - \lambda I)\mathbf{v}, (A - \lambda I)^2\mathbf{v}, \ldots Of generalised eigenvectors. A diagonalisable matrix has all Jordan blocks of size 1×11 \times 1.

Problem. Find the Jordan normal form of

A=(3103)A = \begin{pmatrix} 3 & 1 \\ 0 & 3 \end{pmatrix}
Solution

The characteristic polynomial is det⁡(A−λI)=(3−λ)2\det(A - \lambda I) = (3 - \lambda)^2 So λ=3\lambda = 3 is the Only eigenvalue with algebraic multiplicity 2.

A−3I=(0100)A - 3I = \begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix} Which has rank 1, so the geometric Multiplicity is dim⁡(ker⁡(A−3I))=2−1=1\dim(\ker(A - 3I)) = 2 - 1 = 1.

Since the geometric multiplicity (1) is less than the algebraic multiplicity (2), AA is not Diagonalisable. The Jordan form has one block of size 2:

J=(3103)J = \begin{pmatrix} 3 & 1 \\ 0 & 3 \end{pmatrix}

(In this case, AA is already in Jordan form.) ■\blacksquare

5.6 Spectral Theorem for Real Symmetric Matrices

Section titled “5.6 Spectral Theorem for Real Symmetric Matrices”

Theorem 5.7 (Spectral Theorem). If A∈Mn×n(R)A \in \mathcal{M}_{n \times n}(\mathbb{R}) is symmetric (A=ATA = A^T), then:

  1. All eigenvalues of AA are real.
  2. AA has nn linearly independent orthonormal eigenvectors.
  3. AA is orthogonally diagonalisable: A=QDQTA = QDQ^T where QQ is orthogonal (QTQ=IQ^TQ = I).

Proof. We prove (1) and then (2) and (3) by induction on nn.

(1) Let λ∈C\lambda \in \mathbb{C} be an eigenvalue with eigenvector v∈Cn\mathbf{v} \in \mathbb{C}^n v≠0\mathbf{v} \neq \mathbf{0}. Then

v‾TAv=v‾T(λv)=λv‾Tv\overline{\mathbf{v}}^T A \mathbf{v} = \overline{\mathbf{v}}^T (\lambda \mathbf{v}) = \lambda \overline{\mathbf{v}}^T \mathbf{v}

Since A=ATA = A^T and AA has real entries, A‾=A=AT\overline{A} = A = A^T So

v‾TAv=(Av‾)Tv=(Av‾)Tv=(λ‾ v‾)Tv=λ‾ v‾Tv\overline{\mathbf{v}}^T A \mathbf{v} = (A\overline{\mathbf{v}})^T \mathbf{v} = (\overline{A\mathbf{v}})^T \mathbf{v} = (\overline{\lambda}\,\overline{\mathbf{v}})^T \mathbf{v} = \overline{\lambda}\,\overline{\mathbf{v}}^T \mathbf{v}

Therefore (λ−λ‾)v‾Tv=0(\lambda - \overline{\lambda})\overline{\mathbf{v}}^T\mathbf{v} = 0. Since v‾Tv>0\overline{\mathbf{v}}^T\mathbf{v} \gt 0 We have λ=λ‾\lambda = \overline{\lambda} So λ∈R\lambda \in \mathbb{R}.

(2) and (3) By induction. For n=1n = 1 the result is trivial. Assume it holds for (n−1)×(n−1)(n-1) \times (n-1) Symmetric matrices. Since all eigenvalues are real, AA has a real eigenvalue λ1\lambda_1 with real Eigenvector v1\mathbf{v}_1. Normalise: q1=v1/∥v1∥\mathbf{q}_1 = \mathbf{v}_1 / \lVert \mathbf{v}_1 \rVert.

Let W=q1⊥={w∈Rn:q1Tw=0}W = \mathbf{q}_1^\perp = \{\mathbf{w} \in \mathbb{R}^n : \mathbf{q}_1^T \mathbf{w} = 0\}. For any w∈W\mathbf{w} \in W:

q1T(Aw)=(Aq1)Tw=(λ1q1)Tw=λ1⋅0=0\mathbf{q}_1^T (A\mathbf{w}) = (A\mathbf{q}_1)^T \mathbf{w} = (\lambda_1 \mathbf{q}_1)^T \mathbf{w} = \lambda_1 \cdot 0 = 0

So Aw∈WA\mathbf{w} \in W. Therefore AA restricts to a symmetric linear map A∣W:W→WA|_W : W \to W on an (n−1)(n-1)-dimensional space. By the inductive hypothesis, WW has an orthonormal basis {q2,…,qn}\{\mathbf{q}_2, \ldots, \mathbf{q}_n\} of eigenvectors of A∣WA|_W.

Then {q1,q2,…,qn}\{\mathbf{q}_1, \mathbf{q}_2, \ldots, \mathbf{q}_n\} is an orthonormal eigenbasis for Rn\mathbb{R}^n And A=QDQTA = QDQ^T with Q=[q1∣⋯∣qn]Q = [\mathbf{q}_1 \mid \cdots \mid \mathbf{q}_n]. ■\blacksquare

Problem. Find the eigenvalues, eigenvectors, and diagonalise

A=(4123)A = \begin{pmatrix} 4 & 1 \\ 2 & 3 \end{pmatrix}
Solution

The characteristic polynomial is

det⁡(A−λI)=det⁡(4−λ123−λ)=(4−λ)(3−λ)−2\det(A - \lambda I) = \det\begin{pmatrix} 4 - \lambda & 1 \\ 2 & 3 - \lambda \end{pmatrix} = (4 - \lambda)(3 - \lambda) - 2 =λ2−7λ+10=(λ−5)(λ−2)= \lambda^2 - 7\lambda + 10 = (\lambda - 5)(\lambda - 2)

So the eigenvalues are λ1=5\lambda_1 = 5 and λ2=2\lambda_2 = 2.

For λ1=5\lambda_1 = 5: Solve (A−5I)v=0(A - 5I)\mathbf{v} = \mathbf{0}.

(−112−2)v=0  ⟹  −v1+v2=0  ⟹  v=t(11)\begin{pmatrix} -1 & 1 \\ 2 & -2 \end{pmatrix}\mathbf{v} = \mathbf{0} \implies -v_1 + v_2 = 0 \implies \mathbf{v} = t\begin{pmatrix} 1 \\ 1 \end{pmatrix}

For λ2=2\lambda_2 = 2: Solve (A−2I)v=0(A - 2I)\mathbf{v} = \mathbf{0}.

(2121)v=0  ⟹  2v1+v2=0  ⟹  v=t(1−2)\begin{pmatrix} 2 & 1 \\ 2 & 1 \end{pmatrix}\mathbf{v} = \mathbf{0} \implies 2v_1 + v_2 = 0 \implies \mathbf{v} = t\begin{pmatrix} 1 \\ -2 \end{pmatrix}

Therefore A=PDP−1A = PDP^{-1} with

P=(111−2),D=(5002)P = \begin{pmatrix} 1 & 1 \\ 1 & -2 \end{pmatrix}, \quad D = \begin{pmatrix} 5 & 0 \\ 0 & 2 \end{pmatrix} P−1=1−3(−2−1−11)=(2/31/31/3−1/3)P^{-1} = \frac{1}{-3}\begin{pmatrix} -2 & -1 \\ -1 & 1 \end{pmatrix} = \begin{pmatrix} 2/3 & 1/3 \\ 1/3 & -1/3 \end{pmatrix}

Verification: PDP−1=(111−2)(5002)(2/31/31/3−1/3)PDP^{-1} = \begin{pmatrix} 1 & 1 \\ 1 & -2 \end{pmatrix}\begin{pmatrix} 5 & 0 \\ 0 & 2 \end{pmatrix}\begin{pmatrix} 2/3 & 1/3 \\ 1/3 & -1/3 \end{pmatrix}

=(525−4)(2/31/31/3−1/3)=(10/3+2/35/3−2/310/3−4/35/3+4/3)=(4123)=A= \begin{pmatrix} 5 & 2 \\ 5 & -4 \end{pmatrix}\begin{pmatrix} 2/3 & 1/3 \\ 1/3 & -1/3 \end{pmatrix} = \begin{pmatrix} 10/3 + 2/3 & 5/3 - 2/3 \\ 10/3 - 4/3 & 5/3 + 4/3 \end{pmatrix} = \begin{pmatrix} 4 & 1 \\ 2 & 3 \end{pmatrix} = A. ■\blacksquare

5.7a Intuition: What Do Eigenvalues and Eigenvectors Mean?

Section titled “5.7a Intuition: What Do Eigenvalues and Eigenvectors Mean?”

An eigenvector of a linear transformation AA is a direction that is preserved by AA: applying AA to an eigenvector only stretches or compresses it, without rotating it. The eigenvalue λ\lambda measures how much the eigenvector is stretched: if λ>1\lambda > 1, the eigenvector is elongated; if 0<λ<10 < \lambda < 1, it is compressed; if λ<0\lambda < 0, it is reflected and then scaled.

Geometrically, the eigenvectors of a matrix reveal the “natural axes” of the transformation. For a 2×22 \times 2 matrix, the eigenvectors define the directions along which the transformation acts as simple scaling. In the eigenbasis, the matrix is diagonal: the transformation is just independent stretching along each axis.

This is why diagonalisation is powerful. Computing A100A^{100} directly requires 100 matrix multiplications, but D100D^{100} is trivial: just raise each diagonal entry to the 100th power. The change-of-basis matrix PP handles the translation between the standard basis and the eigenbasis.

Physical examples:

  • In mechanics, the eigenvectors of the inertia tensor are the principal axes of rotation.
  • In vibration analysis, the eigenvectors of the stiffness matrix are the normal modes, and the eigenvalues are the squared natural frequencies.
  • In Google’s PageRank, the eigenvector corresponding to eigenvalue 1 of the web matrix gives the importance ranking of all pages.

Why eigenvalues matter for stability: In the ODE x′=Ax\mathbf{x}' = A\mathbf{x}, the solution is x(t)=eAtx0\mathbf{x}(t) = e^{At}\mathbf{x}_0. If AA is diagonalisable, this becomes x(t)=PeDtP−1x0\mathbf{x}(t) = P e^{Dt} P^{-1} \mathbf{x}_0. The behavior is governed by eλite^{\lambda_i t}: if all Re(λi)<0\mathrm{Re}(\lambda_i) < 0, the system decays to zero (stable); if any Re(λi)>0\mathrm{Re}(\lambda_i) > 0, the system grows without bound (unstable).

Problem. Diagonalise the matrix

A=(210030003)A = \begin{pmatrix} 2 & 1 & 0 \\ 0 & 3 & 0 \\ 0 & 0 & 3 \end{pmatrix}
Solution

Since AA is upper triangular, the eigenvalues are the diagonal entries: λ1=2\lambda_1 = 2, λ2=3\lambda_2 = 3 (with algebraic multiplicity 2).

For λ1=2\lambda_1 = 2: Solve (A−2I)v=0(A - 2I)\mathbf{v} = \mathbf{0}.

A−2I=(010010001)→(010001000)A - 2I = \begin{pmatrix} 0 & 1 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{pmatrix} \to \begin{pmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ 0 & 0 & 0 \end{pmatrix}

Free variable: x1=tx_1 = t. Eigenvector: v1=(1,0,0)T\mathbf{v}_1 = (1, 0, 0)^T.

For λ2=3\lambda_2 = 3: Solve (A−3I)v=0(A - 3I)\mathbf{v} = \mathbf{0}.

A−3I=(−110000000)A - 3I = \begin{pmatrix} -1 & 1 & 0 \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix}

Free variables: x2=sx_2 = s, x3=tx_3 = t. Then x1=sx_1 = s. Eigenvectors: s(1,1,0)T+t(0,0,1)Ts(1, 1, 0)^T + t(0, 0, 1)^T.

The geometric multiplicity of λ2=3\lambda_2 = 3 is 2, equal to its algebraic multiplicity. Therefore AA is diagonalisable with

P=(110010001),D=(200030003)P = \begin{pmatrix} 1 & 1 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{pmatrix}, \quad D = \begin{pmatrix} 2 & 0 & 0 \\ 0 & 3 & 0 \\ 0 & 0 & 3 \end{pmatrix}

Key observation: When the geometric multiplicity equals the algebraic multiplicity for every eigenvalue, the matrix is diagonalisable. When they differ (as in the Jordan form example in Section 5.5), the matrix is not diagonalisable, and the Jordan normal form is the best alternative. ■\blacksquare

Problem. Use the Cayley—Hamilton theorem to compute A10A^{10} for the same matrix AA above.

Solution

The characteristic polynomial is p(λ)=λ2−7λ+10p(\lambda) = \lambda^2 - 7\lambda + 10 So by Cayley—Hamilton, A2=7A−10IA^2 = 7A - 10I.

To find A10A^{10} Divide λ10\lambda^{10} by p(λ)p(\lambda):

λ10=q(λ)(λ2−7λ+10)+r(λ)\lambda^{10} = q(\lambda)(\lambda^2 - 7\lambda + 10) + r(\lambda)

Where r(λ)=aλ+br(\lambda) = a\lambda + b has degree less than 2. Then A10=r(A)=aA+bIA^{10} = r(A) = aA + bI.

To find aa and bb Evaluate at the eigenvalues:

λ10∣λ=5=510=9765625=5a+b\lambda^{10}\big|_{\lambda=5} = 5^{10} = 9765625 = 5a + b

λ10∣λ=2=210=1024=2a+b\lambda^{10}\big|_{\lambda=2} = 2^{10} = 1024 = 2a + b

Subtracting: 3a=9765625−1024=97646013a = 9765625 - 1024 = 9764601 So a=3254867a = 3254867.

b=1024−2⋅3254867=1024−6509734=−6508710b = 1024 - 2 \cdot 3254867 = 1024 - 6509734 = -6508710.

Therefore A10=3254867⋅A−6508710⋅IA^{10} = 3254867 \cdot A - 6508710 \cdot I. ■\blacksquare