灰色关联度分析

clc, clear
a=[0.83  0.90   0.99      0.92  0.87      0.95
326  295     340  287   310   303
21   38      25   19     27   10
3.2 2.4      2.2   2.0   0.9   1.7
0.20    0.25        0.12   0.33     0.20      0.09
0.15    0.20        0.14   0.09 0.15   0.17
250 180     300    200  150   175
0.23    0.15        0.27   0.30 0.18   0.26
0.87    0.95        0.99       0.89 0.82      0.94];
for i=[1 5:9]    %效益型指标标准化
    a(i,:)=(a(i,:)-min(a(i,:)))/(max(a(i,:))-min(a(i,:)));
end
for i=2:4  %成本型指标标准化
   a(i,:)=(max(a(i,:))-a(i,:))/(max(a(i,:))-min(a(i,:))); 
end
[m,n]=size(a);
cankao=max(a')'  %求参考序列的取值
t=repmat(cankao,[1,n])-a;  %求参考序列与每一个序列的差
mmin=min(min(t));   %计算最小差
mmax=max(max(t));  %计算最大差
rho=0.5; %分辨系数
xishu=(mmin+rho*mmax)./(t+rho*mmax)  %计算灰色关联系数
guanliandu=mean(xishu)   %取等权重，计算关联度
[gsort,ind]=sort(guanliandu,'descend')  %对关联度按照从大到小排序

cankao =

 1
 1
 1
 1
 1
 1
 1
 1
 1

xishu =

0.3333    0.4706    1.0000    0.5333    0.4000    0.6667
0.4046    0.7681    0.3333    1.0000    0.5354    0.6235
0.5600    0.3333    0.4828    0.6087    0.4516    1.0000
0.3333    0.4340    0.4694    0.5111    1.0000    0.5897
0.4800    0.6000    0.3636    1.0000    0.4800    0.3333
0.5238    1.0000    0.4783    0.3333    0.5238    0.6471
0.6000    0.3846    1.0000    0.4286    0.3333    0.3750
0.5172    0.3333    0.7143    1.0000    0.3846    0.6522
0.4146    0.6800    1.0000    0.4595    0.3333    0.6296

guanliandu =

0.4630    0.5560    0.6491    0.6527    0.4936    0.6130

gsort =

0.6527    0.6491    0.6130    0.5560    0.4936    0.4630

ind =

 4     3     6     2     5     1