v2009.01.01 - Convex Optimization
v2009.01.01 - Convex Optimization v2009.01.01 - Convex Optimization
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)
- Page 265 and 266: 4.2. FRAMEWORK 265 For the data giv
- Page 267 and 268: 4.2. FRAMEWORK 267 minimizes an aff
- Page 269 and 270: 4.3. RANK REDUCTION 269 whose rank
- Page 271 and 272: 4.3. RANK REDUCTION 271 and where m
- Page 273 and 274: 4.3. RANK REDUCTION 273 4.3.3 Optim
- Page 275 and 276: 4.3. RANK REDUCTION 275 Initialize:
- Page 277 and 278: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 279 and 280: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 281 and 282: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 283 and 284: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 285 and 286: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 287 and 288: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 289 and 290: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 291 and 292: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 293 and 294: 4.4. RANK-CONSTRAINED SEMIDEFINITE
- Page 295 and 296: 4.5. CONSTRAINING CARDINALITY 295 m
- Page 297 and 298: 4.5. CONSTRAINING CARDINALITY 297 m
- Page 299 and 300: 4.5. CONSTRAINING CARDINALITY 299 a
- Page 301 and 302: 4.5. CONSTRAINING CARDINALITY 301 f
- Page 303 and 304: 4.5. CONSTRAINING CARDINALITY 303 n
- Page 305 and 306: 4.5. CONSTRAINING CARDINALITY 305 W
- Page 307 and 308: 4.5. CONSTRAINING CARDINALITY 307 t
- Page 309 and 310: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 311 and 312: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 313 and 314: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 315: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 319 and 320: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 321 and 322: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 323 and 324: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 325 and 326: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 327 and 328: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 329 and 330: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 331 and 332: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 333 and 334: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 335 and 336: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 337 and 338: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 339 and 340: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 341 and 342: 4.7. CONVEX ITERATION RANK-1 341 fi
- Page 343 and 344: 4.7. CONVEX ITERATION RANK-1 343 Gi
- Page 345 and 346: Chapter 5 Euclidean Distance Matrix
- Page 347 and 348: 5.2. FIRST METRIC PROPERTIES 347 co
- Page 349 and 350: 5.3. ∃ FIFTH EUCLIDEAN METRIC PRO
- Page 351 and 352: 5.3. ∃ FIFTH EUCLIDEAN METRIC PRO
- Page 353 and 354: 5.4. EDM DEFINITION 353 The collect
- Page 355 and 356: 5.4. EDM DEFINITION 355 5.4.2 Gram-
- Page 357 and 358: 5.4. EDM DEFINITION 357 D ∈ EDM N
- Page 359 and 360: 5.4. EDM DEFINITION 359 5.4.2.2.1 E
- Page 361 and 362: 5.4. EDM DEFINITION 361 ten affine
- Page 363 and 364: 5.4. EDM DEFINITION 363 spheres: Th
- Page 365 and 366: 5.4. EDM DEFINITION 365 By eliminat
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)