AHP层次分析法示例说明.docx
《AHP层次分析法示例说明.docx》由会员分享,可在线阅读,更多相关《AHP层次分析法示例说明.docx(23页珍藏版)》请在冰豆网上搜索。
AHP层次分析法示例说明
AHP(层次分析法)示例说明
AHP(层次分析法)示例说明
(TheAnalgticHierarachyProcess----AHP)
一.AHP预备知识
为了更好地理解AHP,需要准备一些矩阵方面的知识,以下知识都可以从《线性代数》中找到。
1.1特征根与特征向量
设Aaijmn为n阶方阵,若存在常数和非零n维向量g(g1,g2,,gn),使得
Agg
(1)
则称,是矩阵A的特征根(或特征值),非零向量g是矩阵A关于特征根的特征向量。
1.2特征根的求法
由
(1)得Agg0AEg0,这是一个n元一次线性齐次方程组,该方程组如果
有非零解,则其充分必要条件为:
系数行列式为零,即
AE0
(2)
称
(2)式为矩阵A的特征方程,它是一个一元n次方程,由线性代数基本定理知,该方程有且只
有n个根。
1.3重量模型
设u1,u2,,un为n个物体,重量分别是
g1,g2,
gn。
但是,我们并不知道物体的重量,只
知两两之间重量比的比值:
aij
gi
gj
设准则C为比较重量,问题是:
已知aij(1i,j
n),在准则C下对元素u1,u2,
un排序,也就是按其重量大小排序已知。
g1
g1
g1
g1
g2
gn
Aaijnm
g2
g2
g2
g1
g2
gn
gn
gn
gn
g1
g2
gn
对于以下三个特性:
(1)aij
0
(2)aij
1
aji
(3)aij
ajkaik
aij显然满足
(1)与
(2),但是,(3)式通常不被满足(因为统计或构造这么完整的数据很难)
,
满足
(1)、
(2)的矩阵A为正互反矩阵;满足
(1)、
(2)并且(3)也成立时的矩阵A称为一致性
判断矩阵。
问题是:
已知判断矩阵
A,在准则C下对n个物体排序。
即按重量大小排序。
1–1
AHP(层次分析法)示例说明
如果,aij
gi是,gi,gj是重量的精确值,此时(3)式必定成立,即A是一致性判断矩阵。
gj
令
ggg
2
g
T
1
n
则带入计算,
Agng。
显见n是方阵A的特征根,g是A的与
n对应的特征向量;事实
上此时不难验证:
n是方阵A=(aij)的最大特征根,其余n-1个特征根全为零,而g是A的与最大特征
根n对应的特征向量。
(证明见附录)g的n个分量是物体的相对重量,
因此,可按此对u1,u2,,un
排序。
如果对矩阵A有一个小的扰动,即aij
不再是真实重量的比值,这时显然
A不满足一致性条件,
此时A的最大特征根
max不再是n;因扰动很小,自然
max离n不远,这时
max对应的特征向量虽
然不会是n个物体的真实重量g
g1,g2,
gn
T
,但是,变动也不会太大。
我们设想:
如果扰动
不大,则max离n就不远,此时
max对应的特征向量
g与g差不多,如果
g不改变g的各分量的
大小次序,则g同样给出n个物体u1,u2,
un按重量大小的真实排序。
这样,对不满足一致性的正互反矩阵
A
(aij)nn,我们求其最大特征根
max,再求与max对
应的特征向量g,则可按g对n个物体u1,u2,
un按重量大小排序。
但是,这一番理论有几个疑点:
①当A不满足一致性时,A还有没有最大正的特征根;②既使
A有最大特征根,那么,这个最大特
征根max对应的特征向量的全部分量能否还是正数(重量不可能为负数)?
这两个问题可以用矩阵代数中Perro—Frobineus定理回答。
Perro-Frobineus定理:
正矩阵存在重数为1重的正特征根,其它特征根的模均小于这个正特征根,该正特征根对应的
特征向量可以全部由正分量组成,经“归一化”处理后该特征向量是惟一的。
(证明见itac的ecmp平台文档库中Proof_Of_PF_Theorem.pdf)
Perron定理明白地告诉我们,对正互反矩阵A,既使它不满足一致性,也一定存在最大正的实
特征根,它对应的特征向量的各个分量都可以是正数,并且“归一化”后是惟一的。
但是,我们能
否按这个“归一化”后是惟一的特征向量对n个物体按重量大小排序呢?
或说这个“归一化”后的
特征向量是否会改变扰动前的一致性矩阵A的最大特征根max=n对应的特征向量的各分量间大小
的排序呢?
这个问题太难了,人们简直难于正面明确地回答,而只能给出一个并不是十分令人满意的简接回答。
那就是对判断矩阵Aaij的一致性满意程度进行检验:
我们说过,由于对A不大的扰动,最大特征根离
n不应太远,所以一致性检验自然与
我们可以证明:
只要A的一致性不被满足,那么A的最大特征根
max一定比n大,即max
于正互反矩阵最大特征根随扰动的变大而变大的证明没有找到,忘补充)
n有关。
–n>0。
(对
C.I.
max
n
令
n
1
显然,我们希望C.I.尽量小;但是,C.I.小到什么程度,才能使max与n对应的特征向量“归
一化”后各分量大小次序不被破坏呢?
这仍是一个非常非常困难的问题,可以说,人们难以正面回
答这个问题。
为此,AHP发明者Saaty给出了平均一致性检验值R.I.。
我们重复1000次,对随机判
断矩阵A的最大特征根进行计算后求取算术平均值得到如下平均随机一致性检验指标如下:
1–2
AHP(层次分析法)示例说明
阶数
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
R.I.
0
0
0.52
0.89
1.12
1.26
1.36
1.41
1.46
1.49
1.52
1.54
1.56
1.58
1.59
令
C.I.
C.R.
R.I.
当C.R.0.1时,认为判断矩阵
A的一致性是可以被接受的。
亦即当C.R.
0.1即C.I.0.1R.I.
时,就是说,当给定的判断矩阵A
(aij)的一致性指标C.I.不超过平均随机一致性指标
R.I.的0.1
倍时,认为判断矩阵A
(aij
)的一致性是可以被接受的。
言外之意:
此时的
A的max对应的特征向
量“归一化”后,能给出
n
个物体u1,u2,,un按重量大小的真实排序。
明显看出这个回答不是正
面的,也有些令人难以置信。
但是,这已是目前为止最好的回答了,这也是
AHP理论上不够严谨的
问题。
不过,从应用角度讲,当
C.R.<0.1时,排序的正确性已为所有应用例子所证实。
但是,当
.R.>0.1时,AHP不再适用,这时,只能回头考虑,变更递阶层次结构,或对判断矩阵
A
重新赋值。
C
二.AHP基本步骤
用AHP解决问题,有四个步骤:
1.建立问题的递阶层次结构;
2.构造两两比较判断矩阵;
3.由判断矩阵计算被比较元素相对权重;
4.计算各层元素组合权重,并进行一致性检验。
下面通过一个应用实例说明AHP的每个步骤的实施。
例:
某闹市区一商场附近交通拥挤。
目标G:
改善该街区交通环境。
有三种方案可供选择:
A1:
修天桥或修高架桥;A2:
修地道;A3:
商场搬迁。
选择方案的准则有
5个:
c1:
通车能力;c2:
方便市民;c3:
改造费用;c4:
安全性;c5:
市容美观。
决策步骤:
A.建立问题的递阶层次结构:
1–3
AHP(层次分析法)示例说明
1.目标层最高层:
目标层G:
改变交通环境
2.准则层c1:
通车能力c2:
方便市民c3:
改造费用c4:
安全性c5:
市容美观
3.方案层方案A1方案A2方案A3
递阶层次结构中,每一层的每一个元素均是下一层中每个元素的准则。
B.构造两两比较判断矩阵
构造判断矩阵
A
(aij)nn,在单准则下分别构造,即在
G下对c1c2c3
c4c5,构造判断矩阵;
分别在c1c2
c3c4c5下对A1
A2A3构造判断矩阵。
在单一准则下,如何具体构造两两比较判断矩阵
A
(aij)呢?
即如何具体确定比值
aij呢?
在
AHP中比较常用的是――
1-9比例标度法。
关于1-9
比例标度法的说明:
n个元素u1,u2,
un,两两比较其重要性共要比较
n(n1)次。
第i个元素ui与第j个元素uj
2
重要性之比为aij。
通过使用标度比重,确定
aij,一下是标度值:
aij
1
表示ui
与uj
重量相同,或重要性相同;
aij
3
表示ui
比uj
稍重;
aij
5
表示ui
比uj
明显重;
aij
7
表示ui
比uj
强烈重;
aij
9
表示ui
比uj
极端重;
数2、4、6、8则为上述判断的中值。
两两比较两个元素的重要性,总是在某种准则(准则层比较是以总目标
G为准则,方案层比较,
分别以准则层中各元素为准则)下进行的。
至于为什么取
1-9比例标度,而不取别的,是因为人们
直觉最多只能判断出
9
个等级的差异,再细的差异,人的直觉是分辨不出来的,而两两比较判断矩
阵是领域专家靠感觉去分辨和构造的。
从理论上讲,用
1-15比例标度也未尝不可,只是人的直觉分
辨不出。
对n个物体,两两比较其重要性得判断矩阵
A
(aij)nn,显然aij
满足:
1–4
AHP(层次分析法)示例说明
1
aij0,aij
,aii
1
aji
n
n
1)个判断,所以
A
是正的互反矩阵,且对角线上元素为
1
共计1(
,这样的n阶矩阵可表示
2
为上三角或下三角矩阵。
但
A的元素aij通常不具有传递性,即
:
aijajk
aik
这是由事物的复杂性和人的认识的局限性造成的。
如果式:
aijajk
aik
成立,则称
A是一致性矩阵。
从判断矩阵A出发到导出元素在某种准则
C下按重要性大小的排
序,矩阵A的一致性起着至关重要的作用。
按着1–9比例标度的上述说明,具体构造应用举例的六个准则下的两两比较判断矩阵分别为:
通
方
费
安
市
车
便
用
全
容
G
c1
c2
c3
c4
c5
通车
c1
1
3
5
3
5
方便
c2
1/3
1
3
1
3
费用
c3
1/5
1/3
1
1/3
3
安全
c4
1/3
1
3
1
3
市容
c5
1/5
1/3
1/3
1/3
1
通车
方
能力
A1
A2
A3
便
A1
A2
A3
c1
c2
天桥
天
A1
1
1
5
桥A1
1
3
5
地道
1
1
5
地
1
1
2
/3
1–5
AHP(层次分析法)示例说明
A2
道A2
搬迁
1
搬
1
1
1
A3
/5
1
迁A3
/5
1
/5
/2
费用
安
A3
A1
A2
A3
A1A2
全
c3
c4
天桥
天
1
1
1
4
7
1
A1
桥A1
/2
/3
地道
地
1
1
4
2
1
1
A2
/4
道A2
搬迁
1
搬
1
1
3
1
1
A3
/7
/4
迁A3
市容
c5
A1A2A3
天桥
11
1
A1/2/3
地道
211
A2
搬迁
311
A3
C.由判断矩阵计算被比较元素相对权重对给出的共6个正互反矩阵,分别求:
(1)max
(2)与max对应的特征向量并归一化得排序相对权重向量
(3)每个矩阵求max后,都要进行一致性检验。
例如以c1作准则的判断矩阵为:
1–6
AHP(层次分析法)示例说明
1
1
5
A
1
1
5
1/5
1/5
1
因阶数低,可直接求出最大特征根。
由于
A是一致的,知
max=3,其它的特征根均为0。
下面
来验证这一点:
1
1
5
1
1
5
1
5
|AE|1
1
5
1
1
5
()
1
1/5
1/51
1/51/51
1/5
1
1
5
0
0
()(2
2)
0
0
()2
(
)2
(2)
0
max
3
0
2
3
再例如以准则c2的判断矩阵为:
1
3
5
A
1/3
1
2
1/5
1/2
1
显然A不满足一致性,因为
a12
a23
3
26
a13
5
。
1
3
5
由A1/3
1
2
1/5
1/2
1
1
3
5
1
3
5
2
|A
E|1/31
2
1/31
1
2()
1
1/5
1/21
1/2
1/51/21
1
3
5
0
1/3
(
)(2
2
)
0
1/10
(
)2
1
(
)2(2
)
30
1–7
AHP(层次分析法)示例说明
(
2)
1
2(2
)
30
3
3
2
1
0
30
由于A出现一个小的扰动而不满足一致性,此时不能再有
max=3,而是max>3。
这时,通常用
迭代算法(乘幂法)求解出
max与对应的特征向量,关于
乘幂法可baidu一下。
D.计算各层元素组合权重,并进行一致性检验
(1)设准则层元素C相对于总目标G的排序权重向量为:
a1(a11,a12,,a1m)T(本例中m=5)
(2)方案层各方案
A对准则层各元素
j的排序向量为:
bj2
(b12j,b22
j,
bn2j)
(j
1,2,
m)
(本例中n=3,m=5)
令B2
(b12,b22,,bm2)
(
m=5)
则方案层的n(n=3)个方案相对于总目标的组合权重向量为a2
B2a1:
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
*
最后得到的*,*,*T就是方案A、B、C在总目标G下的排序向量。
(3)对于递阶层次组合判断的一致性检验
我们要逐层计算
C.I.,设得到准则层针对目标层的计算结果为:
C.I.1,R.I.1,C.R.1
方案层针对目标层的相应指标为:
C.I.2
C.I.12,C.I.22,,C.I.2m
a1
R.I.2
R.I.12,R.I.22,
R.I.2m
a1
本例中m=5,则C.R.2C.R.1C.I.2(为什么使用加法)
R.I.2
上面C.I.2i
和R.I.2i
分别是方案层针对准则层的第
i个准则下判断矩阵的一致性指标和平均
随机一致性指标。
当
..
0.1
时,认为递阶层次在
2层水平上整个判断有满意的一致性。
CR2
1–8
AHP(层次分析法)示例说明
附录:
关于(n是方阵A=(aij)的最大特征根,其余n-1个特征根全为零,而g是A的与最大特征根n对应的特征向量)的证明:
证明1:
对于一致性正互反举证:
W1
W1
W1
W1
W2
Wn
W2
W2
W2
AW1
W2
Wn
Wn
Wn
Wn
W1
W2
Wn
很容易看出,每行成比例,因此矩阵的秩
=1,非零特征根有
1个。
并且∑I=∑aii,因此n=∑I=
证明2:
设两两比较相对重量的精确测度为:
W1
W1
W1
W1
W2
Wn
A
W2
W2
W2
W1
W2
Wn
Wn
Wn
Wn
W1
W2
Wn
则特征方程|A
E|0,有一重实根
n及n–1重0根。
证明:
W1
W1
W1
W1
W2
Wn
fn()|AE|
W2
W2
W2
W1
W2
Wn
Wn
Wn
Wn
W1
W2
Wn
W1
W1
W1
W1
W2
Wn
W2
W2
W2
()fn1()
W1
W2
Wn
Wn
Wn
Wn
W1
W2
Wn
B
fn1
1–9
AHP(层次分析法)示例说明
W1
W1
W1
W1
W2
Wn
0
0
B
W1
n1
n1
0
0
W2
0
0
fn
n1
n2
fn22
n1
2fn2
fn
n
2