function??[prindeta]=pca（imSpace）
%主成分分析（PCA）程序。輸入為系數
psi?=?mean（double（imSpace‘））‘;?%求各指標均值
[mn]=size（imSpace）;
for?i?=?1?:?m
????zeroMeanSpace（i:）?=?double（imSpace（i:））?-?psi（i）;%去均值
end
L?=?zeroMeanSpace*?zeroMeanSpace‘;??%計算離差陣
S=L/（n-1）;%計算協方差陣
%求自相關矩陣
r=zeros（mm）;
for?i=1:m
????for?j=1:m
????????r（ij）=L（ij）/sqrt（L（ii）*L（jj））;
????end
end
[eigVecs?eigVals]?=?eig（r）;???%?eig是求解特征值和特征向量???
diagonal=diag（eigVals）;
[diagonal?index]?=?so