Exercise 12.7 - PCA via successive deflation

Answers

The matrix:

𝐈 𝐯 𝐯 T

is unknown as the projection matrix that maps an arbitrary vector into the subspace that is orthogonal to 𝐯 , since:

𝐯 T ( 𝐈 𝐯 𝐯 T ) 𝐱 = 𝐯 T 𝐱 𝐯 T 𝐱 .

The only requirement is that 𝐯 is a unit vector: 𝐯 T 𝐯 = 1 .

For generalization to the projection matrix to the supplementary of a multi-dimensional space spanned by a set of orthogonal and unit bases { 𝐯 1 , , 𝐯 M } (otherwise using the Schmidt procedure first), we only have to use:

𝐈 m = 1 M 𝐯 m 𝐯 m T .

For question (a), we have:

1 N 𝐗 ~ 𝐗 ~ T = 1 N ( 𝐈 𝐯 1 𝐯 1 T ) 𝐗 𝐗 T ( 𝐈 𝐯 1 𝐯 1 T ) T = 1 N 𝐗 𝐗 T 1 N 𝐯 1 𝐯 1 T 𝐗 𝐗 T 𝐯 1 𝐯 1 T = 1 N 𝐗 𝐗 T 𝐯 1 𝐯 1 T λ 1 𝐯 1 𝐯 1 T = 1 N 𝐗 𝐗 T λ 1 𝐯 1 𝐯 1 T .

For question (b), let 𝐂 ~ = 1 N 𝐗 ~ 𝐗 ~ T then its principal eigenvector 𝐮 satisfies:

𝐂 ~ 𝐮 = λ 𝐮 .

To solve this equation, expand 𝐮 onto the eigenvectors of 𝐂 :

𝐮 = d = 1 D f d 𝐯 d ,

then:

𝐂 ~ 𝐮 = d = 1 D f d ( 𝐂 λ 1 𝐯 1 𝐯 1 T ) 𝐯 d = d = 2 D f d λ d 𝐯 d .

Then the eigenequation becomes:

f 1 λ = 0 , d = 2 , , D , f d λ d = f d λ ,

s.t.:

d = 1 D f d 2 = 1 .

To solve for this system, we have:

f 1 = 0 , λ = λ d , d 2 f d = 1 , f d = 1 , d d .

Hence the solution is d = 2 , this finishes the proof.

For question (c), the procedure is roughly:

1.
[ λ , 𝐮 ] = f ( 𝐂 ) ;
2.
Save λ and 𝐮 ;
3.
𝐂 = λ 𝐮 𝐮 T ;
4.
Repeat.
User profile picture
2021-03-24 13:42
Comments