10.03.2015 Views

v2009.01.01 - Convex Optimization

v2009.01.01 - Convex Optimization

v2009.01.01 - Convex Optimization

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

534 CHAPTER 7. PROXIMITY PROBLEMS<br />

verifiable by observing conic dependencies (2.10.3) among the aggregate of<br />

halfspace-description normals. The cone membership constraint in (1311a)<br />

therefore inherently insures existence of a symmetric hollow matrix. <br />

<strong>Optimization</strong> (1311b) would be a Procrustes problem (C.4) were it<br />

not for the hollowness constraint; it is, instead, a minimization over the<br />

intersection of the nonconvex manifold of orthogonal matrices with another<br />

nonconvex set in variable R specified by the hollowness constraint.<br />

We solve problem (1311b) by a method introduced in4.6.0.0.2: Define<br />

R = [r 1 · · · r N+1 ]∈ R N+1×N+1 and make the assignment<br />

⎡<br />

G = ⎢<br />

⎣<br />

r 1<br />

.<br />

r N+1<br />

1<br />

⎤<br />

[ r1 T · · · rN+1 T 1] ⎥<br />

∈ S (N+1)2 +1<br />

⎦<br />

⎡<br />

⎤ ⎡<br />

R 11 · · · R 1,N+1 r 1<br />

.<br />

...<br />

.<br />

∆<br />

= ⎢<br />

⎣ R1,N+1 T ⎥=<br />

⎢<br />

R N+1,N+1 r N+1 ⎦ ⎣<br />

r1 T · · · rN+1 T 1<br />

(1314)<br />

⎤<br />

r 1 r1 T · · · r 1 rN+1 T r 1<br />

.<br />

...<br />

.<br />

r N+1 r1 T r N+1 rN+1 T ⎥<br />

r N+1 ⎦<br />

r1 T · · · rN+1 T 1<br />

where R ∆ ij = r i rj T ∈ R N+1×N+1 and Υ ⋆ ii ∈ R . Since R Υ ⋆ R T = N+1 ∑<br />

ΥiiR ⋆ ii then<br />

problem (1311b) is equivalently expressed:<br />

∥ ∥∥∥ N+1<br />

2<br />

∑<br />

minimize Υ ⋆<br />

R ij , r i<br />

iiR ii − Λ<br />

∥<br />

i=1<br />

F<br />

subject to trR ii = 1, i=1... N+1<br />

trR ij ⎡<br />

= 0,<br />

⎤<br />

i

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

Saved successfully!

Ooh no, something went wrong!