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

650 APPENDIX E. PROJECTION We get P =AA † so this projection matrix must be symmetric. Then for any matrix A∈ R m×n , symmetric idempotent P projects a given vector x in R m orthogonally on R(A). Under either condition (1783) or (1784), the projection Px is unique minimum-distance; for subspaces, perpendicularity and minimum-distance conditions are equivalent. E.3.1 Four subspaces We summarize the orthogonal projectors projecting on the four fundamental subspaces: for A∈ R m×n A † A : R n on R(A † A) = R(A T ) AA † : R m on R(AA † ) = R(A) I −A † A : R n on R(I −A † A) = N(A) I −AA † : R m on R(I −AA † ) = N(A T ) (1786) For completeness: E.7 (1778) N(A † A) = N(A) N(AA † ) = N(A T ) N(I −A † A) = R(A T ) N(I −AA † ) = R(A) (1787) E.3.2 Orthogonal characterization Any symmetric projector P 2 =P ∈ S m projecting on nontrivial R(Q) can be defined by the orthonormality condition Q T Q = I . When skinny matrix Q∈ R m×k has orthonormal columns, then Q † = Q T by the Moore-Penrose conditions. Hence, any P having an orthonormal decomposition (E.3.4) where [287,3.3] (1493) P = QQ T , Q T Q = I (1788) R(P ) = R(Q) , N(P ) = N(Q T ) (1789) E.7 Proof is by singular value decomposition (A.6.2): N(A † A)⊆ N(A) is obvious. Conversely, suppose A † Ax=0. Then x T A † Ax=x T QQ T x=‖Q T x‖ 2 =0 where A=UΣQ T is the subcompact singular value decomposition. Because R(Q)= R(A T ) , then x∈N(A) which implies N(A † A)⊇ N(A). ∴ N(A † A)= N(A).

E.3. SYMMETRIC IDEMPOTENT MATRICES 651 is an orthogonal projector projecting on R(Q) having, for Px∈ R(Q) (confer (1774)) Px − x ⊥ R(Q) in R m (1790) From (1788), orthogonal projector P is obviously positive semidefinite (A.3.1.0.6); necessarily, P T = P , P † = P , ‖P ‖ 2 = 1, P ≽ 0 (1791) and ‖Px‖ = ‖QQ T x‖ = ‖Q T x‖ because ‖Qy‖ = ‖y‖ ∀y ∈ R k . All orthogonal projectors are therefore nonexpansive because √ 〈Px, x〉 = ‖Px‖ = ‖Q T x‖ ≤ ‖x‖ ∀x∈ R m (1792) the Bessel inequality, [92] [197] with equality when x∈ R(Q). From the diagonalization of idempotent matrices (1768) on page 644 P = SΦS T = m∑ φ i s i s T i = i=1 k∑ ≤ m i=1 s i s T i (1793) orthogonal projection of point x on R(P ) has expression like an orthogonal expansion [92,4.10] where Px = QQ T x = k∑ s T i xs i (1794) i=1 Q = S(:,1:k) = [ s 1 · · · s k ] ∈ R m×k (1795) and where the s i [sic] are orthonormal eigenvectors of symmetric idempotent P . When the domain is restricted to range of P , say x=Qξ for ξ ∈ R k , then x = Px = QQ T Qξ = Qξ and expansion is unique. Otherwise, any component of x in N(Q T ) will be annihilated. E.3.2.0.1 Theorem. Symmetric rank/trace. (confer (1781) (1359)) P T = P , P 2 = P ⇔ rankP = trP = ‖P ‖ 2 F and rank(I − P ) = tr(I − P ) = ‖I − P ‖ 2 F (1796) ⋄

E.3. SYMMETRIC IDEMPOTENT MATRICES 651<br />

is an orthogonal projector projecting on R(Q) having, for Px∈ R(Q)<br />

(confer (1774))<br />

Px − x ⊥ R(Q) in R m (1790)<br />

From (1788), orthogonal projector P is obviously positive semidefinite<br />

(A.3.1.0.6); necessarily,<br />

P T = P , P † = P , ‖P ‖ 2 = 1, P ≽ 0 (1791)<br />

and ‖Px‖ = ‖QQ T x‖ = ‖Q T x‖ because ‖Qy‖ = ‖y‖ ∀y ∈ R k . All orthogonal<br />

projectors are therefore nonexpansive because<br />

√<br />

〈Px, x〉 = ‖Px‖ = ‖Q T x‖ ≤ ‖x‖ ∀x∈ R m (1792)<br />

the Bessel inequality, [92] [197] with equality when x∈ R(Q).<br />

From the diagonalization of idempotent matrices (1768) on page 644<br />

P = SΦS T =<br />

m∑<br />

φ i s i s T i =<br />

i=1<br />

k∑<br />

≤ m<br />

i=1<br />

s i s T i (1793)<br />

orthogonal projection of point x on R(P ) has expression like an orthogonal<br />

expansion [92,4.10]<br />

where<br />

Px = QQ T x =<br />

k∑<br />

s T i xs i (1794)<br />

i=1<br />

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

]<br />

∈ R<br />

m×k<br />

(1795)<br />

and where the s i [sic] are orthonormal eigenvectors of symmetric<br />

idempotent P . When the domain is restricted to range of P , say x=Qξ for<br />

ξ ∈ R k , then x = Px = QQ T Qξ = Qξ and expansion is unique. Otherwise,<br />

any component of x in N(Q T ) will be annihilated.<br />

E.3.2.0.1 Theorem. Symmetric rank/trace. (confer (1781) (1359))<br />

P T = P , P 2 = P<br />

⇔<br />

rankP = trP = ‖P ‖ 2 F and rank(I − P ) = tr(I − P ) = ‖I − P ‖ 2 F<br />

(1796)<br />

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

Saved successfully!

Ooh no, something went wrong!