function?[eigvector_PCA?eigvalue_PCA?meanData?new_X]??=?PCA（BpcaRatio）;
%輸入圖像按列存放；輸出按列存放
[NP]?=?size（B）;
%********subtract?mean
meanData=mean（B‘）;
B=B-（ones（P1）*meanData）‘;

%********Find?eigenvectors?vi?of?B‘B?（PxP）
[VD]?=?eig?（1/（P-1）*（B‘*B））;???%scale?factor?gives?eigvals?correct
????????????????????????????????%magnitude?for?large?cov?mat?
????????????????????????????????%（assuming?sample?cov）
????????????????????????????????%（assuming?sample?cov）
%********Sort?eigenvecto