06.09.2021 Views

Linear Algebra, Theory And Applications, 2012a

Linear Algebra, Theory And Applications, 2012a

Linear Algebra, Theory And Applications, 2012a

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.

13.6. POLAR DECOMPOSITIONS 323<br />

Proof: (F ∗ F ) ∗ = F ∗ F and so by Theorem 13.3.3, there is an orthonormal basis of<br />

eigenvectors, {v 1 , ··· , v n } such that<br />

It is also clear that λ i ≥ 0 because<br />

Let<br />

F ∗ F v i = λ i v i ,F ∗ F =<br />

n∑<br />

λ i v i ⊗ v i .<br />

i=1<br />

λ i (v i , v i )=(F ∗ F v i , v i )=(F v i ,Fv i ) ≥ 0.<br />

U ≡<br />

n∑<br />

i=1<br />

λ 1/2<br />

i v i ⊗ v i .<br />

{<br />

Then U 2 = F ∗ F, U = U ∗ , and the eigenvalues of U,<br />

λ 1/2<br />

i<br />

} n<br />

i=1<br />

are all non negative.<br />

Let {Ux 1 , ··· ,Ux r } be an orthonormal basis for U (X) . By the Gram Schmidt procedure<br />

there exists an extension to an orthonormal basis for X,<br />

{Ux 1 , ··· ,Ux r , y r+1 , ··· , y n } .<br />

Next note that {F x 1 , ··· ,Fx r } is also an orthonormal set of vectors in Y because<br />

(F x k ,Fx j )=(F ∗ F x k , x j )= ( U 2 x k , x j<br />

)<br />

=(Uxk ,Ux j )=δ jk .<br />

By the Gram Schmidt procedure, there exists an extension of {F x 1 , ··· ,Fx r } to an orthonormal<br />

basis for Y,<br />

{F x 1 , ··· ,Fx r , z r+1 , ··· , z m } .<br />

Since m ≥ n, there are at least as many z k as there are y k .Nowforx ∈ X, since<br />

{Ux 1 , ··· ,Ux r , y r+1 , ··· , y n }<br />

is an orthonormal basis for X, there exist unique scalars<br />

such that<br />

Define<br />

Thus<br />

x =<br />

Rx ≡<br />

|Rx| 2 =<br />

c 1 , ··· ,c r ,d r+1 , ··· ,d n<br />

r∑<br />

c k Ux k +<br />

k=1<br />

r∑<br />

c k F x k +<br />

k=1<br />

r∑<br />

|c k | 2 +<br />

k=1<br />

Therefore, by Lemma 13.6.1 R ∗ R = I.<br />

Then also there exist scalars b k such that<br />

Ux =<br />

n∑<br />

k=r+1<br />

n∑<br />

k=r+1<br />

n∑<br />

k=r+1<br />

d k y k<br />

d k z k (13.18)<br />

|d k | 2 = |x| 2 .<br />

r∑<br />

b k Ux k (13.19)<br />

k=1

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

Saved successfully!

Ooh no, something went wrong!