评价类模型（一）层次分析法

1、问法:确定评价指标，形成评价体系------>打分法解决评价问题------>需要填好如下表格

 同一颜色的单元格的和为1，它们表示针对某一因素所占的权重或得分。

2、解决评价类问题要解决的三个问题：

① 我们评价的目标是什么？

② 我们为了达到这个目标有哪几种可选的方案？

③ 评价的准则或者说指标是什么？（我们根据什么东西来评价好坏）

3、重要程度表

4、判断矩阵、一致矩阵

判断矩阵： 且满足  的方阵。（也称为正互反矩阵）

一致矩阵：满足  的判断矩阵（正互反矩阵）。

                   特点：各行（列）之间成倍数关系

 5、一致性检验

---

 层次分析法的步骤：

1、分析系统中各个因素之间的关系，构造如下的递阶层次结构。

 2、将每一层次的各元素与上一层次中的某一准则两两比较，构造判断矩阵。

 判断矩阵O-C

 3、对判断矩阵进行一致性检验，并计算被比较元素对于该准则的相对权重。

（1）若判断矩阵是一致矩阵：（不需要进行一致性检验）

进行归一化处理（按列算）

（2）若判断矩阵不是一致矩阵，但是通过了一致性检验：

① 算术平均法求权重

Step 1：将判断矩阵按照列进行归一化处理（每个元素除以其所在列的和）

Step 2：将归一化的矩阵按行求和得到向量

Step 3：将上述向量中的每个元素除以n

 （实际应用中要将n替换成已知数）

② 几何平均法求权重

Step 1：将判断矩阵按行相乘得到一个列向量

Step 2：将上述向量的每个分量开n次方

Step 3：对列向量进行归一化

 （实际应用中要将n替换成已知数）

③ 特征值法求权重

Step 1：求出判断矩阵的最大特征值及其对应的特征向量

Step 2：对特征向量进行归一化

总结：在做题过程中我们发现，在以往的类似题目中用层次分析法时，大都只选择了上述方法其一来求权重，但是上述三种方法求出的权重结果或多或少都有一定的偏差。因此，为了保证结果的稳健性，在本次解题中我们采用三种方法分别求出权重再取其平均值，再根据权重矩阵计算个方案的得分，最后进行排序和综合分析。这样可以尽量减小误差，使得到的结论更加全面有效。

（3）若判断矩阵不是一致矩阵，且没有通过一致性检验：

对判断矩阵进行修正（往一致矩阵上调整）

4、根据权重计算得分并进行排序（用excel，注意用F4锁定单元格）

---

注意

平均随机一致性指标RI的表格中n最多是15，所以评价的决策层不能太多，太多的话n会很大，判断矩阵和一致矩阵差异可能会很大。

如果决策层中指标的数据是已知的------>TOPSIS

---

简单例子

任务：填完下表

先分别填好以下几张表：

 

 

若一致性检验不通过，则修改判断矩阵：

 如果判断矩阵是一致矩阵，则直接归一化求权重:

若一致性检验通过，则三种方法求权重：

 

 

 

得到结果如下：

 继续上述步骤：

 

 总结如下表：

 最后分别计算得分并排序。

---

 代码

clear;clc
J=input('J=');%输入判断矩阵J
%先检验输入的矩阵J是否正确，即J是否是判断矩阵以及是否适合使用层次分析法
flag=0;%判断矩阵J正确无误，符合要求，计flag为0
[m,n]=size(J);%矩阵J有m行n列
if m~=n||m<=1%如果矩阵J不是方阵或者维数小于等于1，则计flag为1
    flag=1;
end
if flag==0
    [r,r]=size(J);%flag==0说明判断矩阵J是一个方阵，所以其行数等于列数记为r
    if r>15%平均随机一致性指标RI的表格中n最多是15,所以若J的维数超过了15，则计flag为2
        flag=2;
    end
end
if flag==0
     if sum(sum(J<=0))>0%如果矩阵J中有元素小于0，则计flag为3
        flag=3;
    end
end
if flag==0
    if sum(sum(J'.*J~=ones(r)))>0%如果矩阵J中aij*aji~=1，则计flag为4
        flag=4;
    end
end
%下面进行一致性检验
[L,Z]=eig(J);%求判断矩阵J的特征值及其对应的特征向量
MaxZ=max(Z(:));%求出判断矩阵J的最大特征值
CI=(MaxZ-r)/(r-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];%平均随机一致性指标
CR=CI/RI(r)%一致性比例
if CR>=0.10%如果一致性比例大于等于0.10，则判断矩阵J没通过一致性检验，需要进行修改，此时计flag为5
    flag=5;
end
%若（修改后）判断矩阵J通过了一致性检验，则下面采用三种方法分别求权重
if flag==0
    %方法一：使用算术平均法求权重
    J1=sum(J);%判断矩阵J按列求和
    J2=repmat(J1,r,1);%将J1复制r*1块
    J3=J./J2;%将判断矩阵J按照列归一化处理
    Q1=(sum(J3,2)./r)%使用算术平均法求出的权重
end
if flag==0
    %方法二：使用几何平均法求权重
    J4=prod(J,2);%将判断矩阵J中的元素按行相乘
    J5=J4.^(1/r);%对J4的所有分量开r次方
    Q2=J5./sum(J5)%使用几何平均法求出的权重
end
if flag==0
    %方法三：使用特征值法求权重
    [h,l]=find(Z==MaxZ,1);%找出特征值矩阵中第一个最大特征值的行列
    Q3=L(:,l)./sum(L(:,l))%使用特征值法求出的权重
end
if flag==0
    Q=(Q1+Q2+Q3)/3%最终的权重结果
end