Linear Algebra
9.18 min read

Motivation for SVD

The spectral theorem gives a beautiful decomposition A=UΛUA = U\Lambda U^\top for symmetric matrices. But what about a general m×nm \times n matrix AA that might not even be square? We need a generalization — something that reveals the geometric structure of any linear transformation, not just symmetric ones.

The key insight: even if AA itself is not symmetric, the products AAA^\top A (an n×nn \times n matrix) and AAAA^\top (an m×mm \times m matrix) are always symmetric and PSD. We can apply the spectral theorem to each of them. The Singular Value Decomposition (SVD) packages this information into a single clean factorization A=UΣVA = U \Sigma V^\top that works for any matrix.

Geometrically, the SVD says every linear map xAx\mathbf{x} \mapsto A\mathbf{x} is a composition of: (1) a rotation/reflection VV^\top in the domain, (2) a scaling Σ\Sigma along coordinate axes, and (3) a rotation/reflection UU in the codomain. This is the most fundamental geometric description of any linear transformation.

Formal View

Remark 9.1 — Need for SVD
The spectral theorem applies only to symmetric matrices. For a general m×nm \times n matrix AA, the products AAA^\top A (size n×nn \times n) and AAAA^\top (size m×mm \times m) are both symmetric and PSD. SVD packages both spectral decompositions into one unified factorization A=UΣVA = U\Sigma V^\top.

Why This Matters

SVD is the Swiss army knife of linear algebra — it solves least squares, compresses data, finds low-rank structure, and enables PCA.

  • Image compression: keep only the top singular values and vectors to get a compact representation.
  • Natural language processing: latent semantic analysis uses SVD to find topics in text.
  • Recommendation systems: matrix factorization via SVD identifies latent factors in user-item matrices.

Quiz

Question 1

The spectral theorem applies to all matrices, not just symmetric ones.

Question 2

For any m×nm \times n matrix AA, the product AAA^\top A is:

Common Mistakes

  • Thinking SVD is just for square matrices — it works for any m×nm \times n matrix, even non-square.
  • Confusing SVD with eigendecomposition — eigendecomposition applies to square matrices, SVD applies universally.