Back to Understanding Linear Algebra

Finding orthogonal bases

Gram-Schmidt converts any basis into an orthonormal one, and recording the steps as it runs produces the QR factorization.

Understanding Linear Algebra, Finding orthogonal bases

Status
FlagsLoading your progress…

10 exercises

Activity 6.4.2

Suppose that WW is a three-dimensional subspace of R4\real^4 with basis: We can see that this basis is not orthogonal by noting that v1⋅v2=8\vvec_1\cdot\vvec_2 = 8. Our goal is to create an orthogonal basis w1\wvec_1, w2\wvec_2, and w3\wvec_3 for WW.

To begin, we declare that w1=v1\wvec_1=\vvec_1, and we call W1W_1 the line defined by w1\wvec_1.

  1. Find the vector v^2\vhat_2 that is the orthogonal projection of v2\vvec_2 onto W1W_1, the line defined by w1\wvec_1.

  2. Form the vector w2=v2−v^2\wvec_2 = \vvec_2-\vhat_2 and verify that it is orthogonal to w1\wvec_1.

  3. Expl…

2.5

Activity 6.4.4

Suppose that AA is the 4×34\times3 matrix whose columns are These vectors form a basis for WW, the subspace of R4\real^4 that we encountered in this activity. Since these vectors are the columns of AA, we have Col(A)=W\col(A) = W.

  1. When we implemented Gram-Schmidt, we first found an orthogonal basis w1\wvec_1, w2\wvec_2, and w3\wvec_3 using Use these expressions to write v1\vvec_1, v2\vvec_2, and v3\vvec_3 as linear combinations of w1\wvec_1, w2\wvec_2, and w3\wvec_3.

  2. We next normalized the orthogonal basis…

2.5

Exercise 1

Suppose that a subspace WW of R3\real^3 has a basis formed by

  1. Find an orthogonal basis for WW.

  2. Find an orthonormal basis for WW.

  3. Find the matrix PP that projects vectors orthogonally onto WW.

  4. Find the orthogonal projection of [34−2]\threevec34{-2} onto WW.

2.5

Exercise 2

Find the QRQR factorization of A=[47−2444]A=\begin{bmatrix} 4 & 7 \\ -2 & 4 \\ 4 & 4 \end{bmatrix}.

2.5

Exercise 3

Consider the basis of R3\real^3 given by the vectors

  1. Apply the Gram-Schmit orthogonalization algorithm to find an orthonormal basis u1\uvec_1, u2\uvec_2, u3\uvec_3 for R3\real^3.

  2. If AA is the 3×33\times3 whose columns are v1\vvec_1, v2\vvec_2, and v3\vvec_3, find the QRQR factorization of AA.

  3. Suppose that we want to solve the equation Ax=b=[−917]A\xvec=\bvec = \threevec{-9}17, which we can rewrite as QRx=bQR\xvec = \bvec.

    1. If we set y=Rx\yvec=R\xvec, the equation QRx=bQR\xvec=\bvec becomes Qy=bQ\yvec=\bvec. Explain how to solve the equation Qy=bQ\yvec=\bvec in a computationally efficient manner.

    2. Ex…

2.5

Exercise 4

Consider the vectors and the subspace WW of R5\real^5 that they span.

  1. Find an orthonormal basis for WW.

  2. Find the 5×55\times5 matrix that projects vectors orthogonally onto WW.

  3. Find b^\bhat, the orthogonal projection of b=[−83−128−4]\bvec=\fivevec{-8}3{-12}8{-4} onto WW.

  4. Express b^\bhat as a linear combination of v1\vvec_1, v2\vvec_2, and v3\vvec_3.

2.5

Exercise 5

Consider the set of vectors

  1. What happens when we apply the Gram-Schmit orthogonalization algorithm?

  2. Why does the algorithm fail to produce an orthogonal basis for R3\real^3?

2.5

Exercise 6

Suppose that AA is a matrix with linearly independent columns and having the factorization A=QRA=QR. Determine whether the following statements are true or false and explain your thinking.

  1. It follows that R=QTAR=Q^TA.

  2. The matrix RR is invertible.

  3. The product QTQQ^TQ projects vectors orthogonally onto Col(A)\col(A).

  4. The columns of QQ are an orthogonal basis for Col(A)\col(A).

  5. The orthogonal complement Col(A)⊥=Nul(QT)\col(A)^\perp = \nul(Q^T).

2.5

Exercise 7

Suppose we have the QRQR factorization A=QRA=QR, where AA is a 7×47\times 4 matrix.

  1. What is the shape of the product QQTQQ^T? Explain the significance of this product.

  2. What is the shape of the product QTQQ^TQ? Explain the significance of this product.

  3. What is the shape of the matrix RR?

  4. If RR is a diagonal matrix, what can you say about the columns of AA?

2.5

Exercise 8

Suppose we have the QRQR factorization A=QRA=QR where the columns of AA are a1,a2,…,an\avec_1,\avec_2,\ldots,\avec_n and the columns of RR are r1,r2,…,rn\rvec_1,\rvec_2,\ldots,\rvec_n.

  1. How can the matrix product ATAA^TA be expressed in terms of dot products?

  2. How can the matrix product RTRR^TR be expressed in terms of dot products?

  3. Explain why ATA=RTRA^TA=R^TR.

  4. Explain why the dot products ai⋅aj=ri⋅rj\avec_i\cdot\avec_j = \rvec_i\cdot\rvec_j.

2.5

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