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

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)

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)

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

Saved successfully!

Ooh no, something went wrong!