Skip to main content

Recitation 10B

Problem 1: Let y=[481],u1=[2/31/32/3],u2=[2/32/31/3]y=\begin{bmatrix}4\\8\\1\end{bmatrix}, u_1=\begin{bmatrix}2/3\\1/3\\2/3\end{bmatrix}, u_2=\begin{bmatrix}-2/3\\2/3\\1/3\end{bmatrix}, and W=Span{u1,u2}W=\mathrm{Span}\{u_1, u_2\}. Let U=[u1u2]U=\begin{bmatrix}u_1 & u_2\end{bmatrix}. Compute UTUU^TU and UUTUU^T. Compute projWy\mathrm{proj}_Wy and (UUT)y(UU^T)y.

Solution: Compute UTU=[u1Tu2T][u1u2]=[u1Tu1u1Tu2u2Tu1u2Tu2]=[1001]U^TU=\begin{bmatrix}u_1^T\\u_2^T\end{bmatrix}\begin{bmatrix}u1 & u_2\end{bmatrix}=\begin{bmatrix}u_1^Tu_1 & u_1^Tu_2 \\u_2^Tu_1 & u_2^Tu_2\end{bmatrix}=\begin{bmatrix}1 & 0 \\ 0 & 1\end{bmatrix}. This implies that u1,u2u_1, u_2 form an orthornomal basis of WW. On the other hand, UUT=[2/32/31/32/32/31/3][2/31/32/32/32/31/3]=19[822254245]UU^T = \begin{bmatrix}2/3 & -2/3\\ 1/3 & 2/3\\ 2/3 & 1/3\end{bmatrix}\begin{bmatrix}2/3 & 1/3 & 2/3 \\ -2/3 & 2/3 & 1/3\end{bmatrix}=\frac{1}{9}\begin{bmatrix}8 & -2 & 2 \\ -2 & 5 & 4 \\ 2 & 4 & 5\end{bmatrix}. Because u1,u2u_1, u_2 form an orthonormal basis, the projection of yy onto WW is given by (yu1)u1+(yu2)u2=[245](y\cdot u_1)u_1 + (y\cdot u_2)u_2 = \begin{bmatrix}2 \\ 4 \\ 5\end{bmatrix}. Another way to compute the projection is to multiply UUTUU^T with yy. In other words, (UUT)y(UU^T)y gives the same vector [245]\begin{bmatrix}2 \\ 4 \\ 5\end{bmatrix}.

Problem 2: Find an orthogonal basis for the column space of the matrix [125114143147121]\begin{bmatrix}1 & 2 & 5\\-1 & 1 & -4\\-1 & 4 & -3\\1 & -4 & 7\\1 & 2 & 1\end{bmatrix}.

Solution: Denote its column vectors by x1,x2,x3x_1, x_2, x_3. Using Gram-Schimidt, we get v1=x1v_1 = x_1, v2=x2x2v1v1v1v1=[30333][10111]v_2 = x_2 - \frac{x_2\cdot v_1}{v_1\cdot v_1}v_1 = \begin{bmatrix}3 \\ 0 \\ 3 \\ -3 \\ 3\end{bmatrix}\sim \begin{bmatrix}1 \\ 0 \\ 1 \\ -1 \\ 1\end{bmatrix} and v3=x3x3v1v1v1v1x3v2v2v2v2=[20222][10111]v_3 = x_3 - \frac{x_3\cdot v_1}{v_1\cdot v_1}v_1 - \frac{x_3\cdot v_2}{v_2\cdot v_2}v_2 = \begin{bmatrix}2 \\ 0 \\ 2 \\ 2 \\ -2\end{bmatrix}\sim \begin{bmatrix}1 \\ 0 \\ 1 \\ 1 \\ -1\end{bmatrix}.