数学建模--（1.1）层次分析法（AHP）
层次分析法
层次分析法（analytic hierarchy process），简称AHP。是建模比赛中比较基础的模型之一，其主要解决评价类的问题。如选择哪种方案最好，哪位员工表现最好等。
它是一个较为 主观 的评价方法，其在赋权得到权重向量的时候，主观因素占比很大。因而在建模比赛中，常常与客观方法得到的权重向量方法进行综合，而得出一个综合的权重向量，进行后续操作。
数学建模模型介绍：
目录
层次分析法
1.建模步骤
2.模型实现
3.模型拓展
4.模型总结
5.代码实现
1.建模步骤
层次分析法进行建模，大致分为以下四步：
- 分析系统中各因素之间的关系，建立系统的递阶层次结构。
- 对于同一层次的个元素关于上一层次中某一准则的重要性两两比较，构造两两比较矩阵（判断矩阵）。
- 由判断矩阵计算被比较元素对于该准则的相对权重，并进行一致性检验（检验通过权重才能用）。
- 填充权重矩阵，根据矩阵计算得分，得出结果。
接下来，我们根据例题讲解，并在相应出进行解释。
2.模型实现
例：小明同学想出去旅游。在查阅了网上的攻略后，他初步选择了苏杭、北戴河和桂林三地之一作为目标景点。请你确定评价指标、形成评价体系来为小明同学选择最佳的方案。
（1）.分析系统中各因素之间的关系，建立系统的递阶层次结构。
首先，需要明确以下三个问题：
- 我们评价的 目标 是什么？ 答：为小明同学选择最佳的旅游景点。
- 评价的 准则 或者说指标是什么？（我们根据什么东西来评价好坏）答：景色、花费、居住、饮食、交通。
- 我们为了达到这个目标有哪几种可选的 方案 ？ 答：三种，分别是去苏杭、去北戴河和去桂林。
根据以上问题，建立层次结构图，旅游地选择层次结构图如下：
（2）. 对于同一层次的个元素关于上一层次中某一准则的重要性两两比较，构造两两比较矩阵（判断矩阵）。
首先，我们先来看看层次分析法最终要得出的结果是什么样子的：
实际上的建模结果就是要填满权重矩阵，即这个表格：
- 其中，蓝色一列代表景色、花费、居住、饮食以及交通的权重，加和为1。（实际上就是准则层关于上一层目标层的重要性）
- 然后同一颜色每一横行，就是三种方案相对于准则层的重要性。如：橙色一行 代表的就是苏杭、北戴河以及桂林关于景色的权重，以此类推。
如何填满这个表格，就需要用判断矩阵得出，这也是构造判断矩阵的意义！
然后，我们看一下如何构造这个判断矩阵，要构造几个？
由上文可知得到这个判断矩阵实际上要分别得出准则层关于目标层的一组权重向量，方案层关于准则层的五组权重向量，实际上我们就需要构造出一个准则层关于目标层的判断矩阵以及五个方案层关于准则层的矩阵，一共六个判断矩阵。（这里采用分治的思想）最终在经过权重计算每组得出一组权重向量，填到相应的表格中。构造的6个判断矩阵如下：
接下来，我们看一下，每个位置应该怎么填。
注意：这个位置不是随便填的，因为影响因子占比很大，有可能第一天我们看重景色，把景色权重写的占比大一些，第二天我们看重饮食了，就把饮食占比大一些，因而常常考虑不周全，而使得不易定量化。所以需要两两比较得出判断矩阵，而两两比较得出重要性填到矩阵中。重要程度如下表：
根据以上这个表格，我们人为的进行填充，得到了下面这个判断矩阵：（实际情况下都是专家填的，但是比赛中大都是我们自己填的，最好有一些理论的依据支撑）
观察一下：上面这个判断矩阵有如下特点：
- aij 表示的意义是，与指标j相比， i的重要程度。
- 当i= j时，两个指标相同，因此同等重要记为1，这就解释了主对角线元素为1
- aj >0 且满足 aij × aji = 1 (我们称满足这一条件的矩阵为正互反矩阵)
(3). 由判断矩阵计算被比较元素对于该准则的相对权重，并进行一致性检验（检验通过权重才能用）。
思考一个问题，拿方案层关于景色的矩阵说明，假设我们填写的判断矩阵是这个样子：
假设：苏杭 = A、北戴河 = B、桂林 = C，
那么由矩阵可以看出，苏杭比北戴河景色好一点 A > B，苏杭和桂林景色一样好 A = C，北戴河比桂林景色好一点 B > C，出现了 矛盾！
这里就不得不提出一个概念叫做 一致矩阵，它在正互反矩阵性质的基础上没有以上的矛盾，可以说：一致矩阵是正互反矩阵的特例。
将上面的矩阵进行改良，得到一致矩阵：
它比正互反矩阵多出两个性质：
- aij = i的重要程度 / j的重要程度，ajk = j的重要程度 / k的重要程度，aik = i的重要程度 / k的重要程度 = aij × ajk。
- 矩阵各行（各列）之间成倍数关系。
我们进行构造矩阵大多是正互反矩阵，难免会出现矛盾，即不容易构造出一致性矩阵，但是我们可以向一致性矩阵靠拢，只要这个差距（CR）不超过一个范围（0.1）那么这个判断矩阵也是可以使用的。这个判断差距的过程叫做 一致性检验。
接下来正式讲解一下一致性检验。
证明过程（只需要了解，可以直接看步骤）：
一致性检验原理： 检验我们构造的判断矩阵和一致矩阵是否有太大的差别。
引理： A 为 n 阶方阵，且 A 的秩 r(A) = 1，则 A 有一个特征值为 A 的迹 tr(A)，其余特征值均为0。因为一致矩阵的各行成比例且不是零矩阵，所以一致矩阵的秩一定为1。
由引理可知： 一致矩阵有一个特征值为 n，其余特征值均为 0。
另外，我们很容易可以得到，特征值为 n 时，对应的特征向量刚好为：
一致性检验的步骤：
第一步：计算一致性指标 CI
第二步：查找对应的平均随机一致性指标 RI
注：在实际运用中，n很少超过10，如果指标的个数大于10，则可考虑建立二级指标体系，或使用我们以后要学习的模糊综合评价模型。此外，RI 可以直接查表使用即可。
第三步：计算一致性比例CR
如果 CR < 0.1, 则可认为判断矩阵的一致性可以接受；否则需要对判断矩阵进行修正。
当然本道例题上面的6个矩阵已经通过了一致性检验，然后接下来要根据判断矩阵计算权重。
(4). 填充权重矩阵，根据矩阵计算得分，得出结果。
计算权重的方法有三种：算数平均法求权重、几何平均法求权重以及特征值法求权重。
一般情况下：第三种特征值法求权重是最常用的，但是建议可以综合三种方法来求得一个综合的权重向量。
下面拿下面这个判断矩阵进行说明：
方法1：算术平均法求权重
第一步： 将判断矩阵按照列归一化（每个元素除以其所在列的和，如1/(1+0.5+0.2)=0.5882）
第二步： 将归一化的列相加（按行求和）
第三步： 将相加后得到的向量中的每个元素除以 n 即可得到权重向量
以表达式进行解释：
方法2：几何平均法求权重
第一步： 将A的元素按照行相乘得到一个新的列向量
第二步： 将新的向量的每个分量开n次方
第三步： 对该列向量进行归一化即可得到权重向量
以表达式进行解释：
求得权重结果如图：
方法3：特征值法求权重
上文在证明一致性检验的时候，提到过如下内容：
由引理可知： 一致矩阵有一个特征值为 n，其余特征值均为 0。
另外，我们很容易可以得到，特征值为 n 时，对应的特征向量刚好为：
那么我们我们直接可以将特征向量归一化即可求得特征向量。
求的结果为：
最大特征值为 ：n=3.0055,
一致性比例 CR = 0.0053 对应的 特征向量:[-0.8902,-0.4132,-0.1918]
归一化后得到权重向量 : [0.5954,0.2764,0.1283]
我们将三种方法求得的权重向量如下图所示：
但在实际建模中建议综合三种方法求得的权重得到一个综合的权重向量更具有说服力！
填充权重矩阵
这里只拿第三种的结果填充权重矩阵
计算得分，得出最终结果：
苏杭得分：指标权重×苏杭与其他两种方案中的权重，即前两列相乘：
0.5954×0.2636＋0.0819×0.4758＋0.4286×0.0538＋0.6337×0.0981＋0.1667×0.1087=0.299
同理：北戴河得分为0.245，桂林得分为0.455。
因此最佳的旅游景点是桂林。
3.模型拓展
- 评价的决策层不能太多，太多的话n会很大，判断矩阵和一致矩阵差异可能会很大。因为平均随机一致性指标 RI 的表格中 n 最多是15，因此应该根据实际情况选择是否应用此方法。
- 如果决策层中指标的数据是已知的，那么层次分析法不容易将这些已知数据应用在其中。如拿上面的例题举例：如果已知景色 、花费、居住、饮食以及交通在三个旅游景点的一些数据，那么如何将这些数据转化为构造判断矩阵的依据，只能为其提供一定的文字说明，而不容易将数据应用到其中。
- 在实际建模中，判断矩阵的数值都是人为填的，具有一定的主观性存在，这时应该搜寻相应的数据让人信服，不能空口无凭。
- 如果说只想拿到的决策因素的权重向量，那大可不必这么麻烦，在第一步递阶层次结构的时候，只需要目标层和准则层即可，构造判断矩阵也只需要构造出一个，并进行检验，检验通过了，差不多就拿到了权重向量。
4.模型总结
步骤总结：
- 分析系统中各因素之间的关系，建立系统的递阶层次结构。
- 对于同一层次的个元素关于上一层次中某一准则的重要性两两比较，构造两两比较矩阵（判断矩阵）。
- 由判断矩阵计算被比较元素对于该准则的相对权重，并进行一致性检验（检验通过权重才能用）。
- 填充权重矩阵，根据矩阵计算得分，得出结果。
5.代码实现
在下面的代码中，我们先计算了权重，然后再进行一致性检验，这是为了顺应计算过程，事实上在逻辑上是说不过去的。
因此在写论文过程中如果用到了层次分析法，一定要先对判断矩阵进行一致性检验。
而且要说明的是，只有非一致性矩阵的判断矩阵才需要进行一致性检验。
如果你的判断矩阵本身就是一个一致矩阵，那么就没有必要进行一致性检验。
输入判断矩阵：
clear;clc
disp('请输入判断矩阵A： ')
% A = input('判断矩阵A=')
A =[1 1 4 1/3 3;
1 1 4 1/3 3;
1/4 1/4 1 1/3 1/2;
3 3 3 1 3;
1/3 1/3 2 1/3 1]方法一：算术平均法求权重：
% 第一步：将判断矩阵按照列归一化（每一个元素除以其所在列的和）
Sum_A = sum(A)
[n,n] = size(A) % 也可以写成n = size(A,1)
% 因为我们的判断矩阵A是一个方阵，所以这里的r和c相同，我们可以就用同一个字母n表示
SUM_A = repmat(Sum_A,n,1) %repeat matrix的缩写
% 另外一种替代的方法如下：
SUM_A = [];
for i = 1:n %循环哦，这一行后面不能加冒号（和Python不同），这里表示循环n次
SUM_A = [SUM_A; Sum_A]
end
clc;A
SUM_A
Stand_A = A ./ SUM_A
% 这里我们直接将两个矩阵对应的元素相除即可
% 第二步：将归一化的各列相加(按行求和)
sum(Stand_A,2)
% 第三步：将相加后得到的向量中每个元素除以n即可得到权重向量
disp('算术平均法求权重的结果为：');
disp(sum(Stand_A,2) / n)
% 首先对标准化后的矩阵按照行求和，得到一个列向量
% 然后再将这个列向量的每个元素同时除以n即可（注意这里也可以用./哦）
方法二：几何平均法求权重：
% 第一步：将A的元素按照行相乘得到一个新的列向量
clc;A
Prduct_A = prod(A,2)
% prod函数和sum函数类似，一个用于乘，一个用于加 dim = 2 维度是行
% 第二步：将新的向量的每个分量开n次方
Prduct_n_A = Prduct_A .^ (1/n)
% 这里对每个元素进行乘方操作，因此要加.号哦。 ^符号表示乘方哦 这里是开n次方，所以我们等价求1/n次方
% 第三步：对该列向量进行归一化即可得到权重向量
% 将这个列向量中的每一个元素除以这一个向量的和即可
disp('几何平均法求权重的结果为：');
disp(Prduct_n_A ./ sum(Prduct_n_A))
方法三：特征值法求权重：
% 第一步：求出矩阵A的最大特征值以及其对应的特征向量
clc
[V,D] = eig(A) %V是特征向量, D是由特征值构成的对角矩阵（除了对角线元素外，其余位置元素全为0）
Max_eig = max(max(D)) %也可以写成max(D(:))哦~
% 那么怎么找到最大特征值所在的位置了？ 需要用到find函数，它可以用来返回向量或者矩阵中不为0的元素的位置索引。
% 那么问题来了，我们要得到最大特征值的位置，就需要将包含所有特征值的这个对角矩阵D中，不等于最大特征值的位置全变为0
% 这时候可以用到矩阵与常数的大小判断运算
D == Max_eig
[r,c] = find(D == Max_eig , 1)
% 找到D中第一个与最大特征值相等的元素的位置，记录它的行和列。
% 第二步：对求出的特征向量进行归一化即可得到我们的权重
V(:,c)
disp('特征值法求权重的结果为：');
disp( V(:,c) ./ sum(V(:,c)) )
% 我们先根据上面找到的最大特征值的列数c找到对应的特征向量，然后再进行标准化。
计算一致性比例CR：
%% 计算一致性比例CR
clc
CI = (Max_eig - 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]; %注意哦，这里的RI最多支持 n = 15
CR=CI/RI(n);
disp('一致性指标CI=');disp(CI);
disp('一致性比例CR=');disp(CR);
if CR<0.10
disp('因为CR < 0.10，所以该判断矩阵A的一致性可以接受!');
else
disp('注意：CR >= 0.10，因此该判断矩阵A需要进行修改!');
end
% % 注意：代码文件仅供参考，一定不要直接用于自己的数模论文中
% % 国赛对于论文的查重要求非常严格，代码雷同也算作抄袭