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

664 APPENDIX E. PROJECTION E.6.3 Orthogonal projection on a vector The formula for orthogonal projection of vector x on the range of vector y (one-dimensional projection) is basic analytic geometry; [12,3.3] [287,3.2] [318,2.2] [332,1-8] 〈y,x〉 〈y,y〉 y = yT x y T y y = yyT y T y x = ∆ P 1 x (1848) where 〈y,x〉/〈y,y〉 is the coefficient of projection on R(y) . An equivalent description is: Vector P 1 x is the orthogonal projection of vector x on R(P 1 )= R(y). Rank-one matrix P 1 is a projection matrix because P 2 1 =P 1 . The direction of projection is orthogonal because P T 1 = P 1 . E.6.4 P 1 x − x ⊥ R(P 1 ) (1849) Orthogonal projection on a vectorized matrix From (1848), given instead X, Y ∈ R m×n , we have the one-dimensional orthogonal projection of matrix X in isomorphic R mn on the range of vectorized Y : (2.2) 〈Y , X〉 〈Y , Y 〉 Y (1850) where 〈Y , X〉/〈Y , Y 〉 is the coefficient of projection. For orthogonal projection, the term outside the vector inner-products 〈 〉 must be identical to the terms inside in three places. E.6.4.1 Orthogonal projection on dyad There is opportunity for insight when Y is a dyad yz T (B.1): Instead given X ∈ R m×n , y ∈ R m , and z ∈ R n 〈yz T , X〉 〈yz T , yz T 〉 yzT = yT Xz y T y z T z yzT (1851) is the one-dimensional orthogonal projection of X in isomorphic R mn on the range of vectorized yz T . To reveal the obscured symmetric projection

E.6. VECTORIZATION INTERPRETATION, 665 matrices P 1 and P 2 we rewrite (1851): y T Xz y T y z T z yzT = yyT y T y X zzT z T z ∆ = P 1 XP 2 (1852) So for projector dyads, projection (1852) is the orthogonal projection in R mn if and only if projectors P 1 and P 2 are symmetric; E.10 in other words, for orthogonal projection on the range of a vectorized dyad yz T , the term outside the vector inner-products 〈 〉 in (1851) must be identical to the terms inside in three places. When P 1 and P 2 are rank-one symmetric projectors as in (1852), (30) R(vec P 1 XP 2 ) = R(vecyz T ) in R mn (1853) and P 1 XP 2 − X ⊥ yz T in R mn (1854) When y=z then P 1 =P 2 =P T 2 and P 1 XP 1 = 〈P 1 , X〉 P 1 = 〈P 1 , X〉 〈P 1 , P 1 〉 P 1 (1855) meaning, P 1 XP 1 is equivalent to orthogonal projection of matrix X on the range of vectorized projector dyad P 1 . Yet this relationship between matrix product and vector inner-product does not hold for general symmetric projector matrices. E.10 For diagonalizable X ∈ R m×m (A.5), its orthogonal projection in isomorphic R m2 on the range of vectorized yz T ∈ R m×m becomes: P 1 XP 2 = m∑ λ i P 1 s i wi T P 2 i=1 When R(P 1 ) = R(w j ) and R(P 2 ) = R(s j ), the j th dyad term from the diagonalization is isolated but only, in general, to within a scale factor because neither set of left or right eigenvectors is necessarily orthonormal unless X is normal [344,3.2]. Yet when R(P 2 )= R(s k ) , k≠j ∈{1... m}, then P 1 XP 2 =0.

664 APPENDIX E. PROJECTION<br />

E.6.3<br />

Orthogonal projection on a vector<br />

The formula for orthogonal projection of vector x on the range of vector y<br />

(one-dimensional projection) is basic analytic geometry; [12,3.3] [287,3.2]<br />

[318,2.2] [332,1-8]<br />

〈y,x〉<br />

〈y,y〉 y = yT x<br />

y T y y = yyT<br />

y T y x = ∆ P 1 x (1848)<br />

where 〈y,x〉/〈y,y〉 is the coefficient of projection on R(y) . An equivalent<br />

description is: Vector P 1 x is the orthogonal projection of vector x on<br />

R(P 1 )= R(y). Rank-one matrix P 1 is a projection matrix because P 2 1 =P 1 .<br />

The direction of projection is orthogonal<br />

because P T 1 = P 1 .<br />

E.6.4<br />

P 1 x − x ⊥ R(P 1 ) (1849)<br />

Orthogonal projection on a vectorized matrix<br />

From (1848), given instead X, Y ∈ R m×n , we have the one-dimensional<br />

orthogonal projection of matrix X in isomorphic R mn on the range of<br />

vectorized Y : (2.2)<br />

〈Y , X〉<br />

〈Y , Y 〉 Y (1850)<br />

where 〈Y , X〉/〈Y , Y 〉 is the coefficient of projection.<br />

For orthogonal projection, the term outside the vector inner-products 〈 〉<br />

must be identical to the terms inside in three places.<br />

E.6.4.1<br />

Orthogonal projection on dyad<br />

There is opportunity for insight when Y is a dyad yz T (B.1): Instead given<br />

X ∈ R m×n , y ∈ R m , and z ∈ R n<br />

〈yz T , X〉<br />

〈yz T , yz T 〉 yzT = yT Xz<br />

y T y z T z yzT (1851)<br />

is the one-dimensional orthogonal projection of X in isomorphic R mn on<br />

the range of vectorized yz T . To reveal the obscured symmetric projection

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

Saved successfully!

Ooh no, something went wrong!