matlab 基于Topsis的熵权法2

基于Topsis的熵权法2

%%基于Topsis的熵权法代码
%% 熵权法部分
clc,clear;
xij=xlsread('accessory1.xls','A','C2:IH147');
[rows,cols]=size(xij);   % 输入矩阵的大小,rows为对象个数，cols为指标个数
k=1/log(rows);        % 求k
zij=xij./repmat(sum(xij.*xij).^0.5,rows,1);
pij=zij./repmat(sum(zij),rows,1);

lnpij=zeros(rows,cols);  % 初始化lnpij
% 计算lnpij
for i=1:rows
    for j=1:cols
       if pij(i,j)==0
          lnpij(i,j)=0;
       else
          lnpij(i,j)=log(pij(i,j));
       end
   end
end
ej=-k*(sum(pij.*lnpij,1)); % 计算熵值ej
dj=1-ej;
wj=dj/sum(dj);
save wj;
%% Topsis部分
m=rows;
n=cols;
%step1&step2
zjplus=max(zij);
zjminus=min(zij);

Diplus=sum([(zij - repmat(zjplus,m,1)) .^ 2 ] .* repmat(wj,m,1),2) .^ 0.5;
Diminus=sum([(zij - repmat(zjminus,m,1)) .^ 2 ] .* repmat(wj,m,1),2) .^ 0.5;
%未归一化得分
Si=Diminus./(Diplus+Diminus);
%归一化得分
Si_uniformized=Si/sum(Si);

其中，标准化处理有多种方式，但差别不大。