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