10.03.2015 Views

v2009.01.01 - Convex Optimization

v2009.01.01 - Convex Optimization

v2009.01.01 - Convex Optimization

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

5.5. INVARIANCE 383<br />

5.5.1.0.1 Example. Translating geometric center to origin.<br />

We might choose to shift the geometric center α c of an N-point list {x l }<br />

(arranged columnar in X) to the origin; [305] [137]<br />

α = α ∆ c = Xb ∆ c = 1 X1 ∈ P ⊆ A (878)<br />

N<br />

where A represents the list’s affine hull. If we were to associate a point-mass<br />

m l with each of the points x l in the list, then their center of mass<br />

(or gravity) would be ( ∑ x l m l )/ ∑ m l . The geometric center is the same<br />

as the center of mass under the assumption of uniform mass density across<br />

points. [190] The geometric center always lies in the convex hull P of the list;<br />

id est, α c ∈ P because bc T 1=1 and b c ≽ 0 . 5.22 Subtracting the geometric<br />

center from every list member,<br />

X − α c 1 T = X − 1 N X11T = X(I − 1 N 11T ) = XV ∈ R n×N (879)<br />

where V is the geometric centering matrix (821). So we have (confer (798))<br />

D(X) = D(XV ) = δ(V T X T XV )1 T + 1δ(V T X T XV ) T − 2V T X T XV ∈ EDM N<br />

5.5.1.1 Gram-form invariance<br />

Following from (880) and the linear Gram-form EDM operator (810):<br />

(880)<br />

<br />

D(G) = D(V GV ) = δ(V GV )1 T + 1δ(V GV ) T − 2V GV ∈ EDM N (881)<br />

The Gram-form consequently exhibits invariance to translation by a doublet<br />

(B.2) u1 T + 1u T ;<br />

D(G) = D(G − (u1 T + 1u T )) (882)<br />

because, for any u∈ R N , D(u1 T + 1u T )=0. The collection of all such<br />

doublets forms the nullspace (898) to the operator; the translation-invariant<br />

subspace S N⊥<br />

c (1876) of the symmetric matrices S N . This means matrix G<br />

can belong to an expanse more broad than a positive semidefinite cone;<br />

id est, G∈ S N + − S N⊥<br />

c . So explains Gram matrix sufficiency in EDM<br />

definition (810). 5.23<br />

5.22 Any b from α = Xb chosen such that b T 1 = 1, more generally, makes an auxiliary<br />

V -matrix. (B.4.5)<br />

5.23 A constraint G1=0 would prevent excursion into the translation-invariant subspace<br />

(numerical unboundedness).

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

Saved successfully!

Ooh no, something went wrong!