matrix theory
Spectral Radius
You should know: eigenvalues and eigenvectors, matrix norms
Overview
The spectral radius ρ(A) of a matrix A is the largest absolute value among all eigenvalues. It governs the long-run behavior of matrix powers: A^n → 0 iff ρ(A) < 1. The spectral radius satisfies ρ(A) ≤ ||A|| for any matrix norm, and equals the infimum of all consistent norms. It is fundamental to the convergence analysis of iterative methods for linear systems.
Intuition
Raise a matrix to a high power: the dominant behavior is controlled by the eigenvalue with the largest absolute value. If that eigenvalue has |λ| < 1, repeated multiplication shrinks all vectors to zero. If |λ| > 1, the iteration blows up. The spectral radius is the threshold — it tells you whether a matrix iteration will converge, diverge, or oscillate.
Formal Definition
For a square matrix A with eigenvalues λ_1, ..., λ_n (counted with multiplicity):
Maximum absolute value of the eigenvalues of A.
Notation
| Notation | Meaning |
|---|---|
| Spectral radius of A | |
| Spectrum of A: the set of all eigenvalues |
Properties
Power convergence
Gelfand formula
Norm inequality
Normal matrices
Spectral radius of A^k
Applications
PageRank's convergence depends on the spectral radius of the Google matrix being less than 1, ensuring the power iteration converges to the stationary distribution.
Worked Examples
- 1
A is upper triangular, so eigenvalues are the diagonal entries: λ_1=2, λ_2=-3.
- 2
Spectral radius is the largest absolute value.
✓ Answer
ρ(A) = 3.
Practice Problems
Find the spectral radius of A = [[0,1],[-2,3]].
The Jacobi iteration for solving Ax=b uses iteration matrix G = D^{-1}(L+U) where A = D - L - U. State the condition on ρ(G) for convergence, and explain why.
Common Mistakes
The spectral radius equals ||A||_2.
||A||_2 = σ_max (largest singular value), while ρ(A) = max|λ_i| (largest eigenvalue magnitude). They are equal only for normal matrices.
Quiz
Summary
- The spectral radius ρ(A) = max|λ_i| is the largest absolute value among eigenvalues of A.
- A^k → 0 as k → ∞ if and only if ρ(A) < 1.
- Gelfand's formula: ρ(A) = lim_{k→∞} ||A^k||^{1/k} for any consistent norm.
- ρ(A) ≤ ||A|| for any sub-multiplicative norm; equality holds for normal matrices in the spectral norm.
- Controls convergence of iterative linear solvers (Jacobi, Gauss-Seidel) and stability of discrete dynamical systems.
References
- WebsiteWikipedia — Spectral radius
Mathematics