数学建模实验报告1层次分析法.doc

数学建模实验报告一、实验要求柴静的纪录片《穹顶之下》从独立媒体人的角度调查了席卷全国多个省份的雾霾的成因, 提出解决的方法有:关停重污染的钢铁厂、提高汽柴油品质、淘汰排放不达标汽车、提高洗煤率等,请仔细观看该纪录片,根据雾霾的成因,选择你认为治理雾霾确实可行的几个方案,并用 AHP 方法给出这几个主要方案的重要性排序。二、前期准备 1、理解层次分析法( AHP )的原理、作用,掌握其使用方法。 2、观看两遍柴静所拍摄的纪录片《穹顶之下》, 选出我认为可较为有效地治理雾霾的几个方法,初步确定各方法的有效性(即权重)。 3、初步拟定三个方案,每个方案中各个治理方法的权重不同。三、思路& 分析 1、根据纪录片《穹顶之下》和个人的经验判断给出各个记录雾霾的方法对于治理雾霾的判断矩阵,以及三个不同方案对于五大措施的判断矩阵。 2、了解了 AHP 的原理后,不难发现 MATLAB 在其中的作用主要是将判断矩阵转化为因素的权重矩阵。当然矩阵要通过一致性检验,得到的权重才足够可靠。 3、分别得到准则层对目标层、方案层对准则层的权重之后,进行层次总排序及一致性检验。得到组合权向量(方案层对目标层)即可确定适用方案。四、实验过程 1、确定层次结构 2、构造判断矩阵(1) 五大措施对于治理雾霾(准则层对目标层)的判断矩阵(2) 三个方案对于五大措施(方案层对准则层)的判断矩阵 3、层次单排序及一致性检验该部分在 MATLAB 中实现, 每次进行一致性检验和权向量计算时, 步骤相同, 输入、输出参数一致。(虽然输入的矩阵阶数可能不同,但可以不把矩阵阶数作为参数输入, 而通过[n,n]=size (A) 来算得阶数。) 因此考虑将这个部分定义为一个函数 judge , 输入一个矩阵 A, 打印一致性检验结果和权向量计算结果,并返回权向量、一致性指标 CI、平均随机一致性指标 RI 。将此脚本存为 ,在另一脚本 中调用。代码如下: 调试通过后,下面便用此函数进行一致性检验及权向量计算: (1) 准则层对目标层( A 矩阵) (2) 方案层对准则层( BB 矩阵) 代码: 结果: 注: 实际实验时, 一开始构造的五个矩阵中有两个没有通过一致性检验。反复调整后方才通过,考虑到实验报告的简洁性要求,不在此赘述调整判断矩阵的过程。 4、层次总排序及一致性检验原理简述如下: 因此,编写如下代码: 其中, max 是权重最大的方案的序号: 结果如下: 五、实验心得体会 1、函数化、模块化的