v2009.01.01 - Convex Optimization

v2009.01.01 - Convex Optimization v2009.01.01 - Convex Optimization

convexoptimization.com
from convexoptimization.com More from this publisher
10.03.2015 Views

644 APPENDIX E. PROJECTION When matrix A in (1762) is skinny full-rank (A † A = I) or has orthonormal columns (A T A = I), either property leads to a biorthogonal characterization of nonorthogonal projection: E.1.1 Biorthogonal characterization Any nonorthogonal projector P 2 =P ∈ R m×m projecting on nontrivial R(U) can be defined by a biorthogonality condition Q T U =I ; the biorthogonal decomposition of P being (confer (1762)) where E.3 (B.1.1.1) and where generally (confer (1791)) E.4 P = UQ T , Q T U = I (1765) R(P )= R(U) , N(P )= N(Q T ) (1766) P T ≠ P , P † ≠ P , ‖P ‖ 2 ≠ 1, P 0 (1767) and P is not nonexpansive (1792). (⇐) To verify assertion (1765) we observe: because idempotent matrices are diagonalizable (A.5), [176,3.3, prob.3] they must have the form (1438) P = SΦS −1 = m∑ φ i s i wi T = i=1 k∑ ≤ m i=1 s i w T i (1768) that is a sum of k = rankP independent projector dyads (idempotent dyads,B.1.1,E.6.2.1) where φ i ∈ {0, 1} are the eigenvalues of P [344,4.1, thm.4.1] in diagonal matrix Φ∈ R m×m arranged in nonincreasing E.3 Proof. Obviously, R(P ) ⊆ R(U). Because Q T U = I R(P ) = {UQ T x | x ∈ R m } ⊇ {UQ T Uy | y ∈ R k } = R(U) E.4 Orthonormal decomposition (1788) (conferE.3.4) is a special case of biorthogonal decomposition (1765) characterized by (1791). So, these characteristics (1767) are not necessary conditions for biorthogonality.

E.1. IDEMPOTENT MATRICES 645 order, and where s i ,w i ∈ R m are the right- and left-eigenvectors of P , respectively, which are independent and real. E.5 Therefore U ∆ = S(:,1:k) = [ s 1 · · · s k ] ∈ R m×k (1769) is the full-rank matrix S ∈ R m×m having m − k columns truncated (corresponding to 0 eigenvalues), while ⎡ Q T = ∆ S −1 (1:k, :) = ⎣ w T1 . w T k ⎤ ⎦ ∈ R k×m (1770) is matrix S −1 having the corresponding m − k rows truncated. By the 0 eigenvalues theorem (A.7.3.0.1), R(U)= R(P ) , R(Q)= R(P T ) , and R(P ) = span {s i | φ i = 1 ∀i} N(P ) = span {s i | φ i = 0 ∀i} R(P T ) = span {w i | φ i = 1 ∀i} N(P T ) = span {w i | φ i = 0 ∀i} (1771) Thus biorthogonality Q T U =I is a necessary condition for idempotence, and so the collection of nonorthogonal projectors projecting on R(U) is the affine subset P k =UQ T k where Q k = {Q | Q T U = I , Q∈ R m×k }. (⇒) Biorthogonality is a sufficient condition for idempotence; P 2 = k∑ s i wi T i=1 k∑ s j wj T = P (1772) j=1 id est, if the cross-products are annihilated, then P 2 =P . Nonorthogonal projection of biorthogonal expansion, x on R(P ) has expression like a Px = UQ T x = k∑ wi T xs i (1773) i=1 E.5 Eigenvectors of a real matrix corresponding to real eigenvalues must be real. (A.5.0.0.1)

644 APPENDIX E. PROJECTION<br />

When matrix A in (1762) is skinny full-rank (A † A = I) or has orthonormal<br />

columns (A T A = I), either property leads to a biorthogonal characterization<br />

of nonorthogonal projection:<br />

E.1.1<br />

Biorthogonal characterization<br />

Any nonorthogonal projector P 2 =P ∈ R m×m projecting on nontrivial R(U)<br />

can be defined by a biorthogonality condition Q T U =I ; the biorthogonal<br />

decomposition of P being (confer (1762))<br />

where E.3 (B.1.1.1)<br />

and where generally (confer (1791)) E.4<br />

P = UQ T , Q T U = I (1765)<br />

R(P )= R(U) , N(P )= N(Q T ) (1766)<br />

P T ≠ P , P † ≠ P , ‖P ‖ 2 ≠ 1, P 0 (1767)<br />

and P is not nonexpansive (1792).<br />

(⇐) To verify assertion (1765) we observe: because idempotent matrices<br />

are diagonalizable (A.5), [176,3.3, prob.3] they must have the form (1438)<br />

P = SΦS −1 =<br />

m∑<br />

φ i s i wi T =<br />

i=1<br />

k∑<br />

≤ m<br />

i=1<br />

s i w T i (1768)<br />

that is a sum of k = rankP independent projector dyads (idempotent<br />

dyads,B.1.1,E.6.2.1) where φ i ∈ {0, 1} are the eigenvalues of P<br />

[344,4.1, thm.4.1] in diagonal matrix Φ∈ R m×m arranged in nonincreasing<br />

E.3 Proof. Obviously, R(P ) ⊆ R(U). Because Q T U = I<br />

R(P ) = {UQ T x | x ∈ R m }<br />

⊇ {UQ T Uy | y ∈ R k } = R(U)<br />

<br />

E.4 Orthonormal decomposition (1788) (conferE.3.4) is a special case of biorthogonal<br />

decomposition (1765) characterized by (1791). So, these characteristics (1767) are not<br />

necessary conditions for biorthogonality.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!