pca基础.m

matlab 使用源码 主成分分析相关源吗 好东西大家一起用

function main()
%*************主成份分析************
%读入文件数据
%X=load('data.txt');
X=input('请输入数据 ');

%==========方法1:求标准化后的协差矩阵,再求特征根和特征向量=================
%标准化处理
X=data;
[p,n]=size(X);
for j=1:n
    mju(j)=mean(X(:,j));
    sigma(j)=sqrt(cov(X(:,j)));
end
for i=1:p
    for j=1:n
        Y(i,j)=(X(i,j)-mju(j))/sigma(j); 
    end
end
sigmaY=cov(Y);

%求X标准化的协差矩阵的特征根和特征向量
[T,lambda]=eig(sigmaY);

disp('特征根(由小到大):');
disp(lambda);
disp('特征向量:');
%disp(T);

for i=1:n
    for j=1:n
       %addT(i,n-j+1)=-T(i,j);
        addT(i,j)=T(i,j);
    end
end
disp(addT);
%方差贡献率;累计方差贡献率
Xsum=sum(sum(lambda,2),1);
for i=1:n
    fai(i)=lambda(i,i)/Xsum;
    %addfai(n-i+1)=fai(i);
    addfai(i)=fai(i);
end
for i=1:n
    psai(i)= sum(sum(lambda(1:i,1:i),2),1)/Xsum;
end
disp('方差贡献率:');
%disp(fai);
disp(addfai);
%disp('累计方差贡献率:');
%disp(psai);
disp('综合评价:');
Q=Y*addT;
R=addfai*Q';
disp(R);

xlswrite('C:\11111111.xls',R');

        

%+============方法2:求X的相关系数矩阵,再求特征根和特征向量================
%X的标准化的协方差矩阵就是X的相关系数矩阵
R=corrcoef(X);
%求X相关系数矩阵的特征根和特征向量
[TR,lambdaR]=eig(R);
disp('特征根(由小到大):');
disp(lambdaR);
disp('特征向量:');
disp(TR);