v2009.01.01 - Convex Optimization
v2009.01.01 - Convex Optimization v2009.01.01 - Convex Optimization
520 CHAPTER 7. PROXIMITY PROBLEMS 7.2.2.2 Applying trace rank-heuristic to Problem 2 Substituting rank envelope for rank function in Problem 2, for D ∈ EDM N (confer (947)) cenv rank(−V T NDV N ) = cenv rank(−V DV ) ∝ − tr(V DV ) (1271) and for desired affine dimension ρ ≤ N −1 and nonnegative H [sic] we get a convex optimization problem minimize ‖ ◦√ D − H‖ 2 F D subject to − tr(V DV ) ≤ κρ (1272) D ∈ EDM N where κ ∈ R + is a constant determined by cut-and-try. The equivalent semidefinite program makes κ variable: for nonnegative and symmetric H minimize D , Y , κ subject to κρ + 2 tr(V Y V ) [ ] dij y ij ≽ 0 , y ij h 2 ij N ≥ j > i = 1... N −1 (1273) − tr(V DV ) ≤ κρ Y ∈ S N h D ∈ EDM N which is the same as (1264), the problem with no explicit constraint on affine dimension. As the present problem is stated, the desired affine dimension ρ yields to the variable scale factor κ ; ρ is effectively ignored. Yet this result is an illuminant for problem (1264) and it equivalents (all the way back to (1257)): When the given measurement matrix H is nonnegative and symmetric, finding the closest EDM D as in problem (1257), (1260), or (1264) implicitly entails minimization of affine dimension (confer5.8.4,5.14.4). Those non−rank-constrained problems are each inherently equivalent to cenv(rank)-minimization problem (1273), in other words, and their optimal solutions are unique because of the strictly convex objective function in (1257). 7.2.2.3 Rank-heuristic insight Minimization of affine dimension by use of this trace rank-heuristic (1271) tends to find the list configuration of least energy; rather, it tends to optimize
7.2. SECOND PREVALENT PROBLEM: 521 compaction of the reconstruction by minimizing total distance. (823) It is best used where some physical equilibrium implies such an energy minimization; e.g., [305,5]. For this Problem 2, the trace rank-heuristic arose naturally in the objective in terms of V . We observe: V (in contrast to VN T ) spreads energy over all available distances (B.4.2 no.20, confer no.22) although the rank function itself is insensitive to choice of auxiliary matrix. 7.2.2.4 Rank minimization heuristic beyond convex envelope Fazel, Hindi, & Boyd [114] [338] [115] propose a rank heuristic more potent than trace (1269) for problems of rank minimization; rankY ← log det(Y +εI) (1274) the concave surrogate function log det in place of quasiconcave rankY (2.9.2.6.2) when Y ∈ S n + is variable and where ε is a small positive constant. They propose minimization of the surrogate by substituting a sequence comprising infima of a linearized surrogate about the current estimate Y i ; id est, from the first-order Taylor series expansion about Y i on some open interval of ‖Y ‖ 2 (D.1.7) log det(Y + εI) ≈ log det(Y i + εI) + tr ( (Y i + εI) −1 (Y − Y i ) ) (1275) we make the surrogate sequence of infima over bounded convex feasible set C where, for i = 0... arg inf rankY ← lim Y i+1 (1276) Y ∈ C i→∞ Y i+1 = arg inf Y ∈ C tr( (Y i + εI) −1 Y ) (1277) Choosing Y 0 =I , the first step becomes equivalent to finding the infimum of trY ; the trace rank-heuristic (1269). The intuition underlying (1277) is the new term in the argument of trace; specifically, (Y i + εI) −1 weights Y so that relatively small eigenvalues of Y found by the infimum are made even smaller.
- Page 469 and 470: 6.6. CORRESPONDENCE TO PSD CONE S N
- Page 471 and 472: 6.6. CORRESPONDENCE TO PSD CONE S N
- Page 473 and 474: 6.7. VECTORIZATION & PROJECTION INT
- Page 475 and 476: 6.7. VECTORIZATION & PROJECTION INT
- Page 477 and 478: 6.8. DUAL EDM CONE 477 When the Fin
- Page 479 and 480: 6.8. DUAL EDM CONE 479 Proof. First
- Page 481 and 482: 6.8. DUAL EDM CONE 481 EDM 2 = S 2
- Page 483 and 484: 6.8. DUAL EDM CONE 483 whose veraci
- Page 485 and 486: 6.8. DUAL EDM CONE 485 6.8.1.3.1 Ex
- Page 487 and 488: 6.8. DUAL EDM CONE 487 has dual aff
- Page 489 and 490: 6.8. DUAL EDM CONE 489 6.8.1.7 Scho
- Page 491 and 492: 6.9. THEOREM OF THE ALTERNATIVE 491
- Page 493 and 494: 6.10. POSTSCRIPT 493 When D is an E
- Page 495 and 496: Chapter 7 Proximity problems In sum
- Page 497 and 498: 497 project on the subspace, then p
- Page 499 and 500: 499 H S N h 0 EDM N K = S N h ∩ R
- Page 501 and 502: 501 7.0.3 Problem approach Problems
- Page 503 and 504: 7.1. FIRST PREVALENT PROBLEM: 503 f
- Page 505 and 506: 7.1. FIRST PREVALENT PROBLEM: 505 7
- Page 507 and 508: 7.1. FIRST PREVALENT PROBLEM: 507 d
- Page 509 and 510: 7.1. FIRST PREVALENT PROBLEM: 509 7
- Page 511 and 512: 7.1. FIRST PREVALENT PROBLEM: 511 w
- Page 513 and 514: 7.1. FIRST PREVALENT PROBLEM: 513 T
- Page 515 and 516: 7.2. SECOND PREVALENT PROBLEM: 515
- Page 517 and 518: 7.2. SECOND PREVALENT PROBLEM: 517
- Page 519: 7.2. SECOND PREVALENT PROBLEM: 519
- Page 523 and 524: 7.2. SECOND PREVALENT PROBLEM: 523
- Page 525 and 526: 7.3. THIRD PREVALENT PROBLEM: 525 g
- Page 527 and 528: 7.3. THIRD PREVALENT PROBLEM: 527 w
- Page 529 and 530: 7.3. THIRD PREVALENT PROBLEM: 529 7
- Page 531 and 532: 7.3. THIRD PREVALENT PROBLEM: 531 7
- Page 533 and 534: 7.3. THIRD PREVALENT PROBLEM: 533 O
- Page 535 and 536: 7.4. CONCLUSION 535 The rank constr
- Page 537 and 538: Appendix A Linear algebra A.1 Main-
- Page 539 and 540: A.1. MAIN-DIAGONAL δ OPERATOR, λ
- Page 541 and 542: A.1. MAIN-DIAGONAL δ OPERATOR, λ
- Page 543 and 544: A.2. SEMIDEFINITENESS: DOMAIN OF TE
- Page 545 and 546: A.3. PROPER STATEMENTS 545 (AB) T
- Page 547 and 548: A.3. PROPER STATEMENTS 547 A.3.1 Se
- Page 549 and 550: A.3. PROPER STATEMENTS 549 For A di
- Page 551 and 552: A.3. PROPER STATEMENTS 551 Diagonal
- Page 553 and 554: A.3. PROPER STATEMENTS 553 For A,B
- Page 555 and 556: A.3. PROPER STATEMENTS 555 A.3.1.0.
- Page 557 and 558: A.4. SCHUR COMPLEMENT 557 A.4 Schur
- Page 559 and 560: A.4. SCHUR COMPLEMENT 559 A.4.0.0.2
- Page 561 and 562: A.5. EIGEN DECOMPOSITION 561 When B
- Page 563 and 564: A.5. EIGEN DECOMPOSITION 563 dim N(
- Page 565 and 566: A.6. SINGULAR VALUE DECOMPOSITION,
- Page 567 and 568: A.6. SINGULAR VALUE DECOMPOSITION,
- Page 569 and 570: A.7. ZEROS 569 Given symmetric matr
7.2. SECOND PREVALENT PROBLEM: 521<br />
compaction of the reconstruction by minimizing total distance. (823) It is best<br />
used where some physical equilibrium implies such an energy minimization;<br />
e.g., [305,5].<br />
For this Problem 2, the trace rank-heuristic arose naturally in the<br />
objective in terms of V . We observe: V (in contrast to VN T ) spreads energy<br />
over all available distances (B.4.2 no.20, confer no.22) although the rank<br />
function itself is insensitive to choice of auxiliary matrix.<br />
7.2.2.4 Rank minimization heuristic beyond convex envelope<br />
Fazel, Hindi, & Boyd [114] [338] [115] propose a rank heuristic more potent<br />
than trace (1269) for problems of rank minimization;<br />
rankY ← log det(Y +εI) (1274)<br />
the concave surrogate function log det in place of quasiconcave rankY<br />
(2.9.2.6.2) when Y ∈ S n + is variable and where ε is a small positive constant.<br />
They propose minimization of the surrogate by substituting a sequence<br />
comprising infima of a linearized surrogate about the current estimate Y i ;<br />
id est, from the first-order Taylor series expansion about Y i on some open<br />
interval of ‖Y ‖ 2 (D.1.7)<br />
log det(Y + εI) ≈ log det(Y i + εI) + tr ( (Y i + εI) −1 (Y − Y i ) ) (1275)<br />
we make the surrogate sequence of infima over bounded convex feasible set C<br />
where, for i = 0...<br />
arg inf rankY ← lim Y i+1 (1276)<br />
Y ∈ C i→∞<br />
Y i+1 = arg inf<br />
Y ∈ C tr( (Y i + εI) −1 Y ) (1277)<br />
Choosing Y 0 =I , the first step becomes equivalent to finding the infimum of<br />
trY ; the trace rank-heuristic (1269). The intuition underlying (1277) is the<br />
new term in the argument of trace; specifically, (Y i + εI) −1 weights Y so that<br />
relatively small eigenvalues of Y found by the infimum are made even smaller.