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

260 CHAPTER 4. SEMIDEFINITE PROGRAMMING 4.2.3.0.1 Corollary. Optimality and strong duality. [312,3.1] [342,1.3.8] For semidefinite programs (584P) and (584D), assume primal and dual feasible sets A ∩ S n + ⊂ S n and C ∗ ⊂ S n × R m (596) are nonempty. ThenX ⋆ is optimal for (P) S ⋆ , y ⋆ are optimal for (D) the duality gap 〈C,X ⋆ 〉−〈b, y ⋆ 〉 is 0 if and only if i) ∃X ∈ A ∩ int S n + or ∃S , y ∈ rel int C ∗ ii) 〈S ⋆ , X ⋆ 〉 = 0 and ⋄ For symmetric positive semidefinite matrices, requirement ii is equivalent to the complementarity (A.7.4) 〈S ⋆ , X ⋆ 〉 = 0 ⇔ S ⋆ X ⋆ = X ⋆ S ⋆ = 0 (609) Commutativity of diagonalizable matrices is a necessary and sufficient condition [176,1.3.12] for these two optimal symmetric matrices to be simultaneously diagonalizable. Therefore rankX ⋆ + rankS ⋆ ≤ n (610) Proof. To see that, the product of symmetric optimal matrices X ⋆ , S ⋆ ∈ S n must itself be symmetric because of commutativity. (1378) The symmetric product has diagonalization [10, cor.2.11] S ⋆ X ⋆ = X ⋆ S ⋆ = QΛ S ⋆Λ X ⋆Q T = 0 ⇔ Λ X ⋆Λ S ⋆ = 0 (611) where Q is an orthogonal matrix. The product of the nonnegative diagonal Λ matrices can be 0 if their main diagonal zeros are complementary or coincide. Due only to symmetry, rankX ⋆ = rank Λ X ⋆ and rankS ⋆ = rank Λ S ⋆ for these optimal primal and dual solutions. (1363) So, because of the complementarity, the total number of nonzero diagonal entries from both Λ cannot exceed n .

4.2. FRAMEWORK 261 When equality is attained in (610) rankX ⋆ + rankS ⋆ = n (612) there are no coinciding main diagonal zeros in Λ X ⋆Λ S ⋆ , and so we have what is called strict complementarity. 4.11 Logically it follows that a necessary and sufficient condition for strict complementarity of an optimal primal and dual solution is X ⋆ + S ⋆ ≻ 0 (613) 4.2.3.1 solving primal problem via dual The beauty of Corollary 4.2.3.0.1 is its conjugacy; id est, one can solve either the primal or dual problem in (584) and then find a solution to the other via the optimality conditions. When a dual optimal solution is known, for example, a primal optimal solution belongs to the hyperplane {X | 〈S ⋆ , X〉=0}. 4.2.3.1.1 Example. Minimum cardinality Boolean. [81] [31,4.3.4] [300] (confer Example 4.5.1.2.1) Consider finding a minimum cardinality Boolean solution x to the classic linear algebra problem Ax = b given noiseless data A∈ R m×n and b∈ R m ; minimize ‖x‖ 0 x subject to Ax = b x i ∈ {0, 1} , i=1... n (614) where ‖x‖ 0 denotes cardinality of vector x (a.k.a, 0-norm; not a convex function). A minimum cardinality solution answers the question: “Which fewest linear combination of columns in A constructs vector b ?” Cardinality problems have extraordinarily wide appeal, arising in many fields of science and across many disciplines. [276] [185] [144] [145] Yet designing an efficient algorithm to optimize cardinality has proved difficult. In this example, we also constrain the variable to be Boolean. The Boolean constraint forces 4.11 distinct from maximal complementarity (4.1.1.1).

260 CHAPTER 4. SEMIDEFINITE PROGRAMMING<br />

4.2.3.0.1 Corollary. Optimality and strong duality. [312,3.1]<br />

[342,1.3.8] For semidefinite programs (584P) and (584D), assume primal<br />

and dual feasible sets A ∩ S n + ⊂ S n and C ∗ ⊂ S n × R m (596) are nonempty.<br />

ThenX ⋆ is optimal for (P)<br />

S ⋆ , y ⋆ are optimal for (D)<br />

the duality gap 〈C,X ⋆ 〉−〈b, y ⋆ 〉 is 0<br />

if and only if<br />

i) ∃X ∈ A ∩ int S n + or ∃S , y ∈ rel int C ∗<br />

ii) 〈S ⋆ , X ⋆ 〉 = 0<br />

and<br />

⋄<br />

For symmetric positive semidefinite matrices, requirement ii is equivalent<br />

to the complementarity (A.7.4)<br />

〈S ⋆ , X ⋆ 〉 = 0 ⇔ S ⋆ X ⋆ = X ⋆ S ⋆ = 0 (609)<br />

Commutativity of diagonalizable matrices is a necessary and sufficient<br />

condition [176,1.3.12] for these two optimal symmetric matrices to be<br />

simultaneously diagonalizable. Therefore<br />

rankX ⋆ + rankS ⋆ ≤ n (610)<br />

Proof. To see that, the product of symmetric optimal matrices<br />

X ⋆ , S ⋆ ∈ S n must itself be symmetric because of commutativity. (1378) The<br />

symmetric product has diagonalization [10, cor.2.11]<br />

S ⋆ X ⋆ = X ⋆ S ⋆ = QΛ S ⋆Λ X ⋆Q T = 0 ⇔ Λ X ⋆Λ S ⋆ = 0 (611)<br />

where Q is an orthogonal matrix. The product of the nonnegative diagonal Λ<br />

matrices can be 0 if their main diagonal zeros are complementary or coincide.<br />

Due only to symmetry, rankX ⋆ = rank Λ X ⋆ and rankS ⋆ = rank Λ S ⋆ for<br />

these optimal primal and dual solutions. (1363) So, because of the<br />

complementarity, the total number of nonzero diagonal entries from both Λ<br />

cannot exceed n .

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

Saved successfully!

Ooh no, something went wrong!