v2009.01.01 - Convex Optimization
v2009.01.01 - Convex Optimization v2009.01.01 - Convex Optimization
664 APPENDIX E. PROJECTION E.6.3 Orthogonal projection on a vector The formula for orthogonal projection of vector x on the range of vector y (one-dimensional projection) is basic analytic geometry; [12,3.3] [287,3.2] [318,2.2] [332,1-8] 〈y,x〉 〈y,y〉 y = yT x y T y y = yyT y T y x = ∆ P 1 x (1848) where 〈y,x〉/〈y,y〉 is the coefficient of projection on R(y) . An equivalent description is: Vector P 1 x is the orthogonal projection of vector x on R(P 1 )= R(y). Rank-one matrix P 1 is a projection matrix because P 2 1 =P 1 . The direction of projection is orthogonal because P T 1 = P 1 . E.6.4 P 1 x − x ⊥ R(P 1 ) (1849) Orthogonal projection on a vectorized matrix From (1848), given instead X, Y ∈ R m×n , we have the one-dimensional orthogonal projection of matrix X in isomorphic R mn on the range of vectorized Y : (2.2) 〈Y , X〉 〈Y , Y 〉 Y (1850) where 〈Y , X〉/〈Y , Y 〉 is the coefficient of projection. For orthogonal projection, the term outside the vector inner-products 〈 〉 must be identical to the terms inside in three places. E.6.4.1 Orthogonal projection on dyad There is opportunity for insight when Y is a dyad yz T (B.1): Instead given X ∈ R m×n , y ∈ R m , and z ∈ R n 〈yz T , X〉 〈yz T , yz T 〉 yzT = yT Xz y T y z T z yzT (1851) is the one-dimensional orthogonal projection of X in isomorphic R mn on the range of vectorized yz T . To reveal the obscured symmetric projection
E.6. VECTORIZATION INTERPRETATION, 665 matrices P 1 and P 2 we rewrite (1851): y T Xz y T y z T z yzT = yyT y T y X zzT z T z ∆ = P 1 XP 2 (1852) So for projector dyads, projection (1852) is the orthogonal projection in R mn if and only if projectors P 1 and P 2 are symmetric; E.10 in other words, for orthogonal projection on the range of a vectorized dyad yz T , the term outside the vector inner-products 〈 〉 in (1851) must be identical to the terms inside in three places. When P 1 and P 2 are rank-one symmetric projectors as in (1852), (30) R(vec P 1 XP 2 ) = R(vecyz T ) in R mn (1853) and P 1 XP 2 − X ⊥ yz T in R mn (1854) When y=z then P 1 =P 2 =P T 2 and P 1 XP 1 = 〈P 1 , X〉 P 1 = 〈P 1 , X〉 〈P 1 , P 1 〉 P 1 (1855) meaning, P 1 XP 1 is equivalent to orthogonal projection of matrix X on the range of vectorized projector dyad P 1 . Yet this relationship between matrix product and vector inner-product does not hold for general symmetric projector matrices. E.10 For diagonalizable X ∈ R m×m (A.5), its orthogonal projection in isomorphic R m2 on the range of vectorized yz T ∈ R m×m becomes: P 1 XP 2 = m∑ λ i P 1 s i wi T P 2 i=1 When R(P 1 ) = R(w j ) and R(P 2 ) = R(s j ), the j th dyad term from the diagonalization is isolated but only, in general, to within a scale factor because neither set of left or right eigenvectors is necessarily orthonormal unless X is normal [344,3.2]. Yet when R(P 2 )= R(s k ) , k≠j ∈{1... m}, then P 1 XP 2 =0.
- 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 and 652: 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: 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.
- Page 703 and 704: E.10. ALTERNATING PROJECTION 703 E
- Page 705 and 706: Appendix F Notation and a few defin
- Page 707 and 708: 707 a.i. c.i. l.i. w.r.t affinely i
- Page 709 and 710: 709 is or ← → t → 0 + as in
- Page 711 and 712: 711 ∑ π(γ) Ξ Π ∏ ψ(Z) D D
- Page 713 and 714: 713 R m×n Euclidean vector space o
664 APPENDIX E. PROJECTION<br />
E.6.3<br />
Orthogonal projection on a vector<br />
The formula for orthogonal projection of vector x on the range of vector y<br />
(one-dimensional projection) is basic analytic geometry; [12,3.3] [287,3.2]<br />
[318,2.2] [332,1-8]<br />
〈y,x〉<br />
〈y,y〉 y = yT x<br />
y T y y = yyT<br />
y T y x = ∆ P 1 x (1848)<br />
where 〈y,x〉/〈y,y〉 is the coefficient of projection on R(y) . An equivalent<br />
description is: Vector P 1 x is the orthogonal projection of vector x on<br />
R(P 1 )= R(y). Rank-one matrix P 1 is a projection matrix because P 2 1 =P 1 .<br />
The direction of projection is orthogonal<br />
because P T 1 = P 1 .<br />
E.6.4<br />
P 1 x − x ⊥ R(P 1 ) (1849)<br />
Orthogonal projection on a vectorized matrix<br />
From (1848), given instead X, Y ∈ R m×n , we have the one-dimensional<br />
orthogonal projection of matrix X in isomorphic R mn on the range of<br />
vectorized Y : (2.2)<br />
〈Y , X〉<br />
〈Y , Y 〉 Y (1850)<br />
where 〈Y , X〉/〈Y , Y 〉 is the coefficient of projection.<br />
For orthogonal projection, the term outside the vector inner-products 〈 〉<br />
must be identical to the terms inside in three places.<br />
E.6.4.1<br />
Orthogonal projection on dyad<br />
There is opportunity for insight when Y is a dyad yz T (B.1): Instead given<br />
X ∈ R m×n , y ∈ R m , and z ∈ R n<br />
〈yz T , X〉<br />
〈yz T , yz T 〉 yzT = yT Xz<br />
y T y z T z yzT (1851)<br />
is the one-dimensional orthogonal projection of X in isomorphic R mn on<br />
the range of vectorized yz T . To reveal the obscured symmetric projection