v2010.10.26 - Convex Optimization
v2010.10.26 - Convex Optimization v2010.10.26 - Convex Optimization
358 CHAPTER 4. SEMIDEFINITE PROGRAMMINGcolumn-sum of X must also be unity. It is this combination of nonnegativity,sum, and sum square constraints that extracts the permutation matrices:(Figure 105) given nonzero vectors A, BminimizeX∈R n×n , G i ∈ S n+1subject to∑‖A − XB‖ 1 + w n 〈G i , W i 〉[ i=1 ] ⎫Gi (1:n, 1:n) X(:, i) ⎬G i =X(:, i) T ≽ 01⎭ ,trG i = 2i=1... n(794)X T 1 = 1X1 = 1X ≥ 0where w ≈10 positively weights the rank regularization term. Optimalsolutions G ⋆ i are key to finding direction matrices W i for the next iterationof semidefinite programs (794) (795):minimizeW i ∈ S n+1 〈G ⋆ i , W i 〉subject to 0 ≼ W i ≼ ItrW i = n⎫⎪⎬, i=1... n (795)⎪⎭Direction matrices thus found lead toward rank-1 matrices G ⋆ i on subsequentiterations. Constraint on trace of G ⋆ i normalizes the i th column of X ⋆ to unitybecause (confer p.419)[ ] XG ⋆ i =⋆ (:, i) [ X ⋆ (:, i) T 1]1(796)at convergence. Binary-valued X ⋆ column entries result from the furthersum constraint X1=1. Columnar orthogonality is a consequence of thefurther transpose-sum constraint X T 1=1 in conjunction with nonnegativityconstraint X ≥ 0 ; but we leave proof of orthogonality an exercise. Theoptimal objective value is 0 for both semidefinite programs when vectors Aand B are related by permutation. In any case, optimal solution X ⋆ becomesa permutation matrix Ξ .
4.6. CARDINALITY AND RANK CONSTRAINT EXAMPLES 359Because there are n direction matrices W i to find, it can be advantageousto invoke a known closed-form solution for each from page 658. What makesthis combinatorial problem more tractable are relatively small semidefiniteconstraints in (794). (confer (790)) When a permutation A of vector B exists,number of iterations can be as small as 1. But this combinatorial Procrustesproblem can be made even more challenging when vector A has repeatedentries.solution via cardinality constraintNow the idea is to force solution at a vertex of permutation polyhedron (100)by finding a solution of desired sparsity. Because permutation matrix X isn-sparse by assumption, this combinatorial Procrustes problem may insteadbe formulated as a compressed sensing problem with convex iteration oncardinality of vectorized X (4.5.1): given nonzero vectors A, BminimizeX∈R n×n ‖A − XB‖ 1 + w〈X , Y 〉subject to X T 1 = 1X1 = 1X ≥ 0where direction vector Y is an optimal solution tominimizeY ∈R n×n 〈X ⋆ , Y 〉subject to 0 ≤ Y ≤ 11 T Y 1 = n 2 − n(525)(797)each a linear program. In this circumstance, use of closed-form solution fordirection vector Y is discouraged. When vector A is a permutation of B ,both linear programs have objectives that converge to 0. When vectors A andB are permutations and no entries of A are repeated, optimal solution X ⋆can be found as soon as the first iteration.In any case, X ⋆ = Ξ is a permutation matrix.4.6.0.0.4 Exercise. Combinatorial Procrustes constraints.Assume that the objective of semidefinite program (794) is 0 at optimality.Prove that the constraints in program (794) are necessary and sufficientto produce a permutation matrix as optimal solution. Alternatively andequivalently, prove those constraints necessary and sufficient to optimallyproduce a nonnegative orthogonal matrix.
- Page 307 and 308: 4.3. RANK REDUCTION 3074.3.3.0.2 Ex
- Page 309 and 310: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 311 and 312: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 313 and 314: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 315 and 316: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 317 and 318: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 319 and 320: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 321 and 322: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 323 and 324: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 325 and 326: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 327 and 328: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 329 and 330: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 331 and 332: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 333 and 334: 4.5. CONSTRAINING CARDINALITY 333mi
- Page 335 and 336: 4.5. CONSTRAINING CARDINALITY 3350R
- Page 337 and 338: 4.5. CONSTRAINING CARDINALITY 337it
- Page 339 and 340: 4.5. CONSTRAINING CARDINALITY 339m/
- Page 341 and 342: 4.5. CONSTRAINING CARDINALITY 341we
- Page 343 and 344: 4.5. CONSTRAINING CARDINALITY 343fl
- Page 345 and 346: 4.5. CONSTRAINING CARDINALITY 345We
- Page 347 and 348: 4.5. CONSTRAINING CARDINALITY 3474.
- Page 349 and 350: 4.5. CONSTRAINING CARDINALITY 349R
- Page 351 and 352: 4.5. CONSTRAINING CARDINALITY 351pe
- Page 353 and 354: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 355 and 356: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 357: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 361 and 362: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 363 and 364: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 365 and 366: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 367 and 368: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 369 and 370: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 371 and 372: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 373 and 374: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 375 and 376: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 377 and 378: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 379 and 380: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 381 and 382: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 383 and 384: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 385 and 386: 4.7. CONSTRAINING RANK OF INDEFINIT
- Page 387 and 388: 4.7. CONSTRAINING RANK OF INDEFINIT
- Page 389 and 390: 4.7. CONSTRAINING RANK OF INDEFINIT
- Page 391 and 392: 4.8. CONVEX ITERATION RANK-1 391whi
- Page 393 and 394: 4.8. CONVEX ITERATION RANK-1 393the
- Page 395 and 396: Chapter 5Euclidean Distance MatrixT
- Page 397 and 398: 5.2. FIRST METRIC PROPERTIES 397to
- Page 399 and 400: 5.3. ∃ FIFTH EUCLIDEAN METRIC PRO
- Page 401 and 402: 5.3. ∃ FIFTH EUCLIDEAN METRIC PRO
- Page 403 and 404: 5.4. EDM DEFINITION 403The collecti
- Page 405 and 406: 5.4. EDM DEFINITION 4055.4.2 Gram-f
- Page 407 and 408: 5.4. EDM DEFINITION 407We provide a
4.6. CARDINALITY AND RANK CONSTRAINT EXAMPLES 359Because there are n direction matrices W i to find, it can be advantageousto invoke a known closed-form solution for each from page 658. What makesthis combinatorial problem more tractable are relatively small semidefiniteconstraints in (794). (confer (790)) When a permutation A of vector B exists,number of iterations can be as small as 1. But this combinatorial Procrustesproblem can be made even more challenging when vector A has repeatedentries.solution via cardinality constraintNow the idea is to force solution at a vertex of permutation polyhedron (100)by finding a solution of desired sparsity. Because permutation matrix X isn-sparse by assumption, this combinatorial Procrustes problem may insteadbe formulated as a compressed sensing problem with convex iteration oncardinality of vectorized X (4.5.1): given nonzero vectors A, BminimizeX∈R n×n ‖A − XB‖ 1 + w〈X , Y 〉subject to X T 1 = 1X1 = 1X ≥ 0where direction vector Y is an optimal solution tominimizeY ∈R n×n 〈X ⋆ , Y 〉subject to 0 ≤ Y ≤ 11 T Y 1 = n 2 − n(525)(797)each a linear program. In this circumstance, use of closed-form solution fordirection vector Y is discouraged. When vector A is a permutation of B ,both linear programs have objectives that converge to 0. When vectors A andB are permutations and no entries of A are repeated, optimal solution X ⋆can be found as soon as the first iteration.In any case, X ⋆ = Ξ is a permutation matrix.4.6.0.0.4 Exercise. Combinatorial Procrustes constraints.Assume that the objective of semidefinite program (794) is 0 at optimality.Prove that the constraints in program (794) are necessary and sufficientto produce a permutation matrix as optimal solution. Alternatively andequivalently, prove those constraints necessary and sufficient to optimallyproduce a nonnegative orthogonal matrix.