Eigenvalues and Eigenvectors

⚡ Key Takeaways

One-line statement: For a square matrix AA, an eigenvector v≠0\mathbf{v} \neq \mathbf{0} satisfies Av=λvA\mathbf{v} = \lambda\mathbf{v}, meaning AA merely scales v\mathbf{v} by the eigenvalue λ\lambda without changing its direction. Key result in words: Eigenvalues are the roots of the characteristic polynomial det⁡(A−λI)=0\det(A - \lambda I) = 0; eigenvectors span the directions left invariant by the linear map. Where it applies: Diagonalization, differential equations, vibrations, quantum mechanics, PCA, Google's PageRank, stability analysis. Who/when: Roots in 18th-century work by Euler and Lagrange on rotation and quadratic forms; named ("eigen-" = "own/characteristic") by David Hilbert (1904).

💡 The Core Insight

Most vectors get knocked off their original line when you apply a matrix. Eigenvectors are the rare vectors that stay on their own line — the matrix only rescales them. If you can find enough of these invariant directions to build a basis, the matrix becomes diagonal in that basis, and applying it (or its powers, or its exponential) reduces to scalar arithmetic on each axis. Eigen-analysis is therefore the search for the coordinate system in which a linear map is as simple as possible.

Formal Statement

DefinitionEigenvalue and Eigenvector

Let A∈Cn×nA \in \mathbb{C}^{n\times n} be a square matrix (or a linear operator T:V→VT: V \to V on a vector space VV). A scalar λ∈C\lambda \in \mathbb{C} is an eigenvalue of AA if there exists a nonzero vector v∈Cn\mathbf{v} \in \mathbb{C}^n such that

1
Av=λv,v≠0.\quad A\mathbf{v} = \lambda\mathbf{v}, \qquad \mathbf{v} \neq \mathbf{0}.

Derivation hint: This is the defining equation; the restriction v≠0\mathbf{v}\neq\mathbf{0} is essential because A0=λ0A\mathbf{0}=\lambda\mathbf{0} holds trivially for every λ\lambda and would render the notion vacuous. Here AA is the n×nn\times n matrix, λ\lambda is the eigenvalue (a scalar), v\mathbf{v} is the corresponding eigenvector, II is the n×nn\times n identity matrix, and 0\mathbf{0} is the zero vector. The vector v\mathbf{v} is said to belong to λ\lambda.

Rewriting [1][1] as (A−λI)v=0(A - \lambda I)\mathbf{v} = \mathbf{0} shows that v\mathbf{v} lies in the null space of A−λIA - \lambda I. For a nonzero solution to exist, A−λIA - \lambda I must be singular:

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

Derivation hint: A homogeneous system Mv=0M\mathbf{v}=\mathbf{0} has a nontrivial solution iff MM is not invertible, which for a square matrix is equivalent to det⁡M=0\det M = 0.

DefinitionCharacteristic Polynomial, Spectrum, Eigenspace

The function pA(λ)=det⁡(A−λI)p_A(\lambda) = \det(A - \lambda I) is the characteristic polynomial of AA; it is a degree-nn polynomial in λ\lambda. Its roots are precisely the eigenvalues, and the multiset of roots is the spectrum σ(A)\sigma(A). For each eigenvalue λ\lambda, the set

3
Eλ=ker⁡(A−λI)={v∈Cn:Av=λv}\quad E_\lambda = \ker(A - \lambda I) = \{\mathbf{v} \in \mathbb{C}^n : A\mathbf{v} = \lambda\mathbf{v}\}

Derivation hint: The kernel of a linear map is always a subspace (closed under addition and scalar multiplication); EλE_\lambda adds the zero vector to the eigenvectors to form the <strong>eigenspace</strong>. is a subspace called the eigenspace of λ\lambda. Its dimension is the geometric multiplicity g(λ)=dim⁡Eλg(\lambda) = \dim E_\lambda, while the multiplicity of λ\lambda as a root of pAp_A is the algebraic multiplicity a(λ)a(\lambda).

TheoremMultiplicity Inequality and Diagonalizability

For every eigenvalue λ\lambda of AA,

4
1≤g(λ)≤a(λ),∑λ∈σ(A)a(λ)=n.\quad 1 \le g(\lambda) \le a(\lambda), \qquad \sum_{\lambda \in \sigma(A)} a(\lambda) = n.

Derivation hint: Geometric multiplicity is at least 1 because λ\lambda being a root guarantees a nonzero kernel vector; the upper bound follows by extending an eigenbasis of EλE_\lambda and computing the characteristic polynomial in block form. The sum equals nn over C\mathbb{C} by the Fundamental Theorem of Algebra. A matrix AA is diagonalizable (i.e. A=PDP−1A = PDP^{-1} with DD diagonal) if and only if g(λ)=a(λ)g(\lambda) = a(\lambda) for every eigenvalue, equivalently iff AA possesses nn linearly independent eigenvectors.

TheoremTrace and Determinant from Eigenvalues

Counting eigenvalues λ1,…,λn\lambda_1, \dots, \lambda_n with algebraic multiplicity,

5
tr⁡(A)=∑i=1nλi,det⁡(A)=∏i=1nλi.\quad \operatorname{tr}(A) = \sum_{i=1}^{n} \lambda_i, \qquad \det(A) = \prod_{i=1}^{n} \lambda_i.

Derivation hint: Expand pA(λ)=∏i(λi−λ)p_A(\lambda) = \prod_i(\lambda_i - \lambda) and match the coefficient of λn−1\lambda^{n-1} (giving the trace) and the constant term pA(0)=det⁡Ap_A(0)=\det A (giving the product).

Why It Holds / Derivation

ProofEigenvalues are roots of the characteristic polynomial

Step 1. Suppose λ\lambda is an eigenvalue with eigenvector v≠0\mathbf{v}\neq\mathbf{0}. From [1][1], Av−λv=0A\mathbf{v} - \lambda\mathbf{v} = \mathbf{0}, so (A−λI)v=0(A-\lambda I)\mathbf{v} = \mathbf{0}.

Step 2. Thus the homogeneous linear system with coefficient matrix M=A−λIM = A - \lambda I has a nontrivial solution v\mathbf{v}. By the invertible matrix theorem, MM cannot be invertible.

Step 3. A square matrix is non-invertible iff its determinant vanishes. Hence det⁡(A−λI)=0\det(A - \lambda I) = 0, which is equation [2][2].

Step 4. Conversely, if det⁡(A−λI)=0\det(A-\lambda I)=0, then A−λIA-\lambda I is singular, so its kernel contains some v≠0\mathbf{v}\neq\mathbf{0}, which is then an eigenvector for λ\lambda. Therefore the eigenvalues are exactly the roots of pA(λ)=det⁡(A−λI)p_A(\lambda)=\det(A-\lambda I).

Step 5. Expanding the determinant, pAp_A is a polynomial of degree exactly nn with leading term (−1)nλn(-1)^n\lambda^n. By the Fundamental Theorem of Algebra it has nn roots in C\mathbb{C} (with multiplicity), proving every complex square matrix has at least one eigenvalue. ■\blacksquare

∎
ProofEigenvectors for distinct eigenvalues are linearly independent

Step 1. Let v1,…,vk\mathbf{v}_1,\dots,\mathbf{v}_k be eigenvectors for distinct eigenvalues λ1,…,λk\lambda_1,\dots,\lambda_k. Suppose, for contradiction, they are dependent; take a minimal dependent subset, reindexed so c1v1+⋯+cmvm=0c_1\mathbf{v}_1 + \cdots + c_m\mathbf{v}_m = \mathbf{0} with all ci≠0c_i\neq 0 and mm smallest.

Step 2. Apply AA: ∑iciλivi=0\sum_i c_i\lambda_i\mathbf{v}_i = \mathbf{0}. Also multiply the original relation by λm\lambda_m: ∑iciλmvi=0\sum_i c_i\lambda_m\mathbf{v}_i = \mathbf{0}.

Step 3. Subtract: ∑i=1m−1ci(λi−λm)vi=0\sum_{i=1}^{m-1} c_i(\lambda_i - \lambda_m)\mathbf{v}_i = \mathbf{0}. Since the λi\lambda_i are distinct, each (λi−λm)≠0(\lambda_i-\lambda_m)\neq 0, giving a shorter nontrivial dependence — contradicting minimality.

Step 4. Hence no dependence exists; the eigenvectors are linearly independent. A corollary: a matrix with nn distinct eigenvalues is automatically diagonalizable. ■\blacksquare

∎
DH
David Hilbert
1862–1943

Though eigenvalue ideas predate him, it was David Hilbert who, around 1904 in his work on integral equations, coined the German prefix eigen ("own," "characteristic," "proper"). He generalized the finite-dimensional spectral picture to infinite-dimensional function spaces, introducing what we now call Hilbert spaces and laying the spectral foundations that John von Neumann would later use to axiomatize quantum mechanics. Hilbert's insight — that an operator's "own values" govern the resonant frequencies of physical and analytic systems — turned eigenanalysis from a computational trick into a structural principle of modern mathematics and physics.

Building Intuition

🧠 Intuition

Picture a linear map as a deformation of a rubber sheet: vectors get stretched, sheared, and rotated. Almost every arrow you draw will point in a new direction after the map. Eigenvectors are the special arrows that survive the deformation pointing the same way — they lie along the "grain" of the transformation. The eigenvalue tells you what happens along that grain: λ>1\lambda > 1 stretches, 0<λ<10 < \lambda < 1 compresses, λ<0\lambda < 0 flips, and ∣λ∣=1|\lambda| = 1 preserves length. If you find enough independent grains to span the whole space, you can describe the map by saying simply "stretch by λ1\lambda_1 along axis 1, by λ2\lambda_2 along axis 2, …" — that is diagonalization. Complex eigenvalues λ=reiθ\lambda = re^{i\theta} signal that the map rotates by angle θ\theta while scaling by rr; there is then no real invariant line, which is why a 2D rotation has no real eigenvectors.

Figure
A shear matrix scaling its eigenvector while a generic vector turns
x y eigenvector v (λ=2) u Au (turned!) Eigenvector stays on its line; generic u rotates.

Worked Examples

ExampleA symmetric 2×2 matrix

Given: A=(2112)A = \begin{pmatrix} 2 & 1 \\ 1 & 2 \end{pmatrix}. Find its eigenvalues and eigenvectors.

Step 1 — Characteristic polynomial.

pA(λ)=det⁡(2−λ112−λ)=(2−λ)2−1=λ2−4λ+3.p_A(\lambda) = \det\begin{pmatrix} 2-\lambda & 1 \\ 1 & 2-\lambda \end{pmatrix} = (2-\lambda)^2 - 1 = \lambda^2 - 4\lambda + 3.

Step 2 — Roots. λ2−4λ+3=(λ−1)(λ−3)=0⇒λ1=1, λ2=3.\lambda^2 - 4\lambda + 3 = (\lambda-1)(\lambda-3) = 0 \Rightarrow \lambda_1 = 1,\ \lambda_2 = 3. Check: tr⁡=4=1+3\operatorname{tr} = 4 = 1+3 ✓, det⁡=3=1⋅3\det = 3 = 1\cdot 3 ✓.

Step 3 — Eigenvector for λ=1\lambda=1. Solve (A−I)v=0(A - I)\mathbf{v} = \mathbf{0}: (1111)v=0⇒v1+v2=0\begin{pmatrix}1&1\\1&1\end{pmatrix}\mathbf{v}=\mathbf{0} \Rightarrow v_1 + v_2 = 0, so v1=(1,−1)⊤\mathbf{v}_1 = (1, -1)^\top.

Step 4 — Eigenvector for λ=3\lambda=3. (−111−1)v=0⇒v1=v2\begin{pmatrix}-1&1\\1&-1\end{pmatrix}\mathbf{v}=\mathbf{0} \Rightarrow v_1 = v_2, so v2=(1,1)⊤\mathbf{v}_2 = (1,1)^\top. Note the eigenvectors are orthogonal (dot product =0=0), as guaranteed for symmetric matrices.

λ1=1, v1=(1−1);λ2=3, v2=(11)\boxed{\lambda_1 = 1,\ \mathbf{v}_1 = \begin{pmatrix}1\\-1\end{pmatrix}; \qquad \lambda_2 = 3,\ \mathbf{v}_2 = \begin{pmatrix}1\\1\end{pmatrix}}
ExampleComplex eigenvalues from a rotation-scaling matrix

Given: A=(0−110)A = \begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix}, a 90° rotation. Find its eigenvalues.

Step 1 — Characteristic polynomial.

pA(λ)=det⁡(−λ−11−λ)=λ2+1.p_A(\lambda) = \det\begin{pmatrix} -\lambda & -1 \\ 1 & -\lambda \end{pmatrix} = \lambda^2 + 1.

Step 2 — Roots. λ2=−1⇒λ=±i.\lambda^2 = -1 \Rightarrow \lambda = \pm i. No real eigenvalues exist — geometrically obvious, since a rotation leaves no real line fixed.

Step 3 — Eigenvector for λ=i\lambda = i. Solve (A−iI)v=0(A - iI)\mathbf{v}=\mathbf{0}: (−i−11−i)v=0\begin{pmatrix}-i&-1\\1&-i\end{pmatrix}\mathbf{v}=\mathbf{0}. From the second row v1=iv2v_1 = i v_2, giving v=(i,1)⊤\mathbf{v} = (i, 1)^\top.

Step 4 — Interpretation. ∣λ∣=1|\lambda| = 1 confirms a pure rotation (no scaling); the argument arg⁡(i)=90°\arg(i) = 90° matches the rotation angle.

λ=±i,v+i=(i1)\boxed{\lambda = \pm i, \quad \mathbf{v}_{+i} = \begin{pmatrix} i \\ 1 \end{pmatrix}}

Common Misconceptions

⚠️ Common mistake

"Every matrix is diagonalizable." False. A matrix is diagonalizable only when geometric and algebraic multiplicities agree for every eigenvalue (equation [4][4]). Matrices that fail this are called defective. Also beware: an eigenvector is a direction, never unique — any nonzero scalar multiple of v\mathbf{v} is also an eigenvector, so eigenspaces, not single vectors, are the invariant objects. Finally, det⁡(A−λI)=0\det(A-\lambda I)=0 is the characteristic equation; do not confuse it with solving Av=0A\mathbf{v}=\mathbf{0} (that finds the null space, the eigenspace of λ=0\lambda = 0).

CounterexampleA defective matrix

Consider A=(2102)A = \begin{pmatrix} 2 & 1 \\ 0 & 2 \end{pmatrix}. Its characteristic polynomial is (2−λ)2(2-\lambda)^2, so λ=2\lambda = 2 with algebraic multiplicity a=2a = 2. But solving (A−2I)v=(0100)v=0(A - 2I)\mathbf{v} = \begin{pmatrix}0&1\\0&0\end{pmatrix}\mathbf{v} = \mathbf{0} forces v2=0v_2 = 0, giving the single eigenvector direction (1,0)⊤(1,0)^\top. Hence g(2)=1<2=a(2)g(2) = 1 < 2 = a(2): there are not enough eigenvectors to span R2\mathbb{R}^2, so AA is not diagonalizable. The best one can do is the Jordan Canonical Form (2102)\begin{pmatrix}2&1\\0&2\end{pmatrix} — a Jordan block.

Connections

connection

Eigen-analysis is the hub of linear algebra. It underlies Diagonalization (writing A=PDP−1A = PDP^{-1}), generalizes to the Jordan Canonical Form when diagonalization fails, and refines into the Singular Value Decomposition for non-square or non-diagonalizable maps (via eigenvalues of A⊤AA^\top A). For symmetric/Hermitian matrices, the Spectral Theorem guarantees a real spectrum and an orthonormal eigenbasis. The Matrix Exponential eAt=PeDtP−1e^{At} = Pe^{Dt}P^{-1} solves linear Systems of Differential Equations, while Markov Chains and PageRank hinge on the dominant eigenvalue. The characteristic polynomial ties eigenvalues to the Cayley–Hamilton Theorem, and to the Determinant and Trace through equation [5][5].

Real-World Applications

🔬 Where this appears

Historical Development

⏳ History
- 1743 — Euler. While solving systems of linear ODEs, Leonhard Euler encounters "characteristic" exponents, implicitly using eigenvalue ideas.
- 1760s — Lagrange. In studying the secular (long-term) motion of planetary orbits, Joseph-Louis Lagrange analyzes the "secular equation" — the characteristic equation in disguise.
- 1820s — Cauchy. Augustin-Louis Cauchy systematizes the theory of quadratic forms and proves that symmetric matrices have real "characteristic roots," an early Spectral Theorem.
- 1846 — Jacobi. Carl Jacobi develops iterative methods to diagonalize symmetric matrices, foreshadowing numerical eigenvalue algorithms.
- 1858 — Cayley. Arthur Cayley founds matrix algebra and (with Hamilton) the Cayley–Hamilton theorem linking a matrix to its characteristic polynomial.
- 1870s — Jordan & Weierstrass. Camille Jordan introduces the canonical form for defective matrices; Karl Weierstrass develops elementary divisors theory.
- 1904 — Hilbert. David Hilbert coins Eigenwert/*Eigenfunktion* and extends spectral theory to infinite dimensions in his work on integral equations.
- 1925–1932 — Heisenberg, von Neumann. Eigenvalues become the measured quantities of quantum mechanics; von Neumann gives the rigorous spectral theorem for unbounded operators.
- 1960s — Francis & Kublanovskaya. The QR algorithm is invented, making large-scale numerical eigenvalue computation practical and underpinning modern scientific computing.

Subtleties & Limits

open-problem

Numerical and theoretical frontiers. Computing eigenvalues of very large or sparse matrices remains demanding: the QR algorithm is robust but O(n3)O(n^3), motivating ongoing research into Krylov-subspace methods (Lanczos, Arnoldi) and randomized algorithms. Eigenvalue conditioning is delicate — for non-normal matrices, eigenvalues can be wildly sensitive to perturbations, leading to the study of pseudospectra (Trefethen), crucial in fluid-dynamics stability where eigenvalues alone mislead. In infinite dimensions the spectrum may be continuous (no eigenvectors at all, as for the position operator), demanding the full machinery of functional analysis. Inverse eigenvalue problems (reconstruct a matrix/structure from prescribed spectra — "can you hear the shape of a drum?", answered no by Gordon–Webb–Wolpert in 1992) remain a rich, partly open area.

See Also

Diagonalization · Spectral Theorem · Jordan Canonical Form · Singular Value Decomposition · Cayley–Hamilton Theorem · Matrix Exponential · PageRank · Principal Component Analysis · Determinants · Linear Transformations