基于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;