数学建模(层次分析法(AHP法))

层次分析法(AHP法)AnalyticHierarchyProcess引言层次分析法（AHP）是美国运筹学家匹茨堡大学教授萨蒂(T.L.Saaty)于上世纪70年代初，为美国国防部研究“根据各个工业部门对国家福利的贡献大小而进行电力分配”课题时，应用网络系统理论和多目标综合评价方法，提出的一种层次权重决策分析方法。这种方法的特点是在对复杂的决策问题的本质、影响因素及其内在关系等进行深入分析的基础上，利用较少的定量信息使决策的思维过程数学化，从而为多目标、多准则或无结构特性的复杂决策问题提供简便的决策方法。是对难于完全定量的复杂系统作出决策的模型和方法。层次分析法在经济、科技、文化、军事、环境乃至社会发展等方面的管理决策中都有广泛的应用。常用来解决诸如综合评价、选择决策方案、估计和预测、投入量的分配等问题。层次分析法建模一、问题的提出日常生活中有许多决策问题。决策是指在面临多种方案时需要依据一定的标准选择某一种方案。例1某人准备选购一台电冰箱他对市场上的6种不同类型的电冰箱进行了解后，选取一些中间指标进行考察。例如电冰箱的容量、制冷级别、价格、型式、耗电量、外界信誉、售后服务等。然后再考虑各种型号冰箱在上述各中间标准下的优劣排序。借助这种排序，最终作出选购决策。在决策时，由于6种电冰箱对于每个中间标准的优劣排序一般是不一致的，因此，决策者首先要对这7个标准的重要度作一个估计，给出一种排序，然后把6种冰箱分别对每一个标准的排序权重找出来，最后把这些信息数据综合，得到针对总目标即购买电冰箱的排序权重。有了这个权重向量，决策就很容易了。例2旅游假期旅游，是去风光秀丽的苏州，还是去凉爽宜人的北戴河，或者是去山水甲天下的桂林？通常会依据景色、费用、食宿条件、旅途等因素选择去哪个地方。例3择业面临毕业，可能有高校、科研单位、企业等单位可以去选择，一般依据工作环境、工资待遇、发展前途、住房条件等因素择业。例4科研课题的选择由于经费等因素，有时不能同时开展几个课题，一般依据课题的可行性、应用价值、理论价值、被培养人才等因素进行选题。分解建立确定计算判断实际问题层次结构多个因素诸因素的相对重要性权向量综合决策一、层次分析法基本原理二、层次分析法的步骤和方法运用层次分析法构造系统模型时，大体可以分为以下四个步骤：1.建立层次结构模型2.构造判断(成对比较)矩阵3.层次单排序及其一致性检验4.层次总排序及其一致性检验将决策的目标、考虑的因素（决策准则）和决策对象按它们之间的相互关系分为最高层、中间层和最低层，绘出层次结构图。最高层：决策的目的、要解决的问题。最低层：决策时的备选方案。中间层：考虑的因素、决策的准则。对于相邻的两层，称高层为目标层，低层为因素层。1建立层次结构模型一个典型的层次可以用下图表示出来：几点注意1.处于最上面的的层次通常只有一个元素，一般是分析问题的预定目标或理想结果。中间层次一般是准则、子准则。最低一层包括决策的方案。层次之间元素的支配关系不一定是完全的，即可以存在这样的元素，它并不支配下一层次的所有元素。2.层次数与问题的复杂程度和所需要分析的详尽程度有关。每一层次中的元素一般不超过9个，因一层中包含数目过多的元素会给两两比较判断带来困难。3.一个好的层次结构对于解决问题是极为重要的。层次结构建立在决策者对所面临的问题具有全面深入的认识基础上，如果在层次的划分和确定层次之间的支配关系上举棋不定，最好重新分析问题，弄清问题各部分相互之间的关系，以确保建立一个合理的层次结构。目标层O(选择旅游地)P2黄山P1桂林P3北戴河准则层方案层C3居住C1景色C2费用C4饮食C5旅途例1.选择旅游地如何在3个目的地中按照景色、费用、居住条件等因素选择.例2大学毕业生就业选择问题获得大学毕业学位的毕业生，在“双向选择”时，用人单位与毕业生都有各自的选择标准和要求。就毕业生来说选择单位的标准和要求是多方面的，例如：①能发挥自己才干作出较好贡献（即工作岗位适合发挥自己的专长）；②工作收入较好（待遇好）；③生活环境好（大城市、气候等工作条件等）；④单位名声好（声誉等）；⑤工作环境好（人际关系和谐等）⑥发展晋升机会多（如新单位或前景好）等。工作选择可供选择的单位P1’P2，Pn贡献收入发展声誉工作环境生活环境目标层准则层方案层将决策问题分为3个或多个层次：最高层：目标层。表示解决问题的目的，即层次分析要达到的总目标。通常只有一个总目标。中间层：准则层、指标层、…。表示采取某种措施、政策、方案等实现预定总目标所涉及的中间环节；一般又分为准则层、指标层、策略层、约束层等。最低层：方案层。表示将选用的解决问题的各种措施、政策、方案等。通常有几个方案可选。每层有若干元素，层间元素的关系用相连直线表示。建立层次结构模型的思维过程的归纳层次分析法所要解决的问题是关于最低层对最高层的相对权重问题，按此相对权重可以对最低层中的各种方案、措施进行排序，从而在不同的方案中作出选择或形成选择方案的原则。在建立递阶层次结构以后，上下层次之间元素的隶属关系就被确定了。假定上一层次的元素Ck作为准则，对下一层次的元素A1,…,An有支配关系，我们的目的是在准则Ck之下按它们相对重要性赋予A1,…,An相应的权重。2构造判断(成对比较)矩阵比较同一层次中每个因素关于上一层次的同一个因素的相对重要性在确定各层次各因素之间的权重时，如果只是定性的结果，则常常不容易被别人接受，因而Saaty等人提出构造：成对比较矩阵A=(aij)nn，即：1.不把所有因素放在一起比较，而是两两相互比较。2.对此时采用相对尺度，以尽可能减少性质不同的诸因素相互比较的困难，以提高准确度。心理学家认为成对比较的因素不宜超过9个，即每层不要超过9个因素。成对比较矩阵是表示本层所有因素针对上一层某一个因素的相对重要性的比较。判断矩阵的元素aij用Saaty的1—9标度方法给出。判断矩阵元素aij的标度方法标度含义1表示两个因素相比，具有同样重要性3表示两个因素相比，一个因素比另一个因素稍微重要5表示两个因素相比，一个因素比另一个因素明显重要7表示两个因素相比，一个因素比另一个因素强烈重要9表示两个因素相比，一个因素比另一个因素极端重要2，4，6，8上述两相邻判断的中值倒数因素i与j比较的判断aij，则因素j与i比较的判断aji=1/aij对于n个元素A1,…,An来说，通过两两比较，得到成对比较（判断）矩阵A=(aij)nn：其中判断矩阵具有如下性质：（1）aij0；（2）aij=1/aji；（3）aii=1。我们称A为正的互反矩阵。根据性质（2）和（3），事实上，对于n阶判断矩阵仅需对其上（下）三角元素共n(n-1)/2个给出判断即可。11/2433217551/41/711/21/31/31/52111/31/5311AijjiijnnijaaaaA1,0,)(要比较各准则C1,C2,…,Cn对目标O的重要性ijjiaCC:A~成对比较阵选择旅游地目标层O(选择旅游地)准则层C3居住C1景色C2费用C4饮食C5旅途C1C2C3C4C5C1C2C3C4C5稍加分析就发现上述成对比较矩阵有问题旅游问题的成对比较矩阵共有6个（一个5阶，5个3阶）。用权值表示影响程度，先从一个简单的例子看如何确定权值。例如一块石头重量记为1，打碎分成n小块，各块的重量分别记为：w1,w2,…wn则可得成对比较矩阵11222112111nnnn由右面矩阵可以看出，jkkiji3层次单排序及其一致性检验即nji,,2,1,2321132321137,2,4aaaaaaAijkjikaaaijkjikaaa但在例2的成对比较矩阵中，在正互反矩阵A中，若,(A的元素具有传递性)则称A为一致阵。定理：n阶正互反阵A的最大特征根maxn,当且仅当=n时A为一致阵一般地，我们并不要求判断具有这种传递性和一致性，这是由客观事物的复杂性与人的认识的多样性所决定的。但在构造两两判断矩阵时，要求判断大体上的一致是应该的。出现甲比乙极端重要，乙比丙极端重要，而丙又比甲极端重要的判断，一般是违反常识的。一个混乱的经不起推敲的判断矩阵有可能导致决策的失误，而且当判断矩阵过于偏离一致性时，用上述各种方法计算的排序权重作为决策依据，其可靠程度也值得怀疑。因而必须对判断矩阵的一致性进行检验。由于λ(A的特征根)连续的依赖于aij，则λ比n大的越多，A的不一致性越严重。引起的判断误差越大。因而可以用λ-n数值的大小来衡量A的不一致程度。1nnCI定义一致性指标:CI=0，有完全的一致性CI接近于0，有满意的一致性CI越大，不一致越严重1.0RICICR一致性检验：利用一致性指标和一致性比率0.1及随机一致性指标的数值表，对A进行检验的过程。一般，当一致性比率的不一致程度在容许范围之内，有满意的一致性，通过一致性检验。否则要重新构造成对比较矩阵A，对aij加以调整。时，认为A定义一致性比率：RICICR判断矩阵一致性检验的步骤如下：(1)计算一致性指标C.I.：1C.I.maxnn其中n为判断矩阵的阶数；(2)查找平均随机一致性指标R.I.：平均随机一致性指标是多次（500次以上）重复进行随机判断矩阵特征根计算之后取算术平均得到的。龚木森、许树柏1986年得出的1—15阶判断矩阵重复计算1000次的平均随机一致性指标如下：阶数12345678R.I.000.520.891.121.261.361.41阶数9101112131415R.I.1.461.491.521.541.561.581.59(3)计算一致性比例C.R.：当C.R.0.1时，一般认为判断矩阵的一致性是可以接受的。否则应对判断矩阵作适当的修正。R.I.C.I.C.R.“选择旅游地”中准则层对目标的权向量及一致性检验11/2433217551/41/711/21/31/31/52111/31/5311A准则层对目标的成对比较阵最大特征根max=5.0735.07350.01851CI一致性指标随机一致性指标RI=1.12(查表)一致性比率CR=0.018/1.12=0.0160.1通过一致性检验旅游问题的成对比较矩阵共有6个（一个5阶，5个3阶）。问题：两两进行比较后，怎样才能知道，下层各因素对上层某因素的影响程度的排序结果呢？计算单一准则下元素的相对权重这一步是要解决在准则Ck下，n个元素A1,…,An排序权重的计算问题。对于n个元素A1,…,An，通过两两比较得到判断矩阵A，解特征根问题Aw=maxw所得到的w(特征向量)经归一化后作为元素A1,…,An在准则Ck下的排序权重，这种方法称为计算排序向量的特征根法。特征根方法的理论依据是如下的正矩阵的Person定理，它保证了所得到的排序向量的正值性和唯一性：定理设n阶方阵A0，max为A的模最大的特征根，则有(1)max必为正特征根，而且它所对应的特征向量为正向量；(2)A的任何其它特征根恒有||max；(3)max为A的单特征根，因而它所对应的特征向量除差一个常数因子外是唯一的。特征根方法中的最大特征根max和特征向量w，可用Matlab软件直接计算。例如：计算矩阵的最大特征值及相应的特征向量。111411/2112411/211/21531/21/41/41/511/31/3111/3311/3222331相应的Matlab程序如下：A=[1,1,1,4,1,1/2;1,1,2,4,1,1/2;1,1/2,1,5,3,1/2;…1/4,1/4,1/5,1,1/3,1/3;1,1,1/3,3,1,1/3;2,2,2,3,3,1];[x,y]=eig(A);eigenvalue=diag(y);lam