10.03.2015 Views

v2009.01.01 - Convex Optimization

v2009.01.01 - Convex Optimization

v2009.01.01 - Convex Optimization

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

314 CHAPTER 4. SEMIDEFINITE PROGRAMMING<br />

constraints in (711). (confer (707)) When a permutation A of vector B exists,<br />

number of iterations can be as small as 1. But this combinatorial Procrustes<br />

problem can be made even more challenging when vector A has repeated<br />

entries.<br />

solution via cardinality constraint<br />

Now the idea is to force solution at a vertex of permutation polyhedron (91)<br />

by finding a solution of desired sparsity. Because permutation matrix X is<br />

n-sparse by assumption, this combinatorial Procrustes problem may instead<br />

be formulated as a compressed sensing problem with convex iteration on<br />

cardinality of vectorized X (4.5.1): given nonzero vectors A, B<br />

minimize<br />

X∈R n×n ‖A − XB‖ 1 + w〈X , Y 〉<br />

subject to X T 1 = 1<br />

X1 = 1<br />

X ≥ 0<br />

(714)<br />

where direction vector Y is an optimal solution to<br />

minimize<br />

Y ∈R n×n 〈X ⋆ , Y 〉<br />

subject to 0 ≤ Y ≤ 1<br />

1 T Y 1 = n 2 − n<br />

(467)<br />

each a linear program. In this circumstance, use of closed-form solution for<br />

direction vector Y is discouraged. When vector A is a permutation of B ,<br />

both linear programs have objectives that converge to 0. When vectors A and<br />

B are permutations and no entries of A are repeated, optimal solution X ⋆<br />

can be found as soon as the first iteration.<br />

In any case, X ⋆ = Ξ is a permutation matrix.<br />

<br />

4.6.0.0.4 Exercise. Combinatorial Procrustes constraints.<br />

Assume that the objective of semidefinite program (711) is 0 at optimality.<br />

Prove that the constraints in program (711) are necessary and sufficient<br />

to produce a permutation matrix as optimal solution. Alternatively and<br />

equivalently, prove those constraints necessary and sufficient to optimally<br />

produce a nonnegative orthogonal matrix.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!