v2010.10.26 - Convex Optimization

v2010.10.26 - Convex Optimization v2010.10.26 - Convex Optimization

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

194 CHAPTER 2. CONVEX GEOMETRY2.13.9 Formulae finding dual cone2.13.9.1 Pointed K , dual, X skinny-or-square full-rankWe wish to derive expressions for a convex cone and its ordinary dualunder the general assumptions: pointed polyhedral K denoted by its linearlyindependent extreme directions arranged columnar in matrix X such thatThe vertex-description is given:rank(X ∈ R n×N ) = N dim aff K ≤ n (412)K = {Xa | a ≽ 0} ⊆ R n (413)from which a halfspace-description for the dual cone follows directly:By defining a matrixK ∗ = {y ∈ R n | X T y ≽ 0} (414)X ⊥ basis N(X T ) (415)(a columnar basis for the orthogonal complement of R(X)), we can sayaff cone X = aff K = {x | X ⊥T x = 0} (416)meaning K lies in a subspace, perhaps R n . Thus a halfspace-descriptionand a vertex-description 2.77 from (314)K = {x∈ R n | X † x ≽ 0, X ⊥T x = 0} (417)K ∗ = { [X †T X ⊥ −X ⊥ ]b | b ≽ 0 } ⊆ R n (418)These results are summarized for a pointed polyhedral cone, havinglinearly independent generators, and its ordinary dual:Cone Table 1 K K ∗vertex-description X X †T , ±X ⊥halfspace-description X † , X ⊥T X T2.77 These descriptions are not unique. A vertex-description of the dual cone, for example,might use four conically independent generators for a plane (2.10.0.0.1, Figure 49) whenonly three would suffice.

2.13. DUAL CONE & GENERALIZED INEQUALITY 1952.13.9.2 Simplicial caseWhen a convex cone is simplicial (2.12.3), Cone Table 1 simplifies becausethen aff coneX = R n : For square X and assuming simplicial K such thatrank(X ∈ R n×N ) = N dim aff K = n (419)we haveCone Table S K K ∗vertex-description X X †Thalfspace-description X † X TFor example, vertex-description (418) simplifies toK ∗ = {X †T b | b ≽ 0} ⊂ R n (420)Now, because dim R(X)= dim R(X †T ) , (E) the dual cone K ∗ is simplicialwhenever K is.2.13.9.3 Cone membership relations in a subspaceIt is obvious by definition (297) of ordinary dual cone K ∗ , in ambient vectorspace R , that its determination instead in subspace S ⊆ R is identicalto its intersection with S ; id est, assuming closed convex cone K ⊆ S andK ∗ ⊆ R(K ∗ were ambient S) ≡ (K ∗ in ambient R) ∩ S (421)because{y ∈ S | 〈y , x〉 ≥ 0 for all x ∈ K} = {y ∈ R | 〈y , x〉 ≥ 0 for all x ∈ K} ∩ S(422)From this, a constrained membership relation for the ordinary dual coneK ∗ ⊆ R , assuming x,y ∈ S and closed convex cone K ⊆ Sy ∈ K ∗ ∩ S ⇔ 〈y , x〉 ≥ 0 for all x ∈ K (423)By closure in subspace S we have conjugation (2.13.1.1):x ∈ K ⇔ 〈y , x〉 ≥ 0 for all y ∈ K ∗ ∩ S (424)

194 CHAPTER 2. CONVEX GEOMETRY2.13.9 Formulae finding dual cone2.13.9.1 Pointed K , dual, X skinny-or-square full-rankWe wish to derive expressions for a convex cone and its ordinary dualunder the general assumptions: pointed polyhedral K denoted by its linearlyindependent extreme directions arranged columnar in matrix X such thatThe vertex-description is given:rank(X ∈ R n×N ) = N dim aff K ≤ n (412)K = {Xa | a ≽ 0} ⊆ R n (413)from which a halfspace-description for the dual cone follows directly:By defining a matrixK ∗ = {y ∈ R n | X T y ≽ 0} (414)X ⊥ basis N(X T ) (415)(a columnar basis for the orthogonal complement of R(X)), we can sayaff cone X = aff K = {x | X ⊥T x = 0} (416)meaning K lies in a subspace, perhaps R n . Thus a halfspace-descriptionand a vertex-description 2.77 from (314)K = {x∈ R n | X † x ≽ 0, X ⊥T x = 0} (417)K ∗ = { [X †T X ⊥ −X ⊥ ]b | b ≽ 0 } ⊆ R n (418)These results are summarized for a pointed polyhedral cone, havinglinearly independent generators, and its ordinary dual:Cone Table 1 K K ∗vertex-description X X †T , ±X ⊥halfspace-description X † , X ⊥T X T2.77 These descriptions are not unique. A vertex-description of the dual cone, for example,might use four conically independent generators for a plane (2.10.0.0.1, Figure 49) whenonly three would suffice.

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

Saved successfully!

Ooh no, something went wrong!