matlab 主成分分析

主成分分析

%主成分分析法
clc,clear;
gj=load('data.txt');
gj=zscore(gj);%数据标准化
r=corrcoef(gj);%计算相关系数矩阵
[x,y,z]=pcacov(r)%y为r的特征值，z为各个成分的贡献率
f=sign(sum(x));%构造元素为+-1的行向量
x=x.*f;%修改特征向量的正负号
num=3;%num为选取的主成分个数
df=gj*x(:,[1:num]);%计算各主成分的得分
tf=df*z(1:num)/100;%计算综合得分
[stf,ind]=sort(tf,'descend');%把得分按照从高到低的次序排列
stf=stf',ind=ind'