matlab 基于Topsis的熵权法

基于Topsis的熵权法

clc,clear;
testData=xlsread('accessory1.xls','C2:IH403');
R=testData;
[rows,cols]=size(R);   % 输入矩阵的大小,rows为对象个数，cols为指标个数
k=1/log(rows);        % 求k
z=R./repmat(sum(R.*R).^0.5,rows,1);
a=sum(z);
p=z./repmat(a,rows,1);

lnpij=zeros(rows,cols);  % 初始化lnfij
% 计算lnfij
for i=1:rows
    for j=1:cols
       if p(i,j)==0
          lnpij(i,j)=0;
       else
          lnpij(i,j)=log(p(i,j));
       end
   end
end

Hj=-k*(sum(p.*lnpij,1)); % 计算熵值Hj
weights=(1-Hj)/(cols-sum(Hj));

cumsum(weights);
save weights;