v2007.09.13 - Convex Optimization
v2007.09.13 - Convex Optimization v2007.09.13 - Convex Optimization
620 APPENDIX E. PROJECTIONE.9.2.2.2 Corollary. Unique projection via dual or normal cone.[73,4.7] (E.10.3.2, confer Theorem E.9.1.0.3) Given point x∈ R n andclosed convex cone K ⊆ R n , the following are equivalent statements:1. point Px is the unique minimum-distance projection of x on K2. Px ∈ K , x − Px ∈ −(K − Px) ∗ = −K ∗ ∩ (Px) ⊥3. Px ∈ K , 〈x − Px, Px〉 = 0, 〈x − Px, y〉 ≤ 0 ∀y ∈ K⋄E.9.2.2.3 Example. Unique projection on nonnegative orthant.(confer (1122)) From Theorem E.9.2.0.1, to project matrix H ∈ R m×n onthe self-dual orthant (2.13.5.1) of nonnegative matrices R m×n+ in isomorphicR mn , the necessary and sufficient conditions are:H ⋆ ≥ 0tr ( (H ⋆ − H) T H ⋆) = 0H ⋆ − H ≥ 0(1795)where the inequalities denote entrywise comparison. The optimal solutionH ⋆ is simply H having all its negative entries zeroed;H ⋆ ij = max{H ij , 0} , i,j∈{1... m} × {1... n} (1796)Now suppose the nonnegative orthant is translated by T ∈ R m×n ; id est,consider R m×n+ + T . Then projection on the translated orthant is [73,4.8]H ⋆ ij = max{H ij , T ij } (1797)E.9.2.2.4 Example. Unique projection on truncated convex cone.Consider the problem of projecting a point x on a closed convex cone thatis artificially bounded; really, a bounded convex polyhedron having a vertexat the origin:minimize ‖x − Ay‖ 2y∈R Nsubject to y ≽ 0(1798)‖y‖ ∞ ≤ 1
E.9. PROJECTION ON CONVEX SET 621where the convex cone has vertex-description (2.12.2.0.1), for A∈ R n×NK = {Ay | y ≽ 0} (1799)and where ‖y‖ ∞ ≤ 1 is the artificial bound. This is a convex optimizationproblem having no known closed-form solution, in general. It arises, forexample, in the fitting of hearing aids designed around a programmablegraphic equalizer (a filter bank whose only adjustable parameters are gainper band each bounded above by unity). [66] The problem is equivalent to aSchur-form semidefinite program (3.1.7.2)minimizey∈R N , t∈Rsubject tot[tI x − Ay(x − Ay) T t]≽ 0(1800)0 ≼ y ≼ 1E.9.3nonexpansivityE.9.3.0.1 Theorem. Nonexpansivity. [124,2] [73,5.3]When C ⊂ R n is an arbitrary closed convex set, projector P projecting on Cis nonexpansive in the sense: for any vectors x,y ∈ R n‖Px − Py‖ ≤ ‖x − y‖ (1801)with equality when x −Px = y −Py . E.16⋄Proof. [40]‖x − y‖ 2 = ‖Px − Py‖ 2 + ‖(I − P )x − (I − P )y‖ 2+ 2〈x − Px, Px − Py〉 + 2〈y − Py , Py − Px〉(1802)Nonnegativity of the last two terms follows directly from the uniqueminimum-distance projection theorem (E.9.1.0.2).E.16 This condition for equality corrects an error in [56] (where the norm is applied to eachside of the condition given here) easily revealed by counter-example.
- Page 570 and 571: 570 APPENDIX D. MATRIX CALCULUSD.2
- Page 572 and 573: 572 APPENDIX D. MATRIX CALCULUSalge
- Page 574 and 575: 574 APPENDIX D. MATRIX CALCULUStrac
- Page 576 and 577: 576 APPENDIX D. MATRIX CALCULUSD.2.
- Page 578 and 579: 578 APPENDIX D. MATRIX CALCULUS
- Page 580 and 581: 580 APPENDIX E. PROJECTIONThe follo
- Page 582 and 583: 582 APPENDIX E. PROJECTIONFor matri
- Page 584 and 585: 584 APPENDIX E. PROJECTION(⇐) To
- Page 586 and 587: 586 APPENDIX E. PROJECTIONNonorthog
- Page 588 and 589: 588 APPENDIX E. PROJECTIONE.2.0.0.1
- Page 590 and 591: 590 APPENDIX E. PROJECTIONE.3.2Orth
- Page 592 and 593: 592 APPENDIX E. PROJECTIONE.3.5Unif
- Page 594 and 595: 594 APPENDIX E. PROJECTIONE.4 Algeb
- Page 596 and 597: 596 APPENDIX E. PROJECTIONa ∗ 2K
- Page 598 and 599: 598 APPENDIX E. PROJECTIONwhere Y =
- Page 600 and 601: 600 APPENDIX E. PROJECTION(B.4.2).
- Page 602 and 603: 602 APPENDIX E. PROJECTIONis a nono
- Page 604 and 605: 604 APPENDIX E. PROJECTIONE.6.4.1Or
- Page 606 and 607: 606 APPENDIX E. PROJECTIONq i q T i
- Page 608 and 609: 608 APPENDIX E. PROJECTIONThe test
- Page 610 and 611: 610 APPENDIX E. PROJECTIONPerpendic
- Page 612 and 613: 612 APPENDIX E. PROJECTIONE.8 Range
- Page 614 and 615: 614 APPENDIX E. PROJECTIONAs for su
- Page 616 and 617: 616 APPENDIX E. PROJECTIONWith refe
- Page 618 and 619: 618 APPENDIX E. PROJECTIONProjectio
- Page 622 and 623: 622 APPENDIX E. PROJECTIONThe foreg
- Page 624 and 625: 624 APPENDIX E. PROJECTION❇❇❇
- Page 626 and 627: 626 APPENDIX E. PROJECTIONE.10 Alte
- Page 628 and 629: 628 APPENDIX E. PROJECTIONbH 1H 2P
- Page 630 and 631: 630 APPENDIX E. PROJECTIONa(a){y |
- Page 632 and 633: 632 APPENDIX E. PROJECTION(a feasib
- Page 634 and 635: 634 APPENDIX E. PROJECTIONwhile, th
- Page 636 and 637: 636 APPENDIX E. PROJECTIONE.10.2.1.
- Page 638 and 639: 638 APPENDIX E. PROJECTION10 0dist(
- Page 640 and 641: 640 APPENDIX E. PROJECTIONE.10.3.1D
- Page 642 and 643: 642 APPENDIX E. PROJECTIONE 3K ⊥
- Page 644 and 645: 644 APPENDIX E. PROJECTION
- Page 646 and 647: 646 APPENDIX F. MATLAB PROGRAMSif n
- Page 648 and 649: 648 APPENDIX F. MATLAB PROGRAMSend%
- Page 650 and 651: 650 APPENDIX F. MATLAB PROGRAMSF.1.
- Page 652 and 653: 652 APPENDIX F. MATLAB PROGRAMScoun
- Page 654 and 655: 654 APPENDIX F. MATLAB PROGRAMSF.3
- Page 656 and 657: 656 APPENDIX F. MATLAB PROGRAMSF.3.
- Page 658 and 659: 658 APPENDIX F. MATLAB PROGRAMS% so
- Page 660 and 661: 660 APPENDIX F. MATLAB PROGRAMS% tr
- Page 662 and 663: 662 APPENDIX F. MATLAB PROGRAMSF.4.
- Page 664 and 665: 664 APPENDIX F. MATLAB PROGRAMSbrea
- Page 666 and 667: 666 APPENDIX F. MATLAB PROGRAMSwhil
- Page 668 and 669: 668 APPENDIX F. MATLAB PROGRAMSF.7
E.9. PROJECTION ON CONVEX SET 621where the convex cone has vertex-description (2.12.2.0.1), for A∈ R n×NK = {Ay | y ≽ 0} (1799)and where ‖y‖ ∞ ≤ 1 is the artificial bound. This is a convex optimizationproblem having no known closed-form solution, in general. It arises, forexample, in the fitting of hearing aids designed around a programmablegraphic equalizer (a filter bank whose only adjustable parameters are gainper band each bounded above by unity). [66] The problem is equivalent to aSchur-form semidefinite program (3.1.7.2)minimizey∈R N , t∈Rsubject tot[tI x − Ay(x − Ay) T t]≽ 0(1800)0 ≼ y ≼ 1E.9.3nonexpansivityE.9.3.0.1 Theorem. Nonexpansivity. [124,2] [73,5.3]When C ⊂ R n is an arbitrary closed convex set, projector P projecting on Cis nonexpansive in the sense: for any vectors x,y ∈ R n‖Px − Py‖ ≤ ‖x − y‖ (1801)with equality when x −Px = y −Py . E.16⋄Proof. [40]‖x − y‖ 2 = ‖Px − Py‖ 2 + ‖(I − P )x − (I − P )y‖ 2+ 2〈x − Px, Px − Py〉 + 2〈y − Py , Py − Px〉(1802)Nonnegativity of the last two terms follows directly from the uniqueminimum-distance projection theorem (E.9.1.0.2).E.16 This condition for equality corrects an error in [56] (where the norm is applied to eachside of the condition given here) easily revealed by counter-example.