综合评价分析法：熵权法matlab实现+层次分析法matlab实现

综合评价分析法：熵权法matlab实现+层次分析法matlab实现

层次分析法matlab实现

clc;clear%层次分析法matlab实现

A=[1    	 1/9	 1/2	 1/4	 1/6
9    	1    	6    	5    	3    
2    	 1/6	1    	 1/2	 1/3
4    	 1/5	2    	1    	 1/5
6    	 1/3	3    	5    	1    
]%判断矩阵，9级度量法
[m,n]=size(A);
[V,D]=eig(A);%求得特征向量和特征值
            %求出最大特征值和它所对应的特征向量
tempNum1=D(1,1);
pos=1;
for h=1:n
    if D(h,h)>tempNum1
        tempNum1=D(h,h);
        pos=h;
    end
end    
w=abs(V(:,pos));
w=w/sum(w);
t=D(pos,pos);
disp('权重w=');disp(w);disp('最大特征根t=');disp(t);
         %以下是一致性检验
CI=(t-n)/(n-1);RI=[0 0 0.58 0.89 1.12 1.26 1.36 1.41 1.46 1.49 1.52 1.54 1.56 1.58 1.59 1.60 1.61 1.615 1.62 1.63];
CR=CI/RI(n);
if CR<0.10
    disp('此矩阵的一致性可以接受!');
    disp('CI=');disp(CI);
    disp('CR=');disp(CR);
else disp('此矩阵的一致性验证失败，请重新进行评分!');
end

熵权法matlab实现

clc,clear%熵权法matlab程序
R=[];%输入矩阵R,列代表评价指标评分，行代表同类别竞争对手
[m,n] = size(R);
R1 = zeros(m,n);
for i = 1:n
    R1(:,i) = (R(:,i)-min(R(:,i)))/(max(R(:,i))-min(R(:,i)));%矩阵标准化
end
ex=sum(R1);

R2=zeros(m,n);
for i=1:m
    for j=1:n
    R2(i,j)=R1(i,j)/ex(j);%%求得pij
    end
end

R3=R2.*log(R2);%求得pij*lnpij
R4=R3;
R4(find(isnan(R4)==1)) = 0;
ex=sum(R4);%求得sum (pij*lnpij)
ex1=-1/log(m)*ex;%每列信息熵
ex2=(1-ex1)/(n-sum(ex1));%每列权重

for i=1:m
R5(i)=sum(R(i,:).*ex2);%计算每个竞争对手得分
end
R5