v2009.01.01 - Convex Optimization

v2009.01.01 - Convex Optimization v2009.01.01 - Convex Optimization

convexoptimization.com
from convexoptimization.com More from this publisher
10.03.2015 Views

316 CHAPTER 4. SEMIDEFINITE PROGRAMMING Assign a rank-1 matrix of variables; symmetric variable matrix X and solution vector x : [ ] [ ] [ ] x [ x G = T 1] X x ∆ xx = 1 x T = T x 1 x T ∈ S N+1 (720) 1 Then design an equivalent semidefinite feasibility problem to find a Boolean solution to Ax ≼ b : find X∈S N x ∈ R N subject to Ax ≼ b G = [ X x x T 1 rankG = 1 (G ≽ 0) δ(X) = 1 ] (721) where x ⋆ i ∈ {−1, 1}, i=1... N . The two variables X and x are made dependent via their assignment to rank-1 matrix G . By (1505), an optimal rank-1 matrix G ⋆ must take the form (720). As before, we regularize the rank constraint by introducing a direction matrix Y into the objective: minimize 〈G, Y 〉 X∈S N , x∈R N subject to Ax ≼ b [ X x G = x T 1 δ(X) = 1 ] ≽ 0 (722) Solution of this semidefinite program is iterated with calculation of the direction matrix Y from semidefinite program (700). At convergence, in the sense (671), convex problem (722) becomes equivalent to nonconvex Boolean problem (719). Direction matrix Y can be an orthogonal projector having closed-form expression, by (1581a), although convex iteration is not a projection method. Given randomized data A and b for a large problem, we find that stalling becomes likely (convergence of the iteration to a positive fixed point

4.6. CARDINALITY AND RANK CONSTRAINT EXAMPLES 317 〈G ⋆ , Y 〉). To overcome this behavior, we introduce a heuristic into the implementation on Wıκımization that momentarily reverses direction of search (like (701)) upon stall detection. We find that rate of convergence can be sped significantly by detecting stalls early. 4.6.0.0.7 Example. Variable-vector normalization. Suppose, within some convex optimization problem, we want vector variables x, y ∈ R N constrained by a nonconvex equality: x‖y‖ = y (723) id est, ‖x‖ = 1 and x points in the same direction as y≠0 ; e.g., minimize f(x, y) x , y subject to (x, y) ∈ C x‖y‖ = y (724) where f is some convex function and C is some convex set. We can realize the nonconvex equality by constraining rank and adding a regularization term to the objective. Make the assignment: ⎡ ⎤ ⎡ x [x T y T 1] G = ⎣ y ⎦ = ⎣ 1 X Z x Z Y y x T y T 1 ⎤ ⎡ ⎦= ∆ ⎣ xx T xy T x yx T yy T y x T y T 1 ⎤ ⎦∈ S 2N+1 (725) where X , Y ∈ S N , also Z ∈ S N [sic] . Any rank-1 solution must take the form of (725). (B.1) The problem statement equivalent to (724) is then written minimize X , Y , Z , x , y subject to (x, y) ∈ C ⎡ f(x, y) + ‖X − Y ‖ F G = ⎣ rankG = 1 (G ≽ 0) tr(X) = 1 δ(Z) ≽ 0 X Z x Z Y y x T y T 1 ⎤ ⎦ (726)

4.6. CARDINALITY AND RANK CONSTRAINT EXAMPLES 317<br />

〈G ⋆ , Y 〉). To overcome this behavior, we introduce a heuristic into the<br />

implementation on Wıκımization that momentarily reverses direction of<br />

search (like (701)) upon stall detection. We find that rate of convergence<br />

can be sped significantly by detecting stalls early.<br />

<br />

4.6.0.0.7 Example. Variable-vector normalization.<br />

Suppose, within some convex optimization problem, we want vector variables<br />

x, y ∈ R N constrained by a nonconvex equality:<br />

x‖y‖ = y (723)<br />

id est, ‖x‖ = 1 and x points in the same direction as y≠0 ; e.g.,<br />

minimize f(x, y)<br />

x , y<br />

subject to (x, y) ∈ C<br />

x‖y‖ = y<br />

(724)<br />

where f is some convex function and C is some convex set. We can realize<br />

the nonconvex equality by constraining rank and adding a regularization<br />

term to the objective. Make the assignment:<br />

⎡ ⎤ ⎡<br />

x [x T y T 1]<br />

G = ⎣ y ⎦ = ⎣<br />

1<br />

X Z x<br />

Z Y y<br />

x T y T 1<br />

⎤<br />

⎡<br />

⎦=<br />

∆ ⎣<br />

xx T xy T x<br />

yx T yy T y<br />

x T y T 1<br />

⎤<br />

⎦∈ S 2N+1 (725)<br />

where X , Y ∈ S N , also Z ∈ S N [sic] . Any rank-1 solution must take the<br />

form of (725). (B.1) The problem statement equivalent to (724) is then<br />

written<br />

minimize<br />

X , Y , Z , x , y<br />

subject to (x, y) ∈ C<br />

⎡<br />

f(x, y) + ‖X − Y ‖ F<br />

G = ⎣<br />

rankG = 1<br />

(G ≽ 0)<br />

tr(X) = 1<br />

δ(Z) ≽ 0<br />

X Z x<br />

Z Y y<br />

x T y T 1<br />

⎤<br />

⎦<br />

(726)

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

Saved successfully!

Ooh no, something went wrong!