weixie.m

来自「matlab 主成分分析法程序」· M 代码 · 共 25 行

clc
data;          % 采集数据
Y=zscore(x)     %Z_score法标准化数据
r=corrcoef(Y)   %标准化后数据的相关矩阵
[D,V]=eig(r)    %相关矩阵的特征向量和特征值
a=fliplr(D(:,6:8))    %将特征向量的后3列翻转（后3列对应最大的四个特征值的特征向量）
Z=Y*a           %主成分决策矩阵
u=[0.4802,0.3856,0.1341]       %主成分指标权重
U=repmat(u,6,1).*Z             %主成分加权决策矩阵
Aplus=max(U)    %正理想方案
Amins=min(U)    %负理想方案

%以下求综合主成分指标qi
%求(8)式中分子矩阵
A=U-repmat(Amins,6,1);
B=A.*A;
C=sum(B');
fenzi=C'
E=U-repmat(Aplus,6,1);
%求(8)式中分母矩阵对应的前一项
F=E.*E;
G=sum(F');
fenmu=G'

q=fenzi./(fenzi+fenmu)          %求综合主成分指标qi