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]

410 CHAPTER 5. EUCLIDEAN DISTANCE MATRIX<br />

In terms of V N , the unit simplex (265) in R N−1 has an equivalent<br />

representation:<br />

S = {s ∈ R N−1 | √ 2V N s ≽ −e 1 } (987)<br />

where e 1 is as in (814). Incidental to the EDM assertion, shifting the<br />

unit-simplex domain in R N−1 translates the polyhedron P in R n . Indeed,<br />

there is a map from vertices of the unit simplex to members of the list<br />

generating P ;<br />

p : R N−1<br />

⎛⎧<br />

⎪⎨<br />

p<br />

⎜<br />

⎝<br />

⎪⎩<br />

−β<br />

e 1 − β<br />

e 2 − β<br />

.<br />

e N−1 − β<br />

⎫⎞<br />

⎪⎬<br />

⎟<br />

⎠<br />

⎪⎭<br />

→<br />

=<br />

⎧<br />

⎪⎨<br />

⎪⎩<br />

R n<br />

x 1 − α<br />

x 2 − α<br />

x 3 − α<br />

.<br />

x N − α<br />

⎫<br />

⎪⎬<br />

⎪⎭<br />

(988)<br />

5.9.1.0.5 Proof. EDM assertion.<br />

(⇒) We demonstrate that if D is an EDM, then each distance-square ‖p(y)‖ 2<br />

described by (985) corresponds to a point p in some embedded polyhedron<br />

P − α . Assume D is indeed an EDM; id est, D can be made from some list<br />

X of N unknown points in Euclidean space R n ; D = D(X) for X ∈ R n×N<br />

as in (798). Since D is translation invariant (5.5.1), we may shift the affine<br />

hull A of those unknown points to the origin as in (928). Then take any<br />

point p in their convex hull (78);<br />

P − α = {p = (X − Xb1 T )a | a T 1 = 1, a ≽ 0} (989)<br />

where α = Xb ∈ A ⇔ b T 1 = 1. Solutions to a T 1 = 1 are: 5.43<br />

a ∈<br />

{e 1 + √ }<br />

2V N s | s ∈ R N−1<br />

(990)<br />

where e 1 is as in (814). Similarly, b = e 1 + √ 2V N β .<br />

P − α = {p = X(I − (e 1 + √ 2V N β)1 T )(e 1 + √ 2V N s) | √ 2V N s ≽ −e 1 }<br />

= {p = X √ 2V N (s − β) | √ 2V N s ≽ −e 1 } (991)<br />

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<br />

hyperplane through the origin orthogonal to 1. Thus the solutions {a} constitute a<br />

hyperplane orthogonal to the vector 1, and offset from the origin in R N by any particular<br />

solution; in this case, a = e 1 .

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

Saved successfully!

Ooh no, something went wrong!