正理想解 Z 是一个方案集 A 中并不存在的虚拟的最佳方案,它的每个属性值都是决策 矩阵中该属性的最好值;而负理想解 Z 则是虚拟的最差方案,它的每个属性值都是决策矩阵中该属性的最差值。
在 n 维空间中,将方案集 A 中的各备选方案a i 与正理想解 Z 和负理想解 Z 的距离进行比较,既靠近正理想解又远离负理想解的方案就是方案集 A 中的最佳方 ∑x 设由决策人给定各属性的权重向量为 w = (w 1, w 2 , , w n ) ,则步骤三,确定正理想解 Z 和负理想解 Z 。
设正理想解Z的第 j个属性值为 z j ,负理想解 Z 第 j 个属性值为z j ,则⎧⎪max z ij , 正理想解 z j = ⎨综合评价与决策方法及其计算机软件实现评价方法大体上可分为两类,其主要区别在确定权重的方法上。
一类是主观赋权法,多 数采取综合咨询评分确定权重,如综合指数法、模糊综合评判法、层次分析法、功效系数法 等。
另一类是客观赋权,根据各指标间相关关系或各指标值变异程度来确定权数,如主成分 分析法、因子分析法、理想解法(也称 TOPSIS 法)等。
目前国内外综合评价方法有数十种 之多,其中主要使用的评价方法有主成分分析法、因子分析、TOPSIS 、秩和比法、灰色关 联、熵权法、层次分析法、模糊评价法、灰色理论法、物元分析法、聚类分析法、价值工程 法、神经网络法等。
1.理想解法目前已有许多解决多属性决策的排序法, 如理想点法、简单线性加权法、加权平方和 法、主成分分析法、功效系数法、可能满意度法、交叉增援矩阵法等。
本节介绍多属性决策 问题的理想解法,理想解法亦称为 TOPSIS 法, 是一种有效的多指标评价方法。
这种方法通 过构造评价问题的正理想解和负理想解, 即各指标的最优解和最劣解, 并用靠近正理想解 和远离负理想解的程度, 通过计算每个方案到理想方案的相对贴近度来对方案进行排序, 从而选出最优方案。
1.1 方法和原理设 多 属 性 决 策 方 案 集 为A = {a 1 , a 2 , , a m } , 衡 量 方 案 优 劣 的 属 性 向 量 为X = {x 1, , x n },这时方案集 A 中的每个方案 a i ( i = 1, , m )的 n 个属性值构成的向量 是 X i = (x i 1, , x in ) ,它作为 n 维空间中的一个点,能唯一地表征方案 a i 。
**案;并可以据此排定方案集 A 中各备选方案的优先序。
用理想解法求解多属性决策问题的概念简单,只要在属性空间定义适当的距离测度就能 计算备选方案与理想解。
TOPSIS 法所用的是欧氏距离。
至于既用正理想解又用负理想解是 因为在仅仅使用正理想解时有时会出现某两个备选方案与正理想解的距离相同的情况,为了 区分这两个方案的优劣,引入负理想解并计算这两个方案与负理想解的距离,与正理想解的 距离相同的方案离负理想解远者为优。
1.2 TOPSIS 法的算法步骤TOPSIS 法的具体算法如下。
步骤一,用向量规划化的方法求得规范决策矩阵。
设多属性决策问题的决策矩阵X = (x ij )m ⨯n ,规范化决策矩阵 Y = ( y ij )m ⨯n ,则y ij = x ijm i =12ij , i = 1,2, , m , j = 1,2, , n (1) 步骤二,构成加权规范阵Z = (z ij )m ⨯n 。
Tz ij = w j ⋅ x ij , i = 1,2, , m , j = 1,2, , n* 0* * 0 0(2)* i ⎪⎩min z ij , j 为效益型属性 j 为成本型属性, j = 1,2, , n(3)⎧⎪max z ij , 负理想解 z 0j = ⎨∑ ( z ∑ ( z 步骤六,按 C 由大到小排列方案的优劣次序。
j ji⎪⎩min z ij ,j 为成本型属性 j 为效益型属性, j = 1,2, , n(4)步骤四,计算各方案到正理想解与负理想解的距离。
备选方案a i 到正理想解的距离为d i * = nj =1ij - z *j )2 ,i = 1,2, , m(5)备选方案a i 到负理想解的距离为d i 0 = nj =1ij- z 0j )2 , i = 1,2, , m(6)步骤五,计算各方案的排队指标值(即综合评价指数)。
C i * = d i 0 /(d i 0 + d *j ) , i = 1,2, , m(7)* i1.3 示例例 1 研究生院试评估。
为了客观地评价我国研究生教育的实际状况和各研究生院的教学质量,国务院学位委员 会办公室组织过一次研究生院的评估。
为了取得经验,先选 5 所研究生院,收集有关数据资 料进行了试评估,表 1 是所给出的部分数据。
表 1 研究生院试评估的部分数据i1 2 3 4 5j人均专著 x 1 (本/人)0.1 0.2 0.4 0.9 1.2生师比 x 25 6 7 10 2科研经费 x 3 (万元/年)5000 6000 7000 10000 400 逾期毕业率 x 4(%) 4.7 5.6 6.7 2.3 1.8解:第一步,数据预处理数据的预处理又称属性值的规范化。
属性值具有多种类型,包括效益型、成本型和区间型等。
这三种属性,效益型属性越大 越好,成本型属性越小越好,区间型属性是在某个区间最佳。
在进行决策时,一般要进行属性值的规范化,主要有如下三个作用:① 属性值有多种 类型,上述三种属性放在同一个表中不便于直接从数值大小判断方案的优劣,因此需要对数 据进行预处理,使得表中任一属性下性能越优的方案变换后的属性值越大。
② 非量纲化, 多属性决策与评估的困难之一是属性间的不可公度性,即在属性值表中的每一列数具有不同 的单位(量纲)。
即使对同一属性,采用不同的计量单位,表中的数值也就不同。
在用各种 多属性决策方法进行分析评价时,需要排除量纲的选用对决策或评估结果的影响,这就是非 量纲化。
③归一化,属性值表中不同指标的属性值的数值大小差别很大,为了直观,更为了 便于采用各种多属性决策与评估方法进行评价,需要把属性值表中的数值归一化,即把表中 数值均变换到[0,1]区间上。
此外,还可在属性规范时用非线形变换或其它办法,来解决或部分解决某些目标的达到 程度与属性值之间的非线性关系,以及目标间的不完全补偿性。
常用的属性规范化方法有以 下几种。
(1)线性变换原始的决策矩阵为X = (x ij ) m ⨯n ,变换后的决策矩阵记为 Y = ( y ij ) m ⨯n , i = 1, , m ,j = 1, , n 。
设 x max 是决策矩阵第 j 列中的最大值,x min 是决策矩阵第 j 列中的最小值。
若j 为效益型属性,则jj 设给定的最优属性区间为[ x j , x j ] , x j 为无法容忍下限, x j 为无法容忍上限,则⎪1, 若x j ≤ x ij ≤ x j⎪1 - ( x ij - x j ) /( x j - x j ), 若x j < x ij ≤ x j ⎩0,设研究生院的生师比最佳区间为 [5,6] , x 2 = 2 , x 2 = 12 。
表 1 的属性 2 的数据处理y ij = x ij / y max 采用上式进行属性规范化时,经过变换的最差属性值不一定为 0,最佳属性值为 1。
若 j 为成本型属性,则y ij = 1- x ij / x max(8)(9)采用上式进行属性规范时,经过变换的最佳属性值不一定为 1,最差属性值为 0。
(2)标准 0-1 变换为了使每个属性变换后的最优值为 1 且最差值为 0,可以进行标准 0-1 变换。
对效益 型属性 j ,令对成本型属性 j ,令y ij =y ij =x ij - x m j i n x m j a x - x m j i n x m j a x - x ijx m j a x - x m j i n(10)(11)(3)区间型属性的变换有些属性既非效益型又非成本型,如生师比。
显然这种属性不能采用前面介绍的两种 方法处理。
0 * ' "⎧1 - ( x 0j - x ij ) /( x 0j - x 'j ), 若x 'j ≤ x ij < x 0j ⎪ 0 * y ij = ⎨* " * * "⎪其它(12) 变换后的属性值 y ij 与原属性值 x ij 之间的函数图形为一般梯形。
当属性值最优区间的上下限相等时,最优区间退化为一个点时,函数图形退化为三角形。
' "见表 2。
表 2 表 2 的属性 2 的数据处理计算的 function mainqujian=[5,6]; lb=2; ub=12; x=[5 6 7 10 2]';y=guifanhua(qujian,lb,ub,x)function y=guifanhua(qujian,lb,ub,x);y=(1-(qujian(1)-x)./(qujian(1)-lb)).*(x>=lb & x<qujian(1))+(x>=qujian(1) &x<=qujian(2))+(1-(x-qujian(2))./(ub-qujian(2))).*(x>qujian(2) & x<=ub);∑x ∑=i 1 ∑=i 1 1 ( x ij - x j )2 , j = 1,2, , n 。
(4)向量规范化无论成本型属性还是效益型属性,向量规范化均用下式进行变换y ij = x ijn i =12 ij(13)这种变换也是线性的,但是它与前面介绍的几种变换不同,从变换后属性值的大小上无法分辨属性值的优劣。
它的最大特点是,规范化后,各方案的同一属性值的平方和为 1,因此常 用于计算各方案与某种虚拟方案(如理想点或负理想点)的欧氏距离的场合。
(5)标准化处理在实际问题中,不同变量的测量单位往往是不一样的。
为了消除变量的量纲效应,使每 个变量都具有同等的表现力,数据分析中常对数据进行标准化处理,即y ij = x ij - x js j,i = 1,2, , m , j = 1,2, , n ,(14)其中 x j = 1 mmx ij , s j =1 m - m表 1 中的数据经标准化处理后的结果见表 3。
表 3 表 1 数据经标准化的属性值表i1 2 3 4 5j人均专著 x 1( y 1 )-0.9741 -0.7623 -0.3388 0.7200 1.3553 生师比 x 2 ( y 2 )-0.3430 0 0.3430 1.3720 -1.3720 科研经费 x 3( y 3 )-0.1946 0.0916 0.3777 1.2362 -1.5109 逾期毕业率 x 4( y 4 )0.2274 0.6537 1.1747 -0.9095 -1.1463计算的Matlab 程序如下: x=[0.1 5 5000 4.7 0.2 6 6000 5.6 0.4 7 7000 6.7 0.9 10 10000 2.3 1.2 2 400 1.8]; y=zscore(x)我们首先对表 1 中属性 2 的数据进行最优值为给定区间时的变换。