Back to Understanding Linear Algebra

Singular Value Decompositions

Factoring any rectangular matrix as U S V^T, built from the eigenvectors of A^T A, yielding singular values and bases for all four fundamental subspaces.

Understanding Linear Algebra, Singular Value Decompositions

Status
FlagsLoading your progress…

14 exercises

Activity 7.4.3

In this activity, we will construct the singular value decomposition of A=[10−1111]A=\begin{bmatrix} 1 & 0 & -1 \\ 1 & 1 & 1 \end{bmatrix}. Notice that this matrix is not square so there are no eigenvalues and eigenvectors associated to it.

  1. Construct the Gram matrix G=ATAG=A^TA and find an orthogonal diagonalization of it.

  2. Identify the singular values of AA and the right singular vectors v1\vvec_1, v2\vvec_2, and v3\vvec_3. What is the dimension of these vectors? How many nonzero singular values a…

2.5

Activity 7.4.4

Let's suppose that a matrix AA has a singular value decomposition A=UΣVTA=U\Sigma V^T where

  1. What is the shape of AA; that is, how many rows and columns does AA have?

  2. Suppose we write a three-dimensional vector x\xvec as a linear combination of right singular vectors: We would like to find an expression for AxA\xvec.

    To begin, VTx=[v1⋅xv2⋅xv3⋅x]=[c1c2c3]V^T\xvec = \threevec{\vvec_1\cdot\xvec} {\vvec_2\cdot\xvec} {\vvec_3\cdot\xvec} = \threevec{c_1}{c_2}{c_3}.

    Now ΣVTx=[2000050000000][c1c2c3]=[20c15c200]\Sigma V^T \xvec = \begin{bmatrix} 20 & 0 & 0 \\ 0 & 5 & 0 \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{bmatrix}\threevec{c_1}{c_2}{c_3} = \cfourvec{20c_1}{5c_2}00.

    And finally, Ax=UΣVTx=[u1u2u3u4][20c15c200]=20c1u1+5c2u2A\xvec = U\Sigma V^T\xvec = \begin{bmatrix} \uvec_1 & \uvec_2 & \uvec_3 & \uvec_4 \end{bmatrix} \cfourvec{20c_1}{5c_2}00 = 20c_1\uvec_1 + 5c_2\uvec_2.

    To summarize, we have Ax=20c1u1+5c2u2A\xvec = 20c_1\uvec_1 + 5c_2\uvec_2.

    What condition on c1c_1, c2c_2, and c3c_3 must be satisfied if x\xvec is a s…

2.5

Activity 7.4.5

Suppose we have a singular value decomposition A=UΣVTA = U\Sigma V^T where

  1. What is the shape of AA? What is rank⁡(A)\rank(A)?

  2. Identify bases for Col(A)\col(A) and Col(AT)\col(A^T).

  3. Explain why

  4. Explain why

  5. If A=UΣVTA = U\Sigma V^T, explain why A=UrΣrVrTA=U_r\Sigma_rV_r^T where the columns of UrU_r are an orthonormal basis for Col(A)\col(A), Σr\Sigma_r is a square, diagonal, invertible matrix, and the columns of VrV_r form an orthonormal basis for Col(AT)\col(A^T).

2.5

Exercise 1

Consider the matrix A=[1210−12]A = \begin{bmatrix} 1 & 2 & 1 \\ 0 & -1 & 2 \\ \end{bmatrix}.

  1. Find the Gram matrix G=ATAG=A^TA and use it to find the singular values and right singular vectors of AA.

  2. Find the left singular vectors.

  3. Form the matrices UU, Σ\Sigma, and VV and verify that A=UΣVTA=U\Sigma V^T.

  4. What is rank⁡(A)\rank(A) and what does this say about Col(A)\col(A)?

  5. Determine an orthonormal basis for Nul(A)\nul(A).

2.5

Exercise 2

Find singular value decompositions for the following matrices:

  1. [000−8]\begin{bmatrix} 0 & 0 \\ 0 & -8 \end{bmatrix}.

  2. [2302]\begin{bmatrix} 2 & 3 \\ 0 & 2 \end{bmatrix}.

  3. [400002]\begin{bmatrix} 4 & 0 & 0 \\ 0 & 0 & 2 \end{bmatrix}

  4. [400002]\begin{bmatrix} 4 & 0 \\ 0 & 0 \\ 0 & 2 \end{bmatrix}

2.5

Exercise 3

Consider the matrix A=[2112]A = \begin{bmatrix} 2 & 1 \\ 1 & 2 \end{bmatrix}.

  1. Find a singular value decomposition of AA and verify that it is also an orthogonal diagonalization of AA.

  2. If AA is a symmetric, positive semidefinite matrix, explain why a singular value decomposition of AA is an orthogonal diagonalization of AA.

2.5

Exercise 4

Suppose that the matrix AA has the singular value decomposition

  1. What are the dimensions of AA?

  2. What is rank⁡(A)\rank(A)?

  3. Find orthonormal bases for Col(A)\col(A), Nul(A)\nul(A), Col(AT)\col(A^T), and Nul(AT)\nul(A^T).

  4. Find the orthogonal projection of b=[102−1]\bvec=\fourvec102{-1} onto Col(A)\col(A).

2.5

Exercise 5

Consider the matrix A=[10−1220−112]A = \begin{bmatrix} 1 & 0 & -1 \\ 2 & 2 & 0 \\ -1 & 1 & 2\\ \end{bmatrix}.

  1. Construct the Gram matrix GG and use it to find the singular values and right singular vectors v1\vvec_1, v2\vvec_2, and v3\vvec_3 of AA. What are the matrices Σ\Sigma and VV in a singular value decomposition?

  2. What is rank⁡(A)\rank(A)?

  3. Find as many left singular vectors uj\uvec_j as you can using the relationship Avj=σjujA\vvec_j=\sigma_j\uvec_j.

  4. Find an orthonormal basis for Nul(AT)\nul(A^T) and use it to construct the matrix UU so that A=UΣVTA=U\Sigma V^T.

  5. State an orthonormal basis for Nul(A)\nul(A) and an…

2.5

Exercise 6

Consider the matrix B=[102−112]B=\begin{bmatrix} 1 & 0 \\ 2 & -1 \\ 1 & 2 \end{bmatrix} and notice that B=ATB=A^T where AA is the matrix in this exercise.

  1. Use your result from this exercise to find a singular value decomposition of B=UΣVTB=U\Sigma V^T.

  2. What is rank⁡(B)\rank(B)? Determine a basis for Col(B)\col(B) and Col(B)⊥\col(B)^\perp.

  3. Suppose that b=[−347]\bvec=\threevec{-3}47. Use the bases you found in the previous part of this exercise to write b=b^+b⊥\bvec=\bhat+\bvec^\perp, where b^\bhat is in Col(B)\col(B) and b⊥\bvec^\perp is in Col(B)⊥\col(B)^\perp.

  4. Find the least-squares approximate solution to the equation Bx=bB\xvec=\bvec.

2.5

Exercise 7

Suppose that AA is a square m×mm\times m matrix with singular value decomposition A=UΣVTA=U\Sigma V^T.

  1. If AA is invertible, find a singular value decomposition of A−1A^{-1}.

  2. What condition on the singular values must hold for AA to be invertible?

  3. How are the singular values of AA and the singular values of A−1A^{-1} related to one another?

  4. How are the right and left singular vectors of AA related to the right and left singular vectors of A−1A^{-1}?

2.5

Exercise 8

  1. If QQ is an orthogonal matrix, remember that QTQ=IQ^TQ=I. Explain why det⁡Q=±1\det Q = \pm 1.

  2. If A=UΣVTA=U\Sigma V^T is a singular value decomposition of a square matrix AA, explain why ∣det⁡A∣|\det A| is the product of the singular values of AA.

  3. What does this say about the singular values of AA if AA is invertible?

2.5

Exercise 9

If AA is a matrix and G=ATAG=A^TA its Gram matrix, remember that

  1. For a general matrix AA, explain why the eigenvalues of GG are nonnegative.

  2. Given a symmetric matrix AA having an eigenvalue λ\lambda, explain why λ2\lambda^2 is an eigenvalue of GG.

  3. If AA is symmetric, explain why the singular values of AA equal the absolute value of its eigenvalues: σj=∣λj∣\sigma_j = |\lambda_j|.

2.5

Exercise 10

Determine whether the following statements are true or false and explain your reasoning.

  1. If A=UΣVTA=U\Sigma V^T is a singular value decomposition of AA, then G=V(ΣTΣ)VTG=V(\Sigma^T\Sigma)V^T is an orthogonal diagonalization of its Gram matrix.

  2. If A=UΣVTA=U\Sigma V^T is a singular value decomposition of a rank 2 matrix AA, then v1\vvec_1 and v2\vvec_2 form an orthonormal basis for the column space Col(A)\col(A).

  3. If AA is a symmetric matrix, then its set of singular values is the same as its set of…

2.5

Exercise 11

Suppose that A=UΣVTA=U\Sigma V^T is a singular value decomposition of the m×nm\times n matrix AA. If σ1,…,σr\sigma_1,\ldots,\sigma_r are the nonzero singular values, the general form of the matrix Σ\Sigma is

  1. If you know that the columns of AA are linearly independent, what more can you say about the form of Σ\Sigma?

  2. If you know that the columns of AA span Rm\real^m, what more can you say about the form of Σ\Sigma?

  3. If you know that the columns of AA are linearly independent and span Rm\real^m, w…

2.5

Adapted from Understanding Linear Algebra © 2023 David Austin. Licensed under CC BY 4.0. Source: Understanding Linear Algebra. Imported at revision 3e6dced.