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.

368 CHAPTER 5. EUCLIDEAN DISTANCE MATRIX<br />

This preceding transformation of trilateration to a semidefinite program<br />

works all the time ((845) holds) despite relaxation (842) because the optimal<br />

solution set is a unique point.<br />

5.4.2.2.5 Proof (sketch). Only the sensor location x 1 is unknown.<br />

The objective function together with the equality constraints make a linear<br />

system of equations in Gram matrix variable G<br />

trG = ‖x 1 ‖ 2 + ‖ˇx 2 ‖ 2 + ‖ˇx 3 ‖ 2 + ‖ˇx 4 ‖ 2<br />

tr(GΦ i1 = ďi1 , i = 2, 3, 4<br />

tr ( ) (846)<br />

Ge i e T i = ‖ˇxi ‖ 2 , i = 2, 3, 4<br />

tr(G(e i e T j + e j e T i )/2) = ˇx T i ˇx j , 2 ≤ i < j = 3, 4<br />

which is invertible:<br />

⎡<br />

⎤<br />

svec(I) T<br />

svec(Φ 21 ) T<br />

svec(Φ 31 ) T<br />

svec(Φ 41 ) T<br />

svec(e 2 e T 2 ) T<br />

svec G =<br />

svec(e 3 e T 3 ) T<br />

svec(e 4 e T 4 ) T<br />

svec ( (e<br />

⎢ 2 e T 3 + e 3 e T 2 )/2 ) T<br />

⎣ svec ( (e 2 e T 4 + e 4 e T 2 )/2 ) T ⎥<br />

svec ( (e 3 e T 4 + e 4 e T 3 )/2 ) ⎦<br />

T<br />

−1⎡<br />

‖x 1 ‖ 2 + ‖ˇx 2 ‖ 2 + ‖ˇx 3 ‖ 2 + ‖ˇx 4 ‖ 2 ⎤<br />

ď 21<br />

ď 31<br />

ď 41<br />

‖ˇx 2 ‖ 2<br />

‖ˇx 3 ‖ 2<br />

‖ˇx 4 ‖ 2<br />

⎢ ˇx T 2ˇx 3<br />

⎥<br />

⎣ ˇx T 2ˇx 4<br />

⎦<br />

ˇx T 3ˇx 4<br />

(847)<br />

That line in the ambient space S 4 of G , claimed on page 367, is traced by<br />

‖x 1 ‖ 2 ∈ R on the right-hand side, as it turns out. One must show this line<br />

to be tangential (2.1.7.2) to S 4 + in order to prove uniqueness. Tangency is<br />

possible for affine dimension 1 or 2 while its occurrence depends completely<br />

on the known measurement data.<br />

<br />

But as soon as significant noise is introduced or whenever distance data is<br />

incomplete, such problems can remain convex although the set of all optimal<br />

solutions generally becomes a convex set bigger than a single point (and still<br />

containing the noiseless solution).

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

Saved successfully!

Ooh no, something went wrong!