《2022年数学建模中的综合评价方法》由会员分享,可在线阅读,更多相关《2022年数学建模中的综合评价方法(20页珍藏版)》请在金锄头文库上搜索。
1、综合评价评价是人类社会中一项经常性的、极重要的认识活动,是决策中的基础性工作。在实际问题的解决过程中,经常遇到有关综合评价问题,如医疗质量的综合评价问题和环境质量的综合评价等。它是根据一个复杂系统同时受到多种因素影响的特点,在综合考察多个有关因素时,依据多个有关指标对复杂系统进行总评价的方法;综合评价的要点:(1)有多个评价指标,这些指标是可测量的或可量化的; (2)有一个或多个评价对象,这些对象可以是人、单位、方案、标书科研成果等;(3)根据多指标信息计算一个综合指标,把多维空间问题简化为一维空间问题中解决,可以依据综合指标值大小对评价对象优劣程度进行排序。综合评价的一般步骤1根据评价目的选
2、择恰当的评价指标,这些指标具有很好的代表性、区别性强,而且往往可以测量,筛选评价指标主要依据专业知识,即根据有关的专业理论和实践,来分析各评价指标对结果的影响,挑选那些代表性、确定性好,有一定区别能力又互相独立的指标组成评价指标体系。2根据评价目的,确定诸评价指标在对某事物评价中的相对重要性,或各指标的权重;3合理确定各单个指标的评价等级及其界限;4根据评价目的,数据特征,选择适当的综合评价方法,并根据已掌握的历史资料,建立综合评价模型;5确定多指标综合评价的等级数量界限,在对同类事物综合评价的应用实践中,对选用的评价模型进行考察,并不断修改补充,使之具有一定的科学性、实用性与先进性,然后推广
3、应用。目前,综合评价有许多不同的方法,如综合指数法、TOPSIS 法、层次分析法、RSR 法、模糊综合评价法、灰色系统法等,这些方法各具特色,各有利弊,由于受多方面因素影响,怎样使评价法更为准确和科学,是人们不断研究的课题。下面仅介绍综合评价的TOPSIS 法、 RSR 法和层次分析法的基本原理及简单的应用。8.1 TOPSIS 法(逼近理想解排序法)Topsis 法是系统工程中有限方案多目标决策分析的一种常用方法。是基于归一化后的原始数据矩阵,找出有限方案中的最优方案和最劣方案(分别用最优向量和最劣向量表示),然后分别计算诸评价对象与最优方案和最劣方案的距离,获得各评价对象与最优方案的相对接
4、近程度,以此作为评价优劣的依据。8.1.1 基本原理TOPSIS 法是 Technique for Order Preference by Similarity to Ideal Solution的缩写, 即逼近于理想解的技术,它是一种多目标决策方法。方法的基本思路是定义决策问题的理想解和负理想解,然后在可行方案中找到一个方案,使其距理想解的距离最近,而距负理想解的距离最远。理想解一般是设想最好的方案,它所对应的各个属性至少达到各个方案中的最好值;负理想解是假定最坏的方案,其对应的各个属性至少不优于各个方案中的最劣值。方案排队的决策规则,是把实际可行解和理想解与负理想解作比较,若某个可行解最靠
5、近理想解,同时又最远离负理想解,则此解是方案集的满意解。8.1.2 距离的测度采用相对接近测度。设决策问题有m 个目标jf(m,j21),n 个可行解)(imiiiZ,Z,ZZ21(n,i21);并设该问题的规范化加权目标的理想解是Z*,其中)(mZ,Z,ZZ21,那么用欧几里得范数作为距离的测度,则从任意可行解iZ到Z的距离为:mjjiji)ZZ(S12i=1 ,n,(8.1) 式中, Zij为第 j 个目标对第i 个方案(解)的规范化加权值。同理,设Z=TmZ,Z,Z)(21为问题的规范化加权目标的负理想解,则任意可行解iZ到负理想解Z之间的距离为:mjjijiZZS12)(i=1 , n
6、 ,(8.2) 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 1 页,共 20 页那么,某一可行解对于理想解的相对接近度定义为:iiiiSSSC0 Ci1 ,i=1,n,(8.3) 于是,若iZ是理想解,则相应的Ci =1;若iZ是负理想解,则相应的Ci =0。iZ愈靠近理想解,Ci愈接近于 1;反之,愈接近负理想解,Ci愈接近于0。那么,可以对Ci进行排队,以求出满意解。8.1.3 TOPSIS 法计算步骤第一步:设某一决策问题,其决策矩阵为A. 由 A 可以构成规范化的决策矩阵Z ,其元素为Zij,且有niijijijffZ12mjni,2
7、,1;,2, 1(8.4) 式中, fij由决策矩阵给出。(8.5) 第二步:构造规范化的加权决策矩阵Z,其元素ZijZij =Wj Ziji=1, ,n; j =1,m(8.6) Wj为第 j 个目标的权。第三步:确定理想解和负理想解。如果决策矩阵Z 中元素 Zij值越大表示方案越好,则m,jZmaxZ,Z,ZZijim2121)((8.7) m,jZminZ,Z,ZZijim2121)((8.8) 第四步:计算每个方案到理想点的距离Si和到负理想点的距离S -i。第五步:按式(8.3)计算 Ci,并按每个方案的相对接近度Ci的大小排序,找出满意解。多目标综合评价排序的方法较多,各有其应用价
8、值。在诸多的评价方法中,TOPSIS 法对原始数据的信息利用最为充分,其结果能精确的反映各评价方案之间的差距,TOPSIS 对数据分布及样本含量,指标多少没有严格的限制,数据计算亦简单易行。不仅适合小样本资料,也适用于多评价对象、多指标的大样本资料。利用TOPSIS 法进行综合评价,可得出良好的可比性评价排序结果。8.1.4应用实例1、TOPSIS 法在医疗质量综合评价中的应用试根据表8.1 数据,采用Topsis 法对某市人民医院19951997 年的医疗质量进行综合评价。表 8.1 某市人民医院19951997 年的医疗质量年度床位周转次数床位周转率(%)平均住院日出入院诊断符合率(%)手
9、术前后诊断符合率(%)三日确诊率( %)治愈好转率(%)病死率(%)危重病人抢救成功率(%)院内感染率(%)1995 20.97 113.81 18.73 99.42 99.80 97.28 96.08 2.57 94.53 4.60 1996 21.41 116.12 18.39 99.32 99.14 97.00 95.65 2.72 95.32 5.99 1997 19.13 102.85 17.44 99.49 99.11 96.20 96.50 2.02 96.22 4.79 在原始数据指标中,平均住院日、病死率、院内感染率三个指标的数值越低越好,这三个指标称为低优指标;其它指标数值
10、越高越好,称为高优指标。是低优指标的可转化为高优指标,其方法为是绝对数低优指标x可使用倒数法(x100) ,是相对数低优指标x,可使用差值法(x1) 。这里,平均住院日采用倒数转化,病死率、院内感染率采用差值转化。转化后数据见表8.2。表 8.2 转化指标值年度床位周床位平均出入院手术前三日治愈病死率危重病院内精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 2 页,共 20 页转次数周转率(%)住院日诊断符合率(%)后诊断符合率(%)确诊率( %)好转率(%)(%) 人抢救成功率(%)感染率(%)1995 20.97 113.81 5.34 99.
11、42 99.80 97.28 96.08 97.43 94.53 95.40 1996 21.41 116.12 5.44 99.32 99.14 97.00 95.65 97.28 95.32 94.01 1997 19.13 102.85 5.73 99.49 99.11 96.20 96.50 97.98 96.22 95.21 根据表 8.2 数据,利用公式(8.4)进行归一化处理,得归一化矩阵值,如表8.3。niijijijffZ12(8.9)例如计算1995 年床位周转次数归一化值,由公式(8.9)得:2221113.1941.2197.2097.20Z509.0其余归一化数值以此
12、类推。表 8.3 归一化矩阵值年度床位周转次数床位周转率平均住院日出入院诊断符合率手术前后诊断符合率三日确诊率治愈好转率病死率危重病人抢救成功率院内感染率1995 0.590 0.592 0.560 0.577 0.580 0.580 0.577 0.577 0.572 0.581 1996 0.602 0.604 0.570 0.577 0.576 0.578 0.575 0.576 0.577 0.572 1997 0.538 0.535 0.601 0.578 0.576 0.574 0.580 0.580 0.583 0.579 由式( 8.7)和式( 8.8)得最优方案和最劣方案:)
13、(mZZZZ,210.581,0.583,0.580,0.580,0.580,0.580,0.578,0.601,0.604,0.602(8.10)(mZZZZ,210.572,0.572,0.576,0.575,0.574,0.576,0.577,0.560,0.535,0.538( 8.11)由式( 8.10) 、 (8.11)和式( 8.1) 、 (8.2)计算各年度D和D,见表 8.4。例如计算1997 年S和S:222579058105350604053806020.S(8.12)094.0222579057205350535053805380.S(8.13)044.0其余各年依次类
14、推。由式( 8.3)计算各年度iC,见表 8.4。例如计算1997 年iC:319.0044.0094.0044.0iC(8.14)其余各年以次类推。表 8.4 不同年度指标值与最优值的相对接近程度及排序结果年份DDiC排序结果1995 0.045 0.078 0.634 2 1996 0.034 0.095 0.736 1 1997 0.094 0.044 0.319 3 由表 8.4 的排序结果可知1996 年医疗质量最好。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页,共 20 页2 TOPSIS 法在环境质量综合评价中的应用实例在环
15、境质量评价中, 把每个样品的监测值和每级的标准值,分别看作TOPSIS 法的决策方案, 由 TOPSIS法可以得到每个样品和每级标准值的Ci值,对 Ci值大小排序, 便可以得到每个样品的综合质量及不同样品间进行综合质量优劣比较。表 8.5 列出所选的参评要素和所确定的评判等级及其代表值表 8.5 某海湾沿岸海水侵染程度分级表参评要素分级 级(无或很轻侵染) 级(轻度侵染 ) 级(较严重侵染) 级(严重侵染 ) 氯离子(mg/l) 100 400 800 2 200 矿化度(mg/l) 500 1 500 2 500 3 500 溴离子(mg/l) 0.25 1.25 2.50 9.00 rHC
16、O3/rCl 1.00 0.31 0.14 0.02 纳 吸 附 比1.40 2.60 4.50 15.50 测得 111和 112水样的各参评要素值如表8.6。表 8.6111和 112水样监测值样品号要素氯离子(mg/l) 矿化度(mg/l) 溴离子(mg/l) rHCO3/rCl 纳 吸 附 比111134.71 542.15 0 0.882 1.576 112152.44 721.18 0.20 1.267 1.366 取海水侵染级标准值和111及 112样品监测值构成TOPSIS 法中的决策矩阵A,那么由式 (8.4)算出 A 的规范化矩阵Z因在制定海水侵染分级标准时,各因子的重要性
17、已隐含在分级标准值中,因此,本文由标准值来确定权重,其计算式如下:niiniiniiISSISSW1)1()1()/(/(8.15) 式中, Wi为i因子的权重;n为标准分级数,在本例中4n;)1(niS为i因子的第1n级标准值;ISi为i因子的第I级标准值。式(8.15)适用于低优指标型因子,在本例中如氯离子、矿化度、溴离子、纳吸附比等,权重计算时用S精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4 页,共 20 页/SI;而对高优指标型因子如rHCO3/rHCl,计算时用S/S。通过计算得权重向量WT=0.1980.119 90.239 80
18、.371 70.076 7由式 (8.6)得加权后的规范化矩阵Z 为由式 (8.7),式 (8.8)得Z=0.1768 0.0898 0.2288 0.2520 0.0719Z-=0.0081 0.0128 0 0.0041 0.0064最后,由式 (8.1),式 (8.2)和式 (8.3)计算iS,iS和 Ci值(表 8.7)。表 8.7 iS,iS和 Ci值表111112iS0.053 5 0.196 2 0.245 5 0.390 5 0.076 7 0.008 7 iS0.355 0 0.263 1 0.209 3 0 0.345 3 0.384 7 iC0.869 0 0.572 8
19、 0.460 2 0 0.818 2 0.977 9 把 Ci排序得C112CC111CCC于是可知: 112样品综合质量优于111样品综合质量,112样品质量优于I 级标准最低界限值,为 I 级; 111样品质量介于I 级和级最低界限值之间,属于级。因此,111 样品为轻度侵染,112样品为无或很轻污染。由监测值也可以知道:111有 4 个因子达到级,1 个因子达到I 级;112有 2个因子达到级(接近 I 级 ),3 个因子达到I 级。因此,本方法评价结果符合客观实际。8.1.5 结论TOPSIS 法是一种多目标决策方法,适用于处理多目标决策问题。本文提出TOPSIS 法应用于环境质量综合
20、评价中,取得较好的效果,与其他方法比较,具有以下优点:1、与环境标准巧妙结合起来,不仅能确定各评价对象所属的级别,还能进行不同评价对象间质量的优劣比较。2、 TOPSIS 法原理简单,能同时进行多个对象评价,计算快捷,结果分辨率高、评价客观,具有较好的合理性和适用性,实用价值较高。TOPSIS 法的缺点是*iC只能反映各评价对象内部的相对接近度,并不能反映与理想的最优方案的相对接近程度。8.2 秩和比法秩和比法是我国统计学家田凤调教授于1988 年提出的一种新的综合评价方法,它是利用秩和比RSR(Rank-sum ratio )进行统计分析的一种方法,该法在医疗卫生等领域的多指标综合评价、统计
21、预测预报、统计质量控制等方面已得到广泛的应用。秩和比是一个内涵较为丰富的综合性指标,它是指行(或列)秩次的平均值,是一个非参数统计量,具有01 连续变量的特征,近年来秩和比统计方法不断完善和充实。8.2.1 分析原理及步骤1、分析原理秩和比是一种将多项指标综合成一个具有01 连续变量特征的统计量,也可看成0100 的计分。多用精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 5 页,共 20 页于现成统计资料的再分析。不论所分析的问题是什么,计算的 RSR 越大越好。 为此, 在编秩时要区分高优指标和低优指标,有时还要引进不分高低的情况。例如,评价预
22、期寿命、受检率、合格率等可视为高优指标;发病率、病死率、超标率为低优指标。在疗效评价中,不变率、微效率等可看作不分高低的指标。指标值相同时应编以平均秩次。秩和比综合评价法基本原理是在一个n 行 m 列矩阵中,通过秩转换,获得无量纲统计量RSR;在此基础上,运用参数统计分析的概念与方法,研究RSR 的分布;以RSR 值对评价对象的优劣直接排序或分档排序,从而对评价对象作出综合评价。2、分析步骤 编秩:将 n 个评价对象的m 个评价指标列成n 行 m 列的原始数据表。编出每个指标各评价对象的秩,其中高优指标从小到大编秩,低优指标从大到小编秩,同一指标数据相同者编平均秩。 计算秩和比( RSR):根
23、据公式mjijinmRRSR1计算,式中i=1,2, n;ijR为第 i 行第 j列元素的秩,最小RSR=1/n,最大 RSR=1。当各评价指标的权重不同时,计算加权秩和比(WRSR),其计算公式为mjijjiRwnwRSR11,Wj 为第 j 个评价指标的权重,Wj=1。通过秩和比 (RSR)值的大小,就可对评价对象进行综合排序,这种利用 RSR 综合指标进行排序的方法称为直接排序。但是在通常情况下还需要对评价对象进行分档,特别是当评价对象很多时,如几十个或几百个评价对象,这时更需要进行分档排序,由此应首先找出RSR 的分布。 计算概率单位(Probit):将RSR(或WRSR)值由小到大排
24、成一列,值相同的作为一组,编制RSR(或 WRSR )频率分布表,列出各组频数f,计算各组累计频数f;确定各组RSR(或 WRSR)的秩次范围R 和平均秩次R;计算累计频率p=AR/n ;将百分率p 转换为概率单位Probit, Probit 为百分率p对应的标准正态离差u 加 5。 计算直线回归方程:以累计频率所对应的概率单位Probit 为自变量,以RSR(或 WRSR )值为因变量,计算直线回归方程,即RSR(WRSR)=a+bProbit。 分档排序:根据标准正态离差分档,分档数目可根据试算结果灵活掌握,最佳分档应该是各档方差一致,相差具有显著性,一般分3-5 档,下面是常用分档数对应
25、的百分位数及概率单位见表8.8。表 8.8 常用分档数及对应概率单位依据各分档情况下概率单位Probit 值,按照回归方程推算所对应的RSR(或 WRSR)估计值对评价对象进行分档排序。具体的分档数根据实际情况决定。8.2.2秩和比法在对某病区护士综合评价中的应用实例某医院对护士考核有4 个指标,它们分别是:业务考核成绩(1x) 、操作考核结果(2x) 、科内测评(3x)和工作量考核(4x) ;下表 8.9 是某病区8 名护士的考核结果:表 8.9 某病区 8 名护士的考核结果待评对象( n)1x2x3x4x精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - -
26、 -第 6 页,共 20 页护士甲86 优- 100 233.9 护士乙92 良98.2 192.9 护士丙88 良99.1 311.1 护士丁72 良95.5 274.9 护士戊70 优97.3 263.6 护士己94 优100 182.3 护士庚84 良91.97 220.6 护士辛50 良91.97 182.0 利用秩和比综合评价法对其进行综合评价。根据秩和比综合评价法的评价步骤,第一步分别对要评价的各项指标进行编秩,由于对护士考核的4 个指标都是高优指标,所以对要评价的各项指标进行编秩如表8.10: 表 8.10 评价的各项指标编秩待评对象( n)1x2x3x4x护士甲86(5)优-(
27、 6)100 ( 7.5)233.9(5)护士乙92(7)良( 3)98.2(5)192.9(3)护士丙88(6)良( 3)99.1(6)311.1(8)护士丁72(3)良( 3)95.5(3)274.9(7)护士戊70(2)优( 7.5)97.3(4)263.6( 6)护士己94(8)优( 7.5)100(7.5)182.3( 2)护士庚84(4)良( 3)91.97(1.5)220.6( 4)护士辛50(1)良( 3)91.97(1.5)182.0( 1)第二步,计算各指标的秩和比(RSR)mjijinmRRSR1其中 m 为指标个数, n 为分组数,ijR为各指标的秩次,RSR 值即为多
28、指标的平均秩次,其值越大越优。各护士 4 项护理考核指标编秩及RSR 值如表 8.11 表 8.11 各护士 4 项护理考核指标编秩及RSR 值待评对象( n)1x2x3x4xRSR 护士甲86(5)优-(6)100 (7.5)233.9(5)0.7344 护士乙92(7)良( 3)98.2(5)192.9(3)0.5313 护士丙88(6)良( 3)99.1(6)311.1(8)0.7188 护士丁72(3)良( 3)95.5(3)274.9(7)0.5000 护士戊70(2)优( 7.5)97.3(4)263.6(6)0.6094 护士己94(8)优( 7.5)100(7.5)182.3(
29、2)0.7813 护士庚84(4)良( 3)91.97(1.5)220.6(4)0.3906 护士辛50(1)良( 3)91.97(1.5)182.0(1)0.2031 如果将 8 名护士进行排序,则可根据 8 名护士的秩和比 (RSR),按由大到小排列就可得到8 名护士由好到差的所有排序;如果要将8 名护士分成几档,则还需继续进行下列工作。第三步,确定RSR 的分布将各指标的RSR 值由小到大进行排列,计算向下累计频率,查百分数与概率单位对照表,求其所对应的概率单位值,见表 8.12 表 8.12 概率单位值RSR f 累积频数R%100)/(nRY 0.2031 1 1 1 12.5 3.
30、8197 0.3906 1 2 2 25.5 4.3255 0.5000 1 3 3 37.5 4.6814 0.5313 1 4 4 50.5 5.0000 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 7 页,共 20 页0.6094 1 5 5 62.5 5.3186 0.7188 1 6 6 75.0 5.6745 0.7344 1 7 7 87.5 6.1503 0.7813 1 8 8 )19.966.8663 其中数据)19 .96是利用%100411)(n估计的。第四步,求回归方程:RSR=A+BY 将概率单位值Y 作为自变量,秩
31、和比RSR 作为因变量 ,经相关和回归分析,因变量RSR 与自变量概率单位值 Y 具有线性相关 (r=0.9528) ,线性回归方程为: RSR=0.1877Y-0.4232, 经 F 检验,F=59.078,P=0.0002,这说明所求线性回归方程具有统计意义。第五步,将8 名护士进行分档,分多少档根据评价对象具体要求确定,如果将8 名护士分为优良差三档,根据统计学家田凤调教授提供的一个分档标准,分档如下表8.13:表 8.13 8 名护士分档表等级Y RRS分档差4 以下0.05 ,方差一致。方差分析: F=43.2921,P中差,均具有显著性意义。在本法编秩中,对于高优指标,最小的指标值
32、编为1,最大的指标值编为n(此点与 RSR 法相同 ),但其余指标值由小到大分别编为1 与 n 之间的线性递增的非整秩次。所编秩次与原指标值之间存在定量的线性对应关系,即原指标值被定量地转换为秩次,而不是简单的等级化,从而避免了秩次化后原指标值定量信息的损失。低优指标的编秩方法相同,但大小方向相反。与 RSR 法比较,非整秩次RSR 法的不足是不能直观地列出秩次,而需经过计算得出,故运算比RSR法多一步。但所增加一点运算换取更准确、更客观的评价结果是值得的。8.3 层次分析法;人们在实际问题中常常会遇到各种各样的决策问题,如旅游地的选取问题,旅游者初次筛选几处旅游地,但每个旅游地的景色、所需费
33、用、居住条件、饮食条件交通等各不相同,根据个人的条件和爱好等如何确定旅游地。再例如,某人准备选购一台电冰箱,他对市场上的6 种不同类型的电冰箱进行了解后,在决定买那一款式是,往往不是直接进行比较,因为存在许多不可比的因素,而是选取一些中间指标进行考察。例如电冰箱的容量、制冷级别、价格、型式、耗电量、外界信誉、售后服务等。然后再考虑各种型号冰箱在上述各中间标准下的优劣排序。借助这种排序,最终作出选购决策。在决策时,由于6 种电冰箱对于每个中间标准的优劣排序一般是不一致的,因此,决策者首先要对这7 个标准的重要度作一个估计,给出一种排序,然后把6 种冰箱分别对每一个标准的排序权重找出来,最后把这些
34、信息数据综合,得到针对总目标即购买电冰箱的排序权重。象这样类似的问题很多,其特点是这类问题所往往涉及到经济、社会、人文等方面的因素。在作比较、判别、评价、决策时这些因素的重要性、影响力或者优先程度往往难以量化,人的主观选择会起着相当重要的作用,这就给用一般的数学方法解决问题带来本质上的困难。层次分析法( analytical hierarchy process ,AHP )是美国匹兹堡大学教授撒泰(A.L.Saaty )于 20 世纪70 年代提出的一种系统分析方法。它综合定性与定量分析,模拟人的决策思维过程,来对多因素复杂系统,特别是难以定量描述的社会系统进行分析。目前, AHP 是分析多目
35、标、多准则的复杂公共管理问题的有力工具。它具有思路清晰、方法简便、适用面广、系统性强等特点,便于普及推广,可成为人们工作和生活中思考问题、解决问题的一种方法。将AHP 引入决策,是决策科学化的一大进步。它最适宜于解决那些难以完全用定量方法进行分析的公共决策问题。应用AHP 解决问题的思路是,首先,把要解决的问题分层次系列化,将问题分解为不同的组成因素,按照因素之间的相互影响和隶属关系将其分层聚类组合,形成一个递阶的、有序的层次结构模型。然后,对模型中每一层次因素的相对重要性,依据人们对客观现实的判断给予定量表示,再利用数学方法确定每一层次全部因素相对重要性次序的权值。最后,通过综合计算各层因素
36、相对重要性的权值,得到最低层(方案层)相当于最高层(总目标)的相当重要性次序的组合权值,以此作为评价和选择方案的依据。AHP 将人们的思维过程和主观判断数学化,不仅简化了系统分析与计算工作,而且有助于决策者保持其思维过程和决策原则的一致性,对于那些难以全部量化处理的复杂的问题,能得到比较满意的决策结果。因此,它在能源政策分析、产业结构研究、科技成果评价、发展战略规划、人才考核评价以及发展目标分析等许多方面得到广泛的应用。下面介绍层次分析法的基本原理、步骤、计算方法、及其应用。8.3.1 层次分析的基本原理为了说明AHP 的基本原理,首先分析下面这个简单的事实。假定我们已知n 只西瓜的每只西瓜的
37、重量分别为1w,2w,nw且总和为1,即11niiw。把这些西瓜两两比较(相除) ,很容易得到表示n 只西瓜相对重量关系的比较矩阵(以后称之为判断矩阵):nnijnnnnnnaWWWWWWWWWWWWWWWWWW212221212111(8.16)精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 10 页,共 20 页显然iia=1,jiijaa1,jkikijaaa,nkji,2 ,1,对于矩阵nnija,如果满足关系jkikijaaa(nkji,2, 1,),则称矩阵具有完全一致性。可以证明具有完全一致性的矩阵A=nnija有以下性质:1)A 的
38、转置亦是一致阵;2)矩阵 A 的最大特征根nmax,其余特征根均为零。3)设Tnuuuu),(21是 A 对应max的特征向量,则jiijuua,nji,2 ,1,。若记nnnnnnWWWWWWWWWWWWWWWWWWA212221212111,nWWWW21,则矩阵A是完全一致的矩阵,且有AW =nnnnnnWWWWWWWWWWWWWWWWWW212221212111nWWW21=nnWnWnW21= nW ( 8.17)即 n 是 n 只西瓜相对重量关系的判断矩阵A 的一个特征根,每只西瓜的重量对应于矩阵A 特征根为 n 的特征向量W 的各个分量。很自然,我们会提出一个相反的问题,如果事先
39、不知道每只西瓜的重量,也没有衡器去称量,我们如果能设法得到判断矩阵A(比较每两只西瓜的重量是容易的),能否导出每只西瓜的重量呢?显然是可以的,在判断矩阵具有完全一致的条件下,我们可以通过解特征值问题AW=maxW 求出正规化特征向量(即假设西瓜总重量为1),从而得到n 只西瓜的相对重量。同样,对于复杂的社会公共管理问题,通过建立层次分析结构模型,构造出判断矩阵,利用特征值方法即可确定各种方案和措施的重要性排序权值,以供决策者参考。对于 AHP ,判断矩阵的一致性是十分重要的。此时矩阵的最大特征根nmax,其余特征根均为零。在一般情况下,可以证明判断矩阵的最大特征根为单根,且nmax。当判断矩阵
40、具有满意的一致性时,最大的矩阵的特征值为n,其余特征根接近于0,这时,基于AHP 得出的结论才基本合理。但由于客观事物的复杂性和人们认识上的多样性,要求判断矩阵都具有完全一致性是不可能的,但我们要求一定程度上的一致,因此对构造的判断矩阵需要进行一致性检验。8.3.2 层次分析法的计算步骤一、建立层次结构模型运用AHP 进行系统分析,首先要将所包含的因素分组,每一组作为一个层次,把问题条理化、层次化,构造层次分析的结构模型。这些层次大体上可分为3 类1、最高层:在这一层次中只有一个元素,一般是分析问题的预定目标或理想结果,因此又称目标层;2、中间层:这一层次包括了为实现目标所涉及的中间环节,它可
41、由若干个层次组成,包括所需要考虑的准则,子准则,因此又称为准则层;3、最底层:表示为实现目标可供选择的各种措施、决策、方案等,因此又称为措施层或方案层。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 11 页,共 20 页层次分析结构中各项称为此结构模型中的元素,这里要注意,层次之间的支配关系不一定是完全的,即可以有元素(非底层元素)并不支配下一层次的所有元素而只支配其中部分元素。这种自上而下的支配关系所形成的层次结构,我们称之为递阶层次结构。递阶层次结构中的层次数与问题的复杂程度及分析的详尽程度有关,一般可不受限制。为了避免由于支配的元素过多而给
42、两两比较判断带来困难,每层次中各元素所支配的元素一般地不要超过9 个,若多于9 个时,可将该层次再划分为若干子层。例如,大学毕业的选择问题,毕业生需要从收入、社会地位及发展机会方面考虑是否留校工作、读研究生、到某公司或当公务员,这些关系可以将其划分为如图8.1 所示的层次结构模型。图 8.1 再如,国家综合实力比较的层次结构模型如图6 .2:图 6 .2 图中,最高层表示解决问题的目的,即应用AHP 所要达到的目标;中间层表示采用某种措施和政策来实现预定目标所涉及的中间环节,一般又分为策略层、约束层、准则层等;最低层表示解决问题的措施或政策(即方案)。然后,用连线表明上一层因素与下一层的联系。
43、如果某个因素与下一层所有因素均有联系,那么称这个因素与下一层存在完全层次关系。有时存在不完全层次关系,即某个因素只与下一层次的部分因素有联系。层次之间可以建立子层次。子层次从属于主层次的某个因素。它的因素与下一层次的因素有联系,但不形成独立层次,层次结构模型往往有结构模型表示。二、构造判断矩阵任何系统分析都以一定的信息为基础。AHP 的信息基础主要是人们对每一层次各因素的相对重要性给出的判断,这些判断用数值表示出来,写成矩阵形式就是判断矩阵。判断矩阵是AHP 工作的出发点,构造判断矩阵是AHP 的关键一步。当上、下层之间关系被确定之后,需确定与上层某元素(目标A 或某个准则Z)相联系的下层各元
44、素在上层元素Z 之中所占的比重。假定 A 层中因素Ak 与下一层次中因素B1,B2, Bn 有联系,则我们构造的判断矩阵如表8.16 所示。表 8.16 判断距阵Ak B1 B2 Bn B1 B2 Bn b11 b21 bn1 b12 b22 bn2 b1n b2n bnn 表 8.16 中, bij 是对于 Ak 而言, Bi 对 Bj 的相对重要性的数值表示,判断矩阵表示针对上一层次某因素而言,本层次与之有关的各因素之间的相对重要性。填写判断矩阵的方法是:向填写人(专家)反复询问:针对判断矩阵的准则,其中两个元素两两比较哪个重要,重要多少。对重要性程度Saaty等人提出用1-9 尺度赋值,
45、见下表8.17 表 8.17 重要性标度含义表重要性标度含义精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 12 页,共 20 页1 表示两个元素相比,具有同等重要性3 表示两个元素相比,前者比后者稍重要5 表示两个元素相比,前者比后者明显重要7 表示两个元素相比,前者比后者强烈重要9 表示两个元素相比,前者比后者极端重要2,4, 6,8 表示上述判断的中间值倒数若元素i与元素 j 的重要性之比为ijb, 则元素 j 与元素i的重要性之比为jib=ijb1设填写后的判断矩阵为nnijbB,则判断矩阵具有如下性质:(1) ijb0,(2) jib=i
46、jb1,(3) iib=1 .,.2, 1ni根据上面性质,判断矩阵具有对称性,因此在填写时,通常先填写iib=1 部分,然后再仅需判断及填写上三角形或下三角形的n(n-1)/2 个元素就可以了。在特殊情况下,判断矩阵可以具有传递性,即满足等式:ikjkijbbb,当上式对判断矩阵所有元素都成立时,则该判断矩阵为一致性矩阵。采用 19 的比例标度的依据是:(1)心理学的实验表明,大多数人对不同事物在相同属性上差别的分辨能力在59 级之间,采用19 的标度反映了大多数人的判断能力;(2)大量的社会调查表明,19的比例标度早已为人们所熟悉和采用;(3)科学考察和实践表明,19 的比例标度已完全能区
47、分引起人们感觉差别的事物的各种属性。因此目前在层次分析法的应用中,大多数都采用尺度。当然,关于不同尺度的讨论一直存在着。三、层次单排序所谓层次单排序是指根据判断矩阵计算对于上一层某因素而言本层次与之有联系的因素的重要性次序的权值。它是本层次所有因素相对上一层而言的重要性进行排序的基础。层次单排序可以归结为计算判断矩阵的特征根和特征向量问题,即对判断矩阵B,计算满足BW =maxW (8. 18)的特征根与特征向量。式中,max为 B 的最大特征根;W 为对应于max的正规化特征向量;W 的分量iw即是相应因素单排序的权值。为了检验矩阵的一致性,需要计算它的一致性指标CI,CI 的定义为CI =
48、1maxnn(8.19)显然,当判断矩阵具有完全一致性时,CI=0。nmax越大, CI 越大,判断矩阵的一致性越差。注意到矩阵 B 的 n 个特征值之和恰好等于n, 所以 CI 相当于除max外其余n-1 个特征根的平均值。为了检验判断矩阵是否具有满意的一致性,需要找出衡量矩阵B 的一致性指标CI 的标准, Saaty引入了随机一致性指标表 8.18。表 8.18 19 矩阵的平均随机一致性指标阶数1 2 3 4 5 6 7 8 9 RI 0.00 0.00 0.58 0.90 1.12 1.24 1.32 1.41 1.45 对于 1 阶、 2 阶判断矩阵, RI 只是形式上的,按照我们对
49、判断矩阵所下的定义,1 阶、 2 阶判断矩阵总是完全一致的。当阶数大于2 时,判断矩阵的一致性指标CI,与同阶平均随机一致性的指标RI 之比RICI称为判断矩阵的随机一致性比率,记为CR。当 CR=RICI0.01 时,判断矩阵具有满意的一致性,否则就需对判断矩阵进行调整。四、层次总排序精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 13 页,共 20 页利用同一层次中所有层次单排序的结果,就可以计算针对上一层次而言本层次所有因素重要性的权值,这就是层次总排序。层次总排序需要从上到下逐层顺序进行,设已算出第k-1 层上 n 个元素相对于总目标的排序
50、为Tknkkwww),()1()1(1)1(,第 k 层kn个元素对于第1k层上第 j 个元素为准则的单排序向量Tkjnkjkjkjkuuuu),()()(2)(1)(.,.2 ,1njknk, 2, 1其中不受第j 个元素支配的元素权重取零,于是可得到nnk阶矩阵)(kU=)()()(2)(1,knkkuuu=)()(21)(2)(22)(21)(1)(12)(11knnknnknkkknkkkkkuuuuuuuuu其中)(kU中的第j列为第 k 层kn个元素对于第1k层上第 j 个元素为准则的单排序向量。记第 k 层上各元素对总目标的总排序为:Tknkkwww),()()(1)(则)(kw
51、)( kU)1(kw)()(21)(2)(22)(21)(1)(12)(11knnknnknkkknkkkkkuuuuuuuuu)1()1(2)1(1knkkwww= njkjkjnnjkjkjnjkjkjwuwuwuk1)1()(1)1()(21)1()(1即有njkjkijkiwuw1)1()()(,kni, 2, 1五、一致性检验为评价层次总排序的计算结果的一致性如何,需要计算与单排序类似的检验量。由高层向下,逐层进行检验。设第k 层中某些因素对k-1 层第 j 个元素单排序的一致性指标为)(kjCI,平均随机一致性指标为)(kjRI,(k 层中与 k-1 层的第 j 个元素无关时,不必
52、考虑),那么第k 层的总排序的一致性比率为:kknjkjkjnjkjkjkRIwCIwCR1)()1(1)()1()(同样当)( kCR 0.10 时,我们认为层次总排序的计算结果具有满意的一致性。8.3.3 层次分析法的应用层次分析法在T .L.Saaty 正式提出以来,由于它在处理复杂的决策问题上的实用性和有效性,很快就在世界范围内得到普遍的重视和广泛的应用,目前精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 14 页,共 20 页它的应用已遍及经济计划和管理、能源政策和分配、行为科学、军事指挥、运输、农业、教育、医疗、环境等领域。从处理的类型
53、看,主要是决策、评价、分析、预测等,这个方法在20 世纪 80 年代初引入我国,很快为广大的数学工作者和有关领域的科技人员所接受,得到了成功的应用。1 旅游地的选择问题:某人准备假期旅游,初次筛选了桂林、黄山和北戴河三处旅游地,但每个旅游地的景色、所需费用、居住条件、饮食条件交通等各不相同,如何在3 个旅游地中按照景色、费用、居住条件、饮食和路途 6 个因素选择一个最佳的旅游地。根据层次分析的基本思想,可分以下几步进行处理:将选择旅游地的决策问题分解为三个层次,最上层为目标层,即选择旅游地,最下层为方案层,有P1(桂林)、 P2(黄山)、 P3(北戴河)三个供选择的地点,中间层为准则层,有C1
54、(景色)、 C2(费用)、C3(居住)、 C4(饮食)、 C5(旅途) 5 个准则,各层间的联系用相连的直线表示如图8.3 所示。目标层O(选择旅游地 ) 准则层C1 景色C2 费用C3 居住C4 饮食C5 旅途方案层P1 桂林P2 黄山P3 北戴河图 8.3 相对于总目标而言,根据旅游者自己的喜好,给出5 个准则之间的相对重要性,利用Saaty 等人提出用 1-9 尺度赋值,构造准则层对目标的成对比较阵构造判断矩阵1135/13/11125/13/13/12/117/14/1557123342/11A这里,判断矩阵是不一致,如231321aaa,为了计算对于上一层选择旅游地而言本层次的5 个
55、准则的重要性次序的权值和判断矩阵是否具有满意的一致性,利用MATLAB软件可求出矩阵A 的最大特征根=5.073 及对应于 的正规化特征向量Twwww),()2(5)2(2)2(1)2(=(0.263,0.475,0.055,0.090,0.110T)。一致性指标018.0155073.5CI随机一致性指标RI=1.12 (n=5, 查表 ),一致性比率CR=0.018/1.12=0.0160.1 通过一致性检验,所以判断矩阵具有满意的一致性。对于上一层选择旅游地而言本层次的5 个准则的重要性次序的权值Twwww),()2(5)2(2)2(1)2(=(0.263,0.475,0.055,0.0
56、90,0.110T)。同样求第3 层(方案 )对第 2 层每一元素 (准则 )的权向量:方案层对C1(景色 )的成对比较阵12/15/1212/15211B,精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 15 页,共 20 页最大特征根,006.31对应于,1的正规化特征向量)3(1u=(0.595,0.276,0.128T)。一致性指标003.0133006.31CI随机一致性指标1RI=0.58 (n=3) ,一致性比率CR=0.003/0.58=0.00520.1 ,通过一致性检验,所以判断矩阵具有满意的一致性。对于上一层景色而言本层次的3
57、 个方案的重要性次序的权值)3(1u(0.595,0.276,0.128T)。方案层对C2(费用 )的成对比较阵1383/1138/13/112B,最大特征根,002.32对应于,2的正规化特征向量)3(2u=(0.082,0.236,0.682T)。一致性指标001. 0133002. 32CI,随机一致性指标2RI=0.58 (n=3) ,一致性比率CR=0.001/0.58=0.00170.1, 通过一致性检验,所以判断矩阵具有满意的一致性。对于上一层费用而言本层次的3 个方案的重要性次序的权值)3(2u(0.082,0.236,0.682T)。方案层对C3(居住 )的成对比较阵1313
58、13123113B最大特征根, 33对应于,3的正规化特征向量)3(3u=(0.429,0.429,0.142T)。一致性指标0133.33CI,随机一致性指标RI=0.58 (n=3) ,一致性比率3RI=0/0.58=00.1, 通过一致性检验,所以判断矩阵具有满意的一致性。对于上一层居住而言本层次的3 个方案的重要性次序的权值)3(3u(0.429,0.429,0.142T)。方案层对C4(饮食 )的成对比较阵114111314314B,最大特征根,009.34对应于,4的正规化特征向量)3(4u=(0.634,0.192,0.174T)。一致性指标0045.0133009. 34CI随
59、机一致性指标4RI=0.58 (n=3) ,一致性比率CR=0.0045/0.58=0.00780.1, 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 16 页,共 20 页通过一致性检验,所以判断矩阵具有满意的一致性。对于上一层饮食而言本层次的3 个方案的重要性次序的权值)3(4u=(0.634,0.192,0.174T)。方案层对C5(旅途 )的成对比较阵144411141115B最大特征根, 35对应于,5的正规化特征向量)3(5u=(0.167,0.167,0.667T)。一致性指标013335CI,随机一致性指标5RI=0.58 (n=
60、3) ,一致性比率CR=0/0.58=00.1, 通过一致性检验,所以判断矩阵具有满意的一致性。对于上一层旅途而言本层次的3 个方案的重要性次序的权值)3(5u=(0.167,0.167,0.667T)。)3(U=)()3(5)3(2)3(1,uuu667.0174.0142.0682.0128.0167.0192.0429.0236.0276.0167.0634.0429.0082.0595.0第 3 层(方案层)对第1 层(目标层)的组合权向量)2()3()3(wWw=109.0098.0054.0476.0264.0667.0174.0142.0682.0128.0167.0192.04
61、29.0236.0276.0167.0634.0429.0082.0595.0455.0245.0299.0即选取桂林、黄山、北戴河的权重分别为0.299、0.245 和 0.455. 组合一致性检验记)0,0 0 5. 0, 0,0 0 1. 0,0 0 3.0(),(54321CICICICICICI)58.0,58. 0,58.0 ,58. 0,58.0(),(54321RIRIRIRIRIRI)3(CI=)2(wCI= 00176.0109.0098.0054.0476.0264.00005.00001.0003.0)3(RI=)2(wRI=58.0109.0098.0054.0476
62、.0264.058.058.058.058.058.0)3(CR=)3()3(RICI=58.000176.0=0.0031.0, 通过一致性检验,我们认为层次总排序的计算结果具有满意的一致性,所以旅游地的选取次序为北戴精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 17 页,共 20 页河、桂林、黄山,它们的权重分别为0.455、0.299 和 0.245. 2、一笔留成利润利用的综合评价背景:某企业有一笔留成利润要由领导决定其用途,总目标是希望能促进工厂更进一步发展。可供选择的方案有:作为奖金发给职工;扩建食堂、托儿所等福利设施;开办职工业余学
63、校进行职工培训;建设图书馆或俱乐部等娱乐设施;引进新设备进行技术改造。衡量这些方案(措施)可从以下三方面着眼:是否调动了职工的生产积极性;是否提高了企业的技术水平;是否改善了职工的物质文化生活状况。现在要对上述五种方案进行优劣性评价,或者说是按优劣顺序把这五种方案排列起来,以便领导从中选择一种方案付诸实施。建立层次结构模型:我们应用AHP 对此问题进行分析后,可建立如图8.4 所示的层次模型。图 8.4 首先根据准则层C 各因素相对于目标层的相对重要性构造比较矩阵并进行计算,所得判断矩阵B 及相应计算结果如下:131331531511B判 断 矩 阵B的 最 大 特 征 值max= 3.038
64、5 , 最 大 特 征 值 相 应 的 正 规 化 特 征 向 量Twwww),()2(3)2(2)2(1)2(=(0.1042,0.6372,0.2583T),一致性指标CI =0.0193,随机一致性指标RI =0.58(3n),所以一致性比率CR =RICI=0.0332 0.10 ,通过一致性检验,所以判断矩阵具有满意的一致性。所以对于目标层而言准则层次的3 个准则的重要性次序的权值Twwww),()2(3)2(2)2(1)2(=(0.1042,0.6372,0.2583T)同样求第3 层(方案 )对第 2 层每一元素 (准则 )的权向量:(1)判断距阵1B(相对于调动生产积极性准则而
65、言,各方案之间的相对重要性比较)及计算结果131315171312214132113151523131745311B,最大特征值max= 5.126 ,对应的正规化特征向量)3(1u=(0.491,0.232,0.092, 0.138,0.046T)。CI=0.032 RI=1.12 CR=0.0280.10 通过一致性检验。(2)判断距阵2B(相对于提高技术水平准则而言,各方案之间的相对重要性比较)及计算结果精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 18 页,共 20 页13315311513351751317112B。最大特征值max=
66、4.117 ,对应的正规化特征向量)3(2u=(0.055,0.564,0.118, 0.263T)。CI =0.039 ,RI =0.90 , CR =0.0430.10 ,通过一致性检验。(3)判断距阵3B(相对于改善职工生活状况准则而言,各方案之间的相对重要性比较)及计算结果113131113131331133113B最大特征值max= 4,对应的正规化特征向量)3(3u=(0.406,0.406,0.094, 0.094T)。CI =0 ,RI =0.90 , CR =00.10 ,通过一致性检验。第三层对第一层总排序计算结果见表8.19:表 8.19 第三层对第一层总排序计算结果C1
67、 C2 C3 层次 P 总排序权值方案排序0.104 0.637 0.285 0.491 0.232 0.092 0.138 0.046 0 0.055 0.564 0.118 0.263 0.406 0.406 0.094 0.094 0 0.157 0.164 0.393 0.113 0.172 4 3 1 5 2 CI = 0.028 ,RI = 0.9231 ,CR = 0.03050.10,通过一致性检验。计算结果表明,为合理使用企业利润留成,对于该企业来说,所提出的五种方案优先次序为:P3 办职工业余和短训班进行职工培训,权值为0.393;P5 引进设备,进行企业技术改造,权值为0
68、.172;P2 扩建职工住宅、食堂、托儿所等集体福利设施,权值为0.164;P1 作为奖金发给职工,权值为0.157;P4 建立图书馆、职工俱乐部等文化设施,权值为0.113。企业领导可根据上述排序结果进行决策。需要注意的是,不同的人对不同企业中的不同情况有不同的判断,用不同的判断值计算的排序结果也不一样。所以应当其请那些对所处理的问题有专门研究的人来作判断,因为他们对所处理的问题和周围的环境了解得比较透彻,能得到合理的判断和正确的排序结果。8.3.4层次分析法的优点和局限性1 系统性层次分析法把研究对象作为一个系统,按照分解、比较判断、综合的思维方式进行决策,成为继机理分析、统计分析之后发展
69、起来的系统分析的重要工具。2 实用性层次分析法把定性和定量方法结合起来,能处理许多用传统的最优化技术无法着手的实际问题,应用范围很广,同时,这种方法使得决策者与决策分析者能够相互沟通,决策者甚至可以直接应用它,这就增精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 19 页,共 20 页加了决策的有效性。3 简洁性具有中等文化程度的人即可以了解层次分析法的基本原理并掌握该法的基本步骤,计算也非常简便,并且所得结果简单明确,容易被决策者了解和掌握。以上三点体现了层次分析法的优点,该法的局限性主要表现在以下几个方面:第一只能从原有的方案中优选一个出来,没
70、有办法得出更好的新方案。第二该法中的比较、判断以及结果的计算过程都是粗糙的,不适用于精度较高的问题。第三从建立层次结构模型到给出成对比较矩阵,人主观因素对整个过程的影响很大,这就使得结果难以让所有的决策者接受。当然采取专家群体判断的办法是克服这个缺点的一种途径。层次分析法是一种定性分析与定量分析相结合的系统分析方法,其解决问题的思路是:首先,把要解决的问题分层系列化。然后,对模型中每一层因素的相对重要性,依据人们对客观现实的判断给予定量表示,再利用数学方法确定每一层次全部因素相对重要性次序的权值。最后,通过综合计算各层因素相对重要性的权值,得到最低层(方案层)相对于最高层(总目标)的相对重要性次序的组合权值,以次作为评价和选择方案的依据。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 20 页,共 20 页