v2009.01.01 - Convex Optimization
v2009.01.01 - Convex Optimization v2009.01.01 - Convex Optimization
652 APPENDIX E. PROJECTION Proof. We take as given Theorem E.2.0.0.1 establishing idempotence. We have left only to show trP =‖P ‖ 2 F ⇒ P T = P , established in [344,7.1]. E.3.3 Summary, symmetric idempotent In summary, orthogonal projector P is a linear operator defined [173,A.3.1] by idempotence and symmetry, and characterized by positive semidefiniteness and nonexpansivity. The algebraic complement (E.2) to R(P ) becomes the orthogonal complement R(I − P ) ; id est, R(P ) ⊥ R(I − P ). E.3.4 Orthonormal decomposition When Z =0 in the general nonorthogonal projector A(A † + BZ T ) (1762), an orthogonal projector results (for any matrix A) characterized principally by idempotence and symmetry. Any real orthogonal projector may, in fact, be represented by an orthonormal decomposition such as (1788). [178,1, prob.42] To verify that assertion for the four fundamental subspaces (1786), we need only to express A by subcompact singular value decomposition (A.6.2): From pseudoinverse (1463) of A = UΣQ T ∈ R m×n AA † = UΣΣ † U T = UU T , I − AA † = I − UU T = U ⊥ U ⊥T , A † A = QΣ † ΣQ T = QQ T I − A † A = I − QQ T = Q ⊥ Q ⊥T (1797) where U ⊥ ∈ R m×m−rank A holds columnar an orthonormal basis for the orthogonal complement of R(U) , and likewise for Q ⊥ ∈ R n×n−rank A . Existence of an orthonormal decomposition is sufficient to establish idempotence and symmetry of an orthogonal projector (1788). E.3.5 Unifying trait of all projectors: direction Relation (1797) shows: orthogonal projectors simultaneously possess a biorthogonal decomposition (conferE.1.1) (for example, AA † for skinny-or-square A full-rank) and an orthonormal decomposition (UU T whence Px = UU T x).
E.3. SYMMETRIC IDEMPOTENT MATRICES 653 E.3.5.1 orthogonal projector, orthonormal decomposition Consider orthogonal expansion of x∈ R(U) : n∑ x = UU T x = u i u T i x (1798) a sum of one-dimensional orthogonal projections (E.6.3), where i=1 U ∆ = [u 1 · · · u n ] and U T U = I (1799) and where the subspace projector has two expressions, (1797) AA † ∆ = UU T (1800) where A ∈ R m×n has rank n . The direction of projection of x on u j for some j ∈{1... n} , for example, is orthogonal to u j but parallel to a vector in the span of all the remaining vectors constituting the columns of U ; u T j(u j u T j x − x) = 0 u j u T j x − x = u j u T j x − UU T x ∈ R({u i |i=1... n, i≠j}) (1801) E.3.5.2 orthogonal projector, biorthogonal decomposition We get a similar result for the biorthogonal expansion of x∈ R(A). Define A ∆ = [a 1 a 2 · · · a n ] ∈ R m×n (1802) and the rows of the pseudoinverse E.8 ⎡ ⎤ a ∗T 1 A † = ∆ a ∗T ⎢ ⎣ 2 ⎥ . ⎦ ∈ Rn×m (1803) a ∗T n under biorthogonality condition A † A=I . In biorthogonal expansion (2.13.8) n∑ x = AA † x = a i a ∗T i x (1804) E.8 Notation * in this context connotes extreme direction of a dual cone; e.g., (361) or Example E.5.0.0.3. i=1
- 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 and 650: E.3. SYMMETRIC IDEMPOTENT MATRICES
- Page 651: 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
- Page 701 and 702: E.10. ALTERNATING PROJECTION 701 E.
E.3. SYMMETRIC IDEMPOTENT MATRICES 653<br />
E.3.5.1<br />
orthogonal projector, orthonormal decomposition<br />
Consider orthogonal expansion of x∈ R(U) :<br />
n∑<br />
x = UU T x = u i u T i x (1798)<br />
a sum of one-dimensional orthogonal projections (E.6.3), where<br />
i=1<br />
U ∆ = [u 1 · · · u n ] and U T U = I (1799)<br />
and where the subspace projector has two expressions, (1797)<br />
AA † ∆ = UU T (1800)<br />
where A ∈ R m×n has rank n . The direction of projection of x on u j for<br />
some j ∈{1... n} , for example, is orthogonal to u j but parallel to a vector<br />
in the span of all the remaining vectors constituting the columns of U ;<br />
u T j(u j u T j x − x) = 0<br />
u j u T j x − x = u j u T j x − UU T x ∈ R({u i |i=1... n, i≠j})<br />
(1801)<br />
E.3.5.2<br />
orthogonal projector, biorthogonal decomposition<br />
We get a similar result for the biorthogonal expansion of x∈ R(A). Define<br />
A ∆ = [a 1 a 2 · · · a n ] ∈ R m×n (1802)<br />
and the rows of the pseudoinverse E.8<br />
⎡ ⎤<br />
a ∗T<br />
1<br />
A † =<br />
∆ a ∗T<br />
⎢<br />
⎣<br />
2 ⎥<br />
. ⎦ ∈ Rn×m (1803)<br />
a ∗T<br />
n<br />
under biorthogonality condition A † A=I . In biorthogonal expansion<br />
(2.13.8)<br />
n∑<br />
x = AA † x = a i a ∗T<br />
i x (1804)<br />
E.8 Notation * in this context connotes extreme direction of a dual cone; e.g., (361) or<br />
Example E.5.0.0.3.<br />
i=1