v2009.01.01 - Convex Optimization
v2009.01.01 - Convex Optimization v2009.01.01 - Convex Optimization
692 APPENDIX E. PROJECTION E.10.1 Distance and existence Existence of a fixed point is established: E.10.1.0.1 Theorem. Distance. [69] Given any two closed convex sets C 1 and C 2 in R n , then P 1 b∈ C 1 is a fixed point of the projection product P 1 P 2 if and only if P 1 b is a point of C 1 nearest C 2 . ⋄ Proof. (⇒) Given fixed point a =P 1 P 2 a∈ C 1 with b =P ∆ 2 a∈ C 2 in tandem so that a =P 1 b , then by the unique minimum-distance projection theorem (E.9.1.0.2) (1931) (b − a) T (u − a) ≤ 0 ∀u∈ C 1 (a − b) T (v − b) ≤ 0 ⇔ ∀v ∈ C 2 ‖a − b‖ ≤ ‖u − v‖ ∀u∈ C 1 and ∀v ∈ C 2 by Schwarz inequality ‖〈x,y〉‖ ≤ ‖x‖ ‖y‖ [197] [266]. (⇐) Suppose a∈ C 1 and ‖a − P 2 a‖ ≤ ‖u − P 2 u‖ ∀u∈ C 1 . Now suppose we choose u =P 1 P 2 a . Then ‖u − P 2 u‖ = ‖P 1 P 2 a − P 2 P 1 P 2 a‖ ≤ ‖a − P 2 a‖ ⇔ a = P 1 P 2 a (1932) Thus a = P 1 b (with b =P 2 a∈ C 2 ) is a fixed point in C 1 of the projection product P 1 P 2 . E.19 E.10.2 Feasibility and convergence The set of all fixed points of any nonexpansive mapping is a closed convex set. [131, lem.3.4] [27,1] The projection product P 1 P 2 is nonexpansive by Theorem E.9.3.0.1 because, for any vectors x,a∈ R n ‖P 1 P 2 x − P 1 P 2 a‖ ≤ ‖P 2 x − P 2 a‖ ≤ ‖x − a‖ (1933) If the intersection of two closed convex sets C 1 ∩ C 2 is empty, then the iterates converge to a point of minimum distance, a fixed point of the projection product. Otherwise, convergence is to some fixed point in their intersection E.19 Point b=P 2 a can be shown, similarly, to be a fixed point of the product P 2 P 1 .
E.10. ALTERNATING PROJECTION 693 (a feasible point) whose existence is guaranteed by virtue of the fact that each and every point in the convex intersection is in one-to-one correspondence with fixed points of the nonexpansive projection product. Bauschke & Borwein [27,2] argue that any sequence monotonic in the sense of Fejér is convergent. E.20 E.10.2.0.1 Definition. Fejér monotonicity. [233] Given closed convex set C ≠ ∅ , then a sequence x i ∈ R n , i=0, 1, 2..., is monotonic in the sense of Fejér with respect to C iff ‖x i+1 − c‖ ≤ ‖x i − c‖ for all i≥0 and each and every c ∈ C (1934) Given x 0 ∆ = b , if we express each iteration of alternating projection by x i+1 = P 1 P 2 x i , i=0, 1, 2... (1935) and define any fixed point a =P 1 P 2 a , then sequence x i is Fejér monotone with respect to fixed point a because ‖P 1 P 2 x i − a‖ ≤ ‖x i − a‖ ∀i ≥ 0 (1936) by nonexpansivity. The nonincreasing sequence ‖P 1 P 2 x i − a‖ is bounded below hence convergent because any bounded monotonic sequence in R is convergent; [222,1.2] [37,1.1] P 1 P 2 x i+1 = P 1 P 2 x i = x i+1 . Sequence x i therefore converges to some fixed point. If the intersection C 1 ∩ C 2 is nonempty, convergence is to some point there by the distance theorem. Otherwise, x i converges to a point in C 1 of minimum distance to C 2 . E.10.2.0.2 Example. Hyperplane/orthant intersection. Find a feasible point (1929) belonging to the nonempty intersection of two convex sets: given A∈ R m×n , β ∈ R(A) C 1 ∩ C 2 = R n + ∩ A = {y | y ≽ 0} ∩ {y | Ay = β} ⊂ R n (1937) the nonnegative orthant with affine subset A an intersection of hyperplanes. Projection of an iterate x i ∈ R n on A is calculated P 2 x i = x i − A T (AA T ) −1 (Ax i − β) (1828) E.20 Other authors prove convergence by different means; e.g., [149] [55]. △
- 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 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: E.10. ALTERNATING PROJECTION 691 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
- Page 715 and 716: 715 H − H + ∂H ∂H ∂H −
- Page 717 and 718: 717 O O sort-index matrix order of
- Page 719 and 720: (x,y) angle between vectors x and y
- Page 721 and 722: Bibliography [1] Suliman Al-Homidan
- Page 723 and 724: BIBLIOGRAPHY 723 [24] Alexander I.
- Page 725 and 726: BIBLIOGRAPHY 725 [52] Stephen Boyd,
- Page 727 and 728: BIBLIOGRAPHY 727 [78] Frank Critchl
- Page 729 and 730: BIBLIOGRAPHY 729 [105] Richard L. D
- Page 731 and 732: BIBLIOGRAPHY 731 [132] Michel X. Go
- Page 733 and 734: BIBLIOGRAPHY 733 [162] T. Herrmann,
- Page 735 and 736: BIBLIOGRAPHY 735 [191] Mark Kahrs a
- Page 737 and 738: BIBLIOGRAPHY 737 [220] K. V. Mardia
- Page 739 and 740: BIBLIOGRAPHY 739 [250] Pythagoras P
- Page 741 and 742: BIBLIOGRAPHY 741 [277] Anthony Man-
692 APPENDIX E. PROJECTION<br />
E.10.1<br />
Distance and existence<br />
Existence of a fixed point is established:<br />
E.10.1.0.1 Theorem. Distance. [69]<br />
Given any two closed convex sets C 1 and C 2 in R n , then P 1 b∈ C 1 is a fixed<br />
point of the projection product P 1 P 2 if and only if P 1 b is a point of C 1<br />
nearest C 2 .<br />
⋄<br />
Proof. (⇒) Given fixed point a =P 1 P 2 a∈ C 1 with b =P ∆ 2 a∈ C 2 in<br />
tandem so that a =P 1 b , then by the unique minimum-distance projection<br />
theorem (E.9.1.0.2)<br />
(1931)<br />
(b − a) T (u − a) ≤ 0 ∀u∈ C 1<br />
(a − b) T (v − b) ≤ 0<br />
⇔<br />
∀v ∈ C 2<br />
‖a − b‖ ≤ ‖u − v‖ ∀u∈ C 1 and ∀v ∈ C 2<br />
by Schwarz inequality ‖〈x,y〉‖ ≤ ‖x‖ ‖y‖ [197] [266].<br />
(⇐) Suppose a∈ C 1 and ‖a − P 2 a‖ ≤ ‖u − P 2 u‖ ∀u∈ C 1 . Now suppose we<br />
choose u =P 1 P 2 a . Then<br />
‖u − P 2 u‖ = ‖P 1 P 2 a − P 2 P 1 P 2 a‖ ≤ ‖a − P 2 a‖ ⇔ a = P 1 P 2 a (1932)<br />
Thus a = P 1 b (with b =P 2 a∈ C 2 ) is a fixed point in C 1 of the projection<br />
product P 1 P 2 . E.19<br />
<br />
E.10.2<br />
Feasibility and convergence<br />
The set of all fixed points of any nonexpansive mapping is a closed convex<br />
set. [131, lem.3.4] [27,1] The projection product P 1 P 2 is nonexpansive by<br />
Theorem E.9.3.0.1 because, for any vectors x,a∈ R n<br />
‖P 1 P 2 x − P 1 P 2 a‖ ≤ ‖P 2 x − P 2 a‖ ≤ ‖x − a‖ (1933)<br />
If the intersection of two closed convex sets C 1 ∩ C 2 is empty, then the iterates<br />
converge to a point of minimum distance, a fixed point of the projection<br />
product. Otherwise, convergence is to some fixed point in their intersection<br />
E.19 Point b=P 2 a can be shown, similarly, to be a fixed point of the product P 2 P 1 .