The Singular Value Decomposition
Prof. Walter Gander ETH Zurich Decenber 12, 2008
Contents
1 The Singular Value Decomposition 1
2 Applications of the SVD 3 2.1 Condition Numbers ...... 4 2.2 Normal Equations and Condition ...... 7 2.3 The Pseudoinverse ...... 7 2.4 Fundamental Subspaces ...... 8 2.5 General Solution of the Linear Least Squares Problem ...... 9 2.6 Fitting Lines ...... 10 2.7 Fitting Ellipses ...... 11
3 Fitting Hyperplanes–Collinearity Test 13
4 Total Least Squares 15
5 Bibliography 18
1 The Singular Value Decomposition
The singular value decomposition (SVD) of a matrix A is very useful in the context of least squares problems. It also very helpful for analyzing properties of a matrix. With the SVD one x-rays a matrix!
Theorem 1.1 (The Singular Value Decomposition, SVD). Let A be an (m×n) matrix with m ≥ n. Then there exist orthogonal matrices U (m × m) and V (n × n) and a diagonal matrix Σ = diag(σ1, . . . , σn)(m × n) with σ1 ≥ σ2 ≥ ... ≥ σn ≥ 0, such that
A = UΣV T holds. If σr > 0 is the smallest singular value greater than zero then the matrix A has rank r.
Definition 1.1 The column vectors of U = [u1,..., um] are called the left singular vectors and similarly V = [v1,..., vn] are the right singular vectors. The values σi are called the singular values of A.
Proof The 2-norm of A is defined by kAk2 = maxkxk=1 kAxk. Thus there exists a vector x with kxk = 1 such that z = Ax, kzk = kAk2 =: σ. Let y := z/kzk. We have obtained Ax = σy with kxk = kyk = 1. Now we construct vectors orthogonal to x and similarly to y and form the matrices V = [x,V1] and U = [y,U1]. This could e.g. be done in Matlab by the command > V = orth([x, rand(n, n-1)]), U = orth([y, rand(n, n-1)])
1 Now T T T T T y y Ax y AV1 σ w A1 = U AV = T A [x,V1] = T T = U1 U1 Ax U1 AV1 0 B T T T T T because y Ax = y σy = σy y = σ and U1 Ax = σU1 y = 0 since U1 ⊥ y. T T We claim that w := y AV1 = 0. In order to prove that we compute σ σ2 + kwk2 A = 1 w Bw and conclude from that equation that 2 σ 2 2 2 2 2 2 2 A1 = (σ + kwk )) + kBwk ≥ (σ + kwk ) . w
T Now since V and U are orthogonal kA1k2 = kU AV k2 = kAk2 = σ holds and