v2009.01.01 - Convex Optimization
v2009.01.01 - Convex Optimization v2009.01.01 - Convex Optimization
306 CHAPTER 4. SEMIDEFINITE PROGRAMMING a cardinality-1 constraint in its nonnegative equivalent (464), then this cardinality-1 recursive reconstruction algorithm continues to hold for a signed variable as in (460). This cardinality-1 reconstruction algorithm also holds more generally when affine subset A has any higher dimension n−m . But we leave proof of this curious result pending because it is practically useless; by that we mean, it is more efficient to search over the columns of matrix A for a cardinality-1 solution known a priori to exist. 4.5.2 constraining cardinality of signed variable Now consider a feasibility problem equivalent to the classical problem from linear algebra Ax = b , but with an upper bound k on cardinality ‖x‖ 0 : for vector b∈R(A) find x ∈ R n subject to Ax = b ‖x‖ 0 ≤ k (692) where ‖x‖ 0 ≤ k means vector x has at most k nonzero entries; such a vector is presumed existent in the feasible set. Convex iteration (4.5.1) works with a nonnegative variable; absolute value |x| is therefore needed here. We propose that nonconvex problem (692) can be equivalently written as a sequence of convex problems that move the cardinality constraint to the objective: minimize x∈R n 〈|x| , y〉 subject to Ax = b ≡ minimize 〈t , y + ε1〉 x∈R n , t∈R n subject to Ax = b −t ≼ x ≼ t (693) minimize y∈R n 〈t ⋆ , y + ε1〉 subject to 0 ≼ y ≼ 1 y T 1 = n − k (467) where ε is a relatively small positive constant. This sequence is iterated until a direction vector y is found that makes |x ⋆ | T y ⋆ vanish. The term 〈t , ε1〉 in (693) is necessary to determine absolute value |x ⋆ | = t ⋆ (3.1.3) because vector y can have zero-valued entries. By initializing y to (1−ε)1,
4.5. CONSTRAINING CARDINALITY 307 the first iteration of problem (693) is a 1-norm problem (456); id est, minimize 〈t , 1〉 x∈R n , t∈R n subject to Ax = b −t ≼ x ≼ t ≡ minimize x∈R n ‖x‖ 1 subject to Ax = b (460) Subsequent iterations of problem (693) engaging cardinality term 〈t , y〉 can be interpreted as corrections to this 1-norm problem leading to a 0-norm solution; vector y can be interpreted as a direction of search. 4.5.2.1 local convergence As before (4.5.1.2), convex iteration (693) (467) always converges to a locally optimal solution. 4.5.2.2 simple variations on a signed variable Several useful equivalents to linear programs (693) (467) are easily devised, but their geometrical interpretation is not as apparent: e.g., equivalent in the limit ε→0 + minimize 〈t , y〉 x∈R n , t∈R n subject to Ax = b (694) −t ≼ x ≼ t minimize y∈R n 〈|x ⋆ |, y〉 subject to 0 ≼ y ≼ 1 y T 1 = n − k (467) We get another equivalent to linear programs (693) (467), in the limit, by interpreting problem (460) as infimum to a vertex-description of the 1-norm ball (Figure 61, Example 3.1.3.0.1, confer (459)): minimize a∈R 2n 〈a , y〉 subject to [A −A ]a = b a ≽ 0 minimize y∈R 2n 〈a ⋆ , y〉 subject to 0 ≼ y ≼ 1 y T 1 = 2n − k (467) (695)
- Page 255 and 256: 4.2. FRAMEWORK 255 sets are closed
- Page 257 and 258: 4.2. FRAMEWORK 257 4.2.1.1.3 Exampl
- Page 259 and 260: 4.2. FRAMEWORK 259 4.2.2 Duals The
- Page 261 and 262: 4.2. FRAMEWORK 261 When equality is
- Page 263 and 264: 4.2. FRAMEWORK 263 The pseudoinvers
- 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: 4.5. CONSTRAINING CARDINALITY 305 W
- 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 and 316: 4.6. CARDINALITY AND RANK CONSTRAIN
- Page 317 and 318: 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-
4.5. CONSTRAINING CARDINALITY 307<br />
the first iteration of problem (693) is a 1-norm problem (456); id est,<br />
minimize 〈t , 1〉<br />
x∈R n , t∈R n<br />
subject to Ax = b<br />
−t ≼ x ≼ t<br />
≡<br />
minimize<br />
x∈R n ‖x‖ 1<br />
subject to Ax = b<br />
(460)<br />
Subsequent iterations of problem (693) engaging cardinality term 〈t , y〉 can<br />
be interpreted as corrections to this 1-norm problem leading to a 0-norm<br />
solution; vector y can be interpreted as a direction of search.<br />
4.5.2.1 local convergence<br />
As before (4.5.1.2), convex iteration (693) (467) always converges to a locally<br />
optimal solution.<br />
4.5.2.2 simple variations on a signed variable<br />
Several useful equivalents to linear programs (693) (467) are easily devised,<br />
but their geometrical interpretation is not as apparent: e.g., equivalent in<br />
the limit ε→0 +<br />
minimize 〈t , y〉<br />
x∈R n , t∈R n<br />
subject to Ax = b<br />
(694)<br />
−t ≼ x ≼ t<br />
minimize<br />
y∈R n 〈|x ⋆ |, y〉<br />
subject to 0 ≼ y ≼ 1<br />
y T 1 = n − k<br />
(467)<br />
We get another equivalent to linear programs (693) (467), in the limit, by<br />
interpreting problem (460) as infimum to a vertex-description of the 1-norm<br />
ball (Figure 61, Example 3.1.3.0.1, confer (459)):<br />
minimize<br />
a∈R 2n 〈a , y〉<br />
subject to [A −A ]a = b<br />
a ≽ 0<br />
minimize<br />
y∈R 2n 〈a ⋆ , y〉<br />
subject to 0 ≼ y ≼ 1<br />
y T 1 = 2n − k<br />
(467)<br />
(695)