Understanding Singular Values
A singular value measures the magnitude by which a matrix stretches or compresses vectors along principal axes. For a matrix A, singular values are derived from the eigenvalues of ATA (or A*A for complex matrices, where * denotes conjugate transpose).
Mathematically, if λ is an eigenvalue of ATA, then σ = √λ is a singular value of A. Since ATA is positive semi-definite, all eigenvalues are non-negative, guaranteeing real, non-negative singular values.
The largest singular value equals the operator norm of the matrix under the Euclidean metric—a critical quantity in numerical analysis. Smaller singular values reveal directions where the matrix has minimal influence, crucial for understanding rank-deficient or nearly singular systems.
Computing Singular Values
To find singular values of matrix A:
1. Compute ATA
2. Find eigenvalues λ₁, λ₂, …, λₙ of ATA
3. Singular values: σᵢ = √λᵢ
For a 2×2 matrix:
σ₁,₂ = √[(trace(ATA) ± √(trace²(ATA) − 4·det(ATA))) ÷ 2]
σᵢ— The i-th singular valueλᵢ— The i-th eigenvalue of A^T·AA^T— Transpose of matrix A
Singular Values vs. Eigenvalues
Though related, singular values and eigenvalues serve different purposes:
- Applicability: Every matrix has singular values. Only square matrices possess eigenvalues.
- Sign: Singular values are always real and non-negative. Eigenvalues may be negative or complex.
- Interpretation: Singular values describe how the matrix acts on arbitrary vectors; eigenvalues describe invariant directions where vectors are merely scaled.
- Symmetric matrices: For a symmetric positive semi-definite matrix, singular values equal the absolute values of eigenvalues, but they coincide only when all eigenvalues are non-negative.
This distinction explains why singular value decomposition (SVD) is more universally applicable than eigenvalue decomposition in machine learning, statistics, and image processing.
Special Cases: Diagonal and Symmetric Matrices
Diagonal matrices offer the simplest case: singular values are simply the absolute values of the diagonal entries. This makes them numerically stable and computationally trivial.
Symmetric matrices have an elegant property: their singular values equal the absolute values of their eigenvalues. If the matrix is also positive definite or positive semi-definite, singular values and eigenvalues are identical, making computation straightforward.
These cases are rarely encountered in practice, but recognising them can accelerate mental calculation and build intuition for how matrix structure influences singular values across larger, less structured examples.
Common Pitfalls and Practical Tips
Singular value computation requires care to avoid numerical errors and conceptual misunderstandings.
- Don't confuse eigenvalues with singular values — Singular values come from <strong>A</strong><sup>T</sup><strong>A</strong>, not <strong>A</strong> itself. Extracting eigenvalues of the original matrix will give wrong results unless the matrix is symmetric. Always form the Gram matrix first.
- Watch for rank-deficiency — Zero or near-zero singular values indicate linear dependence or near-dependence among rows or columns. In numerical calculations, a threshold is often used to distinguish true zeros from rounding errors—a value like 1e-10 times the largest singular value.
- Condition number reveals numerical stability — The ratio of the largest to smallest non-zero singular value (the condition number) predicts how sensitive computations involving the matrix are to perturbations. Large condition numbers warn of potential numerical instability in solving linear systems or least-squares problems.
- Rectangular matrices are common — Unlike eigenvalue decomposition, SVD handles non-square matrices seamlessly. Singular values remain real and non-negative regardless of matrix shape, making SVD a universal tool across data science and engineering applications.