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.

670 APPENDIX E. PROJECTION<br />

As for all subspace projectors, range of the projector is the subspace on which<br />

projection is made: {P 1 Y P 2 | Y ∈ R m×p }. Altogether, for projectors P 1 and<br />

P 2 of any rank, this means projection P 1 XP 2 is unique minimum-distance,<br />

orthogonal<br />

P 1 XP 2 − X ⊥ {P 1 Y P 2 | Y ∈ R m×p } in R mp (1870)<br />

and P 1 and P 2 must each be symmetric (confer (1852)) to attain the infimum.<br />

E.7.2.0.1 Proof. Minimum Frobenius norm (1869).<br />

Defining P ∆ = A 1 (A † 1 + B 1 Z T 1 ) ,<br />

inf ‖X − A 1 (A † 1 + B 1 Z1 T )X(A †T<br />

2 + Z 2 B2 T )A T 2 ‖ 2 F<br />

B 1 , B 2<br />

= inf ‖X − PX(A †T<br />

2 + Z 2 B2 T )A T 2 ‖ 2 F<br />

B 1 , B 2<br />

(<br />

)<br />

= inf tr (X T − A 2 (A † 2 + B 2 Z2 T )X T P T )(X − PX(A †T<br />

2 + Z 2 B2 T )A T 2 )<br />

B 1 , B 2<br />

(<br />

= inf tr X T X −X T PX(A †T<br />

2 +Z 2 B2 T )A T 2 −A 2 (A †<br />

B 1 , B 2<br />

2+B 2 Z2 T )X T P T X<br />

)<br />

+A 2 (A † 2+B 2 Z2 T )X T P T PX(A †T<br />

2 +Z 2 B2 T )A T 2<br />

(1871)<br />

Necessary conditions for a global minimum are ∇ B1 =0 and ∇ B2 =0. Terms<br />

not containing B 2 in (1871) will vanish from gradient ∇ B2 ; (D.2.3)<br />

(<br />

∇ B2 tr −X T PXZ 2 B2A T T 2 −A 2 B 2 Z2X T T P T X+A 2 A † 2X T P T PXZ 2 B2A T T 2<br />

+A 2 B 2 Z T 2X T P T PXA †T<br />

2 A T 2+A 2 B 2 Z T 2X T P T PXZ 2 B T 2A T 2<br />

= −2A T 2X T PXZ 2 + 2A T 2A 2 A † 2X T P T PXZ 2 +<br />

)<br />

2A T 2A 2 B 2 Z2X T T P T PXZ 2<br />

= A T 2<br />

(−X T + A 2 A † 2X T P T + A 2 B 2 Z2X T T P T PXZ 2<br />

(1872)<br />

= 0 ⇔<br />

R(B 1 )⊆ N(A 1 ) and R(B 2 )⊆ N(A 2 )<br />

(or Z 2 = 0) because A T = A T AA † . Symmetry requirement (1868) is implicit.<br />

Were instead P T = ∆ (A †T<br />

2 + Z 2 B2 T )A T 2 and the gradient with respect to B 1<br />

observed, then similar results are obtained. The projector is unique.<br />

Perpendicularity (1870) establishes uniqueness [92,4.9] of projection P 1 XP 2<br />

on a matrix subspace. The minimum-distance projector is the orthogonal<br />

projector, and vice versa.<br />

<br />

)

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

Saved successfully!

Ooh no, something went wrong!