Eigenvalues & EigenvectorsTopic #25 of 30

Applications of Eigenvalues

Matrix powers, differential equations, Markov chains, and principal component analysis.

Overview

Eigenvalues and eigenvectors have countless applications across mathematics, science, and engineering. They simplify complex problems by revealing the fundamental modes of a system.

Matrix Powers

If AA is diagonalizable with A=PDP1A = PDP^{-1}:

An=PDnP1A^n = PD^n P^{-1}

where:

Dn=[λ1n000λ2n000λnn]D^n = \begin{bmatrix} \lambda_1^n & 0 & \cdots & 0 \\ 0 & \lambda_2^n & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & \lambda_n^n \end{bmatrix}

Example

A=[4213],P=[2111],D=[5002]A = \begin{bmatrix} 4 & 2 \\ 1 & 3 \end{bmatrix}, \quad P = \begin{bmatrix} 2 & 1 \\ 1 & -1 \end{bmatrix}, \quad D = \begin{bmatrix} 5 & 0 \\ 0 & 2 \end{bmatrix}

A3=PD3P1 where D3=[125008]A^3 = PD^3 P^{-1} \text{ where } D^3 = \begin{bmatrix} 125 & 0 \\ 0 & 8 \end{bmatrix}

Fibonacci Numbers

The Fibonacci sequence Fn+2=Fn+1+FnF_{n+2} = F_{n+1} + F_n can be expressed as:

[Fn+1Fn]=[1110]n[10]\begin{bmatrix} F_{n+1} \\ F_n \end{bmatrix} = \begin{bmatrix} 1 & 1 \\ 1 & 0 \end{bmatrix}^n \begin{bmatrix} 1 \\ 0 \end{bmatrix}

Eigenvalues of [1110]\begin{bmatrix} 1 & 1 \\ 1 & 0 \end{bmatrix}:

λ1=1+521.618 (golden ratio φ)\lambda_1 = \frac{1+\sqrt{5}}{2} \approx 1.618 \text{ (golden ratio } \varphi)

λ2=1520.618\lambda_2 = \frac{1-\sqrt{5}}{2} \approx -0.618

Closed form (Binet's formula):

Fn=φnψn5F_n = \frac{\varphi^n - \psi^n}{\sqrt{5}}

Markov Chains

For transition matrix PP:

  • λ=1\lambda = 1 is always an eigenvalue
  • Stationary distribution is eigenvector for λ=1\lambda = 1
  • Long-term behavior determined by eigenvalue 1

Example

P=[0.80.30.20.7](transition probabilities)P = \begin{bmatrix} 0.8 & 0.3 \\ 0.2 & 0.7 \end{bmatrix} \quad \text{(transition probabilities)}

Stationary: solve (PI)π=0(P - I)\boldsymbol{\pi} = \mathbf{0}

π=[0.60.4]\boldsymbol{\pi} = \begin{bmatrix} 0.6 \\ 0.4 \end{bmatrix} (60% state 1, 40% state 2)

Differential Equations

For dxdt=Ax\frac{d\mathbf{x}}{dt} = A\mathbf{x}:

Solution:

x(t)=c1eλ1tv1+c2eλ2tv2++cneλntvn\mathbf{x}(t) = c_1 e^{\lambda_1 t}\mathbf{v}_1 + c_2 e^{\lambda_2 t}\mathbf{v}_2 + \cdots + c_n e^{\lambda_n t}\mathbf{v}_n

Each term corresponds to an eigenmode.

Stability Analysis

  • All eigenvalues have negative real parts → stable (solutions decay)
  • Any eigenvalue has positive real part → unstable (solutions grow)
  • Eigenvalues on imaginary axis → oscillatory

Principal Component Analysis (PCA)

Given data covariance matrix CC:

  • Eigenvalues = variance in each principal direction
  • Eigenvectors = principal directions
  • Sort by eigenvalue for dimensionality reduction

Data variance along vi=λi\text{Data variance along } \mathbf{v}_i = \lambda_i

Total variance=λ1+λ2++λn\text{Total variance} = \lambda_1 + \lambda_2 + \cdots + \lambda_n

Vibration Analysis

For mechanical systems:

  • Eigenvalues → natural frequencies (ω2=λ\omega^2 = \lambda)
  • Eigenvectors → mode shapes

Mass-spring system: Mx¨=KxM\ddot{\mathbf{x}} = -K\mathbf{x}

Eigenvalue problem: Kv=ω2MvK\mathbf{v} = \omega^2 M\mathbf{v}

Image Compression

Singular Value Decomposition (SVD) uses eigenvalues:

A=UΣVTA = U\Sigma V^T

Keep largest singular values (related to eigenvalues of ATAA^T A) for compression.

Google PageRank

Web page importance via eigenvector of link matrix:

  • Pages are nodes
  • Links define transition probabilities
  • Dominant eigenvector gives page rankings

Quantum Mechanics

  • Eigenvalues = measurable quantities (energy, momentum)
  • Eigenvectors = quantum states
  • Schrödinger equation is an eigenvalue problem

H^ψ=Eψ\hat{H}\psi = E\psi

Summary of Applications

FieldEigenvalue Interpretation
MechanicsNatural frequencies
CircuitsResonant frequencies
StatisticsVariance components
Graph TheoryConnectivity, PageRank
Quantum PhysicsEnergy levels
Control TheorySystem stability
Machine LearningFeature importance