层次分析Matlab数学建模

1、层次分析法层次分析法（ The analytic hierarchy process）简称 AHP，在 20 世纪 70 年代中期由美国运筹学家托马斯·塞蒂（）正式提出。它是一种定性和定量相结合的、系统化、层次化的分析方法。 由于它在处理复杂的决策问题上的实用性和有效性，很快在世界范围得到重视。它的应用已遍及经济计划和管理、能源政策和分配、行为科学、军事指挥、运输、农业、教育、人才、医疗和环境等领域。层次分析法的基本思路与人对一个复杂的决策问题的思维、判断过程大体上是一样的。不妨用假期旅游为例：假如有3 个旅游胜地A、 B、 C 供你选择，你会根据诸如景色、费用和居住、饮食、 旅途条

2、件等一些准则去反复比较这3 个候选地点首先，你会确定这些准则在你的心目中各占多大比重，如果你经济宽绰、醉心旅游，自然分别看重景色条件，而平素俭朴或手头拮据的人则会优先考虑费用，中老年旅游者还会对居住、饮食等条件寄以较大关注。其次， 你会就每一个准则将3 个地点进行对比，譬如 A 景色最好， B 次之； B 费用最低，C 次之； C 居住等条件较好等等。最后，你要将这两个层次的比较判断进行综合，在A、 B、C 中确定哪个作为最佳地点。层次分析法的基本步骤1、建立层次结构模型。在深入分析实际问题的基础上，将有关的各个因素按照不同属性自上而下地分解成若干层次，同一层的诸因素从属于上一层的因素或对上层

3、因素有影响，同时又支配下一层的因素或受到下层因素的作用。 最上层为目标层， 通常只有 1 个因素， 最下层通常为方案或对象层， 中间可以有一个或几个层次， 通常为准则或指标层。 当准则过多时(譬如多于 9 个 )应进一步分解出子准则层。2、构造成对比较阵。从层次结构模型的第2 层开始，对于从属于(或影响 )上一层每个因素的同一层诸因素，用成对比较法和1 9 比较尺度构造成对比较阵，直到最下层。3、计算权向量并做一致性检验。对于每一个成对比较阵计算最大特征根及对应特征向量，利用一致性指标、随机一致性指标和一致性比率做一致性检验。若检验通过， 特征向量(归一化后 )即为权向量：若不通过，需重新构造

4、成对比较阵。4、计算组合权向量并做组合一致性检验。计算最下层对目标的组合权向量，并根据公式做组合一致性检验， 若检验通过， 则可按照组合权向量表示的结果进行决策， 否则需要重新考虑模型或重新构造那些一致性比率较大的成对比较阵。层次分析法的优点运用层次分析法有很多优点， 其中最重要的一点就是简单明了。 层次分析法不仅适用于存在不确定性和主观信息的情况， 还允许以合乎逻辑的方式运用经验、 洞察力和直觉。 也许层次分析法最大的优点是提出了层次本身， 它使得买方能够认真地考虑和衡量指标的相对重要性。建立层次结构模型将问题包含的因素分层：最高层（解决问题的目的）；中间层（实现总目标而采取的各精选文库种措

5、施、必须考虑的准则等。也可称策略层、约束层、准则层等）；最低层（用于解决问题的各种措施、方案等） 。把各种所要考虑的因素放在适当的层次内。用层次结构图清晰地表达这些因素的关系。例 1 购物模型某一个顾客选购电视机时，对市场正在出售的四种电视机考虑了八项准则作为评估依据，建立层次分析模型如下：例 2 选拔干部模型对三个干部候选人y1、 y2 、y3，按选拔干部的五个标准：品德、才能、资历、年龄和群众关系，构成如下层次分析模型：假设有三个干部候选人y1、 y2 、 y3，按选拔干部的五个标准：品德，才能，资历，年龄和群众关系，构成如下层次分析模型Image:选拔干部层次分析模型.gif-2精选文库

6、构造成对比较矩阵比较第i 个元素与第j 个元素相对上一层某个因素的重要性时，使用数量化的相对权重 aij 来描述。设共有n 个元素参与比较，则称为成对比较矩阵。成对比较矩阵中aij 的取值可参考 Satty的提议，按下述标度进行赋值。 aij 在1-9 及其倒数中间取值。aij = 1，元素 i与元素 j对上一层次因素的重要性相同；aij = 3，元素 i比元素 j略重要；aij = 5，元素 i比元素 j重要；aij = 7， 元素i 比元素j 重要得多；aij = 9，元素 i比元素 j的极其重要；aij = 2n， n=1,2,3,4，元素i 与 j 的重要性介于aij = 2n - 1

7、 与 aij = 2n + 1 之间；， n=1,2,.,9 ， 当且仅当 aji = n。成对比较矩阵的特点：。（备注：当i=j 时候， aij = 1）对例 2 ， 选拔干部考虑5 个条件：品德x1，才能 x2，资历 x3，年龄 x4，群众关系x5。某决策人用成对比较法，得到成对比较阵如下：a14 = 5 表示品德与年龄重要性之比为5，即决策人认为品德比年龄重要。作一致性检验从理论上分析得到：如果A 是完全一致的成对比较矩阵，应该有但实际上在构造成对比较矩阵时要求满足上述众多等式是不可能的。因此退而要求成对比较矩阵有一定的一致性，即可以允许成对比较矩阵存在一定程度的不一致性。由分析可知，对

8、完全一致的成对比较矩阵，其绝对值最大的特征值等于该矩阵的维数。对成对比较矩阵 的一致性要求， 转化为要求： 的绝对值最大的特征值和该矩阵的维数相差不大。-3精选文库检验成对比较矩阵A 一致性的步骤如下：计算衡量一个成对比较矩阵A （n>1 阶方阵）不一致程度的指标CI：RI是 这 样 得 到 的 : 对 于 固 定 的n, 随 机 构 造 成 对 比 较 阵A,其 中aij是 从1,2, ,9,1/2,1/3, ,1/9 中随机抽取的 . 这样的 A 是不一致的 , 取充分大的子样得到A 的最大特征值的平均值n123456789RI000.58 0.90 1.12 1.24 1.32 1

9、.41 1.45注解 :从有关资料查出检验成对比较矩阵A 一致性的标准RI：RI 称为平均随机一致性指标，它只与矩阵阶数n 有关。按下面公式计算成对比较阵A 的随机一致性比率CR：。判断方法如下：当 CR<0.1 时，判定成对比较阵A 具有满意的一致性，或其不一致程度是可以接受的；否则就调整成对比较矩阵A，直到达到满意的一致性为止。例如对例2 的矩阵计算得到,查得 RI=1.12，。这说明A 不是一致阵，但A 具有满意的一致性，A 的不一致程度是可接受的。此时 A 的最大特征值对应的特征向量为 U=(-0.8409,-0.4658,-0.0951,-0.1733,-0.1920) 。 这

10、个向量也是问题所需要的。 通常要将该向量标准化： 使得它的各分量都大于零， 各分量之和等于 1。该特征向量标准化后变成 U = (0.475,0.263,0.051,0.103,0.126)Z 。经过标准化后这个向量称为权向量。这里它反映了决策者选拔干部时，视品德条件最重要，其次是才能，再次-4精选文库是群众关系，年龄因素，最后才是资历。各因素的相对重要性由权向量U 的各分量所确定。求 A 的特征值的方法，可以用 MATLAB 语句求 A 的特征值 : Y,D =eig（ A）， D 为成对比较阵 的特征值， Y 的列为相应特征向量。在实践中，可采用下述方法计算对成对比较阵A = (aij)

11、的最大特征值max(A)和相应特征向量的近似值。定义可以近似地看作A 的对应于最大特征值的特征向量。计算可以近似看作A 的最大特征值。实践中可以由来判断矩阵A 的一致性。层次总排序及决策现在来完整地解决例2 的问题，要从三个候选人y1,y2,y3 中选一个总体上最适合上述五个条件的候选人。对此，对三个候选人y = y1,y2,y3 分别比较他们的品德(x1)，才能 (x2)，资历 (x3)，年龄 (x4)，群众关系 (x5)。先成对比较三个候选人的品德，得成对比较阵经计算， B1 的权向量 x1(Y) = (0.082,0.244,0.674)z故 B1 的不一致程度可接受。 x1(Y)可以直

12、观地视为各候选人在品德方面的得分。类似地，分别比较三个候选人的才能，资历，年龄，群众关系得成对比较阵-5精选文库通过计算知，相应的权向量为它们可分别视为各候选人的才能分，资历分，年龄分和群众关系分。经检验知 B2,B3,B4,B5 的不一致程度均可接受。最后计算各候选人的总得分。y1 的总得分从计算公式可知， y1 的总得分 (y1)实际上是 y1 各条件得分 x1(y1) , x2(y1) ,., x5(y1) ,的加权平均 , 权就是各条件的重要性。同理可得 y2,Y3 的得分为 z(y2) = 0.243,z(y3) = 0.4520.4570.2630.0510.1030.126总得分

13、Y10.0820.6060.4290.6360.1670.305Y20.2440.2650.4290.1850.1670.243Y30.6740.1290.1430.1790.6670.452即排名 :Y3 > Y1 > Y2比较后可得：候选人y3 是第一干部人选。层次分析法的用途举例例如，某人准备选购一台电冰箱，他对市场上的6 种不同类型的电冰箱进行了解后，在决定买那一款式时，往往不是直接拿电冰箱整体进行比较，因为存在许多不可比的因素，而是选取一些中间指标进行考察。例如电冰箱的容量、制冷级别、价格、型号、耗电量、外界信誉、 售后服务等。 然后再考虑各种型号冰箱在上述各中间标准下的

14、优劣排序。借助这种排-6精选文库序，最终作出选购决策。 在决策时， 由于 6 种电冰箱对于每个中间标准的优劣排序一般是不一致的，因此，决策者首先要对这7 个标准的重要度作一个估计，给出一种排序，然后把6种冰箱分别对每一个标准的排序权重找出来，最后把这些信息数据综合， 得到针对总目标即购买电冰箱的排序权重。有了这个权重向量，决策就很容易了。层次分析法应用的程序运用 AHP 法进行决策时，需要经历以下4 个步骤：1、建立系统的递阶层次结构；2、构造两两比较判断矩阵； （正互反矩阵）3、针对某一个标准，计算各备选元素的权重；4、计算当前一层元素关于总目标的排序权重。5、进行一致性检验。应用层次分析法

15、的注意事项如果所选的要素不合理，其含义混淆不清，或要素间的关系不正确，都会降低AHP 法的结果质量，甚至导致AHP 法决策失败。为保证递阶层次结构的合理性，需把握以下原则：1、分解简化问题时把握主要因素，不漏不多；2、注意相比较元素之间的强度关系，相差太悬殊的要素不能在同一层次比较。层次分析法应用实例1、建立递阶层次结构；2、构造两两比较判断矩阵；（正互反矩阵）对各指标之间进行两两对比之后，然后按 9 分位比率排定各评价指标的相对优劣顺序，依次构造出评价指标的判断矩阵。3、针对某一个标准，计算各备选元素的权重；关于判断矩阵权重计算的方法有两种，即几何平均法（根法） 和规范列平均法（和法）。（

16、1）几何平均法（根法）-7精选文库计算判断矩阵A 各行各个元素mi 的乘积；计算 mi 的 n 次方根；对向量进行归一化处理；该向量即为所求权重向量。（ 2）规范列平均法（和法）计算判断矩阵 A 各行各个元素 mi 的和；将 A 的各行元素的和进行归一化；该向量即为所求权重向量。计算矩阵 A 的最大特征值 ?max对于任意的i=1,2, ,n, 式中为向量AW 的第 i 个元素（ 4）一致性检验构造好判断矩阵后，需要根据判断矩阵计算针对某一准则层各元素的相对权重，并进行一致性检验。虽然在构造判断矩阵A 时并不要求判断具有一致性，但判断偏离一致性过大也是不允许的。因此需要对判断矩阵A 进行一致性

17、检验。clcclear%修改对比矩阵、一致性检验就可以a=1,1,1,4,1,1/21,1,2,4,1,1/21,1/2,1,5,3,1/21/4,1/4,1/5,1,1/3,1/31,1,1/3,3,1,12,2,2,3,3,1;x,y=eig(a);eigenvalue=diag(y);lamda=eigenvalue(1);ci1=(lamda-6)/ 5;cr1=ci1/1.24w1=x(:,1)/sum(x(:,1)b1=1,1/4,1/2;4,1,3;2,1/3,1;x,y=eig(b1);eigenvalue=diag(y);lamda=eigenvalue(1);ci21=(l

18、amda-3)/ 2;cr21=ci21/0.58-8精选文库w21=x(:,1)/sum(x(:,1)b2=11/41/5;411/2;521;x,y=eig(b2);eigenvalue=diag(y);lamda=eigenvalue(1);ci22=(lamda-3)/ 2;cr22=ci22/0.58w22=x(:,1)/sum(x(:,1)b3=131/3;1/311/7;371;x,y=eig(b3);eigenvalue=diag(y);lamda=eigenvalue(1);ci23=(lamda-3)/ 2;cr23=ci23/0.58w23=x(:,1)/sum(x(:,1)b4=11/35;317;1/51/71;x,y=eig(b4);eigenvalue=diag(y);lamda=eigenvalue(1);ci24=(lamda-3)/ 2;cr24=ci24/0.58w24=x(:,1)/sum(x(:,1)b5=117;117;1/71/71;x,y=eig(b5);eigenvalue=diag(y);lamda=eigenvalue(2