主成分分析的matlab程序

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