熵权法 The entropy weight method(EWM)

 1 % 初始数据矩阵R
 2 R=xlsread("C:\Users\PC\Desktop\数据.xlsx");
 3 % 输入矩阵的大小,rows为对象个数，cols为指标个数
 4 [rows,cols]=size(R);   
 5 k=1/log(rows);         % 求k
 6 Rmin = min(R);
 7 Rmax = max(R);
 8 A = max(R) - min(R);
 9 y = R - repmat(Rmin,2033,1);
10 for j = 1 : size(y,2)
11      y(:,j) = y(:,j)/A(j);
12 end
13 % 求Y（i,j）
14 S = sum(y,1);
15 Y = zeros(rows,cols); 
16 for i = 1 : size(Y,2)
17     Y(:,i) = y(:,i)/S(i);
18 end
19 % 初始化lnYij
20 lnYij=zeros(rows,cols);  
21 % 计算lnYij
22 for i=1:rows
23     for j=1:cols
24         if Y(i,j)==0
25             lnYij(i,j)=0;
26         else
27             lnYij(i,j)=log(Y(i,j));
28         end
29     end
30 end
31 % 计算熵值Hj
32 ej=-k*(sum(Y.*lnYij,1)); 
33 weights=(1-ej)/(cols-sum(ej));
34 F = zeros(rows,cols);
35 for k = 1 : size(R,2)
36      F(:,k) = weights(k)*y(:,k);
37 end
38 format long
39 % F即为对变量进行熵权法客观赋权后，计算获得的综合评分
40 F = sum(F,2) ;