2021-02-10

层次分析法之购置个人电脑的决策

摘要：摘要是论文内容不加注释和评论的简短陈述，其作用是使读者不阅读论文全文即能获得必要的信息。
一般说来，摘要应包含以下五个方面的内容：
①研究的主要问题；
②建立的什么模型；
③用的什么求解方法；
④主要结果（简单、主要的）；
⑤自我评价和推广。
摘要中尽量不要有数学表达式。

关键字：*** *** ***

一 问题的重述
数学建模竞赛要求解决给定的问题，所以一般应以“问题的重述”开始。
此部分的目的是要吸引读者读下去，所以文字不可冗长，内容选择不要过于分散、琐碎，措辞要精练。
这部分的内容是将原问题进行整理，将已知和问题明确化即可。
注意：
在写这部分的内容时，绝对不可照抄原题！
应为：在仔细理解了问题的基础上，用自己的语言重新将问题描述一篇。应尽量简短，没有必要像原题一样面面俱到。
问题一：你要购置一台个人电脑，考虑电脑的功能价格等因素，如何做出决策。

二 问题的分析
2.1 对问题一的分析

购置电脑需考虑四方面：电脑预算，个人需求，电脑售后，电脑体验感。我们可以利用层次分析法计算出这几款电脑的得分即最终权重，其中将电脑预算，个人需求，电脑售后，电脑体验感作为四个评价指标。

三 模型假设

1）假设有P1,P2,P3三款电脑供你选择，你会根据电脑预算，个人需求，电脑售后，电脑体验感这些准则来反复比较这三款电脑。
2）假设比较三款电脑的同时，你的四个评价指标是恒定的。

四 符号说明

五 模型的建立与求解
5.1 问题一购置一台适合的电脑

5.1.1建立层次结构模型

层次结构模型分为三层，最上层为目标层O，即购置一台最合适的电脑；最下层为方案层，即提供选择的三款电脑P1,P2,P3；中间层为准则层，包括电脑预算C1，个人需求C2，电脑售后C3，电脑体验感C4四个指标：

5.1.2模型求解

①构造判断矩阵O-C:将准则层中的三个元素C1,C2,C3,C4两两比较，得出成对比较矩阵。

求解O-C判断矩阵的特征值得λmax =4.0566 ，且由特征值法计算出的权重向量ω =[0.3183，0.4869，0.1221，0.0727],由公式，计算得到CI= 0.0189 , CR=0.0212 < 0.1 ,通过一致性检验。

②构造判断矩阵C1-P,C2-P,C3-P,C4-P。


5.1.3模型求解结果与分析

根据5.1.2结合公式
，我们计算出P层每个方案所占的最终权重。将最终权重汇总成如下表格：
三款电脑的权重（得分）排序为P2>P1>P3，P2款电脑得分最高，所以在模型假设成立的条件下优先购买P2款电脑。分析：①上述判断矩阵的CR值均小于0.1，说明判断矩阵具有满意的一致性，可将其特征向量用于权重的定量描述。②从准则层对于目标层的权重向量来看，个人需求占的权重最大，电脑体验感占的权重最小，说明购买电脑更看重的是个人需求，电脑体验感对于购买电脑来说重要性相对较小。③从方案层对于目标层的权重向量来看，P2款电脑要优于其他两款电脑。

5.1.4模型的检验

稳健性检验。本文是采用特征值法计算权重，而不同的计算方法可能导致结果有所偏差，为了保证结果的稳健性，分别采用算数平均法，几何平均法再次计算最终权重，并进行比较和综合分析，这样避免了采用单一方法所产生的偏差，得出的结论更全面，更有效。

可以看出，三种方法计算出的最终权重几乎一样，所以说明我们的模型具有稳健性，有效性。
六 模型的评价与推广
这一部分应包括：
①你的模型完成了什么工作？达到了什么目的？得出了什么规律？
②你的建模方法是否有创造性？为今后的工作提供了什么思路？结果有什么理论或实际用途？
③模型中有何不足之处？有何改进建议？
④模型中有何遗留未解决的问题？以及解决这些问题可能的关键点和方向。

这一部分一定要有！

七 我们的建议
关于…的建议

八 参考文献
引用别人的成果或其他公开的资料(包括网上查到的资料)必须按照规定的参考文献的表述方式在正文引用处和参考文献中均明确列出。正文引用处用方括号标示参考文献的编号，如[1][3]等；引用书籍还必须指出页码。参考文献按正文中的引用次序列出，其中

书籍的表述方式为：
[编号] 作者，书名，出版地：出版社，出版年。

参考文献中期刊杂志论文的表述方式为：
[编号] 作者，论文名，杂志名，卷期号：起止页码，出版年。

参考文献中网上资源的表述方式为：
[编号] 作者，资源标题，网址，访问时间（年月日）。

九 附录
不便于编入正文的资料都收集在这里。
可以包括：
①某一问题的详细证明或求解过程；
②比较大的流程图；
③计算机的源程序及结果（或者使用傻瓜式操作软件的截图，例如SPSS）；
④较繁杂的图表或计算结果（一般结果只要不超过A4一页，尽量都放在正文中）。

层次分析法代码：
disp('请输入判断矩阵A')
A=input('A=');
[n,n] = size(A);
[V,D] = eig(A);
Max_eig = max(max(D));
[r,c]=find(D == Max_eig , 1);
disp('特征值法求权重的结果为：');
disp( V(:,c) ./ sum(V(:,c)) )
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];  
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
算术平均法求权重代码：
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))