数学建模——熵权法

clc;clear;
x=xlsread('熵权法.xlsx');
ind=ones(size(x,2),1);%正向写1，负向写2
[n,m]=size(x);%n个样本，m个指标，数据的归一化处理
for i=1:m
    if ind(i)==1%正向指标归一化
        X(:,i)=guiyi(x(:,i),1.0,0.002,1);
    else%负向指标归一化
        X(:,i)=guiyi(x(:,i),2.0,0.002,1);
    end
end
%计算第j个指标下，第i个样本占该指标的比重p(i,j)
for i=1:n
    for j=1:m
        p(i,j)=X(i,j)/sum(X(:,j));
    end
end
%计算第j个指标的熵值e(j)
k=1/log(n);
for j=1:m
    e(j)=-k*sum(p(:,j).*log(p(:,j)));
end
g=ones(1,m)-e;%计算信息熵冗余度
w=g./sum(g);%求权值w
s=X*w';%求综合得分

function y=guiyi(x,type,ymin,ymax)
%实现正向或负向指标归一化，返回归一化后的数据矩阵
%x为原始数据矩阵，一行代表一个样本，每列对应一个指标
%type设定正向指标1，负向指标2
%ymin，ymax为归一化的区间端点
[n,m]=size(x);
y=zeros(n,m);
xmin=min(x);
xmax=max(x);
switch type
    case 1
        for j=1:m
            y(:,j)=(ymax-ymin)*(x(:,j)-xmin(j))/(xmax(j)-xmin(j))+ymin;
        end
    case 2
         for j=1:m
            y(:,j)=(ymax-ymin)*(xmax(j)-x(:,j))/(xmax(j)-xmin(j))+ymin;
         end
end