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

652 APPENDIX E. PROJECTION Proof. We take as given Theorem E.2.0.0.1 establishing idempotence. We have left only to show trP =‖P ‖ 2 F ⇒ P T = P , established in [344,7.1]. E.3.3 Summary, symmetric idempotent In summary, orthogonal projector P is a linear operator defined [173,A.3.1] by idempotence and symmetry, and characterized by positive semidefiniteness and nonexpansivity. The algebraic complement (E.2) to R(P ) becomes the orthogonal complement R(I − P ) ; id est, R(P ) ⊥ R(I − P ). E.3.4 Orthonormal decomposition When Z =0 in the general nonorthogonal projector A(A † + BZ T ) (1762), an orthogonal projector results (for any matrix A) characterized principally by idempotence and symmetry. Any real orthogonal projector may, in fact, be represented by an orthonormal decomposition such as (1788). [178,1, prob.42] To verify that assertion for the four fundamental subspaces (1786), we need only to express A by subcompact singular value decomposition (A.6.2): From pseudoinverse (1463) of A = UΣQ T ∈ R m×n AA † = UΣΣ † U T = UU T , I − AA † = I − UU T = U ⊥ U ⊥T , A † A = QΣ † ΣQ T = QQ T I − A † A = I − QQ T = Q ⊥ Q ⊥T (1797) where U ⊥ ∈ R m×m−rank A holds columnar an orthonormal basis for the orthogonal complement of R(U) , and likewise for Q ⊥ ∈ R n×n−rank A . Existence of an orthonormal decomposition is sufficient to establish idempotence and symmetry of an orthogonal projector (1788). E.3.5 Unifying trait of all projectors: direction Relation (1797) shows: orthogonal projectors simultaneously possess a biorthogonal decomposition (conferE.1.1) (for example, AA † for skinny-or-square A full-rank) and an orthonormal decomposition (UU T whence Px = UU T x).

E.3. SYMMETRIC IDEMPOTENT MATRICES 653 E.3.5.1 orthogonal projector, orthonormal decomposition Consider orthogonal expansion of x∈ R(U) : n∑ x = UU T x = u i u T i x (1798) a sum of one-dimensional orthogonal projections (E.6.3), where i=1 U ∆ = [u 1 · · · u n ] and U T U = I (1799) and where the subspace projector has two expressions, (1797) AA † ∆ = UU T (1800) where A ∈ R m×n has rank n . The direction of projection of x on u j for some j ∈{1... n} , for example, is orthogonal to u j but parallel to a vector in the span of all the remaining vectors constituting the columns of U ; u T j(u j u T j x − x) = 0 u j u T j x − x = u j u T j x − UU T x ∈ R({u i |i=1... n, i≠j}) (1801) E.3.5.2 orthogonal projector, biorthogonal decomposition We get a similar result for the biorthogonal expansion of x∈ R(A). Define A ∆ = [a 1 a 2 · · · a n ] ∈ R m×n (1802) and the rows of the pseudoinverse E.8 ⎡ ⎤ a ∗T 1 A † = ∆ a ∗T ⎢ ⎣ 2 ⎥ . ⎦ ∈ Rn×m (1803) a ∗T n under biorthogonality condition A † A=I . In biorthogonal expansion (2.13.8) n∑ x = AA † x = a i a ∗T i x (1804) E.8 Notation * in this context connotes extreme direction of a dual cone; e.g., (361) or Example E.5.0.0.3. i=1

E.3. SYMMETRIC IDEMPOTENT MATRICES 653<br />

E.3.5.1<br />

orthogonal projector, orthonormal decomposition<br />

Consider orthogonal expansion of x∈ R(U) :<br />

n∑<br />

x = UU T x = u i u T i x (1798)<br />

a sum of one-dimensional orthogonal projections (E.6.3), where<br />

i=1<br />

U ∆ = [u 1 · · · u n ] and U T U = I (1799)<br />

and where the subspace projector has two expressions, (1797)<br />

AA † ∆ = UU T (1800)<br />

where A ∈ R m×n has rank n . The direction of projection of x on u j for<br />

some j ∈{1... n} , for example, is orthogonal to u j but parallel to a vector<br />

in the span of all the remaining vectors constituting the columns of U ;<br />

u T j(u j u T j x − x) = 0<br />

u j u T j x − x = u j u T j x − UU T x ∈ R({u i |i=1... n, i≠j})<br />

(1801)<br />

E.3.5.2<br />

orthogonal projector, biorthogonal decomposition<br />

We get a similar result for the biorthogonal expansion of x∈ R(A). Define<br />

A ∆ = [a 1 a 2 · · · a n ] ∈ R m×n (1802)<br />

and the rows of the pseudoinverse E.8<br />

⎡ ⎤<br />

a ∗T<br />

1<br />

A † =<br />

∆ a ∗T<br />

⎢<br />

⎣<br />

2 ⎥<br />

. ⎦ ∈ Rn×m (1803)<br />

a ∗T<br />

n<br />

under biorthogonality condition A † A=I . In biorthogonal expansion<br />

(2.13.8)<br />

n∑<br />

x = AA † x = a i a ∗T<br />

i x (1804)<br />

E.8 Notation * in this context connotes extreme direction of a dual cone; e.g., (361) or<br />

Example E.5.0.0.3.<br />

i=1

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

Saved successfully!

Ooh no, something went wrong!