层次分析法建模
层次分析法(AHP-Analytic Hierachy process)---- 多目标决策方法
70 年代由美国运筹学家T·L·Satty提出的,是一种定性与定量分析相结合的多目标决策分析方法论。吸收利用行为科学的特点,是将决策者的经验判断给予量化,对目标(因素)结构复杂而且缺乏必要的数据情况下,采用此方法较为实用,是一种系统科学中,常用的一种系统分析方法,因而成为系统分析的数学工具之一。
传统的常用的研究自然科学和社会科学的方法有:
机理分析方法:利用经典的数学工具分析观察的因果关系;
统计分析方法:利用大量观测数据寻求统计规律,用随机数学方法描述(自然现象、
社会现象)现象的规律。
基本内容:(1)多目标决策问题举例AHP建模方法
(2)AHP建模方法基本步骤
(3)AHP建模方法基本算法
(3)AHP建模方法理论算法应用的若干问题。
参考书:1、姜启源,数学模型(第二版,第9章;第三版,第8章),高等教育出版社
2、程理民等,运筹学模型与方法教程,(第10章),清华大学出版社
3、《运筹学》编写组,运筹学(修订版),第11章,第7节,清华大学出版社
一、问题举例:
A.大学毕业生就业选择问题
获得大学毕业学位的毕业生,“双向选择”时,用人单位与毕业生都有各自的选择标准和要求。就毕业生来说选择单位的标准和要求是多方面的,例如:
①能发挥自己的才干为国家作出较好贡献(即工作岗位适合发挥专长);
②工作收入较好(待遇好);
③生活环境好(大城市、气候等工作条件等);
④单位名声好(声誉-Reputation);
⑤工作环境好(人际关系和谐等)
⑥发展晋升(promote, promotion)机会多(如新单位或单位发展有后劲)等。
问题:现在有多个用人单位可供他选择,因此,他面临多种选择和决策,问题是他将如何作出决策和选择?——或者说他将用什么方法将可供选择的工作单位排序?
工作选择
贡献收入发展声誉工作环境生活环境
B.假期旅游地点选择
暑假有3个旅游胜地可供选择。例如:1P :苏州杭州,2P 北戴河,3P 桂林,到底到哪个
地方去旅游最好?要作出决策和选择。为此,要把三个旅游地的特点,例如:①景色;②费用;③居住;④环境;⑤旅途条件等作一些比较——建立一个决策的准则,最后综合评判确定出一个可选择的最优方案。
目标层
准则层
方案层
C .资源开发的综合判断
7种金属可供开发,开发后对国家贡献可以通过两两比较得到,决定对哪种资源先开发,效用最用。
二、问题分析:
例如旅游地选择问题:一般说来,此决策问题可按如下步骤进行:
(S1)将决策解分解为三个层次,即:
目标层:(选择旅游地) 准则层:(景色、费用、居住、饮食、旅途等5个准则)
方案层:(有1P ,2P ,3P 三个选择地点)
并用直线连接各层次。
(S2)互相比较各准则对目标的权重,各方案对每一个准则的权重。这些权限重在人的思维过
程中常是定性的。
例如:经济好,身体好的人:会将景色好作为第一选择;
中老年人:会将居住、饮食好作为第一选择; 经济不好的人:会把费用低作为第一选择。
而层次分析方法则应给出确定权重的定量分析方法。
(S3)将方案后对准则层的权重,及准则后对目标层的权重进行综合。 (S4)最终得出方案层对目标层的权重,从而作出决策。
以上步骤和方法即是AHP 的决策分析方法。
三、确定各层次互相比较的方法——成对比较矩阵和权向量
在确定各层次各因素之间的权重时,如果只是定性的结果,则常常不容易被别人接受,因而Santy 等人提出:一致矩阵法..... 即:1. 不把所有因素放在一起比较,而是两两相互比较
2. 对此时采用相对尺度,以尽可能减少性质不同的诸因素相互比较的困难,提高准确度。
因素比较方法 —— 成对比较矩阵法:
目的是,要比较某一层n 个因素n C C C , ,,21 对上一层因素O 的影响(例如:旅游决策解中,比较景色等5个准则在选择旅游地这个目标中的重要性)。
采用的方法是:每次取两个因素i C 和j C 比较其对目标因素O 的影响,并用ij a 表示,全部比较的结果用成对比较矩阵表示,即:
)1( 1
,0 ,)(
ij ij ij
ji ij nxn ij a a a a a a A 或 (1) 由于上述成对比较矩阵有特点: ji
ij ij ij a a a a A 1 ,0 , )(
故可称A 为正互反矩阵:显然,由 ji
ij a a 1
,即:1 ji ij a a ,故有:1 ji a
例如:在旅游决策问题中:
2112 a =(费用)(景色)21C C 表示: 2O 1O 21的重要性为(费用)对目标
的重要性为景色)对目标(C C
故:),费用重要性为即景色重要性为21(2
112 a
14413 a = (居住条件)(景色)31C C 表示: 1O C 4O (3
1的重要性为(居住条件)对目标的重要性为景色)对目标C 即:景色为4,居住为1。
17723 a =(居住条件)(费用)32
C C 表示: 1O C 7O (32的重要性为(居住条件)对目标
的重要性为费用)对目标C
即:费用重要性为7,居住重要性为1。
因此有成对比较矩阵:
113
5
13
1112513131211714
1
553374121
21
A ??问题:稍加分析就发现上述成对比较矩阵的问题: ① 即存在有各元素的不一致性,例如:
既然:4
1114a ;22113313113212112
a a C C a C C a 所以应该有:188412
1
3
12
31213223
C C C C a a C C a
而不应为矩阵A 中的1
723 a
②成对比较矩阵比较的次数要求太 ,因:n 个元素比较次数为:!
2)
1(2
n n C n 次, 因此,问题是:如何改造成对比较矩阵,使由其能确定诸因素n C C , ,1 对上层因素O 的权重?
对此Saoty 提出了:在成对比较出现不一致情况下,计算各因素n C C , ,1 对因素(上层因素)O 的权重方法,并确定了这种不一致的容许误差范围。
为此,先看成对比较矩阵的完全一致性——成对比较完全一致性
四:一致性矩阵
Def :设有正互反成对比较矩阵:
1 a , , 1 , , 1 1nn 221122222212211121121111n n n n n n j i
ij n
n n
n W W W W a W W a W W a W W a W W a W W a W W a W W a W W a A
(4) 除满足:(i )正互反性:即
)1 ( 1
ji ij ji
ij ij a a a a a 或 而且还满足:(ii )一致性:即
n 2, 1,j i,
h
a h a a k a a a a j i kj i j i ij A A A n
n
n n n n n W W W W W W W W W W W
W W W W W W W A
2
1
2221212111
321W W W W n nW nW nW W W
W W W W W W W W W W W W W A n n n n n n n
212121
121
110)( W nI A n W W W W 21
n
n M M M , , ,21 n
W W W , , ,21
n W W W W 21A
称特征根法,求权向量的方法量权向量,此种用特征向为即对上层因素O的权重,,C ,,C C ,就表示诸因素=W=则归一化后的特征向量,=:重量向量 为特征根的特征向量为
以的特征根为n 21 1W W W W ,121 i n W W W n n A
n W W W W 21n M M M , , ,21 A A A A W ij a ij a A i W A W W W W n A n max A
A W A W k k
k lim W W k k
lim A A k k
lim A A A *A A A A *A A A A
W W W e A e e A k T k
k lim
111
e W A n n n A I C A A A
W
W )(-及 W max W A =nW AW m ax A W A A A
n
k n k kk n k n
n k
k
n a 1
1
1
1
1
=A 1
ii a n max A n max A n
max
n
a ii
k h
max
max
max
max k k
n
1
1
max
max
n n n
I C k k
n max 0 I C A I
C A A A
10 I C A I R A I R I R
① I R I
C R C
I R A A I C I R n A )(j i a ij
91ij a ji a ji
ij a a 11 ii a 然后再计算A 的一致性指标I C ,因此A 是非常不一致的,此时,I C 值相当大.
② 如此构造相当多的A ,再用它们的I C 平均值作为随机一致性指标。
③ Satty 对于不同的1( n n ~11),用100~500个样本A 计算出上表所列出的随机一致性
指标I R 作为修正值表。
3. 一致性检验指标的定义——一致性比率R C 。
由随机性检验指标R C 可知:
当2 ,1 n 时,0 I R ,这是因为1, 2阶正互反阵总是一致阵。
对于3 n 的成对比较阵A ,将它的一致性指标I C 与同阶(指n 相同)的随机一致性指标I R 之比称为一致性比率——简称一致性指标, 即有: 一致性检验指标的定义——一致性比率 定义:I R I C R C
: I
R I
C R C 当:10 I
R I
C R C 时,认为主观判断矩阵A 的不一致程度在容许范围之内,
可用其特征向量作为权向量。否则,对主观判断矩阵A 重新进行成对比较,构重新的主观判断矩阵A 。 注:上式10 I
R I
C R C 的选取是带有一定主观信度的。
六、标度——比较尺度解:
在构造正互反矩阵时,当比较两个可能是有不同性质的因素i C 和j C 对于上层因素O 的影响时,采用什么样的相对刻度较好,即ij a 的元素的值在(1~9)或(1~9
1)或更多
的数字,Satty 提出用1~9尺度最好,即ij a 取值为1~9或其互反数1~9
1,心理学家也
提出:人们区分信息等级的极限解能力为7±2。可见对n n 阶矩阵,只需作出2
)
1( n n 个判断值即可
注:以上比较的标度Satty 曾用过多种标度比较层,得到的结论认为:1~9尺度不仅在较简单的尺度中最好,而且比较的结果并不劣于较为复杂的尺度。Satty 曾用的比较尺度为:
① 1~3, 1~5, 1~6,…, 1~11,以及
② )1.0( d ~)9.0( d ,其中4 ,3 ,2 ,1 d ③ p
1~P
9,其中 5 ,4 ,3 ,2 P …
等共27种比较尺度,对放在不同距离处的光源亮度进行比较判断,并构造出成对比较矩阵,计算出权向量。同时把计算出来的这些权向量与按照物理学中光强度定律和其他物理知识得到的实际权向量进行对比。结果也发现1~9的比较标度不仅简单,而效果也较好(至少不比其他更复杂的尺度差)
因而用1~9的标度来构造成对比较矩阵的元素较合适。
七、组合权向量的计算——层次总排序的权向量的计算
层次分析法的基本思想:
(1) 计算出下一层每个元素对上一层每个元素的权向量W
def :层次总排序,计算同一层次所有元素对最高层相对重要性的排序权值。
当然要先:①构造下一层每个元素对上一次每个元素的成对比较矩阵 ②计算出成对比较矩阵的特征向量(和法,根法,幂法) ③由特征向量求出最大特征根m ax (由和法,根法,幂法求得) ④用最大特征根m ax 用方式 1
max
n n
I C 及I
R R
C R C
对成对比较矩阵进行一致性检,并通过。
(2) 并把下层每个元素对上层每个元素的权向量按列排成以下表格形式:例,假定:上
层A 有m 个元素,m A A A , , ,21 ,且其层次总排序权向量为m a a a , , ,21 ,
下层B 有n 个元素n B B B , , ,21 ,则按j B 对 i A 个元素的单排序权向量的列向量为ij b ,即有:
注:①若下层元素k B 与上层元素j A 无关系时,取0 kj b ②总排序权向量各分量的计算公式:),,1(1
n i b
a W ij
m
j j i
(3) 对层次总排序进行一致性检验:从高层到低层逐层进行,如果
如果B 层次某些元素对j A 单的排序的一致性指标为j CI ,相应的平均随机一致性指标为
j RI ,则B 层总排序随机一致性比率为:
m
j j
j
m
j j
j
RI a
CI
a R C 1
1
当10 CR 时,认为层次总排序里有满意的一致性,否则应重新调整判断矩阵的元素取值。
八、层次分析法的基本步骤:
(S1)建立层次结构模型
将有关因素按照属性自上而下地分解成若干层次: 同一层各因素从属于上一层因素,或对上层因素有影响,同时又支配下一层的因素或受到下层因素的影响。
最上层为目标层(一般只有一个因素),最下层为方案层或对象层/决策层,中间可以有1个或几个层次,通常为准则层或指标层。
当准则层元素过多(例如多于9个)时,应进一步分解出子准则层。
(S2)构造成对比较矩阵,以层次结构模型的第2层开始,对于从属于(或影响及)上一层每个因素的同一层诸因素,用成对比较法和1~9比较尺度构造成对比较矩阵,直到最下层。
(S3)计算(每个成对比较矩阵的)权向量并作一致性检验
① 对每一个成对比较矩阵计算最大特征根m ax 及对应的特征向量(和法、根法、幂法等)
n W W W 1
② 利用一致性指标I C ,随机一致性指标R C 和一致性比率作一致性检验
I R I C CR ③ 若通过检验(即1.0 R C ,或1.0 I C )则将上层出权向量
n W W W 1归一化之后
作为(j B 到j A )的权向量(即单排序权向量) ④ 若1.0 R C 不成立,则需重新构造成对比较矩阵
(S4)计算组合权向量并作组合一致性检验——即层次总排序
① 利用单层权向量的权值m j W W W n j , ,11
构组合权向量表:并计算出特征根,组
合特征向量,一致性
② 若通过一致性检验,则可按照组合权向量 n W W 1的表示结果进行决策(
n W W 1中i W 中最大者的最优),即:
T
n i W W W W W ,,:max *1
③ 若未能通过检验,则需重新考虑模型或重新构造那些一致性比率,CR 较大的成对比较
矩阵
九、特征根的近似求法(实用算法)
层次分析法的基本思路是计算上层每个元素对下一层次各元素的权向量(即最大特征根
m ax 对应的特征向量
n W W W 1),以及组合权向量及一致性检验问题。
计算判断矩阵最大特征根和对应阵向量,并不需要追求较高的精确度,这是因为判断矩阵本身有相当的误差范围。而且优先排序的数值也是定性概念的表达,故从应用性来考虑也希望使用较为简单的近似算法。常用的有以下求特征根的近似求法:“和法”、“根法”、“幂法”,具体如下:
1.“和法”求最大特征根和对应特征向量(近似解)
(S1)将矩阵nxm ij a A )( 的每一列向量的归一化得: n i ij
ij ij a a W 1
~
(利用数据验证即为:每个
位置的数除以该列的合计)
(S2)对ij W ~
按行求和得: n
j ij i W W 1
~~
(S3)将i W ~
归一化,即有:
n
i i
i
i W W W 1
~
~,则有特征向量:
n W W W 1~ (S4)计算与特征向量
n W W W 1对应的最大特征根m ax 的近似值: n i i i
W AW n 1max )(1
此方法:实际上是将A 的列向量归一化后取平均值作为A 的特征向量。
解释: 当A 为一致矩阵时,它的每一列向量都是特征向量
可以在A 的不一致性不严重时,取A 的列向量(归一化后)的平均值作为近似特征向
量是合理的(有依据的)。
2.“根法”求最大特征根特征向量近似值:
步骤与“和法”相同,只是在(S2)时:对归一化后的列向量按行“求和”改为按行“求
积”再取n 次方根,即:n
n
j ij i W W 1
1~~
。
即有具体步骤:
(S1)将矩阵min )(ij a A 的每一列向量归一化得:
n
i ij
ij
ij a
a W 1
~
(S2)对归一化以后的列向量各元素:
n
i ij
ij
ij a
a W 1
~
按行“求和”并开n 次方根得:n
n
j ij i W W 1
1~~
(S3)再将i W ~归一化得:
n i n n j ij n
n
j ij n i i i i W W W W W 11
11
11~~~~
得到特征向量近似值:
n W W W W
21 (S4)计算最大特征根:
i
i
W A n )(1max 作为最大特征根的近似值。 注:“根法”是将“和法”中求列向量的算术平均值改为求几何平均值。
3.“幂法”求最大特征根:
(S1)任取n 维归一化初始向量)
0(W
(S2)计算 ,2 ,1 ,0 ,~)()
1( k AW W
k k
(S3))1(~ k W 归一化,即令:
n
i k i
k k W W W 1
)
1()1()1(~~
(S4)对预先给定的 ,当) , ,2 ,1()()1(n i W W k i k i 时,)
1( k W
即为所求的特征向
量;否则返回(S2) (S5)计算最大特征根, n i k i
k i W W n 1)()1(max
~
1 以上用幂法求最大特征根m ax 对应特征向量的迭代方法,其收敛性由TH1(教材P325)中
的3)W e
A e e
A k T k
k lim ,其中
11 e ,W 是对应m ax 的归一化向量 特征。 (证明:可以将A 化为标准形证明)保证。)
0(W
任意选取,也可以取由“根法”、“和法”得
到的
n W W W 21
注:在以上求特征根和特向量的方法中“和法”最简单。
例:在旅游问题中,求目标层到准则层的成对比较矩阵为A 的特征向量和最大特征根:
准则层:
方案层:
1135131112513131211
7
1
415571233421
1A =
113
2
.0333.01122
.0333.0333.05.01143
.025.055712
33
45.01
利用“和法”求A 的特征向量
n W W 1和特征根max
(S1)将
nxn
ij
W A 的元素按列归一化得:
097.0095
.0176.0098.0085.0097.0095.0118.0098.0085.0032.0048.0059.0070.0064.0484.0476.0411.0489.0510.029.0286.0235.0245.0265.0~nxn
ij
W A
333
.1011333.0535.10115.0531732174043.22.02.0143.015.0917.3333.0333.025.0215
4322
各列归一化的分母
(S2)将
nxn W A ij ~中元素ij W ~
按行求和得各行元素之和: n
j ij i W W 1
~~
(S3)再将上述矩阵向量归一化得到特征向量近似值,
W
W A i ~511.0493.0273.037.2312.1~
102.0099.0055.0474.0262.0511.0493.0273.037.2312.1999.41~
1n i i i W W 特征向量
其中
999.4)511.0493.0273.037.2312.1(~
5
1
i W
(S4)计算与特征向量相对应最大特征根(的近似值)
51
5414313212111max
511W W a W W a W W a W W a W W ia W A n n
j i i j n j i i j n j i i j n j i i j i j n i i
i
0802
.5401.255
1
373.598.4960.4038.505.551
102.0548.0099.0493.0055.0273.0474.0388.2262.0323.151102.0102.0099.0165.0095.0087.0 099
.0102.0099.011.0095.0087.0055.0034.00495.0055.0068.0066.0 474
.05.0495.0385.0474.0524.0262.0306.0297.022.0237.0263.051102.0102.0099.0055.0474.0262.0 11 3 0.2 333.0099.0102.0099.0055.0474.0262.0 11 2 0.2 337.0 055.0102.0099.0055.0474.0262.0 0.333 0.5 0.1431 25.0474.0102.0099.0055.0474.0262.05 5 7 12262.0102.0099.0055.0474.0262.033 4 0.5 151
故有最大特征根
102.0099.0055.0474.0262.0 , 0802.5max
对A 一致性检验指标:02.04
0802
.0450802.51
max
n n
CI
1.0018.012
.102
.012.1
CR RI 故通过检验。
十、应用实例
对前面旅游问题进行决策
目标层:
准则层:
决策层:
已知:①目标A 对准则5 4, 3, 2, 1,i i B 的权重向量为:
T
W 102.0 099.0 055.0 474.0 262.0 (由前面已算出)
,并已通过一致性检验。
②准则54321 , , , ,B B B B B 相对于321 , ,P P P 的成对比较矩阵为
1B 对321 , ,P P P 作用的成对比较矩阵为:
15121
521223332
312322
211312
111Y Y b b b
b b b b b b B 同样2B 对321 , ,P P P 作用的成对比较矩阵为:
13
83113
8131
11B
131
131133
3Y Y B
114111314314B
144411
1411
15B 解:
对以上每个比较矩阵都可计算出最大特征根m ax 及对象的特征向量W (即权重向量),并进行一致性检验:CR RI CI
以1B 为例用“和法”求出1B 的特征根m ax 及对立的特征向量1W
15.02.0215.0521
1B
(S1)对1B 按列归一化得:
125.0143.0118.025.0286.0294.0625.0571.0588.0~
1ij W B
(S2)对按列归一化反向量再按行求和:
n
j ij W W 1386.083.0784.1~
(S3)对W ~
按行归一化得到特征向量W n i i
i W W 1
~~
129.0277.0595.0386.083.0784.1386.0386.083.0784.183.0386.083.0784.1784.1
(S4)计算特征根)
(m ax 1
B
15.02.0215
.052
1
B 1
1
1
max
i i
i
W W B n
007.3022.93
1
3007.3015.331129.0387.0277.0833.0595.0794.131129.0129.0139.0119.0277.0258.0277.0298.0595.0645.0554.0595.031129.0129.0277.0595.01 5.0 2.0277.0129.0277.0595.02 1 5.0595.0129.0277.0595.05 2 131)(max
1
B 一致性检验:
1
.0006.058
.00035.058
.01.00035.02
007
.0133007.31max
RI CI CR RI n m
CI
故通过检验,既成对矩阵1B 可以接受。
同样步骤对5432 , , , B B B B ,对 , , , ,54321P P P P P 的影响 用特征向量
5432 , , ,B B B B W
W
W
W
表示
最大特征根用:
54
32max max max max , , ,B B B B 表示
并分别计算一致性检验指标: 5432 B B B B CI CI CI CI
58.0 58.0 58.03 RI
)()()()
(5432 B B B B CR CR CR RI
CI
CR
列表如下: