层次分析法步骤解析—根法、和法、幂法

层次分析法(AHP)AHP(AnalyticHierarchyProcess)方法，是由20世纪70年代由美国著名运筹学学家T.L.Satty提出的。它是指将决策问题的有关元素分解成目标、准则、方案等层次，在此基础上进行定性分析和定量分析的一种决策方法。这一方法的特点，是在对复杂决策问题的本质、影响因素及其内在关系等进行深入分析之后，构建一个层次结构模型，然后利用较少的定量信息，把决策的思维过程数学化，从而为求解多准则或无结构特性的复杂决策问题提供了一种简便的决策方法。AHP十分适用于具有定性的，或定性定量兼有的决策分析。这是一种十分有效的系统分析和科学决策方法，现在已广泛地应用在企业信用评级、经济管理规划、能源开发利用与资源分析、城市产业规划、企业管理、人才预测、科研管理、交通运输、水资源分析利用等方面。一、递阶层次结构的建立一般来说，可以将层次分为三种类型：（1）最高层：只包含一个元素，表示决策分析的总目标，因此也称为总目标层。（2）中间层：包含若干层元素，表示实现总目标所涉及的各子目标，包含各种准则、约束、策略等，因此也称为目标层。（3）最低层：表示实现各决策目标的可行方案、措施等，也称为方案层。典型的递阶层次结构如下：m一个好的递阶层次结构对解决问题极为重要，因此在建立递阶层次结构时，应注意到：（1）从上到下顺序地存在支配关系，用直线段（作用线）表示上一层次因素与下一层次因素之间的关系，同一层次及不相邻元素之间不存在支配关系。（2）整个结构不受层次限制。（3）最高层只有一个因素，每个因素所支配元素一般不超过9个，元素过多可进一步分层。（4）对某些具有子层次结构可引入虚元素，使之成为典型递阶层次结构。二、构造比较判断矩阵设有m个目标（方案或元素），根据某一准则，将这m个目标两两进行比较，把第i个目标（i=1,2,…,m）对第j个目标的相对重要性记为aij，(j=1,2,…,m)，这样构造的m阶矩阵用于求解各个目标关于某准则的优先权重，成为权重解析判断矩阵，准则1准则2准则3准则m1子准则1子准则2子准则3子准则m2方案1方案2方案3方案n总目标简称判断矩阵，记作A=(aij)m×m。Satty于1980年根据一般人的认知习惯和判断能力给出了属性间相对重要性等级表（见表4-4）。利用该表取aij的值，称为1-9标度方法。表4-4目标重要性判断矩阵A中元素的取值相对重要性定义说明1同等重要两个目标同样重要3略微重要由经验或判断，认为一个目标比另一个略微重要5相当重要由经验或判断，认为一个目标比另一个重要7明显重要深感一个目标比另一个重要，且这种重要性已有实践证明9绝对重要强烈地感到一个目标比另一个重要得多2,4,6,8两个相邻判断的中间值需要折中时采用若决策者能够准确估计aij（i,j,k=1,2,…,m），则有：aij=1/ajiaij=aik·akjaii=1定义4-1设A=(aij)m×m，A0，（即aij0;i,j=1,2,…,m），如果满足条件（1）aii=1（i=1,2,…,m）；（2）aij=1/aji（i,j=1,2,…,m），则称矩阵A为互反正矩阵。定义4-2设A=(aij)m×m，A0，如果满足条件aij=aik·akj（i,j,k=1,2,…,m）则称矩阵A为一致性矩阵。定理4-1对于任何一个m阶互反正矩阵A，均有max≥m，其中max是矩阵A的最大特征值。定理4-2m阶互反正矩阵A为一致性矩阵的充分必要条件是A的最大特征根为m。三、单准则下的排序层次分析法的信息基础是比较判断矩阵。由于每个准则都支配下一层若干因素，这样对于每一个准则及它所支配的因素都可以得到一个比较判断矩阵。因此根据比较判断矩阵如何求得各因素w1,w2,…,wm对于准则A的相对排序权重的过程称为单准则下的排序。这里设A=(aij)m×m，A0。（一）本征向量法利用AW=W求出所有的值，其中max为的最大值，求出max对应的特征向量W*，然后把特征向量W*规一化为向量W，则W=[w1,w2,…wm]T为各个目标的权重。求需要解m次方程，当m≥3时，计算比较麻烦，可以利用matlab来求解。（二）判断矩阵的近似解法判断矩阵是决策者主观判断的定量描述，求解判断矩阵不要求过高的精度。这里，介绍三种近似计算方法：根法、和法及幂法。幂法适于在计算机上运算。1、根法（1）A中每行元素连乘并开m次方，得到向量Tm),...,,(**2*1*其中，mmjijiaw1*（2）对W*作归一化处理，得到权重向量W=(w1,w2,…wm)T，其中miiii**/（3）对A中每列元素求和，得到向量S=(s1,s2,…sm)，其中sj=miija1（4）计算max的值，SWwsimii1max=miiiwAWm1)(12、和法（1）将A的元素按列作归一化处理，得矩阵Q=(qij)m×m。其中，mkkjijijaaq1/（2）将Q的元素按行相加，得向量Tm),...,,(21。其中，mjijiq1（3）对向量作归一化处理，得权重向量W=(w1,w2,…wm)T，其中mkkiiw1/（4）求出最大特征值miiiwAWm1max)(13、幂法幂法是一种逐步迭代的方法，经过若干次迭代计算，按照规定的精度，求出判断矩阵A的最大特征值及其对应的特征向量。定理3设矩阵A=(aij)m×m，A0，则CWeAeeAkTkklim，其中，W是A的最大特征值对应的的特征向量，C为常数，向量e=(1,1,…,1)T。幂法的计算步骤是：①任取初始正向量X(0)=(x1(0),x2(0),…,xm(0))T，计算0)0()0()0()0(0/},{maxmXYxXmii②迭代计算，对于k=0,1,2,…计算1)1()1()1()1(1)()1(/},{,maxkkkkiikkkkmXYxXmAYX③精度检查。当kkmm1时，转入步骤④；否则，令k=k+1，转入步骤②。④求最大特征值和对应的特征向量，将Y(k+1)归一化，即1max1)1()1(,/kmikikmyYW例判断矩阵125A=1/2171/51/71用幂法计算A的最大特征值max及其对应额特征向量。精度=0.0001。解：取初始向量X(0)=(1,1,1)T，迭代过程见下表kX(k)Y(k)011111118.00008.50001.34290.941210.158023.73122．57660.489110.69060.131133.03672.10830.426810.69430.141543.09612.18480.440610.70570.142353.12292.20180.443110.70500.141963.11952.19830.442610.70470.141973.11892.19800.442610.70470.141983.11892.19800.442610.70470.1419由上表看出，当k=7时，|m8-m7|=|3.1189-3.1189|=00.0001,迭代终止，得到max=3.1189，W=（0.5415,0.3816,0.0769）T四、单准则下的一致性检验由于客观事物的复杂性，会使我们的判断带有主观性和片面性，完全要求每次比较判断的思维标准一致是不太可能的。因此在我们构造比较判断矩阵时，我们并不要求n(n-1)/2次比较全部一致。但这可能出现甲与乙相比明显重要，乙与丙相比极端重要，丙与甲相比明显重要，这种比较判断会出现严重不一致的情况。我们虽然不要求判断具有一致性，但一个混乱的，经不起推敲的比较判断矩阵有可能导致决策的失误，所以我们希望在判断时应大体一致。而上述计算权重的方法，当判断矩阵过于偏离一致性时，其可靠程度也就值得怀疑了。因此，对于每一层次作单准则排序时，均需要作一致性的检验。一致性指标（ConsistencyIndex,CI）:1maxmmCI随机指标（RandomIndex,RI）一致性比率（ConsistencyRate,CR）:CR=CI/RI当CR取0.1时，最大特征值'max=CI·(m-1)+m=0.1·RI·(m-1)+m表4-5随机指标RI，'max取值表m123456789RI000.580.901.121.241.321.411.45'max3.1164.275.456.627.798.9910.16表中当n=1,2时，RI=0，这是因为1,2阶判断矩阵总是一致的。当n≥3时，若CR0.1即max'max，认为比较判断矩阵的一致性可以接受，否则应对判断矩阵作适当的修正，直到max小于'max通过一致性检验时，求得的W才有效。五、层次总排序计算同一层次中所有元素对最高层（总目标）的相对重要性标度（又称权重向量）称为层次总排序。1、层次总排序的步骤为：（1）计算同一层次所有因素对最高层相对重要性的权重向量，这一过程是自上而下逐层进行；（2）设已计算出第k-1层上有nk-1个元素相对总目标的权重向量为w(k-1)=(w1(k-1),w2(k-1),…,wn(k-1)(k-1))T（3）第k层有个nk个元素，他们对于上一层次（第k-1层）的某个元素j的单准则权重向量为pj(k)=(w1j(k),w2j(k),…,wnkj)(k))T（对于与k-1层第j个元素无支配关系的对应wij取值为0）；（4）第k层相对总目标的权重向量为wk=(p1(k),p2(k),…pk-1(k),)w(k-1)2、层次总排序的一致性检验人们在对各层元素作比较时，尽管每一层中所用的比较尺度基本一致，但各层之间仍可能有所差异，而这种差异将随着层次总排序的逐渐计算而累加起来，因此需要从模型的总体上来检验这种差异尺度的累积是否显著，检验的过程称为层次总排序的一致性检验。第k层的一致性检验指标CIk=(CI1(k-1),CI2(k-1),…,CInK(k-1))w(k-1)RIk=(RI1(k-1),RI2(k-1),…,RInK(k-1))w(k-1)CRk=CRk-1+CIk/RIk(3≤k≤n)当CRk0.1,可认为评价模型在第k层水平上整个达到局部满意一致性。六、递阶层次结构权重解析过程1、树状结构目标体系目标可分为多个层次，每个下层目标都隶属于一个而且只隶属一个上层目标，下层目标是对上层目标的具体说明。对于树状结构的目标体系，需由上而下逐步确定权重，即由树干向树梢，求树杈各枝相对于树杈的权重。2、网状结构目标体系网状结构的目标也分为多个层次，每个下层目标隶属于某几个上层目标（至少有一个下层目标隶属于不止一个上层目标）。七、AHP方法的基本步骤层次分析法大体分为以下六个步骤：1、明确问题；建立层次结构；2、两两比较，建立判断矩阵；2、层次单排序及其一致性检验；3、层次总排序及其一致性检验；4、根据分析计算结果，考虑相应的决策。例4-4（1）建立方案评价的递阶层次结构模型决策目标A商品价格B质量C企业财务状况D售后服务水平E候选人X候选人Y候选人Z（2）构造比较判断矩阵设以A为比较准则，B、C、D、E的两两比较判断矩阵为A，类似地构造矩阵B、C、D、E。ABCDEB1234C1/2123D1/31/213/2E1/41/32/31BXYZX13/24/5Y2/318/15Z5/415/81CXYZX13/21/2Y2/311/2Z221DXYZX132Y1/313/4Z1/24/31EXYZX13/21/2Y2/311/3Z231对于上述各比较判断矩阵，用Matlab数学软件求出其最大特征值及其对应的特征向量，将特征向量经归一化后，即可得相应的层次单排序的相对重要性权重向量，以及一致性指标CI和一致性比率CR，列表如下：层次单排序及一致性计算结果矩阵层次单排序的权重向量max'max)1/()(maxmmCIRICR=CI/RIA(0.4694,0