1、第九章第九章 层次分析层次分析 The Analytic Hierarchy Process (AHP)第九章 层次分析在管理中,人们常常需要对一些情况作出决策:例如企业的决策者要决定购置哪种设备,上马什么产品;经理要从若干求职者中决定录用哪些人员;地区、部门官员要对人口、交通、经济、环境等领域的发展规划作出决策。在日常生活中也常会遇到,在多种类不同特征的商品中选购。报考学校选择志愿。毕业时选择工作岗位等。这一系列的问题,单纯靠构造一个数学模型来求解的方法往往行不通,而用完全主观的定夺也常常表现为举棋不定,而最终选择不理想,甚至不满意的决策方案。面对这样的问题,运筹学者开始了对人们思维决策过程
2、进行分析、研究。第九章 层次分析美国运筹学家,T.L.Saaty等人在九十年代提出了一种能有效处理这类问题的实用方法,称之为层次分析法(AHP法)T.L.Saaty等曾把它用于电力工业计划,运输业研究,美国高等教育事业1985-2000展望,1985年世界石油价格预测等方面。第九章 层次分析这种方法的特征:定性与定量相结合,把人们的思维过程层次化,数量化。AHP法作为一种决策方法是在1982年11月召开的中美能源、资源、环境学术会议上,有Saaty学生H.Gholamnezhad首先向中国介绍的。以后层次分析法在中国得到很大的发展,很快应用到能源系统分析,城市规划,经济管理科研成果评价的许多领
3、域。第九章 层次分析9.1 层次分析法的基本步骤 运用AHP法进行决策时,大体可以分为4个步骤进行: (1) 分析系统中各个因素的关系,建立系统的递阶层次结构; (2) 对同一层次的各元素关于上一层次中某一准则的重要性进行两两比较,构造两两比较判断矩阵;第九章 层次分析 (3) 由判断矩阵计算被比较元素对于该准则的相对权重; (4) 计算各层元素对系统目标的合成权重,并进行排序。第九章 层次分析 一、建立层次分析的结构模型: 用AHP分析问题,首先要把问题条理化、层次化,构造层次分析的结构模型。这些层次大体上可分为3类: 1、最高层:在这一层次中只有一个元素,一般是分析问题的预定目标或理想结果
4、,因此又称目标层;第九章 层次分析 2、中间层:这一层次包括了为实现目标所涉及的中间环节,它可由若干个层次组成,包括所需要考虑的准则,子准则,因此又称为准则层; 3、最底层:表示为实现目标可供选择的各种措施、决策、方案等,因此又称为措施层或方案层。 层次分析结构中各项称为此结构模型中的元素。第九章 层次分析决策目标准则1方案1准则m1准则2子准则1方案2子准则2方案mr子准则m2第九章 层次分析 注:层次之间的支配关系不一定是完全的,即可以有元素(非底层元素)并不支配下一层次的所有元素而只支配其中部分元素。这种自上而下的支配关系所形成的层次结构,我们称之为递阶层次结构。 递阶层次结构中的层次数
5、与问题的复杂程度及分析的详尽程度有关,一般可不受限制。 第九章 层次分析 为了避免由于支配的元素过多而给两两比较判断带来困难,每层次中各元素所支配的元素一般地不要超过9个,若多于9个时,可将该层次再划分为若干子层。 例1、某顾客选购电冰箱时,对市场上正在出售的四种电冰箱考虑6项准则作为评价依据,得到如下层次分析模型:第九章 层次分析目标层:准则层:方案层:信 誉 T 1A型 式 T 2B价 格 T 3C容 量 T 4D制 冷 级 别 T 5耗 电 量 T 6选 购 电 冰 箱第九章 层次分析 例2、选择科研课题: 某研究单位现有3个科研课题,限于人力物力,只能承担其中一个课题,如何选择? 考虑
6、下列因素: 成果的贡献大小,对人材培养的作用,课题可行性。 在成果贡献方面考察:应用价值及科学第九章 层次分析 意义(理论价值,对某科技领域的推动作用); 在课题可行性方面考虑:难易程度(难易程度与自身的科技力量的一致性),研究周期(预计需要花费的时间),财政支持(所需经费,设备及经费来源,有关单位支持情况等)。第九章 层次分析目标层合理选择科研课题A成果贡献B1人才培养B2课题可行性B3课题D1课题D2课题D3应用价值 c1科学意义 c2难易程度 c3研究周期 c4财政支持 c5第九章 层次分析方案层准则层 例3、设某港务局要改善一条河道的过河运输条件,为此需要确定是否要建立桥梁或隧道以代替
7、现有轮渡。 此问题中过河方式的确定取决于过河方式的效益与代价(即成本)。通常我们用费效比(效益/代价)作为选择方案的标准。为此构造以下两个层次分析的结构模型。第九章 层次分析准则层过河的效益A经济效益B1社会效益B2环境效益B3桥梁D1隧道D2渡船D3收入 c2岸间商业 c3节省时间c1当地商业c4建筑就业c5安全可靠c6交往沟通c7自豪感c8舒 适c9进出方便c10美 化c11第九章 层次分析方案层目标层过河的代价A经济代价B1社会代价B2环境代价B3桥梁D1投入资金c1操作维护c2冲击渡船业c3冲击生活方式c4交通拥挤 c5居民搬迁 c6汽车排废物 c7对水的污染 c8对生态的破坏c9隧道
8、D2渡船D3第九章 层次分析目标层准则层方案层 二、构造判断矩阵: 上、下层之间关系被确定之后,需确定与上层某元素Z(目标A或某个准则Z)相联系的下层元素(x1,x2,xn)各在上层元素Z之中所占的比重。 方法:每次取2个元素,如xi,xj,以aij表示 xi 和 xj 对Z的影响之比。这里得到的A=(aij)nn称为两两比较的判断矩阵。第九章 层次分析 Saaty建议用19及其倒数做为标度来确定aij的值,19比例标度的含义: xi比xj强(重要)的程度 xi/ xj 相等 稍强 强 很强 绝对强 aij 1 2 3 4 5 6 7 8 9 19标度的理由:两两比较的心理习惯, 显然,判断矩
9、阵A的元素有如下特征: 第九章 层次分析 1 aij0 2aji=1/aij 3aii=1 我们称判断矩阵A为正互反矩阵。 第九章 层次分析 例如在例2中,准则层B对目标层作因素两两比较,并可建立下面判断矩阵: B1:B2为3 B1:B3为1 认为人才培养比另二项稍重要,另二项差不多相同重要。第九章 层次分析判断矩阵 B1 B2 B3 B1 1 3 1 A= B2 1/3 1 1/3 B3 1 3 1第九章 层次分析三、单一准则下元素相对排序权重计算及判断矩阵一致性检验:1、单一准则下元素排序: 求判断矩阵A的最大特征值max及标准化(归一化)的特征向量W。W的向量为同一层次中相应元素对于上一
10、层次中某个因素相对重要性的排序权重。有wi0,i, 。 n1ii1w第九章 层次分析 在构造判断矩阵时,各层元素间两两比较时,aij应有某种传递性质,即若甲比乙重要,乙比丙重要,合理地应有甲比丙更重要,在数值上表示为aijajk=aik 即 若xi与xj相比aij=3,xj与xk相比ajk=2,那么有传递性的判断应xj与xk相比,ajk=6 。第九章 层次分析 2、判断矩阵的一致性概念:、判断矩阵的一致性概念: 判断矩阵是各元素均为正数的矩阵这种正矩阵有下列重要性质。第九章 层次分析定理设n阶方阵A为正矩阵, max为A的最大模特征值,u =(u1,u2,un)T为max的相应特征向量。 、m
11、ax 0,ui 0,i =1,2,n、max是单特征根;(因此 u 除差一常数因子外是唯一的)、A的任何其它特征值,有max| |。第九章 层次分析定义:若正互反矩阵A满足aijajk=aik i ,j ,k =1,2,n 则称A为一致阵。一致阵的重要性质:设A是一致阵, 1A的转置亦是一致阵; aij=1/aji ,aij=1 ,i ,j=1,2,n; 由定义 aijajk=aik 则显然第九章 层次分析 2A的每一行均为任意指定的另一行的正数倍,从而A的秩为1。(即只有一个非零特征值,其余n-1个为0特征值); 考虑第行元素ai1,ai2,ain 对于第k行元素ak1,ak2,akn j=
12、1,2,n, aij=aikakj 即第行各元素分别为第k行各元素的aik倍。 第九章 层次分析 3A的最大特征根max= n,其余特征根皆为零; 4设u=(u1,u2,un)T是A对应max的特征向量,则aij=ui /uj i ,j =1, 2, , n 容易验证:对于n及向量u=(u1,u2,un)T 若aij=ui /uj ij 则 Au=nu (i, )又由定理1及性质2可知 max=n,u满足4 第九章 层次分析injijnjijnuuua 115若A为判断矩阵,那么A对应于max =n 的标准化(归一化)特征向量 u=(u1, u2,un)T 就是一组排序权向量。 (归一化 )由
13、性质4即知。1.2 进一步地有如下定理 定理2、n阶正互反矩阵A=(aij)nn是一致阵的充分必要条件为max=n nii1u1第九章 层次分析Proof : “必要性”即是上面性质3已证 “充分性”设A的最大特征值为max,相应特征向量u=(u1,un)T Au= max u 分量形式:对 i =1,2,n 由定理1知ui0 ,于是max= 注意aij=1,max-1= aij uj /ui n1jimaxjijuua nij1jijijuua/nij1j第九章 层次分析求和(把i=1,n的各式相加):nmax-n= aij uj /ui 注意 aji=1/aij 整理上式得:nmax-n=
14、 (aij uj /ui +1/ aij uj /ui )( )nij1jn1i nijn1i11第九章 层次分析*( )式末端=n2-n=n(n-1)注意:当x0时 x+(1/x)2当且仅当x=1时等号成立 。于是:aij ( uj /ui )+ (1/ aij)( uj /ui) 2( )式右端 2 = 2(n-1)+(n-2)+2+1=n(n-1)=左端 当且仅当 aij (uj /ui)=1时等号成立第九章 层次分析* 11n1inij* aij ( uj /ui )即aijajk=(ui /uj) (uj /uk)= uj /uk=ajk故A是一致阵。 由于客观事物的复杂性与人的认识
15、的多样性,我们得到的判断矩阵常常不具有传递性和一致性,但应该要求这些判断大体是一致的。 当判断矩阵过于偏离一致性时,它的可靠性值得怀疑,为此需对判断矩阵进行一致性检验。第九章 层次分析一致性检验步骤:、计算一致性指标C.I.=(max-n)/(n-1) (ConsisTeney Index)、查找相应的平均随机一致性指标R.I.(Random Index) 115阶正互反矩阵计算1000次得到的平均随 机一致性指标: 矩阵阶数 1 2 3 4 5 6 7 8 R.I. 0 0 0.52 0.89 1.12 1.26 1.36 1.41第九章 层次分析矩阵阶数 9 10 11 12 13 14
16、15 R.I. 1.46 1.49 1.52 1.54 1.56 1.58 1.59计算:R.I.=(max-n)/(n-1), max为m次判断矩阵max的平均值。max产生方法:取定阶数n,随机构造正互反矩阵=(ij)nn ,ij在1, 2, , 9, 1/2, 1/3, , 1/9这17个数中随机抽取,第九章 层次分析(只需取n(n-1)/2个,对角元为1,其余按正互反性得到)取充分大的子样计算所有的最大特征值,然后求平均即为max 。、计算一致性比率C.R. (consistency ratio) C.R.= C.I./R.I.当C.R.0.1时 认为判断矩阵的一致性是可接受的。当C.
17、R. 0.1时 应修正判断矩阵。第九章 层次分析例如 对前面矩阵 1 3 1 A= 1/3 1 1/3 1 3 1 计算出 max=3归一化向量u=(3/7,1/7,3/7)T C.I.=(max-3)/(3-1)=0C.R.=0 是一致阵。第九章 层次分析例: 1 2 5 A= 1/2 1 7 1/5 1/7 1 计算出 max=3.1189,u=(0.5415,0.3816,0.0761)T C.I.=(3.1189-3)/(3-1)=0.05945 查表得R.I.=0.52 C.R.=0.05945/0.52=0.11430.1,应修正判断矩阵第九章 层次分析 四、计算各层元素对目标层的
18、总排序权重: 层次总排序过程:计算同一层次所有因素对于最高层(总目标)相对重要性的排序权值。 从最高层到底层逐层进行: 设已算出第k-1层上nk-1个元素相对于总目标的排序为 w(k-1)=(w1(k-1),w2(k-1),w n (k-1)T第九章 层次分析K-1 第k层nk个元素对于第k-1层上第j个元素为准则的单排序向量 uj(k)=(u1j(k),u2j(k),un j(k)T j=1,2,nk-1 其中不受第j个元素支配的元素权重取零,于是可得到nknk-1阶矩阵 u11(k) u12(k) u1n (k) U(k)= u21(k) u22(k) u2 n (k) un 1(k) u
19、n 2(k) un n (k)第九章 层次分析kkkkk-1k-1k-1第k层上各元素对总目标的总排序w(k)为: w(k)=(w1(k),w2(k),wn (k)T w(k)=U(k)w(k-1)分量形式:wi(k)= uij(k)wj(k-1) i=1,2,n于是可得到公式:w(k)=U(k)U(k-1) U(3)w(2)w(2)为第二层上元素对目标的排序(即是单层排序)第九章 层次分析k11knj 各层总排序的一致性检验: 由高层向下,逐层进行检验,设第k层 中某些因素对k-1层第j个元素单排序的 一致性指标为C.I.j(k),平均随机一致性指标为R.I.j(k),(k层中与k-1层的第
20、j个元素无关时,不必考虑),那么第k层的总排序的一致性比率为: C.R.(k)= / )()1(.kjn1jkjIRwk 第九章 层次分析)()1(.kjn1jkjICwk 当C.R.(k)。幂法一致产生,使得,其它0 第九章 层次分析 幂法是处理这类矩阵求最大特征值及特征向量的一个简单而有效的方法。 幂法原理:设n阶矩阵A的特征值为1, 2,n有如下性质: 1 2 3 n 有n个线性无关的特征向量u1,u2,un x(0) Rn,则可表示为 x(0)= iui第九章 层次分析n1i利用迭代公式 x(k+1)=Ax(k) k=0,1,得到点列x(0),x(1),x(2),显然,x(k+1)=A
21、(k)x0 =Ak iui = iAkui = i ikui = ikiui+ i(i / 1)kui n1in1in2i第九章 层次分析n1i由于 i/ 1 0,且最大分量,且最大分量=max xi ,0,1 i n=y=(1/) xx=Ay,=max xi 1 i n - ?停停 特征值特征值 特征向量特征向量xNOyes此法当矩阵一致性较好时,收敛很快。在实用上常用下面的一些更为简单的方法(仅对近似一致性矩阵适应)。2、方根法:步骤: 、求Mi=( aij)1/n i=1,2,n 、标准化(归一化):Wi=Mi / Mj 、max=(1/n) (AW)i /Wi n1jn1jn1i第九章
22、 层次分析EX.EX. 1 3 1 M1= =1.4422 A= 1/3 1 1/3 M2= =0.4807 1 3 1 M3= =1.4422 w1=0.4286 归一化: w2=0.1428 w3=0.4286 Aw=(1.2856,0.4285,1.2856)Tmax=2.9999 313133/ 113/ 13131第九章 层次分析3、和积法:步骤:、求(每列归一化) bij=aij / akj i,j=1,2n、行求和Mi= bij i= 1,2,n再归一化:Wi=Mi / Mj i= 1,2,n 、max=(1/n) (AW)i/Wi n1kn1jn1jn1i第九章 层次分析例:
23、1 3 1 3/7 3/7 3/7 M1=9/7 A= 1/3 1 1/3 B= 1/7 1/7 1/7 M2=3/7 1 3 1 3/7 3/7 3/7 M3=9/7 Mj=3 w2=1/7 Aw=(9/7,3/7,9/7)T w3=3/7 max=3显然,当A是一致阵时, max=n,对归一化的waij=wi/wj 第九章 层次分析 w1=3/731j 方根法:Mi=( aij)1/n=Wi/S S=( Wj)1/n i=1,2,n 归一化后w即为(w1,w2,wn)T max=(1/n) (Aw)i / wi (Aw=nw) =n2/n =nn1jn1i第九章 层次分析n1j和积法: a
24、kj= wk/wjbij=aij/ akj=wi / wk Mi= bij=(nwi)/ wk 归一化后w即为 (w1,w2,wn)T同理max=n当A近似一致阵时,这些量是近似的。例: 1 2 5 A= 1/2 1 3 1/5 1/3 1 n1j第九章 层次分析nk 1nk 1nk 1nk 1nk 1用幂法:取x(0)=(1,0,0)T k x1 x2 x3 y1 y2 y3 0 1 0 0 1 1 0 0 1 1 0.5 0.2 1 1 0.5 0.2 2 3 1.6 0.5667 3 1 .5333 .1889 3 3.0111 1.6 0.5667 3.0111 1 .5314 .18
25、82 4 3.0037 1.5959 0.5653 3.0037 1 .5313 .1882 5 3.0037 1.5959 0.5653 3.0037第九章 层次分析max=3.0037 C.I.= (max-3)/(3-1)=0.00185C.R.=C.I./R.I.=0.00185/0.52=0.0036满足一致阵要求 。u=(3.0037,1.5959,0.5653)T归一化得:w=(0.5816,0.3090,0.1094)T第九章 层次分析用方根法:、M1= = =2.1544 M2= =1.1447 M3= =0.4055、归一化:M1+ M2+ M3=3.7046 w1=2.1
26、544/3.7046=0.5815 w2=1.1447/3.7046=0.3090 w3=0.4055/3.07046=0.1095 w=(0.5815,0.3090,0.1095)T第九章 层次分析5*2*110332/33315/1、Aw=(1.7470,0.9283,0.3388)T 1 1.7470 0.9283 0.3288 3 0.5815 0.3090 0.1095 3.0037第九章 层次分析max=+=用和积法: 、 1 2 5 0.5882 0.6 0.5556 A= 1/2 1 3 B= 0.2941 0.3 0.3333 1/5 1/3 1 0.1177 0.1 0.1
27、111、行求和M=(1.7438,0.9274,0.3288)T M1+M2+M3=3 归一化:w=(0.5813,0.3091,0.1096)T第九章 层次分析列归一化、Aw=(1.7475,0.9286,0.3289)T 1 1.7475 0.9286 0.3289 3 0.5813 0.3091 0.1096 3.0038第九章 层次分析max=+=二、残缺判断与群组决策: 1、残缺判断及处理方法: 应用AHP进行决策时,每个准则应有一个判断矩阵,需进行 n(n-1)/2 次两两比较 (判断矩阵的上或下三角)。 当层次很多,因素复杂时,判断量很大,可能出现某个参与决策的专家对某些判断缺少
28、把握,或不想发表意见,使判断矩阵残缺。第九章 层次分析 可接受的残缺判断矩阵 若任一残缺元素都可通过已给出的元素间接获得的残缺判断矩阵。 根据一致性的条件:间接获得的元素指,若aij缺少可由aij=aikakj或更一般地aij=aik ak k ak k ak j得到。第九章 层次分析11232s可接受的残缺矩阵的排序向量计算 常用的有特征根方法,对数最小二乘法及最小偏差法等。 特征根法:设A对应max的特征向量w =(w1,w2,wn)T 由一致性条件知 aij = wi /wj,特征根法即把缺少的的元素用wi /wj来替代。第九章 层次分析设原判断矩阵A=(aij)nn构造辅助矩阵 C =
29、 (cij)nn使 cij= aij , aij0 wi /wj , aij=0例:设 1 2 0 A A= 12 1 2 是可接受的残缺矩阵 0 12 1第九章 层次分析辅助矩阵 1 2 w1w3 C= 1/2 1 2 w1w3 1/2 1解特征根问题:cw = maxw展开:左=(2w1+2w2,1/2w1+w2+2w3,1/2w2+2w3)T = max(w1, w2, w3)T解得:解得: maxmax=3 =3 w =(0.5714=(0.5714,0.28570.2857,,0.1429,0.1429)T T第九章 层次分析可以看出:C 的特征值问题等价于 2 2 0 = 1/2
30、1 2 0 1/2 2 的特征值问题(AwAw= maxw与Cw= maxw相同)第九章 层次分析故只需求下列矩阵的特征值及特征向量 = (aij )nn aij 当当aij0,ij aij = 0 当当aij=0 mi+1 当当ij时,时,mi为第为第i行中残缺元素的个数行中残缺元素的个数求解 w = maxw 可得不完整信息下的排列向量第九章 层次分析 (3)一致性检验: max-n C.I.= (n-1)-( ) 当C.R.=C.I./R.I.0.1时 认为有满意的 一致性。 nii/nm1第九章 层次分析 2. 2.群组决策:群组决策: 为使决策科学化、民主化,一个复杂系统通常是由多个
31、决策者(专家)或决策部门参与决策的。群组决策问题是指采取一定的方法以使决策者的决策综合成一个较合理的结果的过程。 第九章 层次分析应做好如下工作:重视并做好专家咨询工作; 、合理选择咨询对象;(专长及熟悉的领域) 、创造适合于咨询工作的良好环境;(介绍AHP方法,提供信息,独立思考) 、正确的咨询方法;(通过咨询确定递阶层次结构,设计好表格) 第九章 层次分析、及时分析专家咨询信息,必要时要进行反馈及多轮次咨询群组决策综合分析方法:两类方法: 、将各专家的判断矩阵综合,得到综合判断矩阵,再计算排序。第九章 层次分析 、先求各专家判断矩阵的排序向量,再综合成群组排序向量。 设S个专家的判断矩阵:
32、 Ak=(aij (k) k=1,2, S 分别求出它们各自的排序向量 wk=(w1(k) ,w2(k) wn(k)T (k)第九章 层次分析再记平均综合向量为w =(w1,w2, ,wn)T方法1.加权几何平均综合排序向量法: 计算 wj= wj / (归一化) 其中,k k 0且 k为第k个决策者的权重. niiw1kkk11j=1,2j=1,2.n.n第九章 层次分析ssjjjjwwww )()()()()2()1(21 对可采用性的考察:计算wj的标准差:j=其相应于新的总体判断矩阵A=(aij) (aij=wi /wj)的总体标准差:第九章 层次分析(K)2 2 skjj-wws1
33、)()1/(1ij=个体标准差:(k)= 当总体标准差满足要求时,这组群组判断可采用,当个体标准差(k)满足要求时,认为第k个决策者的决策可通过,否则将信息反馈给有关专家,供修改时参考。第九章 层次分析 skijij-aas1 )()1/(1(K)2 2 njjkj-wwn12)( )()1/(1方法2.加权算术平均综合向量法:计算 W=1Wj(1)+ 2Wj(2)+sWj(s) k0,可类似地根据,式判断可采用性。第九章 层次分析 s1kk1 9.3应用举例一、一、某工厂有一笔企业留成利润,要决 定如何使用。 供选择方案: 作奖金,集体福利设施,引入设备技术 建立如下层次分析模型:第九章 层
34、次分析目标层:准则层C:方案层P:合理使用留成利润 A改善职工生活条件C3提高技术水平C2调动职工积极性C1引进设备技术P3福利P2奖金P1第九章 层次分析A-C判断矩阵: A C1 C2 C3 w(2) C1 1 1/5 1/3 0.105 C2 5 1 3 0.637 C3 3 1/3 1 0.258 max=3.038 归一化特征向量w(2) C.I.=0.019 C.R.=0.032760.1 满意的一致性第九章 层次分析C1-P: C1 P1 P2 U1(3) P1 1 1/3 0.75 P2 3 1 0.75 max=2 C.I.=0 第九章 层次分析C2-P: C2 P2 P3
35、U2(3) P2 1 1/5 0.167 P3 5 1 0.833max=2 C.I.=0第九章 层次分析C3-P: C3 P1 P3 U3(3) P1 1 2 0.667 P2 1/2 1 0.333 max=2 C.I.=0 第九章 层次分析 0.25 0 0.667 U(3)= 0.75 0.167 0.333 0 0.833 0w(3)=U(3)w(2)=(0.198,0.27,0.531)T得到P3优于P1又优于P2,从分配上可以用53.1%来引进新设备,新技术;用19.8%来发奖金;用29.1%来改善福利。第九章 层次分析二、二、层次分析法对于下面几种情况的优化问题特别适用:问题中
36、除可计量的量外,还存在不可计量的量时,可用AHP通过对不可计量的量与可计量的量的相对比较,而获得相对的量测;当优化问题的结构难以事先确定,而在很大程度上取决于决策者的经验时;第九章 层次分析各变量不独立,有内部相关性时;目标与约束,约束与约束之间紧密联系时;多目标问题;第九章 层次分析 在用AHP法解决优化问题时,常用的有两种方式:当模型中涉及不可计量的量时,用AHP法的比例标度来确定目标函数,约束函数的权重(系数)直接采用AHP模型 AHP法有广泛的应用前景,可以用来决定其它方面的一些问题。下面举一个解决优化问题的例子。第九章 层次分析例:最佳食品搭配问题! 假设某人有3种食品可供选择:肉,
37、面包,蔬菜它们所含营养成分及单价如下表:食品 维生素A 维生素B2 热量 单价 搭配量 (国际 (毫克/克) (千卡/克) (元/克) 单价/克) 肉 0.3527 0.0021 2.86 0.0055 X1面包 0 0.0006 2.76 0.0012 X2蔬菜 25.0 0.002 0.25 0.0014 X3第九章 层次分析该人体重55公斤,每天对各种营养的最小需求为:维生素A:7500 国际单位维生素B2:1.6338 毫克热量:2050 千卡问题:应如何搭配食品?(自然的想法是:使在保证营养的情况下支出最小)第九章 层次分析 容易建立如下线性规划模型:min Z=0.0055 x1+
38、0.0012 x2+0.0014 x3s.t. 0.3527 x1+25.0 x37500 0.0021 x1+0.0006 x2+0.002x31.6338 2.86 x1+2.76 x2+0.25 x32050 x1,x2,x30利用单纯形法可得解利用单纯形法可得解 x*=(0, 689.44, 610.67)T z*1.67第九章 层次分析 即,不吃肉,面包689.44克,蔬菜610.67克,每日支出1.67元。显然这个最优方案是行不通的,它没有考虑本人对食品的偏好。我们可根据偏好加约束: x1140, x2450, x3不限 得到线性规划解: x*=(245.44, 450.00 42
39、4.19)T Z*=2.48元 第九章 层次分析其次,在这里各营养成分被看成同样重要,起决定因素的是支出。但实际上,营养价值与支出都需考虑,只是地位(权重)不同。这样无法建立目标函数。下面用层次分析法来处理问题:层次结构: 第九章 层次分析每日需求 R支出 C 营养 N维生素 A 维生素 B2维生素 Q肉 me面包 br蔬菜 ve第九章 层次分析对于一个中等收入的人,满足营养要求比支出更重要。于是: R N C w(2) N 1 3 0.75 C 1/3 1 0.25 max=2 C.I.=0第九章 层次分析 N A B2 Q w1(3) A 1 1 2 0.4 B2 1 1 2 0.4 Q
40、1/2 1/2 1 0.2 max=3 C.I.=0第九章 层次分析 0.4 0w(3) = 0.4 0 0.25 0.2 0 0.25 =(0.3, 0.3, 0.15, 0.25)T 0 1最底层(方案层)对准则层的单排列权重,只需对题目给的数据归一化即可。 由于要支出最小价格倒数,价格倒数归一:( 181.818,833.333,714.286 )T 于是得到第九章 层次分析 A B2 Q C(价格) me 0.0139 0.4468 0.4872 0.1057U(4) br 0.0000 0.1277 0.4702 0.4819 ve 0.9861 0.4255 0.0426 0.43
41、10合成权重w(4) = U(4)w(3) = (0.24, 0.23, 0.53)T第九章 层次分析设 x1=0.24k, x2=0.23k, x3=0.53k则 min Z = 0.002338k s.t. 13.3346k 7500 0.0017k 1.6338 1.4537k 2050 k 0 解得:k = 1410.20 变为第九章 层次分析x1=338.45克,x2=324.35克,x3=749.41克Z=3.30元 满足条件此时各营养成分含量如下:维生素A:18804.52国际单位维生素B2:2.400毫克热量Q:2050.01千卡若认为总支出太大,可适当降低第二层中营养的权重
42、。第九章 层次分析若改为 R N C w(2) N 1 1 0.5 C 1 1 0.5 max=2 C.I.=0 第九章 层次分析 其余不变: 0.4 0 w(3) = 0.4 0 0.5 = (0.2,0.2,0.1, 0.5)T 0.2, 0.2 0.5 0 1 .0319 .4468 .4872 .1081 0.2 w(4) =U(4)w(3) = .0000 .1277 .4702 .4819 0.2 .9861 .4255 .0426 .4130 0.1 =(0.193, 0.314, 0.493) T第九章 层次分析 类似上面可解得:设 x1=0.193k, x2=0.314k,
43、x3=0.493k min Z=0.0021285k 则 s.t. 12.3931k7500 0.0016k1.6338 1.54187k2050 0.193k140, 0.314k450, k0、变为第九章 层次分析得解 k=1329.56于是 x1=256.61克, x2=419.48克, x3=655.47克, Z=2.83元即每日肉256.61克,面包419.48克,蔬菜655.47克,总支出2.83元第九章 层次分析各营养成分含量如下: 维生素A:16479.33国际单位 维生素B2:2.100毫克 热量 Q:2050.01千卡第九章 层次分析 下图为日支出对于营养权重变化的灵敏度曲线。它们基本上位于线性规划、的可行解目标值(支出为1.673.80元)范围内。10515支出(元)1.671, 3.800.50.750.251第九章 层次分析 灵敏度分析: 判断矩阵经常会受到扰动:矩阵的不一致性可看做是一种扰动;决策者做出的决策差异总是存在的(同一决策者,不同时刻,地点对同一问题判断存在差异);计算中亦带来误差。第九章 层次分析 灵敏度分析即指由于求解过程中引入误差,而导致问题解的变化。若变化大,则结果一般是靠不住的。 灵敏度分析方法:矩阵摄动法 第九章 层次分析