1、 层次分析法建模层次分析法建模 人们在日常生活中经常碰到许多决议问题:买一件衬衫,你要在棉、丝、涤纶及花、白、方格之中作出选择;请朋友吃饭,要策划是办家宴或去饭店,是吃中餐还是西餐或自助餐;假期旅游,是去风光绮丽苏杭,还是去迷人北戴河海滨,或者去山水甲天下桂林。假如认为这些日常小事无须作决议问题认真对待话,那么当你面对报考学校,挑选专业,或者选择工作岗位时候,就要慎重考虑、重复比较,尽也许地作出满意得决议了。1第1页第1页 从事各种职业人也经常面对抉择:一个厂长要决定购买哪种设备,上马什么产品;科技人员要选择研究课题;医生要为疑难病症拟定治疗方案;经理要从若干应试者中选拔秘书;各地域各部门官员
2、则要对人口、交通、经济、环境等领域发展规划作出决议。人们在处理上面这些决议问题时候,要考虑原因有多有少,有大有小,但是一个共同特点是他们通常涉及到经济、社会、人文等方面原因。再做比较、判断、评价、决议时,这些原因主要性、影响力或者优先程度往往难以量化,2第2页第2页人主观选择(当然要依据客观实际)会起着相称主要作用,这就给普通数学办法处理问题带来本质上困难。T.L.Saaty等人在七十年代提出了一个能有效地处理这样一类问题实用办法,称层次分析法层次分析法(Analytic Hierarchy Process,简记AHP).这是一个定性和定量相结合、系统化、层次化分析办法。过去研究自然和社会现象
3、主要有机理分析和统计分析两种办法,起着用典型数学工具分析现象因果关系,后者以随机数学为工具,共过大量观测数据寻求统计规律,近年来发展系统分析又是一个办法,而层次分析法就是系统分析数学工具之一。3第3页第3页 下面先简介层次分析法基本环节和应用实例,再讨论该办法在理论、计算以及建模等方面若干问题。4第4页第4页层次分析法基本环节层次分析法基本环节 层次分析法基本思绪与人对一个复杂决议问题思维、判断过程大体上是同样。不妨用前面提到假期旅游为例,假如有P1、P2、P33个旅游胜地供你选择,你会依据诸如景色、费用和居住、饮食、旅途条件等一些准则去重复比较那3个候选地点。首先,你会拟定这些准则在你心目中
4、各占多大比重,假如你经济宽裕、醉心旅游,自然尤其看重景色条件,而平素朴素或手头拮据人则会优先考虑费用,中老年旅游者还会对居住、饮食等条件寄以较大关注。5第5页第5页 另一方面,你会就每一个准则将3个地点进行对比,譬如P1景色最好,P2次之;P2费用最低,P3次之;P3居住等条件比较好等等。最后,你要将这两个比较判断进行综合,在P1、P2、P3中拟定哪个作为最佳地点。上面思维过程能够加工整理成为下几种环节:1、讲决议问题分解为3个层次、最上层为目的层,即选择旅游地,最下层为方案层,有P1、P2、P33个供选择地点,中间层为准则层,有景色、费用、居住、饮食、旅途5个准则,各层间联系涌现联直线表示(
5、图9-1)。6第6页第6页目的层准则层方案层图图9-1 选择旅游地层次结构7第7页第7页 2、通过互相比较拟定各准则对于目的权重,及各方案对于每一准则权重。这些权重在人思维过程中通常是定性,而在层次分析法中则要给出得到权重定量办法。3、将方案层对准则层权重及准则对目的层权重进行综合,最后拟定方案层对目的层权重。在层次分析法中要给出进行综合计算办法。层次分析法将定性分析与定量计算结合起来完毕上述环节,给出决议结果。下面我们来阐明如何比较同一层各原因对上层原因影响,从而拟定它们在上层原因中占权重。8第8页第8页 成对比较距阵和权向量 包括到社会、经济、人文等原因决议问题主要困难在于,这些原因通常不
6、易定量地测量。人们凭自己经验和知识进行判断,当原因较多时给出结果往往是不全方面和不准确,假如只是定性结果,则经常不容易被他人接受。Saaty等人做法,一是不把全部原因放在一起比较,而是两两相互对比,而是对比时采取相对尺度,以尽也许地降低性质不同诸原因相互比较困难,提升准确度。9第9页第9页 假设要比较某一层n个原因C1,C2,,Cn对上层一个原因O影响,如旅游决议问题中比较景色等5个准则在选择旅游地这个目的中主要性。每次取两个原因Ci和Cj,用aij表示Ci和Cj对O影响之比,所有比较结果可用对比比较距阵 A=(aij)nn,aij 0 ,(1)表示。由(1)给出aij特点,A称为正互反矩正互
7、反矩阵阵。显然比由aii=1。如用C1,,C5依次表示景色、费用、饮食、旅游5个准则,设某人用成对比较距阵(正互反阵)为10第10页第10页(2)(2)中a12=,表示景色C1和给用C2对选择旅游地这个目的O主要性之比为1:2;a13=4表示景色C1和居住条件C3之比为4:1;a23=7表示费用C2与居住条件C3之比为7:1。能够看出在此人选择旅游地时,费用原因最主要,景色次之。如何由成对比较阵拟定诸原因C1,Cn对上层原因O权重。11第11页第11页 仔细分析一下(2)式给出成对比较阵A能够发觉,既然C1与C2之比为1:2;C1与C3之比为4:1。那么C2与C3之比由于8:1而不是7:1才干
8、阐明成对比较是一致。但是,n个原因要做 次,所有一致要求是太苛刻了。Saaty等人给出了在成对比较不一致情况下计算各原因C1,,Cn对原因O权重办法,并且拟定了这种不一致允许范围。为了阐明这点我们先当作对比较完全一致。设想把一块大石头O砸成n块小石头C1,,Cn,假如准确地称出它们重量为w1,,wn,在作成对比较时令aij=wi/wj,12第12页第12页那么得到(3)这些比较显然是一致,n块小石头对大石头权重(即在大石头中占比重)可用向量 w=(w1,w2,,wn)T13第13页第13页表示,假如大石头为单位重量,则有显然,A各个列向量与w仅相差一个百分比因子。普通地,假如一个正互反阵A满足
9、 aijajk=aik,i,j,k=1,2,,n (4)则A称为一致性距阵,简称一致一致阵阵。(3)式给出A显然是一致阵。容易证实n阶一致阵A有下列性质。1、A秩为1,A唯一非零特性根为n;2、A任一列(行)向量都是相应于特性根n特性值。14第14页第14页 假如得到成对比较阵是一致阵,像(3)式A,自然应取相应于特性根n。归一化特性向量(即分量之和为1)表示诸原因C1,Cn对上层原因O权重,这个向量称为权向量权向量。假如成对比较阵A不是一致阵,但在不一致允许范围内(下面将阐明如何拟定这个范围),Saaty等人提议用相应于A最大特性根(即作)特性向量(归一化后)作为权向量w,即w满足 Aw=w
10、 (5)直观地看,由于矩阵A特性根和特性向量连续地依赖于矩阵元素aij,因此当aij离一致性要求不远时,A特性根和特性向量也与一致阵相差不大。15第15页第15页 (5)式表示方法称为由成对比较阵求权向量特性根法。求和w简便算法和特性根法更深入意义,以及其它求权向量方法见9.3阶。比较尺度 当比较两个可能含有不同性质原因Ci和Cj对于一个上层原因O影响时,采取什么样尺度aij很好呢?Saaty等人提出用1-9尺度,即aij取值范围是1,2,9及其互反数1,。理由以下。16第16页第16页 1、再进行定性成对比较时,人们头脑中通常有5中明显等级,用1-9尺度能够以便地表示下列。表表9-1 1-9
11、尺度尺度aij含义含义 尺度aij含 义1Ci与Cj得影响相同相同3Ci与Cj得影响稍强稍强5Ci与Cj得影响强强7Ci与Cj得影响明显地强9Ci与Cj得影响绝对地强绝对地强2,4,6,8Ci与Cj得影响之比在上述两个相邻等级之间1,,,Ci与Cj得影响之比为上面aij互反数17第17页第17页2、心理学家认为,进行成对比较原因太多,将超出人判断能力,最多达之72范围。如以9个为限,用1-9尺度表示它们之间差异正适当。3、Saaty曾用1-3,1-5,1-17,(d+0.1)-(d+0.9)(d=1,2,3,4),1p-9p(p=2,3,4,5)等共27中比较尺度,对在不同距离出判断某光源亮度
12、等实例结构成对比较阵,并算出权向量。把这些权向量与按照光强定律等物理知识得到实际权向量进行对比发觉,1-9尺度不但在简朴尺度中最好,而且结果并不劣于较复杂尺度。18第18页第18页 依据上述定理和连续地依赖于aij事实可知比n大得多,A不一致程度越严重,用特性向量作为权向量引起判断误差越大。因而能够用-n数值大小来衡量A不一致程度,Saaty将 定义为一致指标一致指标.CI=0时A为一致阵;CI越大A不一致程度越严重,注意到An个特性根之和等于A对角元素之和(为何?),而A对角元素均为1,因此特性根之和 .(不妨记1=)。由此可知,一致性指标CI相称于处外其余n-1个特性根平均值(取绝对值)。
13、19第19页第19页 为了确定A不一致程度允许范围,需要找出衡量A一致性指标CI标准。Saaty又引入了所谓随机一致性指标RI,计算RI过程是:对于固定n,随机地结构正互反阵A/(它元素aij/(ij)从1-9,1-1/9中取随机值,aji/为aij/互反数,aii/=1),然后计算A/一致性指标CI。能够想象到,A/是非常不一致,它CI相称大。如此结构相称多A/,用他们CI平均值作为随机一致性指标。Saaty对于不同n(=111),用100-500个样本A/算出随机一致性指标RI数值以下。20第20页第20页 表表9-2 随机一致性指标随机一致性指标RI数值数值 n1234567891011
14、RI000.580.901.121.241.321.411.451.491.51表中n-1,2时RI=0,是由于1,2阶正互反阵总是一致阵。对于n3成对比较阵A,将它一致性指标CI与同阶(指n相同)随机一致性指标RI之比成为一致性比率一致性比率CR 当(7)21第21页第21页时认为A不一致程度在允许范围之内,可用其特性向量作为权向量。不然要重新进行成对比较,对A加以调整。顺便指出,(7)式中0.1选取是带有一定主观信度。对于A利用(6),(7)式和表9-2进行检查成为一致性检查一致性检查。对于(2)式给出A能够算出,=5.073,归一化特性向量 w=(0.263,0.475,0.055,0.
15、099,0.110)T.由(6)式22第22页第22页在表9-2中查出RI=1.12.按(7)式计算,于是通过了一致性检查,故上述w可作为权向量。23第23页第23页 组合权向量组合权向量 在旅游决议问题中我们已经得到了第2层(准则层)对第1层(目的层,只有一个原因)权向量,记作w(2)=(w1(2),,w5(2)T(即由(2)式A算出w)。用同样办法结构第3层(方案层,见图9-1)对第2层每一个准则成对比较阵,不妨设他们为24第24页第24页 这里距阵Bk(k=1,,5)中元素bij(k)是方案(旅游地)Pi与Pj对于准则Ck(景色、费用等)优越性比较尺度。由第3 层成对比较阵Bk计算出权向
16、量wk(3),最大特性根k和一致性指标CIk,结果列入下表 表表9-3 旅游决旅游决议问题议问题第第3层计层计算算结结果果k123450.5950.0820.4290.6330.166wk(3)0.2770.2360.4290.1930.1660.1290.6820.1420.1750.6683.0053.00233.0093CIk0.0030.00100.005025第25页第25页不难看出,由于n=3时随机一致性指标RI=0.58(表9-2),因此上面CIk均可通过一致性检查。下面问题是由各准则对目的权向量w(2)和各方案对一准则权向量wk(3)(k=1,,5),计算各方案对目的权向量,称
17、为组合权向量组合权向量,记作w(3)。对于方案P1,它在景色等5个准则中权重用wk(3)第1个分量表示(图9-3中wk(3)第一行),而5各准则对于目的权重又用权向量w(2)表示,因此方案P1在目的中组合权重应为它们相应两两乘积之和,即 0.595*0.263+0.82*0.475+0.429*1.055 +0.633*0.099+0.166*0.110=0.300 26第26页第26页同样能够计算P2,P3 在目的中组合权重为0.246和0.456,于是组合向量w(3)=(0.300,0.246,0.456)T.结果表明方案P3在旅游地选择中占权重近于1/2,远不小于P1,P2,应作为第1选
18、择地点。27第27页第27页 由上述计算可知,对于3个层次决议问题。若第1层只有1个原因,第2、3层分别有n、m个原因,记第2、3层对第1、2层权向量分别为:w(2)=(w1(2),,wn(2)T wk(3)=(wk1(3),,wkn(3)T,k=1,2,n以wk(3)为列向量构成矩阵 W(3)=w1(3),,wn(3)则第3层对第1层组合权向量为 w(3)=W(3)w(2)(8)28第28页第28页 更普通地,若共有s层,则第k层对第1层(设只有1个原因)组合权向量为 w(k)=W(k)w(k),k =3,4,,s (9)其中W(k)是以第k层对第k-1层权向量为列向量构成矩阵。于是最下层(
19、第s层)对最上面一层组合权向量为 w(s)=W(s)W(s-1)W(3)w(2)(10)29第29页第29页 组合一致性检查组合一致性检查 在层次分析整个计算过程中,除了对每个成对比较阵进行一致性检查,以判断每个权向量是否能够应用外,还要进行所谓组合一致性检查,一边拟定组合权向量是否能够作为最后决议依据。组合一致性检查能够逐层进行。若第p层一致性指标为CI1(p),,CIn(p)(n是第p-1层原因数目),随机一致性指标为RI1(p),,RIn(p),定义 CI(p)=CI1(p),CIn(p)w(p-1)(11)RI(p)=RI1(p),RIn(p)w(p-1)(12)30第30页第30页则
20、第p层对第1层组合一致性比率组合一致性比率为 CR(p)=CR(p-1)+CI(p)/RI(p),p=3,4,,s (13)其中CR(2)为由(7)式计算一致性比率。最后,当最下层对最上层组合一致性比率 CR(s)0.1 (14)认为整个层次比较判断通过一致性检查。在旅游决议问题中能够算出CI(3)=0.00176,RI(3)=0.58,前面已有了CR(2)=0.016,于是由(13)得到31第31页第31页 CR(3)=0.016+0.00176/0.58=0.0190.1经过了组合一致性检验,前面得到组合权向量w(3)能够作为最终决议依据.在本节最终,将层次分析法基本步骤归纳以下。1、建立
21、层次结构模型 在深入分析实际问题基础上,将相关各个原因按照不同属性自上而下地分解成若干层次。同一层诸原因从属于上一层原因或对上层原因有影响,同时又支配下一层原因或受到下层原因作用。最上层为目标层,通常只有1个原因,最下层通常为方案层或对象层,中间能够由1个或几个层次,通常为准则或指标层。32第32页第32页当准则层过多时(譬如多于9个)应进一步分解出子准则层。2、结构成对比较阵结构成对比较阵 从层次结构模型第2层开始,对于从属于(或影响及)上一层每个原因同一层诸原因,用成对比较法和1-9比较尺度结构成对比较阵,直到最下层。3、计算权向量并做出一致性检查计算权向量并做出一致性检查 对于每一个成对比较阵计算最大特性根及相应特性向量(计算办法见9.3节),利用一致性指标,随机一致性指标和一致性比率做一致性检查。若检查通过,特性向量(归一化后)及做权向量;若不通过,需重新结构成对比较阵。33第33页第33页 4、计算组合权向量并做组合一致性检查计算组合权向量并做组合一致性检查 利用(10)式计算最下层对目的层组合权向量,并依据(11)(14)式做组合一致性检查。若检查通过,则可按照组合权向量表示结果进行决议,不然需要重新考虑模型或重新结构那些一致性比率CR较大成对比较阵。34第34页第34页