多层次灰色综合评价法及ahp验证

多层次灰色综合评价法及ahp验证

% grey_correlation_appraisal_ahp.m

clear all

clc

%指标数

a1_0=[2421 7409 2732 12188];

a2_0=[1293 4372 1350 4018];

a3_0=[300 0 100 100];

a4_0=[200 190 240 240];

a5_0=[2000 1150 2000 7791];

a6_0=[22 1148 35 931];

a7_0=[0.035 0.13 0.045 -0.088];

a8_0=[4 0 0 3];

a9_0=[50 165 100 220];

a10_0=[1 0 2 0];

%待判数据矩阵

A=[a1_0',a2_0',a3_0',a4_0',a5_0',a6_0',a7_0',a8_0',a9_0',a10_0']'; p=0.6;

for i=1:10

B(i,:)=(A(i,:)-min(A(i,:)))/(max(A(i,:))-min(A(i,:)));

end

%最佳值取每列的最大值（指标的最大值）

for i=1:10

V0(i)=max(B(i,:));

end

for i=1:10

for j=1:4

C(i,j)=abs(B(i,j)-V0(i));

end

end

r_min=min(min(C));

r_max=max(max(C));

% 计算相关系数E

i=1;

for i=1:10

for j=1:4

E(i,j)=(r_min+p*r_max)/(C(i,j)+p*r_max);

end

end

E;

% A的权重向量

Wa =[0.1062 0.2605 0.6333];

% B1的权重向量

Wb1= [0.2198 0.4265 0.0769 0.1648 0.1119];

% B2的权重向量

Wb2=[0.1667 0.8333];

% B3的权重向量

Wb3=[0.2519 0.5889 0.1593];

% B1的指标关联度

Rb1=Wb1*E(1:5,:);

% B2的指标关联度

Rb2=Wb2*E(6:7,:);

% B3的指标关联度

Rb3=Wb3*E(8:10,:);

% A的指标关联度

RA=Wa*[Rb1;Rb2;Rb3];

fprintf('利用多层次灰色综合评价计算结果为：\n'); fprintf('A的指标关联度为[%f %f %f %f]\n',RA);

% 利用层次分析法验证

% W为由层次分析法得到的各指标的权重系数

W =[0.1062*Wb1 0.2605*Wb2 0.6333*Wb3]; RAHP=B'*W';

fprintf('利用层次分析法计算结果为：\n');

fprintf('评价结果大小为[%f %f %f %f]\n',RAHP); % 将结果显示出来

subplot(2,2,1);

plot(RA);

subplot(2,2,3);

bar(RA); %柱状图

subplot(2,2,2);

plot(RAHP);

subplot(2,2,4);

bar(RAHP);