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)

E.1. IDEMPOTENT MATRICES 645<br />

order, and where s i ,w i ∈ R m are the right- and left-eigenvectors of P ,<br />

respectively, which are independent and real. E.5 Therefore<br />

U ∆ = S(:,1:k) = [ s 1 · · · s k<br />

]<br />

∈ R<br />

m×k<br />

(1769)<br />

is the full-rank matrix S ∈ R m×m having m − k columns truncated<br />

(corresponding to 0 eigenvalues), while<br />

⎡<br />

Q T = ∆ S −1 (1:k, :) = ⎣<br />

w T1<br />

.<br />

w T k<br />

⎤<br />

⎦ ∈ R k×m (1770)<br />

is matrix S −1 having the corresponding m − k rows truncated. By the<br />

0 eigenvalues theorem (A.7.3.0.1), R(U)= R(P ) , R(Q)= R(P T ) , and<br />

R(P ) = span {s i | φ i = 1 ∀i}<br />

N(P ) = span {s i | φ i = 0 ∀i}<br />

R(P T ) = span {w i | φ i = 1 ∀i}<br />

N(P T ) = span {w i | φ i = 0 ∀i}<br />

(1771)<br />

Thus biorthogonality Q T U =I is a necessary condition for idempotence, and<br />

so the collection of nonorthogonal projectors projecting on R(U) is the affine<br />

subset P k =UQ T k where Q k = {Q | Q T U = I , Q∈ R m×k }.<br />

(⇒) Biorthogonality is a sufficient condition for idempotence;<br />

P 2 =<br />

k∑<br />

s i wi<br />

T<br />

i=1<br />

k∑<br />

s j wj T = P (1772)<br />

j=1<br />

id est, if the cross-products are annihilated, then P 2 =P .<br />

<br />

Nonorthogonal projection of<br />

biorthogonal expansion,<br />

x on R(P ) has expression like a<br />

Px = UQ T x =<br />

k∑<br />

wi T xs i (1773)<br />

i=1<br />

E.5 Eigenvectors of a real matrix corresponding to real eigenvalues must be real.<br />

(A.5.0.0.1)

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

Saved successfully!

Ooh no, something went wrong!