层次分析法是一种主观赋值的评价方法。层次分析法将与决策有关的元素分解成目标、准则、方案分析。是一种灵活有效的决策方法。
AHP层次分析法建模大致:
1. 建立系统的递阶层次结构(找关系,分层次)。
2. 根据同一层元素对上一层评价的重要性,构建比较矩阵(判断矩阵)。
3. 由判断矩阵计算该评价的相对权重
4. 进行一致性检验。(比对检验表)
5. 填充权重矩阵,根据矩阵计算得分。
求权重方法:
1. 算数平均法求权重
2. 几何平均法求权重
3. 特征值法求权重
一致性检验求权重:
得到权重矩阵后可以计算最大特征根
再依次求CI、RI。RI需要比对表格
一致性比例CR: CR=CI/RI
判断:
CR<0.1时表明判断矩阵的一致性程度在允许的范围内,此时可以开展权向量计算;
若CR>=0.1,说明判断矩阵出现逻辑错误,需要重新构建。
评分:
matlab实现:
clear;clc;
A=[1 2 3 4 5;
1/2 1 2 3 4;
1/3 1/2 1 2 3;
1/4 1/3 1/2 1 2;
1/5 1/4 1/3 1/2 1];
%disp('请输入准则层判断矩阵A(n阶)');
%A=input('A=');
[n,n]=size(A);
Sum_A=sum(A);
SUM_A=repmat(Sum_A,n,1);
Stand_A=A./SUM_A;
disp('算术平均法求权重的结果为:')
disp(sum(Stand_A,2)./n);
Prduct_A=prod(A,2);%包含每一行元素乘积的列向量
Prduct_n_A=Prduct_A.^(1/n);
disp('几何平均法求权重的结果为:')
disp(Prduct_n_A./sum(Prduct_n_A));
[V,D]=eig(A);%求得特征向量和特征值
%V是包含A的特征向量矩阵,D是一个对角矩阵
%求出最大特征值和它所对应的特征向量
tempNum=D(1,1);
pos=1;
for h=1:n
if D(h,h)>tempNum
tempNum=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.52 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