Chapter 14

Linear, Bilinear, and Quadratic Forms

[51/100%]
Show
LaTeX
Problem 14.1

Show that a symmetric bilinear form x′Ayx'Ay (in nn-dimensional vectors xx and yy) can be expressed in terms of the corresponding quadratic form, that is, the quadratic form whose matrix is AA. Do so by verifying that

x′Ay=(1/2)[(x+y)′A(x+y)−x′Ax−y′Ay]. x'Ay = (1/2)\left[(x+y)'A(x+y) - x'Ax - y'Ay\right] .
?
Problem 14.2

Show that corresponding to any quadratic form x′Axx'Ax (in the nn-dimensional vector xx) there exists a unique upper triangular matrix BB such that x′Axx'Ax and x′Bxx'Bx are identically equal, and express the elements of BB in terms of the elements of AA.

?
Problem 14.3

Show, by example, that the sum of two positive semidefinite matrices can be positive definite.

?
Problem 14.4

Show, via an example, that there exist (nonsymmetric) nonsingular positive semidefinite matrices.

?
Problem 14.5

Show, by example, that there exist an n×nn \times n positive semidefinite matrix AA and an n×mn \times m matrix PP (where m<nm < n) such that P′APP'AP is positive definite.

?
Problem 14.6

For an n×nn \times n matrix AA and an n×mn \times m matrix PP, it is the case that (1) if AA is nonnegative definite, then P′APP'AP is nonnegative definite; (2) if AA is nonnegative definite and rk⁡(P)<m\operatorname {rk}\left(P\right) < m, then P′APP'AP is positive semidefinite; and (3) if AA is positive definite and rk⁡(P)=m\operatorname {rk}\left(P\right) = m, then P′APP'AP is positive definite. Convert these results, which are for nonnegative definite (positive definite or positive semidefinite) matrices, into equivalent results for nonpositive definite matrices.

?
Problem 14.7

Let {X1,…,Xr}\{ X_1, \ldots , X_r\} represent a set of matrices from a linear space V\mathcal{V}. And, let A={aij}A = \{ a_{ij}\} represent the r×rr \times r matrix whose ijijth element is ⟨Xi,Xj⟩\left\langle X_i, X_j \right\rangle -- this matrix is referred to as the Gram matrix (or the Gramian) of the set {X1,…,Xr}\{ X_1, \ldots , X_r\} and its determinant is referred to as the Gramian (or the Gram determinant) of {X1,…,Xr}\{ X_1, \ldots , X_r\}.

?
(a)

Show that AA is symmetric and nonnegative definite.

(b)

Show that X1,…,XrX_1, \ldots , X_r are linearly independent if and only if AA is nonsingular.

Problem 14.8

Let A={aij}A = \{ a_{ij}\} represent an n×nn \times n symmetric positive definite matrix, and let B={bij}=A−1B = \{ b_{ij}\} = A^{-1}. Show that, for i=1,…,ni = 1, \ldots , n,

bii≥1/aii, b_{ii} \geq 1/a_{ii} ,

with equality holding if and only if, for all j≠ij \neq i, aij=0a_{ij} = 0.

?
Problem 14.9

Let AA represent an m×nm \times n matrix and DD a diagonal matrix such that A=PDQA = PDQ for some matrix PP of full column rank and some matrix QQ of full row rank. Show that rk⁡(A)\operatorname {rk}\left(A\right) equals the number of nonzero diagonal elements in DD.

?
Problem 14.10

Let AA represent an n×nn \times n symmetric idempotent matrix and VV an n×nn \times n symmetric positive definite matrix. Show that rk⁡(AVA)=tr⁡(A)\operatorname {rk}\left(AVA\right) = \operatorname {tr}(A).

?
Problem 14.11

Show that if an n×nn \times n matrix AA is such that x′Ax≠0x'Ax \neq 0 for every n×1n \times 1 nonnull vector xx, then AA is either positive definite or negative definite.

?
Problem 14.12
?
(a)

Let AA represent an n×nn \times n symmetric matrix of rank rr. Take PP to be an n×nn \times n nonsingular matrix and DD an n×nn \times n diagonal matrix such that A=P′DPA = P'DP -- the existence of such matrices is guaranteed. The number, say mm, of diagonal elements of DD that are positive is called the index of inertia of AA (or of the quadratic form x′Axx'Ax whose matrix is AA). Show that the index of inertia is well-defined in the sense that mm does not vary with the choice of PP or DD. That is, show that, if P1P_1 and P2P_2 are nonsingular matrices and D1D_1 and D2D_2 diagonal matrices such that A=P1′D1P1=P2′D2P2A = P_1'D_1P_1 = P_2'D_2P_2, then D2D_2 contains the same number of positive diagonal elements as D1D_1. Show also that the number of diagonal elements of DD that are negative equals r−mr - m.

(b)

Let AA represent an n×nn \times n symmetric matrix. Show that A=P′diag⁡(Im,−Ir−m,0)PA = P'\operatorname {diag}(I_m, -I_{r-m}, 0)P for some n×nn \times n nonsingular matrix PP and some nonnegative integers mm and rr. Show further that mm equals the index of inertia of the matrix AA and that r=rk⁡(A)r = \operatorname {rk}\left(A\right).

(c)

An n×nn \times n symmetric matrix BB is said to be congruent to an n×nn \times n symmetric matrix AA if there exists an n×nn \times n nonsingular matrix PP such that B=P′APB = P'AP. (If BB is congruent to AA, then clearly AA is congruent to BB.) Show that BB is congruent to AA if and only if BB has the same rank and the same index of inertia as AA. This result is called Sylvester's law of inertia, after James Joseph Sylvester (1814--1897).

(d)

Let AA represent an n×nn \times n symmetric matrix of rank rr and with index of inertia mm. Show that AA is nonnegative definite if and only if m=rm = r and is positive definite if and only if m=r=nm = r = n.

Problem 14.13

Let AA represent an n×nn \times n symmetric nonnegative definite matrix of rank rr. Then, there exists an n×rn \times r matrix BB (of rank rr) such that A=BB′A = BB'. Let XX represent any n×mn \times m matrix (where m≥rm \geq r) such that A=XX′A = XX'.

?
(a)

Show that X=PBXX = P_BX.

(b)

Show that X=(B,0)QX = (B, 0)Q for some orthogonal matrix QQ.

Problem 14.14

Show that if a symmetric matrix AA has a nonnegative definite generalized inverse, then AA is nonnegative definite.

?
Problem 14.15

Suppose that an n×nn \times n matrix AA has an LDU decomposition, say A=LDUA = LDU, and let d1,d2,…,dnd_1, d_2, \ldots , d_n represent the diagonal elements of the diagonal matrix DD. Show that

∣A∣=d1d2⋯dn. \left|A\right| = d_1d_2 \cdots d_n .
?
Problem 14.16
?
(a)

Suppose that an n×nn \times n matrix AA (where n≥2n \geq 2) has a unique LDU decomposition, say A=LDUA = LDU, and let d1,d2,…,dnd_1, d_2, \ldots , d_n represent the first, second, , nnth diagonal elements of DD. Show that di≠0d_i \neq 0 (i=1,2,…,n−1)(i = 1, 2, \ldots , n-1) and that dn≠0d_n \neq 0 if and only if AA is nonsingular.

(b)

Suppose that an n×nn \times n (symmetric) matrix AA (where n≥2n \geq 2) has a unique U′DUU'DU decomposition, say A=U′DUA = U'DU, and let d1,d2,…,dnd_1, d_2, \ldots , d_n represent the first, second, , nnth diagonal elements of DD. Show that di≠0d_i \neq 0 (i=1,2,…,n−1)(i = 1, 2, \ldots , n-1) and that dn≠0d_n \neq 0 if and only if AA is nonsingular.

Problem 14.17

Suppose that an n×nn \times n (symmetric) matrix AA has a unique U′DUU'DU decomposition, say A=U′DUA = U'DU. Use the result of Part (b) of Exercise 14.16 to show that AA has no LDU decompositions other than A=U′DUA = U'DU.

?
Problem 14.18

Show that if a nonsingular matrix has an LDU decomposition, then that decomposition is unique.

?
Problem 14.19

Let AA represent an n×nn \times n matrix (where n≥2n \geq 2). By for instance using the results of Exercises 14.16, 14.17, and 14.18, show that if AA has a unique LDU decomposition or (in the special case where AA is symmetric) a unique U′DUU'DU decomposition, then the leading principal submatrices (of AA) of orders 1,2,…,n−11, 2, \ldots , n-1 are nonsingular and have unique LDU decompositions.

?
Problem 14.20
?
(a)

Let A={aij}A = \{ a_{ij}\} represent an m×nm \times n nonnull matrix of rank rr. Show that there exist an m×mm \times m permutation matrix PP and an n×nn \times n permutation matrix QQ such that

PAQ=(B11B12B21B22), PAQ = \left(\begin{smallmatrix} B_{11} & B_{12} \\ B_{21} & B_{22} \end{smallmatrix}\right) ,

where B11B_{11} is an r×rr \times r nonsingular matrix whose leading principal submatrices (of orders 1,2,…,r−11, 2, \ldots , r-1) are nonsingular.

(b)

Let B=(B11B12B21B22)B = \left(\begin{smallmatrix} B_{11} & B_{12} \\ B_{21} & B_{22} \end{smallmatrix}\right) represent any m×nm \times n nonnull matrix of rank rr such that B11B_{11} is an r×rr \times r nonsingular matrix whose leading principal submatrices (of orders 1,2,…,r−11, 2, \ldots , r-1) are nonsingular. Show that there exists a unique decomposition of BB of the form

B=(L1L2)D(U1,U2), B = \left(\begin{smallmatrix} L_1 \\ L_2 \end{smallmatrix}\right) D (U_1, U_2) ,

where L1L_1 is an r×rr \times r unit lower triangular matrix, U1U_1 is an r×rr \times r unit upper triangular matrix, and DD is an r×rr \times r diagonal matrix. Show further that this decomposition is such that B11=L1DU1B_{11} = L_1DU_1 is the unique LDU decomposition of B11B_{11}, DD is nonsingular, L2=B21U1−1D−1L_2 = B_{21}U_1^{-1}D^{-1}, and U2=D−1L1−1B12U_2 = D^{-1}L_1^{-1}B_{12}.

Problem 14.21

Show, by example, that there exist n×nn \times n (nonsymmetric) positive semidefinite matrices that do not have LDU decompositions.

?
Problem 14.22

Let AA represent an n×nn \times n nonnegative definite (possibly nonsymmetric) matrix that has an LDU decomposition, say A=LDUA = LDU. Show that the diagonal elements of the diagonal matrix DD are nonnegative.

?
Problem 14.23

Let AA represent an m×km \times k matrix of full column rank. And, let A=QRA = QR represent the QR decomposition of AA; that is, let QQ represent the unique m×km \times k matrix whose columns are orthonormal with respect to the usual inner product and let RR represent the unique k×kk \times k upper triangular matrix with positive diagonal elements such that A=QRA = QR. Show that A′A=R′RA'A = R'R (so that A′A=R′RA'A = R'R is the Cholesky decomposition of A′AA'A).

?
Problem 14.24

Let AA represent an m×km \times k matrix of rank rr (where rr is possibly less than kk). Consider the decomposition A=QR1A = QR_1, where QQ is an m×rm \times r matrix with orthonormal columns and R1R_1 is an r×kr \times k submatrix whose rows are the rr nonnull rows of a k×kk \times k upper triangular matrix RR having rr positive diagonal elements and n−rn - r null rows. (Such a decomposition can be obtained by using the results of Exercise 6.4 -- refer to Exercise 6.5.) Generalize the result of Exercise 14.23 by showing that if the inner product with respect to which the columns of QQ are orthonormal is the usual inner product, then A′A=R′RA'A = R'R (so that A′A=R′RA'A = R'R is the Cholesky decomposition of A′AA'A).

?
Problem 14.25

Let A={aij}A = \{ a_{ij}\} represent an n×nn \times n matrix that has an LDU decomposition, say A=LDUA = LDU. And, define G=U−1D−L−1G = U^{-1}D^{-}L^{-1} (which is a generalized inverse of AA).

?
(a)

Show that

G=D−L−1+(I−U)G=U−1D−+G(I−L). G = D^{-}L^{-1} + (I-U)G = U^{-1}D^{-} + G(I-L) .
(b)

For i=1,…,ni = 1, \ldots , n, let did_i represent the iith diagonal element of the diagonal matrix DD; and, for i,j=1,…,ni, j = 1, \ldots , n, let ℓij,uij\ell_{ij}, u_{ij}, and gijg_{ij} represent the ijijth elements of L,UL, U, and GG, respectively. Take D−=diag⁡(d1∗,…,dn∗)D^{-} = \operatorname {diag}(d_1^{*}, \ldots , d_n^{*}), where di∗=1/did_i^{*} = 1/d_i, if di≠0d_i \neq 0, and di∗d_i^{*} is an arbitrary scalar, if di=0d_i = 0. Show that

gii=di∗−∑k=i+1nuikgki=di∗−∑k=i+1ngikℓki(E.1) g_{ii} = d_i^{*} - \sum _{k=i+1}^{n} u_{ik}g_{ki} = d_i^{*} - \sum _{k=i+1}^{n} g_{ik}\ell _{ki} \tag {E.1}

and that

g_{ij} = \begin{cases} -\sum _{k=j+1}^{n} g_{ik}\ell _{kj}, & \text{for } j < i , \tag {E.2a} \\ -\sum _{k=i+1}^{n} u_{ik}g_{kj}, & \text{for } j > i \tag {E.2b} \end{cases}

(where the degenerate sums ∑k=n+1ngikℓki\sum_{k=n+1}^{n} g_{ik}\ell_{ki} and ∑k=n+1nuikgki\sum_{k=n+1}^{n} u_{ik}g_{ki} are to be interpreted as 00).

(c)

Devise a recursive procedure that uses the formulas from Part (b) to generate a generalized inverse of AA.

Problem 14.26

Verify that a principal submatrix of a skew-symmetric matrix is skew-symmetric.

?
Problem 14.27
?
(a)

Show that the sum of skew-symmetric matrices is skew-symmetric.

(b)

Show that the sum A1+A2+⋯+AkA_1+A_2+\cdots +A_k of n×nn \times n nonnegative definite matrices A1,A2,…,AkA_1, A_2, \ldots , A_k is skew-symmetric if and only if A1,A2,…,AkA_1, A_2, \ldots , A_k are skew-symmetric.

(c)

Show that the sum A1+A2+⋯+AkA_1+A_2+\cdots +A_k of n×nn \times n symmetric nonnegative definite matrices A1,A2,…,AkA_1, A_2, \ldots , A_k is a null matrix if and only if A1,A2,…,AkA_1, A_2, \ldots , A_k are null matrices.

Problem 14.28
?
(a)

Let A1,A2,…,AkA_1, A_2, \ldots , A_k represent n×nn \times n nonnegative definite matrices. Show that tr⁡(∑i=1kAi)≥0\operatorname {tr}(\sum_{i=1}^{k} A_i) \geq 0, with equality holding if and only if ∑i=1kAi\sum_{i=1}^{k} A_i is skew-symmetric or equivalently if and only if A1,A2,…,AkA_1, A_2, \ldots , A_k are skew-symmetric.

(b)

Let A1,A2,…,AkA_1, A_2, \ldots , A_k represent n×nn \times n symmetric nonnegative definite matrices. Show that tr⁡(∑i=1kAi)≥0\operatorname {tr}(\sum_{i=1}^{k} A_i) \geq 0, with equality holding if and only if ∑i=1kAi=0\sum_{i=1}^{k} A_i = 0 or equivalently if and only if A1,A2,…,AkA_1, A_2, \ldots , A_k are null matrices.

Remark.
?

That ∑i=1kAi\sum_{i=1}^{k} A_i being skew-symmetric is equivalent to A1,A2,…,AkA_1, A_2, \ldots , A_k being skew-symmetric is the result of Part (b) of Exercise 14.27.

Problem 14.29

Show, via an example, that (for n>1n > 1) there exist n×nn \times n (nonsymmetric) positive definite matrices AA and BB such that tr⁡(AB)<0\operatorname {tr}(AB) < 0.

?
Problem 14.30
?
(a)

Show, via an example, that (for n>1n > 1) there exist n×nn \times n symmetric positive definite matrices AA and BB such that the product ABAB has one or more negative diagonal elements (and hence such that ABAB is not nonnegative definite).

(b)

Show, however, that the product of two n×nn \times n symmetric positive definite matrices cannot be nonpositive definite.

Problem 14.31

Let A={aij}A = \{ a_{ij}\} and B={bij}B = \{ b_{ij}\} represent n×nn \times n matrices, and take CC to be the n×nn \times n matrix whose ijijth element cij=aijbijc_{ij} = a_{ij}b_{ij} is the product of the ijijth elements of AA and BB -- CC is the so-called Hadamard product of AA and BB. Show that if AA is nonnegative definite and BB is symmetric nonnegative definite, then CC is nonnegative definite. Show further that if AA is positive definite and BB is symmetric positive definite, then CC is positive definite.

?
Remark.
?

Taking x=(x1,…,xn)′x = (x_1, \ldots , x_n)' to be an arbitrary n×1n \times 1 vector and F=(f1,…,fn)F = (f_1, \ldots , f_n) to be a matrix such that B=F′FB = F'F, begin by showing that x′Cx=tr⁡(AH)x'Cx = \operatorname {tr}(AH), where H=G′GH = G'G with G=(x1f1,…,xnfn)G = (x_1f_1, \ldots , x_nf_n).

Problem 14.32

Let A1,A2,…,AkA_1, A_2, \ldots , A_k and B1,B2,…,BkB_1, B_2, \ldots , B_k represent n×nn \times n symmetric nonnegative definite matrices. Show that tr⁡(∑i=1kAiBi)≥0\operatorname {tr}(\sum_{i=1}^{k} A_iB_i) \geq 0, with equality holding if and only if, for i=1,2,…,ki = 1, 2, \ldots , k, AiBi=0A_iB_i = 0, thereby generalizing the results of Part (b) of Exercise 14.28.

?
Problem 14.33

Let AA represent a symmetric nonnegative definite matrix that has been partitioned as

A=(TUVW) A = \left(\begin{smallmatrix} T & U \\ V & W \end{smallmatrix}\right)

where TT (and hence WW) is square. Show that VT−UVT^{-}U and UW−VUW^{-}V are symmetric and nonnegative definite.

?
Problem 14.34

Show, via an example, that there exists an (m+n)×(m+n)(m+n) \times (m+n) (nonsymmetric) positive semidefinite matrix AA of the form A=(TUVW)A = \left(\begin{smallmatrix} T & U \\ V & W \end{smallmatrix}\right), where TT is of dimensions m×mm \times m, WW of dimensions n×nn \times n, UU of dimensions m×nm \times n, and VV of dimensions n×mn \times m, for which C(U)⊄C(T)C(U) \not\subset C(T) and/or R(V)⊄R(T)R(V) \not\subset R(T), the expression rk⁡(T)+rk⁡(W−VT−U)\operatorname {rk}\left(T\right) + \operatorname {rk}\left(W-VT^{-}U\right) does not necessarily equal rk⁡(A)\operatorname {rk}\left(A\right), and the formula

(T−+T−UQ−VT−−T−UQ−−Q−VT−Q−),(*) \left(\begin{smallmatrix} T^{-} + T^{-}UQ^{-}VT^{-} & -T^{-}UQ^{-} \\ -Q^{-}VT^{-} & Q^{-} \end{smallmatrix}\right) , \tag {*}

where Q=W−VT−UQ = W - VT^{-}U, does not necessarily give a generalized inverse of AA.

?
Problem 14.35

Show, via an example, that there exists an (m+n)×(m+n)(m+n) \times (m+n) symmetric partitioned matrix AA of the form A=(TUU′W)A = \left(\begin{smallmatrix} T & U \\ U' & W \end{smallmatrix}\right), where TT is of dimensions m×mm \times m, UU of dimensions m×nm \times n, and WW of dimensions n×nn \times n, such that TT is nonnegative definite and (depending on the choice of T−T^{-}) the Schur complement W−U′T−UW - U'T^{-}U of TT relative to T−T^{-} is nonnegative definite, but AA is not nonnegative definite.

?
Problem 14.36

An n×nn \times n matrix A={aij}A = \{ a_{ij}\} is said to be diagonally dominant if, for i=1,2,…,ni = 1, 2, \ldots , n, ∣aii∣>∑j=1(j≠i)n∣aij∣\left|a_{ii}\right| > \sum_{j=1(j \neq i)}^{n} \left|a_{ij}\right|. (In the degenerate special case where n=1n = 1, AA is said to be diagonally dominant if it is nonnull).

?
(a)

Show that a principal submatrix of a diagonally dominant matrix is diagonally dominant.

(b)

Let A={aij}A = \{ a_{ij}\} represent an n×nn \times n diagonally dominant matrix, partition AA as A=(A11ab′ann)A = \left(\begin{smallmatrix} A_{11} & a \\ b' & a_{nn} \end{smallmatrix}\right) [so that A11A_{11} is of dimensions (n−1)×(n−1)(n-1) \times (n-1)], and let C=A11−(1/ann)ab′C = A_{11} - (1/a_{nn})ab' represent the Schur complement of anna_{nn}. Show that CC is diagonally dominant.

(c)

Show that a diagonally dominant matrix is nonsingular.

(d)

Show that a diagonally dominant matrix has a unique LDU decomposition.

(e)

Let A={aij}A = \{ a_{ij}\} represent an n×nn \times n symmetric matrix. Show that if AA is diagonally dominant and if the diagonal elements a11,a22,…,anna_{11}, a_{22}, \ldots , a_{nn} of AA are all positive, then AA is positive definite.

Problem 14.37

Let A={aij}A = \{ a_{ij}\} represent an n×nn \times n symmetric positive definite matrix. Show that det⁡(A)≤∏i=1naii\det (A) \leq \prod_{i=1}^{n} a_{ii}, with equality holding if and only if AA is diagonal.

?
Problem 14.38

Let A=(abcd)A = \left(\begin{smallmatrix} a & b \\ c & d \end{smallmatrix}\right), where a,b,ca, b, c, and dd are scalars.

?
(a)

Show that AA is positive definite if and only if a>0a > 0, d>0d > 0, and ∣b+c∣/2<ad\left|b+c\right|/2 < \sqrt{ad}.

(b)

Show that, in the special case where AA is symmetric (i.e., where c=bc = b), AA is positive definite if and only if a>0a > 0, d>0d > 0, and ∣b∣<ad\left|b\right| < \sqrt{ad}.

Problem 14.39

By, for example, making use of the result of Exercise 14.38, show that if an n×nn \times n matrix A={aij}A = \{ a_{ij}\} is symmetric positive definite, then, for j≠i=1,…,nj \neq i = 1, \ldots , n,

∣aij∣<aiiajj≤max⁡(aii,ajj). \left|a_{ij}\right| < \sqrt{a_{ii}a_{jj}} \leq \max (a_{ii}, a_{jj}) .
?
Problem 14.40

Show, by example, that it is possible for the determinants of both leading principal submatrices of a 2×22 \times 2 symmetric matrix to be nonnegative without the matrix being nonnegative definite and that, for n≥3n \geq 3, it is possible for the determinants of all nn leading principal submatrices of an n×nn \times n symmetric matrix to be nonnegative and for the matrix to be nonsingular without the matrix being nonnegative definite.

?
Problem 14.41

Let V\mathcal{V} represent a subspace of Rn×1\mathbb {R}^{n \times 1} of dimension rr (where r≥1r \geq 1). Take B=(b1,b2,…,br)B = (b_1, b_2, \ldots , b_r) to be any n×rn \times r matrix whose columns b1,b2,…,brb_1, b_2, \ldots , b_r form a basis for V\mathcal{V}, and let LL represent any left inverse of BB. Let gg represent a function that assigns the value x∗yx * y to an arbitrary pair of vectors xx and yy in V\mathcal{V}.

?
(a)

Let ff represent an arbitrary inner product for Rr×1\mathbb {R}^{r \times 1}, and denote by s⋅ts \cdot t the value assigned by ff to an arbitrary pair of rr-dimensional vectors ss and tt. Show that gg is an inner product (for V\mathcal{V}) if and only if there exists an ff such that (for all xx and yy in V\mathcal{V})

x∗y=(Lx)⋅(Ly). x * y = (Lx) \cdot (Ly) .
(b)

Show that gg is an inner product (for V\mathcal{V}) if and only if there exists an r×rr \times r symmetric positive definite matrix WW such that (for all xx and yy in V\mathcal{V})

x∗y=x′L′WLy. x * y = x'L'WLy .
(c)

Show that gg is an inner product (for V\mathcal{V}) if and only if there exists an n×nn \times n symmetric positive definite matrix WW such that (for all xx and yy in V\mathcal{V})

x∗y=x′Wy. x * y = x'Wy .
Problem 14.42

Let V\mathcal{V} represent a linear space of m×nm \times n matrices, and let A⋅BA \cdot B represent the value assigned by a quasi-inner product to any pair of matrices AA and BB (in V\mathcal{V}). Show that the set

U={A∈V:A⋅A=0}, \mathcal{U} = \{ A \in \mathcal{V} : A \cdot A = 0\} ,

which comprises every matrix in V\mathcal{V} with a zero quasi norm, is a linear space.

?
Problem 14.43

Let WW represent an m×mm \times m symmetric positive definite matrix and VV an n×nn \times n symmetric positive definite matrix.

?
(a)

Show that the function that assigns the value tr⁡(A′WBV)\operatorname {tr}(A'WBV) to an arbitrary pair of m×nm \times n matrices AA and BB qualifies as an inner product for the linear space Rm×n\mathbb {R}^{m \times n}.

(b)

Show that the function that assigns the value tr⁡(A′WB)\operatorname {tr}(A'WB) to an arbitrary pair of m×nm \times n matrices AA and BB qualifies as an inner product for Rm×n\mathbb {R}^{m \times n}.

(c)

Show that the function that assigns the value tr⁡(A′WBW)\operatorname {tr}(A'WBW) to an arbitrary pair of m×mm \times m matrices AA and BB qualifies as an inner product for Rm×m\mathbb {R}^{m \times m}.

Problem 14.44

Let AA represent a q×pq \times p matrix, BB a p×np \times n matrix, and CC an m×qm \times q matrix. Show that

?
(a)

CAB(CAB)−C=CCAB(CAB)^{-}C = C if and only if rk⁡(CAB)=rk⁡(C)\operatorname {rk}\left(CAB\right) = \operatorname {rk}\left(C\right), and

(b)

B(CAB)−CAB=BB(CAB)^{-}CAB = B if and only if rk⁡(CAB)=rk⁡(B)\operatorname {rk}\left(CAB\right) = \operatorname {rk}\left(B\right).

Problem 14.45

Let U\mathcal{U} represent a subspace of Rn×1\mathbb {R}^{n \times 1}, let XX represent an n×pn \times p matrix whose columns span U\mathcal{U}, and let WW and VV represent n×nn \times n symmetric positive definite matrices. Show that each of the following two conditions is necessary and sufficient for the projection PX,WyP_{X,W}y of yy on U\mathcal{U} with respect to WW to be the same (for every yy in Rn\mathbb {R}^{n}) as the projection PX,VyP_{X,V}y of yy on U\mathcal{U} with respect to VV:

?
(a)

V=PX,W′VPX,W+(I−PX,W)′V(I−PX,W)V = P_{X,W}'VP_{X,W} + (I-P_{X,W})'V(I-P_{X,W});

(b)

there exist a scalar cc, a p×pp \times p matrix KK, and an n×nn \times n matrix HH such that

V=cW+WXKX′W+(I−PX,W)′H(I−PX,W). V = cW + WXKX'W + (I-P_{X,W})'H(I-P_{X,W}) .
Problem 14.46

Let yy represent an nn-dimensional column vector, let U\mathcal{U} represent a subspace of Rn×1\mathbb {R}^{n \times 1}, and let XX represent an n×pn \times p matrix whose columns span U\mathcal{U}. ``Recall'' that, for any m×nm \times n matrix LL, for any subspace W\mathcal{W} of Rn×1\mathbb {R}^{n \times 1}, and for V={v∈Rm:v=Lx for some x∈W}\mathcal{V} = \{ v \in \mathbb {R}^{m} : v = Lx \text{ for some } x \in \mathcal{W}\},

x⊥HW  ⟺  Lx⊥IV,(*) x \perp _H \mathcal{W} \iff Lx \perp _I \mathcal{V} , \tag {*}

where H=L′LH = L'L (and where xx represents an arbitrary n×1n \times 1 vector), and

C(LZ)=V,(**) C(LZ) = \mathcal{V} , \tag {**}

where ZZ is any n×qn \times q matrix whose columns span W\mathcal{W}. Use results (*) and (**) to show that, for any n×nn \times n symmetric nonnegative definite matrix WW, y⊥WUy \perp_W \mathcal{U} if and only if X′Wy=0X'Wy = 0.

?
Problem 14.47

Let WW represent an n×nn \times n symmetric nonnegative definite matrix.

?
(a)

Show that, for any n×pn \times p matrix XX and any n×qn \times q matrix UU such that C(U)⊂C(X)C(U) \subset C(X),

  1. WPX,WU=WUWP_{X,W}U = WU, and U′WPX,W=U′WU'WP_{X,W} = U'W;

  2. PU,WPX,W=PU,WP_{U,W}P_{X,W} = P_{U,W}, and PX,W′WPU,W=WPX,WPU,W=WPU,WP_{X,W}'WP_{U,W} = WP_{X,W}P_{U,W} = WP_{U,W}.

(b)

Show that, for any n×pn \times p matrix XX and any n×qn \times q matrix UU such that C(U)=C(X)C(U) = C(X),

WPU,W=WPX,W. WP_{U,W} = WP_{X,W} .
Problem 14.48

Let U\mathcal{U} represent a subspace of Rn×1\mathbb {R}^{n \times 1}, let AA represent an n×nn \times n matrix and WW an n×nn \times n symmetric nonnegative definite matrix, and let XX represent any n×pn \times p matrix whose columns span U\mathcal{U}.

?
(a)

Show that AA is a projection matrix for U\mathcal{U} with respect to WW if and only if

A=PX,W+(I−PX,W)XK A = P_{X,W} + (I-P_{X,W})XK

for some p×np \times n matrix KK.

(b)

Show that if AA is a projection matrix for U\mathcal{U} with respect to WW, then WA=WPX,WWA = WP_{X,W}.

Problem 14.49

Let AA represent an n×nn \times n matrix and WW an n×nn \times n symmetric nonnegative definite matrix.

?
(a)

Show (by, e.g., using the results of Exercise 14.48) that if A′WA=WAA'WA = WA [or, equivalently, if (I−A)′WA=0(I-A)'WA = 0], then AA is a projection matrix with respect to WW, and in particular AA is a projection matrix for C(A)C(A) with respect to WW, and, conversely, show that if AA is a projection matrix with respect to WW, then A′WA=WAA'WA = WA.

(b)

Show that if AA is a projection matrix with respect to WW, then in particular AA is a projection matrix for C(A)C(A) with respect to WW.

(c)

Show that AA is a projection matrix with respect to WW if and only if WAWA is symmetric and WA2=WAWA^2 = WA.

Problem 14.50

Let U\mathcal{U} represent a subspace of Rn×1\mathbb {R}^{n \times 1}, let XX represent an n×pn \times p matrix whose columns span U\mathcal{U}, and let WW and VV represent n×nn \times n symmetric nonnegative definite matrices. Show (by, e.g., making use of the result of Exercise 14.46) that each of the following two conditions is necessary and sufficient for every projection of yy on U\mathcal{U} with respect to WW to be a projection (for every yy in Rn\mathbb {R}^{n}) of yy on U\mathcal{U} with respect to VV:

?
(a)

X′VPX,W=X′VX'VP_{X,W} = X'V, or, equivalently, X′V(I−PX,W)=0X'V(I-P_{X,W}) = 0;

(b)

there exists a p×pp \times p matrix QQ such that VX=WXQVX = WXQ, or, equivalently, C(VX)⊂C(WX)C(VX) \subset C(WX).

Problem 14.51

Let XX represent an n×pn \times p matrix and WW an n×nn \times n symmetric nonnegative definite matrix. As in the special case where WW is positive definite, let

CW⊥(X)={y∈Rn×1:y⊥WC(X)}. C_W^{\perp }(X) = \{ y \in \mathbb {R}^{n \times 1} : y \perp _W C(X)\} .
?
(a)

By, for example, making use of the result of Exercise 14.46, show that

CW⊥(X)=N(X′W)=C(I−PX,W). C_W^{\perp }(X) = N(X'W) = C(I-P_{X,W}) .
(b)

Show that

dim⁡[CW⊥(X)]=n−rk⁡(WX)≥n−rk⁡(X)=n−dim⁡[C(X)]. \dim [C_W^{\perp }(X)] = n - \operatorname {rk}\left(WX\right) \geq n - \operatorname {rk}\left(X\right) = n - \dim [C(X)] .
(c)

By, for example, making use of the result of Exercise 14.46, show that, for any solution b∗b^{*} to the linear system X′WXb=X′WyX'WXb = X'Wy (in bb), the vector y−Xb∗y - Xb^{*} is a projection of yy on CW⊥(X)C_W^{\perp }(X) with respect to WW.