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

410 CHAPTER 5. EUCLIDEAN DISTANCE MATRIX In terms of V N , the unit simplex (265) in R N−1 has an equivalent representation: S = {s ∈ R N−1 | √ 2V N s ≽ −e 1 } (987) where e 1 is as in (814). Incidental to the EDM assertion, shifting the unit-simplex domain in R N−1 translates the polyhedron P in R n . Indeed, there is a map from vertices of the unit simplex to members of the list generating P ; p : R N−1 ⎛⎧ ⎪⎨ p ⎜ ⎝ ⎪⎩ −β e 1 − β e 2 − β . e N−1 − β ⎫⎞ ⎪⎬ ⎟ ⎠ ⎪⎭ → = ⎧ ⎪⎨ ⎪⎩ R n x 1 − α x 2 − α x 3 − α . x N − α ⎫ ⎪⎬ ⎪⎭ (988) 5.9.1.0.5 Proof. EDM assertion. (⇒) We demonstrate that if D is an EDM, then each distance-square ‖p(y)‖ 2 described by (985) corresponds to a point p in some embedded polyhedron P − α . Assume D is indeed an EDM; id est, D can be made from some list X of N unknown points in Euclidean space R n ; D = D(X) for X ∈ R n×N as in (798). Since D is translation invariant (5.5.1), we may shift the affine hull A of those unknown points to the origin as in (928). Then take any point p in their convex hull (78); P − α = {p = (X − Xb1 T )a | a T 1 = 1, a ≽ 0} (989) where α = Xb ∈ A ⇔ b T 1 = 1. Solutions to a T 1 = 1 are: 5.43 a ∈ {e 1 + √ } 2V N s | s ∈ R N−1 (990) where e 1 is as in (814). Similarly, b = e 1 + √ 2V N β . P − α = {p = X(I − (e 1 + √ 2V N β)1 T )(e 1 + √ 2V N s) | √ 2V N s ≽ −e 1 } = {p = X √ 2V N (s − β) | √ 2V N s ≽ −e 1 } (991) 5.43 Since R(V N )= N(1 T ) and N(1 T )⊥ R(1) , then over all s∈ R N−1 , V N s is a hyperplane through the origin orthogonal to 1. Thus the solutions {a} constitute a hyperplane orthogonal to the vector 1, and offset from the origin in R N by any particular solution; in this case, a = e 1 .

5.9. BRIDGE: CONVEX POLYHEDRA TO EDMS 411 that describes the domain of p(s) as the unit simplex Making the substitution s − β ← y S = {s | √ 2V N s ≽ −e 1 } ⊂ R N−1 + (987) P − α = {p = X √ 2V N y | y ∈ S − β} (992) Point p belongs to a convex polyhedron P−α embedded in an r-dimensional subspace of R n because the convex hull of any list forms a polyhedron, and because the translated affine hull A − α contains the translated polyhedron P − α (931) and the origin (when α∈A), and because A has dimension r by definition (933). Now, any distance-square from the origin to the polyhedron P − α can be formulated {p T p = ‖p‖ 2 = 2y T V T NX T XV N y | y ∈ S − β} (993) Applying (940) to (993) we get (985). (⇐) To validate the EDM assertion in the reverse direction, we prove: If each distance-square ‖p(y)‖ 2 (985) on the shifted unit-simplex S −β ⊂ R N−1 corresponds to a point p(y) in some embedded polyhedron P − α , then D is an EDM. The r-dimensional subspace A − α ⊆ R n is spanned by p(S − β) = P − α (994) because A − α = aff(P − α) ⊇ P − α (931). So, outside the domain S − β of linear surjection p(y) , the simplex complement \S − β ⊂ R N−1 must contain the domain of the distance-square ‖p(y)‖ 2 = p(y) T p(y) to remaining points in the subspace A − α ; id est, to the polyhedron’s relative exterior \P − α . For ‖p(y)‖ 2 to be nonnegative on the entire subspace A − α , −V T N DV N must be positive semidefinite and is assumed symmetric; 5.44 −V T NDV N ∆ = Θ T pΘ p (995) where 5.45 Θ p ∈ R m×N−1 for some m ≥ r . Because p(S − β) is a convex polyhedron, it is necessarily a set of linear combinations of points from some 5.44 The antisymmetric part ( −V T N DV N − (−V T N DV N) T) /2 is annihilated by ‖p(y)‖ 2 . By the same reasoning, any positive (semi)definite matrix A is generally assumed symmetric because only the symmetric part (A +A T )/2 survives the test y T Ay ≥ 0. [176,7.1] 5.45 A T = A ≽ 0 ⇔ A = R T R for some real matrix R . [287,6.3]

5.9. BRIDGE: CONVEX POLYHEDRA TO EDMS 411<br />

that describes the domain of p(s) as the unit simplex<br />

Making the substitution s − β ← y<br />

S = {s | √ 2V N s ≽ −e 1 } ⊂ R N−1<br />

+ (987)<br />

P − α = {p = X √ 2V N y | y ∈ S − β} (992)<br />

Point p belongs to a convex polyhedron P−α embedded in an r-dimensional<br />

subspace of R n because the convex hull of any list forms a polyhedron, and<br />

because the translated affine hull A − α contains the translated polyhedron<br />

P − α (931) and the origin (when α∈A), and because A has dimension r by<br />

definition (933). Now, any distance-square from the origin to the polyhedron<br />

P − α can be formulated<br />

{p T p = ‖p‖ 2 = 2y T V T NX T XV N y | y ∈ S − β} (993)<br />

Applying (940) to (993) we get (985).<br />

(⇐) To validate the EDM assertion in the reverse direction, we prove: If<br />

each distance-square ‖p(y)‖ 2 (985) on the shifted unit-simplex S −β ⊂ R N−1<br />

corresponds to a point p(y) in some embedded polyhedron P − α , then D<br />

is an EDM. The r-dimensional subspace A − α ⊆ R n is spanned by<br />

p(S − β) = P − α (994)<br />

because A − α = aff(P − α) ⊇ P − α (931). So, outside the domain S − β<br />

of linear surjection p(y) , the simplex complement \S − β ⊂ R N−1 must<br />

contain the domain of the distance-square ‖p(y)‖ 2 = p(y) T p(y) to remaining<br />

points in the subspace A − α ; id est, to the polyhedron’s relative exterior<br />

\P − α . For ‖p(y)‖ 2 to be nonnegative on the entire subspace A − α ,<br />

−V T N DV N must be positive semidefinite and is assumed symmetric; 5.44<br />

−V T NDV N ∆ = Θ T pΘ p (995)<br />

where 5.45 Θ p ∈ R m×N−1 for some m ≥ r . Because p(S − β) is a convex<br />

polyhedron, it is necessarily a set of linear combinations of points from some<br />

5.44 The antisymmetric part ( −V T N DV N − (−V T N DV N) T) /2 is annihilated by ‖p(y)‖ 2 . By<br />

the same reasoning, any positive (semi)definite matrix A is generally assumed symmetric<br />

because only the symmetric part (A +A T )/2 survives the test y T Ay ≥ 0. [176,7.1]<br />

5.45 A T = A ≽ 0 ⇔ A = R T R for some real matrix R . [287,6.3]

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

Saved successfully!

Ooh no, something went wrong!