Chapter 12

Projections and Projection Matrices

[8/100%]
Show
LaTeX
Problem 12.1

Let Y\mathbf{Y} represent a matrix in a linear space V\mathcal{V}, let U\mathcal{U} and W\mathcal{W} represent subspaces of V\mathcal{V}, and take {X1,…,Xs}\{ \mathbf{X}_{1}, \ldots , \mathbf{X}_{s}\} to be a set of matrices that spans U\mathcal{U} and {Z1,…,Zt}\{ \mathbf{Z}_{1}, \ldots , \mathbf{Z}_{t}\} to be a set that spans W\mathcal{W}. Verify that Y⊥U\mathbf{Y} \perp \mathcal{U} if and only if Y⋅Xi=0\mathbf{Y} \cdot \mathbf{X}_{i} = 0 for i=1,…,si = 1, \ldots , s (i.e., that Y\mathbf{Y} is orthogonal to U\mathcal{U} if and only if Y\mathbf{Y} is orthogonal to each of the matrices X1,…,Xs\mathbf{X}_{1}, \ldots , \mathbf{X}_{s}); and, similarly, that U⊥W\mathcal{U} \perp \mathcal{W} if and only if Xi⋅Zj=0\mathbf{X}_{i} \cdot \mathbf{Z}_{j} = 0 for i=1,…,si = 1, \ldots , s and j=1,…,tj = 1, \ldots , t (i.e., that U\mathcal{U} is orthogonal to W\mathcal{W} if and only if each of the matrices X1,…,Xs\mathbf{X}_{1}, \ldots , \mathbf{X}_{s} is orthogonal to each of the matrices Z1,…,Zt\mathbf{Z}_{1}, \ldots , \mathbf{Z}_{t}).

?
Problem 12.2

Let U\mathcal{U} and V\mathcal{V} represent subspaces of Rm×n\mathbb {R}^{m \times n}. Show that if dim⁡(V)>dim⁡(U)\dim (\mathcal{V}) > \dim (\mathcal{U}), then V\mathcal{V} contains a nonnull matrix that is orthogonal to U\mathcal{U}.

?
Problem 12.3

Let U\mathcal{U} represent a subspace of the linear space Rm\mathbb {R}^{m} of all mm-dimensional column vectors. Take M\mathcal{M} to be the subspace of Rm×n\mathbb {R}^{m \times n} defined by W∈M\mathbf{W} \in \mathcal{M} if and only if W=(w1,…,wn)\mathbf{W} = (\mathbf{w}_{1}, \ldots , \mathbf{w}_{n}) for some vectors w1,…,wn\mathbf{w}_{1}, \ldots , \mathbf{w}_{n} in U\mathcal{U}. Let Z\mathbf{Z} represent the projection (with respect to the usual inner product) of an m×nm \times n matrix Y\mathbf{Y} on M\mathcal{M}, and let X\mathbf{X} represent any m×pm \times p matrix whose columns span U\mathcal{U}. Show that Z=XB∗\mathbf{Z} = \mathbf{X}\mathbf{B}^{*} for any solution B∗\mathbf{B}^{*} to the linear system

X′XB=X′Y(in B). \mathbf{X}'\mathbf{X}\mathbf{B} = \mathbf{X}'\mathbf{Y} \quad \text{(in }\mathbf{B}\text{)}.
?
Problem 12.4

The projection (with respect to the usual inner product) of an nn-dimensional column vector y\mathbf{y} on a subspace U\mathcal{U} of Rn\mathbb {R}^{n} in the special case where n=3n = 3, y=(3,−38/5,74/5)′\mathbf{y} = (3, -38/5, 74/5)' and U=sp⁡{x1,x2,x3}\mathcal{U} = \operatorname {sp}\{ \mathbf{x}_{1}, \mathbf{x}_{2}, \mathbf{x}_{3}\}, with

x1=(036),x2=(−224),x3=(−212), \mathbf{x}_{1} = \left(\begin{smallmatrix} 0 \\ 3 \\ 6 \end{smallmatrix}\right), \quad \mathbf{x}_{2} = \left(\begin{smallmatrix} -2 \\ 2 \\ 4 \end{smallmatrix}\right), \quad \mathbf{x}_{3} = \left(\begin{smallmatrix} -2 \\ 1 \\ 2 \end{smallmatrix}\right),

was determined to be the vector (3,22/5,44/5)′(3, 22/5, 44/5)' --- and it was observed that x1\mathbf{x}_{1} and x2\mathbf{x}_{2} are linearly independent and that x3=x2−(1/3)x1\mathbf{x}_{3} = \mathbf{x}_{2} - (1/3)\mathbf{x}_{1}, with the consequence that dim⁡(U)=2\dim (\mathcal{U}) = 2. Recompute the projection of y\mathbf{y} on U\mathcal{U} (in this special case) by taking X\mathbf{X} to be the 3×23 \times 2 matrix

(0−23264) \left(\begin{smallmatrix} 0 & -2 \\ 3 & 2 \\ 6 & 4 \end{smallmatrix}\right)

and carrying out the following two steps: (1) compute the solution to the normal equations X′Xb=X′y\mathbf{X}'\mathbf{X}\mathbf{b} = \mathbf{X}'\mathbf{y}; and (2) postmultiply X\mathbf{X} by the solution you computed in Step (1).

?
Problem 12.5

Let X\mathbf{X} represent any n×pn \times p matrix. If a p×np \times n matrix B∗\mathbf{B}^{*} is a solution to the linear system X′XB=X′\mathbf{X}'\mathbf{X}\mathbf{B} = \mathbf{X}' (in B\mathbf{B}), then B∗\mathbf{B}^{*} is a generalized inverse of X\mathbf{X} and XB∗\mathbf{X}\mathbf{B}^{*} is symmetric. Show that, conversely, if a p×np \times n matrix G\mathbf{G} is a generalized inverse of X\mathbf{X} and if XG\mathbf{X}\mathbf{G} is symmetric, then X′XG=X′\mathbf{X}'\mathbf{X}\mathbf{G} = \mathbf{X}' (i.e., G\mathbf{G} is a solution to X′XB=X′\mathbf{X}'\mathbf{X}\mathbf{B} = \mathbf{X}').

?
Problem 12.6

Using the result of Part (b) of Exercise 9.3 (or otherwise), show that, for any nonnull symmetric matrix A\mathbf{A},

PA=B(TB)−1T, \mathbf{P}_{\mathbf{A}} = \mathbf{B}(\mathbf{T}\mathbf{B})^{-1}\mathbf{T} ,

where B\mathbf{B} is any matrix of full column rank and T\mathbf{T} any matrix of full row rank such that A=BT\mathbf{A} = \mathbf{B}\mathbf{T}. (That TB\mathbf{T}\mathbf{B} is nonsingular follows from the result of Exercise 8.3.)

?
Problem 12.7

Let V\mathcal{V} represent a kk-dimensional subspace of the linear space Rn\mathbb {R}^{n} of all nn-dimensional column vectors. Take X\mathbf{X} to be any n×pn \times p matrix whose columns span V\mathcal{V}, let U\mathcal{U} represent a subspace of V\mathcal{V}, and define A\mathbf{A} to be the projection matrix for U\mathcal{U}. Show (1) that a matrix B\mathbf{B} (of dimensions n×nn \times n) is such that By\mathbf{B}\mathbf{y} is the projection of y\mathbf{y} on U\mathcal{U} for every y∈V\mathbf{y} \in \mathcal{V} if and only if B=A+Z∗\mathbf{B} = \mathbf{A} + \mathbf{Z}_{*} for some solution Z∗\mathbf{Z}_{*} to the homogeneous linear system X′Z=0\mathbf{X}'\mathbf{Z} = \mathbf{0} (in an n×nn \times n matrix Z\mathbf{Z}) and (2) that, unless k=nk = n, there is more than one matrix B\mathbf{B} such that By\mathbf{B}\mathbf{y} is the projection of y\mathbf{y} on U\mathcal{U} for every y∈V\mathbf{y} \in \mathcal{V}.

?
Problem 12.8

Let {A1,…,Ak}\{ \mathbf{A}_{1}, \ldots , \mathbf{A}_{k}\} represent a nonempty linearly independent set of matrices in a linear space V\mathcal{V}. And, define (as in Gram-Schmidt orthogonalization) kk nonnull orthogonal linear combinations, say B1,…,Bk\mathbf{B}_{1}, \ldots , \mathbf{B}_{k}, of A1,…,Ak\mathbf{A}_{1}, \ldots , \mathbf{A}_{k} as follows:

B1=A1,B2=A2−x12B1,⋮Bj=Aj−xj−1,jBj−1−⋯−x1jB1,⋮Bk=Ak−xk−1,kBk−1−⋯−x1kB1, \begin{align} \mathbf{B}_{1} & = \mathbf{A}_{1}, \\ \mathbf{B}_{2} & = \mathbf{A}_{2} - x_{12}\mathbf{B}_{1}, \\ & \vdots \\ \mathbf{B}_{j} & = \mathbf{A}_{j} - x_{j-1,j}\mathbf{B}_{j-1} - \cdots - x_{1j}\mathbf{B}_{1}, \\ & \vdots \\ \mathbf{B}_{k} & = \mathbf{A}_{k} - x_{k-1,k}\mathbf{B}_{k-1} - \cdots - x_{1k}\mathbf{B}_{1}, \end{align}

where (for i<j=1,…,ki < j = 1, \ldots , k)

xij=Aj⋅BiBi⋅Bi. x_{ij} = \frac{\mathbf{A}_{j} \cdot \mathbf{B}_{i}}{\mathbf{B}_{i} \cdot \mathbf{B}_{i}} .

Show that Bj\mathbf{B}_{j} is the (orthogonal) projection of Aj\mathbf{A}_{j} on some subspace Uj\mathcal{U}_{j} (of V\mathcal{V}) and describe Uj\mathcal{U}_{j} (j=2,…,kj = 2, \ldots , k).

?