v2009.01.01 - Convex Optimization
v2009.01.01 - Convex Optimization v2009.01.01 - Convex Optimization
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) ⋄
- Page 599 and 600: C.2. TRACE, SINGULAR AND EIGEN VALU
- Page 601 and 602: C.2. TRACE, SINGULAR AND EIGEN VALU
- Page 603 and 604: C.3. ORTHOGONAL PROCRUSTES PROBLEM
- Page 605 and 606: C.4. TWO-SIDED ORTHOGONAL PROCRUSTE
- Page 607 and 608: C.4. TWO-SIDED ORTHOGONAL PROCRUSTE
- Page 609 and 610: Appendix D Matrix calculus From too
- Page 611 and 612: D.1. DIRECTIONAL DERIVATIVE, TAYLOR
- Page 613 and 614: D.1. DIRECTIONAL DERIVATIVE, TAYLOR
- Page 615 and 616: D.1. DIRECTIONAL DERIVATIVE, TAYLOR
- Page 617 and 618: D.1. DIRECTIONAL DERIVATIVE, TAYLOR
- Page 619 and 620: D.1. DIRECTIONAL DERIVATIVE, TAYLOR
- Page 621 and 622: D.1. DIRECTIONAL DERIVATIVE, TAYLOR
- Page 623 and 624: D.1. DIRECTIONAL DERIVATIVE, TAYLOR
- Page 625 and 626: D.1. DIRECTIONAL DERIVATIVE, TAYLOR
- Page 627 and 628: D.1. DIRECTIONAL DERIVATIVE, TAYLOR
- Page 629 and 630: D.1. DIRECTIONAL DERIVATIVE, TAYLOR
- Page 631 and 632: D.2. TABLES OF GRADIENTS AND DERIVA
- Page 633 and 634: D.2. TABLES OF GRADIENTS AND DERIVA
- Page 635 and 636: D.2. TABLES OF GRADIENTS AND DERIVA
- Page 637 and 638: D.2. TABLES OF GRADIENTS AND DERIVA
- Page 639 and 640: Appendix E Projection For any A∈
- Page 641 and 642: 641 U T = U † for orthonormal (in
- Page 643 and 644: E.1. IDEMPOTENT MATRICES 643 where
- Page 645 and 646: E.1. IDEMPOTENT MATRICES 645 order,
- Page 647 and 648: E.1. IDEMPOTENT MATRICES 647 When t
- Page 649: E.3. SYMMETRIC IDEMPOTENT MATRICES
- Page 653 and 654: E.3. SYMMETRIC IDEMPOTENT MATRICES
- Page 655 and 656: E.5. PROJECTION EXAMPLES 655 E.4.0.
- Page 657 and 658: E.5. PROJECTION EXAMPLES 657 a ∗
- Page 659 and 660: E.5. PROJECTION EXAMPLES 659 E.5.0.
- Page 661 and 662: E.6. VECTORIZATION INTERPRETATION,
- Page 663 and 664: E.6. VECTORIZATION INTERPRETATION,
- Page 665 and 666: E.6. VECTORIZATION INTERPRETATION,
- Page 667 and 668: E.6. VECTORIZATION INTERPRETATION,
- Page 669 and 670: E.7. ON VECTORIZED MATRICES OF HIGH
- Page 671 and 672: E.7. ON VECTORIZED MATRICES OF HIGH
- Page 673 and 674: E.8. RANGE/ROWSPACE INTERPRETATION
- Page 675 and 676: E.9. PROJECTION ON CONVEX SET 675 A
- Page 677 and 678: E.9. PROJECTION ON CONVEX SET 677 W
- Page 679 and 680: E.9. PROJECTION ON CONVEX SET 679 P
- Page 681 and 682: E.9. PROJECTION ON CONVEX SET 681 E
- Page 683 and 684: E.9. PROJECTION ON CONVEX SET 683 T
- Page 685 and 686: E.9. PROJECTION ON CONVEX SET 685
- Page 687 and 688: E.10. ALTERNATING PROJECTION 687 E.
- Page 689 and 690: E.10. ALTERNATING PROJECTION 689 b
- Page 691 and 692: E.10. ALTERNATING PROJECTION 691 a
- Page 693 and 694: E.10. ALTERNATING PROJECTION 693 (a
- Page 695 and 696: E.10. ALTERNATING PROJECTION 695 wh
- Page 697 and 698: E.10. ALTERNATING PROJECTION 697 E.
- Page 699 and 700: E.10. ALTERNATING PROJECTION 699 10
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 />
⋄