层次分析法（原理+代码）

层次分析法-介绍

1.对评价类问题进行打分

例：

2.确定层级结构

3.确定权重指标

4.用“两两比较法”确定指标重要性

例：


矩阵a[i][j]表示i相较于j的重要程度，这个矩阵称为层次分析法中的判断矩阵

5.确定待选方案的最终得分

层次分析法-定义

什么是层次分析法

层次分析法-矩阵的一致性

一致矩阵和不一致矩阵

使用判断矩阵之前要检验矩阵的一致性

一致性检验方法

层次分析法-权重的计算

对于一致性矩阵进行归一化

对于非一致性矩阵用算术平均法求权重

算术平均法求权重论文表述

特征值法求权重（重点）

两种求权重方法的区别和联系

如果想要凑字数，可以两种全用然后求取平均值。

统计权重结果（层次总排序）

第一列的权重计算是通过“两两比较法”确定指标重要性这一部分内容来确定的

得出表格后要进行一致性检验

总结

层次结构图见层次分析-定义处

应用

（1）

（2）

这个论文存在一些问题，需先进行一致性检验然后计算权重，若一致性检验未通过则计算权重没有意义
（3）

层次分析图示要具体，如（1）（3）

层次分析法代码

%% 层次分析法
% 只有非一致矩阵才需要一致性检验
% % 要先进行一致性检验，通过后才能求权重，我们这里放到了最后做的一致性检验
%% 获取判断矩阵
disp('请输入判断矩阵A')
A=input('A=');
[n,n] = size(A);
%% 方法1： 算术平均法求权重
Sum_A = sum(A); %将A矩阵按列求和
SUM_A = repmat(Sum_A,n,1);
Stand_A = A ./ SUM_A;

disp('算术平均法求权重的结果为：');
w1 = sum(Stand_A,2)./n;
disp(w1)

%% 方法2： 特征值法求权重
[V,D] = eig(A);
Max_eig = max(max(D));% 先按列求最大值，得到行向量，再从这个向量里面求最大值
[r,c]=find(D == Max_eig , 1);
disp('特征值法求权重的结果为：');
w2 = V(:,c) ./ sum(V(:,c));
disp(w2)
disp('两种方法的平均权值为：');
disp((w1 + w2) ./ 2);
%% 计算一致性比例CR
CI = (Max_eig - n) / (n-1);
RI=[0 0.0001 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，注意检查！
% 这里n=2时，RI=0，我们为了避免分母为0，将这里的第二个元素改为了很接近0的正数
CR=CI/RI(n);
disp('最大特征值为:');
disp(Max_eig);
disp('一致性指标CI=');disp(CI);
disp('一致性比例CR=');disp(CR);
if CR<0.10
    disp('CR<0.10，该判断矩阵A的一致性可以接受!');
else
    disp('注意：CR >= 0.10，该判断矩阵需要进行修改!');
end