灰色关联分析发

森之张卫东

1. <div>%灰色关联分析
   
2. clc
   
3. clear
   
4. a=[0.83  0.90 0.99   0.92 0.87   0.95
   
5. 326  295  340  287 310   303
   
6. 21  38      25   19  27   10
   
7. 3.2 2.4      2.2   2.0  0.9   1.7
   
8. 0.20 0.25     0.12   0.33  0.20   0.09
   
9. 0.15 0.20     0.14   0.09 0.15   0.17
   
10. 250 180     300    200 150   175
    
11. 0.23 0.15     0.27   0.30 0.18   0.26
    
12. 0.87 0.95     0.99    0.89 0.82   0.94];
    
13. for i=[1 5:9]    %效益型指标标准化
    
14.     a(i,:)=(a(i,:)-min(a(i,:)))/(max(a(i,:))-min(a(i,:)));
    
15. end
    
16. for i=2:4  %成本型指标标准化
    
17.    a(i,:)=(max(a(i,:))-a(i,:))/(max(a(i,:))-min(a(i,:)));
    
18. end
    
19. [m,n]=size(a);
    
20. disp('求参考序列的取值')
    
21. cankao=max(a')'  %求参考序列的取值
    
22. t=repmat(cankao,[1,n])-a;  %求参考序列与每一个序列的差
    
23. mmin=min(min(t));   %计算最小差
    
24. mmax=max(max(t));  %计算最大差
    
25. %----------------------------------------------------
    
26. %分辨系数
    
27. rho=0.5;
    
28. %----------------------------------------------------
    
29. disp('计算灰色关联系数')
    
30. xishu=(mmin+rho*mmax)./(t+rho*mmax)  %计算灰色关联系数
    
31. %---------------------------------------------------------------
    
32. %若用层次分析法得出的权重，只需在 xishu 矩阵上修改即可，
    
33. %用 xishu 里面的数据再点乘以权重矩阵即可。
    
34. %点乘即为对应位置相乘！
    
35. %---------------------------------------------------------------
    
36. disp('取等权重，计算关联度')
    
37. guanliandu=mean(xishu)   %取等权重，计算关联度</div><div><div>disp('对关联度按照从大到小排序')
    
38. [gsort,ind]=sort(guanliandu,'descend')  %对关联度按照从大到小排序</div><div><div><div><div><div><div><div><div><div><div>
    
39. </div></div></div></div></div></div></div></div></div></div></div>

复制代码