4-主成分分析法

clc, clear;
pp = xlsread('E:\a-模組化\第六輪\結果\附件4.一次風機部分監測資料.xls', 2,'L21:P3411');
gj=pp;   %原始資料 
gj=zscore(gj); %資料标準化
r=corrcoef(gj);  %計算相關系數矩陣
%下面利用相關系數矩陣進行主成分分析，x的列為r的特征向量，即主成分的系數
[x,y,z]=pcacov(r) %y為r的特征值，z為各個主成分的貢獻率
f=repmat(sign(sum(x)),size(x,1),1); %構造與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'

mz = [z';cumsum(z')];
mz = mz';