(完整版)层次分析法及matlab程序
圣诞用英语怎么说-江西护理职业学院
层次分析法建模
层次分析法(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)机会多(如新单位或单位发展有后劲)等。
问题:现在有多个用人单位可供他选择,因此,他面临多种选择和决策,问题是他将如何
作出决
策和选择?——或者说他将用什么方法将可供选择的工作单位排序?
工作选择
贡献收入 发展 声誉 工作环境 生活环境
1
B.假期旅游地点选择 <
br>暑假有3个旅游胜地可供选择。例如:
P
1
:苏州杭州,
P
2
北戴河,
P
3
桂林,到底到哪个
地方去旅游最好?要作出决策和选择
。为此,要把三个旅游地的特点,例如:①景色;②费用;
③居住;④环境;⑤旅途条件等作一些比较—
—建立一个决策的准则,最后综合评判确定出一
个可选择的最优方案。
目标层
选择旅游地
景费居饮旅
准则层
色 用 住 食 途
方案层
P
1
P
2
P
3
C.资源开发的综合判断
7种金属可供开发,开发后对国家贡献可以通过两两比较
得到,决定对哪种资源先开发,
效用最用。
对经济发展、贡献U
经济价值开採费 风险费 要求量 战略重要性 交通条件
铁In 铜Co 磷酸盐 钿Ur 铝Al
金Go
可供选择的单位P
1’
P
2 ‘
-----
P
n
二、问题分析:
例如旅游地选择问题:一般说来,此决策问题可按如下步骤进行:
2
(S1)将决策解分解为三个层次,即:
目标层:(选择旅游地)
准则层:(景色、费用、居住、饮食、旅途等5个准则)
方案层:(有
P
1
,
P
2
,
P
3
三个选择地点)
并用直线连接各层次。
(S2)互相比较各准则对目标的权重,各方案对每一个准则的权重。
这些权限重在人的思维过
程中常是定性的。
例如:经济好,身体好的人:会将景色好作为第一选择;
中老年人:会将居住、饮食好作为第一选择;
经济不好的人:会把费用低作为第一选择。
而层次分析方法则应给出确定权重的定量分析方法。
(S3)将方案后对准则层的权重,及准则后对目标层的权重进行综合。
(S4)最终得出方案层对目标层的权重,从而作出决策。
以上步骤和方法即是AHP的决策分析方法。
三、确定各层次互相比较的方法——成对比较
矩阵和权向量
在确定各层次各因素之间
的权重时,如果只是定性的结果,则常常不容易被别人接受,因
而Santy等人提出:一致矩阵法
.....
即:1. 不把所有因素放在一起比较,而是两两相互比较
2.
对此时採用相对尺度,以尽可能减少性质不同的诸因素相互比较的困难,提高准确度。
因素比较方法 —— 成对比较矩阵法:
目的是,要比较某一层
n
个因素<
br>C
1
,C
2
,
, C
n
对上一
层因素O的影响(例如:旅游决策解
中,比较景色等5个准则在选择旅游地这个目标中的重要性)。 <
br>採用的方法是:每次取两个因素
C
i
和
C
j
比较其对
目标因素O的影响,并用
a
ij
表示,全部
比较的结果用成对比较矩阵表示,
即:
A(a
ij
)
nxn
,
a
ij
0, a
ji
1
(或a
ij
a
ij
1)
(1)
a
ij
1
a
ji
由于上述成对比较矩阵有特点:
A(a
ij
)
, a
ij
0,
a
ij
故可称
A
为正互反矩阵:显然,由
a
ij
1
,即:
a
ij
a
ji
1
,故有:
a
ji
1
a
ji
3
例如:在旅游决策问题中:
(O的重要性为1
C
C(景色)
1
景色)对目标
a
12
1
=
1
表示:
2
C(费用)
O的重要性为2
2
2
C(费用)对目标
1,费用重要性为2)
故:
a
12
1
(即景色重要性为
2
O的重要性为4
C<
br>1
(景色)对目标
C(景色)
1
4
a
13
4
=
表示:
1
C(居
住条件)
C(居住条件)对目标O的重要性为1
3
3
即:景色为4
,居住为1。
O的重要性为7
C
2
(费用)对目标
C(
费用)
2
a
23
7
7
=
表示:
1
C(居住条件)
O的重要性为1
3
3
C(居住条件)对目标
即:费用重要性为7,居住重要性为1。
1
2
1
因此有成对比较矩阵:
A
4
1
3
1
3
1
2
1
1
7
1
5
1
5
4
7
1
2
3
33
55
11
23
11
11
??问题:稍加分析就发现上述成对比较矩阵的问题:
①
即存在有各元素的不一致性,例如:
既然:
a
12
C
1
1
C
411
a
21
2;
a
13
1
a
31
C
2
2C
3
1a
13
4
所以应该有:
a23
C
2
a
21
C
1
2
8
8
1
1
C
3
a31
C
3
4
C
1
C
2
而不应为矩阵<
br>A
中的
a
23
7
1
②成对比较矩阵比较的次数要求太 ,因:
n
个元素比较次数为:
C
n
2
n(n
1)
次,
2!因此,问题是:如何改造成对比较矩阵,使由其能确定诸因素
C
1
,
, C
n
对上层因素O的
权重?
对此Saoty提出了:在成对比较出现不一致情况下,计算各因素
C
1
,
,
C
n
对因素(上层
因素)O的权重方法,并确定了这种不一致的容许误差范围。
为此,先看成对比较矩阵的完全一致性——成对比较完全一致性
4
四:一致性矩阵
Def
:
设有正互反成对比较矩阵:
a
W
1
W
11
<
br>W
1 a
12
1
1 ,
, a
W
1
1n
1
W
2
W
n
a
W
2W
2
W
2
21
a
22
1
,
A
W
1
W
,
a
2n
2
W
n
a
W
(4)
i
ij
W
j
a
W
n
W
n
W
n
n1
W
a
n2
a
nn
1
1
W
2W
n
除满足:(i)正互反性:即
a
1
ij
0 a
ij
a
( 或 a
ij
a
ji
1)
ji
而且还满足:(ii)一致性:即
a
a
i
ij
a
a
ah
i
k
a
kj
i
i, j1, 2, n
有点点
错误
j
a
j
h
则称满足上述条件的正互反对称矩阵A为一致性矩阵,
简称一致阵。
一致性矩阵(一致阵)性质:
性质1:
A
的秩
Rank(A)=1显然
A
的唯一非0的特征根为n
性质2:
A
的任一列(行)向量都是对应特征根
n
的特征向量:
即有(特征向量、特征值):
W
1
W
1
WW
W
1
W
<
br>W
1
W
2
W
n
A
22
2
W
W
1
1
W
W
W
2
,则向量
W
W<
br>n
W
2
n
W
n
<
br>
n
W<
br>1
W
2
W
n
W
3<
br>
5
W
1
W
1
W
1
W
2
满足:<
br>AW
W
n
W
n
WW
2
1
即:
(AnI)W0
W
1
W
nW
1
1
W
n
W
nW
<
br>
2
2
nW
W
n
W
nW
W
n
n
n
我的理解:通过A(变换A与W中的
元素
有关)变换将一致W矩阵变成权向量
W(特征向量),如果正互反矩阵W’接近一致矩
阵,同样
的道理变换A可以将W’变成权向量
(这里的权向量与W’稍有不同)
W
1
W
2
W
启发与思
考:既然一致矩阵有以上性质,即
n
个元素
W
1,
W
2,
W
3 ,
…
W
n
构成的向量
W
n
是一致矩阵A的特征向量,则可以把向量W归一化后的向量
,看成是诸元素
W
1,
W
2,
W
3 ,
…
W
n
目标的权向量,因此,可以用求
A
的特征根和特征
向量的办法,求出元素
W
1,
W
2,
W
3 ,
…
W
n
相对于目标O的劝向量。
解释:
一致矩阵即:
n
件物体
M
1
, M
2
,
,M
n
,它们重量分别为
W
1
,
W
2
,
,W
n
,将他们两两比较
W
1
W
2
重量,其比值构成一
致矩阵,若用重量向量
W
右乘
A
,则
W
n
6
A的特征根为n,
W
1
W
2
:重量向量
W=
,则归一化后的特征向量
以n为特征根的特征向量为
W
:
n
W
1
W=
W
i
=1
,就表示诸因素C
1
,C
2
,
,C
n
对上层因素O的权重,即为
W
权向量,此种用特征向量求权向量的方法 称特征根法,
分析:
W
1
W
2
若重量
向量
W
未知时,则可由决策者对物体
M
1
,
M
2
,
,M
n
之间两两相比关系,
W
n
主观作出比值的判
断,或用Delphi(调查法)来确定这些比值,使
A
矩阵(不一定有一致性)
为已
知的,并记此主观判断作出的矩阵为(主观)判断矩阵
A
,并且此
A
(不一致
)在不一致
的容许范围内,再依据:
A
ij
的特征根或和特征向量
W
连续地依赖于
ij
矩阵的元素
a
,即当
a
离一致性
的要求不太远时,
A
的特
征根
i
和特征值(向量)
W
与一致矩阵
A
的特征根
和特征
向量
W
也相差不
大的道理:由特征向量
W
求权向量
W
的方
法即为特征向量法,并由此
引出一致性检查的方法。
问题:Remark
以上讨论的用求特征根来求权向量
W
的方法和思路,在理论上应解决以下问题:
1. 一致阵的性质1是说:一致阵的最大特征根为
n
(即必要条件),但用特征根来
求特征向量
时,应回答充分条件:即正互反矩阵是否存在正的最大特征根和正的特征向量?且如果正互<
br>
7
反矩阵
A
的最大特征根
ma
x
n
时,
A
是否为一致阵?
2. 用主观判断矩阵
A<
br>的特征根
和特征向量
W
连续逼近一致阵
A
的特征根
和特征向量
W
时,即:
由
lim
k
k
k
得到:
limW
k
W
k
即:
limA
k
A
k
是否在理论上有依据。
3.一般情况下,主观判断矩阵
A<
br>在逼近于一致阵
A
的过程中,用与
A
接近的
A*
来代
替
A
,
即有
A*A
,这种近似的替代一致性矩阵
A
的作法,就导致了产生的偏差估计问题,即一致
性检验问题,即要确定一种一致性检验判断指标,由此
指标来确定在什么样的允许范围内,
主观判断矩阵是可以接受的,否则,要
两两比较构造主观判断矩阵。此问题即一致性检
验问题的内容。
以上三个问题:前两个问题由
数学严格比较可获得(见教材P325,定理1、定理2)。第
3个问题:Satty给出一致性指标(
TH1,TH2介绍如下:)
附:
Th1:(教材P326,perronTh 比隆
1970 )对于正矩阵
A
(
A
的所有元素为正数)
(1)
A
的最大特征根是正单根
;
(2)
<
br>对应正特征向量
W
(
W
的所有分量为正数)
1<
br>
k
Ae
1
(3)
lim
Tk
W
其中:
e
为半径向量
,
W
是对应
的归一化特征向量
k
eAe<
br>
1
证明:(3)可以通过将
A
化为标准形证明
Th2:
n
阶正互反阵A的最大特征根
n
;
当
n
时,
A
是一致阵
8
五、一致性检验——一致性指标:
1.一致性检验指标的定义和确定——
CI
(平均值)的定义:
当人们对
复杂事件的各因素,采用两两比较时,所得到的主观判断矩阵
A
,一般不可直接
保证正
互反矩阵
A
就是一致正互反矩阵
A
,因而存在误差(及误差估计问题)。这种
误差,必
然导致特征值和特征向量之间的误差
(
)及
W-W
。此时就导致问题
AW=
max
W
与问
题
AWnW
之间的差别。(上述问题中
max
是主观判断矩阵A
的特征值,
W
是带有偏差的相
对权向量)。这是由判断矩阵不一致性所
引起的。
因此,为了避免误差太大,就要衡量主观判断矩阵
A
的一致性。
因为:
①当主观判断矩阵
A
为一致阵
A
时就有:
=
k
k
1n
1
nn
k
a
k
k
1
n
A
为一致阵时有:
a
ii
1
(a[ii]为对角线k
1k
1
nn
上的值,按照一致性矩阵的理解,它
应该为1)
此时存在唯一的非
max
n
(由一致阵性质1:R
ark(4)=1,
A
有唯一非O最大特征根且
max
n
)
②当主观判断矩阵
A
不是一致矩阵时,此时一般有:
max
n
(Th2)
此时,应有:
max
即:
max
n
k
a
h
max
k
ii
n
(不大理解)
k
max
所以,可以取其平均值作为检验主观判断矩阵的准则,一致性的指标,
即:
C
I
显然:
(1)
当
max
n
时,有:
CI0
,
A
为完全一致性
max
n
n
1
k
max
k
n
1
9
(2)
CI值越大,主观判断矩阵
A
的完全一致性越差,即:
A
偏离
A越远(用特征向量
作为权向量引起的误差越大)
(3) 一般
CI01<
br>,认为主观判断矩阵
A
的一致性可以接受,否则应重新进行两两比
较,构造主观
判断矩阵。
2.随机一致性检验指标——
RI
问题:实际操作时发现:
主观判断矩阵
A
的维数越大,判断的一致性越差,故应放宽对高维矩
阵的一致性要求。
于是引入修正值
RI
来校正一致性检验指标:即定义
RI
的修正值表为:
A
的维数
1 2 3
4 5 6 7 8 9
0.00 0.00
0.58 0.96 1.12 1.24 1.32 1.41 1.45
RI
并定义新的一致性检验指标为:
C
R
C
I
R
I
随机一致性检验指标——
RI
的解释:
为确定
A
的不一致程度的容许范围,需要确定衡量
A
的一致性指示
C
I
的标准。于是Satty
又引入所谓随机一致性指标
RI
,其定义和计
算过程为:
(i
j)
从1~9和1~
1
中随机取值,①
对固定的
n
,随机构造正互反阵
A
,其元素
a
i
j
9
与
a
ji
的互反性,即:
aij
1
且满足
a
ij
a
ji
1
. ,且
a
ii
② 然后再计算
A<
br>
的一致性指标
CI
,因此
A
是非常不一致的,
此时,
CI
值相当大.
③ 如此构造相当多的
A
,再
用它们的
CI
平均值作为随机一致性指标。
④ Satty对于不同的
n
(n1
~11),用100~500个样本
A
计算出上表所列出的随机一
致性
指标
RI
作为修正值表。
3.
一致性检验指标的定义——一致性比率
CR
。
由随机性检验指标
CR
可知:
当
n1,
2
时,
RI0
,这是因为1, 2阶正互反阵总是一致阵。
对于
n3
的成对比较阵
A
,将它的一致性指标
CI
与同阶(指n
相同)的随机一致性指
标
RI
之比称为一致性比率——简称一致性指
标,
即有: 一致性检验指标的定义——一致性比率
10
定义:
C
R
C
IC
I
:
C
R
R
IR
I
C
I
当:
C
R
01
时,认为主观判断矩阵
A
的不一致程度在容许范围之内,R
I
可用其特征向量作为权向量。否则,对主观判断矩阵
A
重
新进行成对比较,构重新的主观
判断矩阵
A
。
注:上式
C
R
C
I
01
的选取是带有一定主观信度的。
R
I
六、标度——比较尺度解:
在构造正互反矩阵时,当比较两
个可能是有不同性质的因素
C
i
和
C
j
对于上层因素O的影响时,採用什么样的相对刻度较好,即
a
ij
的元素的值在(1~9)或(1
~
1
)或更多
9
的数字,Satty提出用1~9尺度最好,即
a<
br>ij
取值为1~9或其互反数1~
1
,心理学家也
9
提出:人
们区分信息等级的极限解能力为
7
±2。可见对
nn
阶矩阵,只需作出判断值即可
标度
a
ij
1
3
5
7
9
2,4,6,8,
倒数1,
定
义
n(n1)
个
2
因素
i
与因素
j
相
同重要
因素
i
比因素
j
稍重要
因素
i
比因素
j
较重要
因素
i
比因素
j
非常重要
因素
i
比因素
j
绝对重要
因素
i
与因素
j
的重要性的比较值介于上述两
个相邻等级之间
因素
j
与
因素
i
比较得到判断值为
a
ij
的互反
数,
111
11111
, , , , , , ,
23456789
a
ji
1
a
ii
1
a
ij
注:以上比较的标度Satt
y曾用过多种标度比较层,得到的结论认为:1~9尺度不仅在较简单
的尺度中最好,而且比较的结果并
不劣于较为复杂的尺度。Satty曾用的比较尺度为:
① 1~3, 1~5, 1~6,…,
1~11,以及
②
(d0.1)
~
(d0.9)
,其中
d1, 2, 3,
4
11
③
1
p
~
9
,其中
P2, 3, 4,
5
…
等共27种比较尺度,对放在不同距离处的光源亮度进行比较判断,并构造出成对比较矩
阵,计
算出权向量。同时把计算出来的这些权向量与按照物理学中光强度定律和其他物理知识得到的实际权向量进行对比。结果也发现1~9的比较标度不仅简单,而效果也较好(至少不比其他更
复杂
的尺度差)
因而用1~9的标度来构造成对比较矩阵的元素较合适。
P
七、组合权向量的计算——层次总排序的权向
量的计算
层次分析法的基本思想:
(1)
计算出下一层每个元素对上一层每个元素的权向量
W
def:层次总排序,计算同一层次所有元素对最高层相对重要性的排序权值。
当然要先:①构造下一层每个元素对上一次每个元素的成对比较矩阵
②计算出成对比较矩阵的特征向量(和法,根法,幂法)
③由特征向量求出最大特征根
max
(由和法,根法,幂法求得)
④用最大特征根
max
用方式
C
I
max
n
n
1
及
C
R
C
R
对成对比较矩
R
I
阵进行一致性检,并通过。
(2) 并把下层每个元素对上层每个元素的权向量按列排成以下
表格形式:例,假定:上
层
A
有
m
个元素,
A
1<
br>, A
2
,
,
A
m
,且其层次总排序权向量为
a
1
, a
2
,
, a
m
,
下层
B
有
n
个元
素
B
1
, B
2
,
,
B
n
,则按
B
j
对
A
i
个元素的单排
序权向量的列向
量为
b
ij
,即有:
层次
A
1
A
1
…
A
1
m
(权向
B
层总是排序权重
量、列向量)
12
a
1
a
2
…
a
m
B
1
B
2
B
n
b
11
b
12
b
1m
W1
a
j
b
1j
W
2
<
br>
a
j
b
2j
j
1
m
j
1
m
m
b
12
b
22
b
2m
b
n1
b
n2
b
nm<
br>
W
n
a
j
b<
br>nj
j
1
max
CI
计算出最大特根(方法:和法、根法、幂法)
一致性检验
C
I
一致性检验比率
max
n
n
1
CI
C
I
j
C
R
R
I
m
a
jCI
3
a
j
m
j
检验
CR01
否?
RI
j
注:①若下层元
素
B
k
与上层元素
A
j
无关系时,取
b
k
j
0
②总排序权向量各分量的计算公式:
W
i
<
br>
ab
j
1
m
jij
(i
1,
,n)
(3)
对层次总排序进行一致性检验:从高层到低层逐层进行,如果
如果
B
层次某些元素对
A
j
单的排序的一致性指标为
CI
j
,相应的平均随机一致
性指标为
RI
j
,则
B
层总排序随机一致性比率为:
C
R
aCI
j
1
m
j
m
j
a
j
1
j
R
I
j
当
CR01
时,认为层次总排序里有满意的一致性,否则应重新调整
判断矩阵的元素取值。
八、层次分析法的基本步骤:
(S1)建立层次结构模型
将有关因素按照属性自上而下地分解成若干层次:
同一层各因素从属于上一层因素,或对上层
因素有影响,同时又支配下一层的因素或受
到下层因素的影响。
13
最上层为目标层(一般只有一个因素),最下层为方案层或对象层决策层,中间可以有
1个 或几个层次,通常为准则层或指标层。
当准则层元素过多(例如多于9个)时,应进一步分解出子准则层。
(S2)构造成对比较矩 阵,以层次结构模型的第2层开始,对于从属于(或影响及)上一
层每个因素的同一层诸因素,用成对比 较法和1~9比较尺度构造成对比较矩阵,直到最下
层。
(S3)计算(每个成对比较矩阵的)权向量并作一致性检验
① 对每一个成对比较矩阵计算 最大特征根
max
及对应的特征向量(和法、根法、幂法等)
W
1
W
W
n
② 利用一致性指标
CI
,随机一致性指标
CR
和一致性比率作一致性检验
CR
C
I
R
I
W
1
③ 若通过检验(即
CR0.1
,或
CI0.1
)则将上层出权向量
W
归一化之后
W
n
作为(
B
j
到
A
j
)的权向量(即单排 序权向量)
④ 若
CR0.1
不成立,则需重新构造成对比较矩阵
(S4)计算组合权向量并作组合一致性检验——即层次总排序
W
1
① 利用单层权向量的权值
W
j
j
1,
, m
构组合权向量表:并计算出特征根,组
W
n
合特征向量,一致性
14
上
单 层
层 重
权
量
向
下层 量
层次
A
1
A
1
…
A
1
m
W
1
W
计算组合权向量
W
n
其中
W<
br>i
m
a
1
a
2
…
a
m
aW
j
1
j
m
ij
B
1
B
2
B
n
最大特征根
max
(i)W
11
W
12
W
1m
W
1
a
j
b
1j
W
2
a
j
b
2j
j
1
m
j
1
m
W
12
W
22
W
2m
W
n1
W
n2
W
nm
W
n
a
j
b<
br>nj
j
1
和法、根法、幂法
一致性检验
CI
CI
j
(j)max
n
n
1
CI0.1
?
一致性随机检验
RI
RI
j
对照表
一致性比率
CR
CI
CR
jRI
m
a
j
CI
j
CR01
?
a
j
m
j
RI
2j
W
1
W
1
② 若通过
一致性检验,则可按照组合权向量
W
的表示结果进
行决策(
W
W
W
n
n
中
W
i
中最大者的最优),即:
W*
maxW:W
i
W
1
,
,W
n
T
③ 若未能通过检验,则需重新考虑模型或重新构造那些一致性比率,
CR
较大
的成对比较
矩阵
九、特征根的近似求法(实用算法)
层次分析法的基本思路是计算上层每个元素对下一层次各元素的权向量(即最大特征根
15
W
1
m
ax
对应的特征向量
W
),以及组合权向
量及一致性检验问题。
W
n
计算判断矩
阵最大特征根和对应阵向量,并不需要追求较高的精确度,这是因为判断矩阵
本身有相当的误差范围。而
且优先排序的数值也是定性概念的表达,故从应用性来考虑也希望
使用较为简单的近似算法。常用的有以
下求特征根的近似求法:“和法”、“根法”、“幂法”,具
体如下:
1.“和法”求最大特征根和对应特征向量(近似解)
a
ij
~
(
S1)将矩阵
A(a
ij
)
nxm
的每一列向量的归一化得:W
ij
n
(利用数据验证即为:每个
a
ij
i
1
位置的数除以该列的合计)
n
~~
~
(S2)对
W
ij
按行求和得:
W
i
W
ij
j
1
(S3)将
W
i
归一化,即有:
W
i
~
~
W
i
W
1
~
,则有特征向量:
W
n
~
W
W
i
n
i
1
W
1
1
n
(AW)
i
(S4)计算与
特征向量
W
对应的最大特征根
max
的近似值:
max
nW
i
1
i
W
n
此
方法:实际上是将
A
的列向量归一化后取平均值作为
A
的特征向量。
解释:
当
A
为一致矩阵时,它的每一列向量都是特征向量
W
可以在
A
的不一致性不严重时,取
A
的列向量(归
一化后)的平均值作为近似特征向
量是合理的(有依据的)。
2.“根法”求最大特征根特征向量近似值:
步骤与“和法”相同,只是在(S2)时:对归
一化后的列向量按行“求和”改为按行“求
~
~
积”
再取
n
次方根,即:
W
i
W
ij
。
j
1
n
1
n
即有具体步骤:
16
(S1)将矩阵
A(a
ij
)
min
的每一列向量归一化得:
W
ij<
br>
~
a
ij
a
i
1
n
ij
(S2)对归一化以后的列向量各元素:
W
ij
~
a
ij
a
i
1
n
ij
n
~
~
按行“求和”并开
n次方根得:
W
i
W
ij
j
1
1
n
<
br>n
~
W
ij
~
W
i
j
1
~
(S3)再将<
br>W
i
归一化得:
W
i
n
1~
n
n
W
i
n
~
W
ij
i
1
i
1
j
1
W
1
W
2
W
得到特征向量近似值:
W
n
(S4)计算最大特征根:
max
<
br>1
n
1
(AW)
i
作为最大特征根的近似值。
nW
i
注:“根法”是将“和法”中求列向量的算术平均值改为求几何平均值。
3.“幂法”求最大特征根:
(S1)任取
n
维归一化初始向量
W
(S2)计算
W
(0)
~
(k
1)
AW
(k)
,
k
0, 1, 2,
~
(k
1)
~
(
k
1)
(k
1)
W
(S
3)
W
归一化,即令:
W
~
(k
1)
W
i
i
1
n
(S4)对预先给定的
,当
W
i
量;否则返回(S2)
(S5)计算最大特征根,
max
(k
1)
W
i
(k)
(i
1, 2,
, n)
时,
W
(k1)
即为所求的特征向
~
1
n
W
i
(k
1)
<
br>(k)
n
i
1
W
i
17
以上用幂法求最大特征根
max
对应特征向量的
迭代方法,其收敛性由TH1(教材P325)中
1
Ae<
br>的3)
lim
Tk
W
,其中
e
,
W
是对应
max
的归一化向量
特征。
k
eAe
1
k
(证明:可以将
A
化为标准形证明)保证。
W
(0)
任意选取,也可以取由“根法”、“和法”得
W
1
到的
W
W
2
W
n
注:在以上求特征根和特向量的方法中“和法”最简单。
例:在旅游问题中,求目标层到准则层的成对比较矩阵为
A
的特征向量和最大特征根:
选择旅游地
准则层:
景费居饮旅
色 用 住 食 途
方案层:
P
1
P
2
P
3
1
1
2
433
<
br>21755
11
10.5433
11
21755
A
4
1
23
10.50.333
<
br>
11
7
11
=
0.250.143
2
3
0.
3330.2211
0311
1
5
1
35
311
.3330.
2
18
1.312<
br>
2.37
AW
~~
i
0.273
W
0.493
0.511
W1
利用“和法”求
A
的特征向量
W
和特征根
max
W
n
(S1)将
A
Wij
nxn
的元素按列归一化得:
0.265
0.2450.2350.2860.29
0.4890.4110.4760.484
A
W
~
ij
<
br>0.510
nxn
0.0640.0700.0590.0480
.032
0.0850.0980.1180.0950.
097
0.0850.0980.1760.0950.097
2
1
2
0.25
0.333
0.333
3.917
2
0.5
1
0.143
0.2
0.2
2.043
3
4
7
1
2
3
17
各列归一化的分母
4
3
5
0.5
1
1
10.5
5
3
5
0.333
1
1
10.333
(S2)将
A
W
~<
br>~
~
n
~
ij
nxn
中元素
W<
br>ij
按行求和得各行元素之和:
W
i
W
ij
j
1
(S3)再将上述矩阵向量归一化得到特征向量近似值,
1.312
0.262
2.37
0.474
W
W
~
i
1
n
0.273
0.
W
4.999
.493
05
5
0.099
特征向量
i
0
i
1
0.511
0.102
5
其中
W
~
i
(1.312
2.37
0
.273
0.493
0.511)
4.999
1
(S4)计算与特征向量相对应最大特征根(的近似值)
max<
br>
1
n
AW
i
n
i
1
W
i
n
nnn
a
1
iaW
2j
W
i
a
3
j
W
i
a
4j
W
i
1j
ii
j
1i
j
1i
j
1
a
5j
W
i
5
W
i
j
1
1
W
2
W
3
W
4
W
5
19
0.262
0.262
0.262
0.4740.4740.4
74
1 0.5 4
33
0.055
21 7 5
5
0.055
0.25 0.1431
0.5 0.333
0.055
0.099
0.099
0.099
0.102
0.102
0.10
2
1
5
0.262
0.474
0.055
0.262
0.474
0.262
0.474
0.337 0.2
2 11
0.055
0.333 0.2
3 11
0.055
0.09
9
0.099
0.102
0.102
0.099
0.102<
br>
<
br>1
0.263
0.237
0.22
0.297
0.3060.524
0.474
0.385
0.495
0.
5
0.262
5
0.474
0.066
0.068
0.055
0.0495
0
.0340.087
0.095
0
0.055
.11
0.099
0.102
0.099
0.087
0.095
0.165
0.09
9
0.102
0.102
1<
br>
1.3232.3880.2730.4930
5
0.
262
0.474
0.055
0.099
.548
0.102
1
5
5.05
5.038
4.960
4.9
8
5.373
1
5
25.401
5.0802
0.262
0.
474
故有最大特征根
max
5.080
2 , W
0.055
0.099
0.102
20
对
A
一致性检验指标:
CI
max
n
n
1
5.0802<
br>
50.0802
0.02
44
RI
1.12
CR
故通过检验。
0.02
0.018
0.1
1.
12
十、应用实例
对前面旅游问题进行决策
目标层:
A
选择旅游地点
0.262 0.474 0.099 0.102
0.055
景费居饮
色 用 住 食
准则层:
B
1
B
2
B
3
B
4
0.595 0.129 0.129
0.277
决策层:
P
1
P
2
P
3
已知:①目标
A
对准则
B
i
i1, 2, 3,
4, 5
的权重向量为:
T
,并已通过一致性检
W
0.262 0.474
0.055 0.099 0.102
(由前面已算出)
旅
途
B
5
验。
②准则
B
1
,
B
2
, B
3
, B
4
,
B
5
相对于
P
1
, P
2
,
P
3
的成对比较矩阵为
B
1
对
P
1
,
P
2
, P
3
作用的成对比较矩阵为:
b
1
1
b
12
B
1
b
21b
22
b
31
b
32
b
13
12
b
23
<
br>
Y
2
1
1
b
33
Y
5
2
5
2
1
同样
B
2
对
P
1
, P
2
, P
3
作用的成对比较矩阵为:
21
1
B
1
<
br>
3
8
1
3
1
3<
br>1
1
8
1
B
1
3
3
Y
1
3
1
1
Y
3
3
3
1
134<
br>
11
B
4
1
11
B
5
1
1
3
1
44
11
4
解:
1
4
1
4
1
对以上每个比较矩阵都可计算出最大特征根
max
及
对象的特征向量
W
(即权重向量),并
进行一致性检验:
CIRI
CR
以
B
1
为例用“和法”求出
B
1
的
特征根
max
及对立的特征向量
W
1
25
1
B
1
0.512
0.20.51
0.5880.57
10.625
~
(S1)对
B
1
按列归一化
得:
B
1
W
ij
0.2940.2860.2
5
0.1180.1430.125
1.784
~
(S2)对按列归一化反向量再
按行求和:
W
W
ij
0.83<
br>
j
1
0.386
<
br>n
~
W
i
~
(S3)对
W
按行归一化得到特
征向量
W
W
n
~
W<
br>i
i
1
1.784
0.5
95
1.784
0.83
0.386
0.83
W
0.277
1.784
0.
83
0.386
0.386
0.129
1.784<
br>
0.83
0.386
1
(S4)计算特征根
max
(B)
max
1
BW
1
B
0.5
i
n
i
1
Wi
1
5
12
0.20.51
2
22
B)
(
max
1
0.595
0.595
0.595
1 2 50.2770.5 1 20.2770.2
0.5 10.277
0.129
0.129
0.129
1
3
0.5950.2770.12
9
1
0.595
0.554
0.645
0.2
98
0.277
0.258
0.119
0.139
0.129
<
br>
3
0.5950.2770.129
1
1.7940.8330.387
3
0.5950.2770.129
11
3.0
15
3.007
3
9.0223.00733
一致性检验:
CI
max
mn
1
RI
0.58
CI0.0035
CR
0.006
0.1
RI0.58
3.0
07
30.007
0.0035
0.1
3
12
故通过检验,既成对矩阵
B
1
可以接受。
同样步骤对
B
2
, B
3
,
B
4
, B
5
,对
P
1
,
P
2
, P
3
, P
4
, P
5
的影响
用特征向量
W
B
2
,
W
B
B
3
,
W
B
4
,
W
B
5
表示
35
24
,
max
,
max
,
max
最大特征根用:
max
表示
B
B
B
并分别计算一致性检
验指标:
CI
RI
B
2
3
CI
B
3
CI
B
4
CI
B
5
0.58 0.58
0.58
CR
列表如下:
权
准则层
值
决策层
(B
2
)
CI
CR
(B
3
)
CR
(B
4
)
CR
(B
5
)
RI
组合权向量
B1
B2 B3 B4 B5
0.262 0.474 0.055
0.099 0.102
W
i
a
j
b
ij
j
1
n
n
P
1
0.595
0.082 0.429 0.633 0.166
W
1
a
j
b
1j
0.299
j
1
23
P
2
0.277 0.236 0.429 0.193 0.166
W
2
a
j
b
2j
0.246
j
1
5
5
P
3
0.129
0.682 0.142 0.175 0.668
W
3
a
j
b
3j
0.456
j
1
max
CI
RI
CR
3.007 3.002 3 3.009 3
0.0035 0.001 0 0.005 0
0.58
0.58 0.58 0.58 0.58
0.006
n
其中
W
1
, W
2
,
W
3
的计算公式为:
W
i
ab
(i
1,
,n)
j
1
j
ij
0.595
0.082
5
W
1
a
j
b
1j
0.262, 0.474, 0.055, 0.099, 0.102
<
br>0.429
j
1
0.633
0.166
0.262
0.595
0.474
0.082.
0.055
0.429
0.099
0.633
0.1
02
0.166
0.156
0.039
<
br>0.024
0.063
0.017
0.299
0.595
0.082
<
br>
0.277 0.236 0.429 0.193 0.166<
br>
0.429
0.633
0.166
W
2
a
j
b
2j
j
1
5
<
br>0.246
0.595
0.082
5
W
3
a
j
b
3j
<
br>
0.129 0.682 0.142 0.175 0.668<
br>
0.429
j
1
<
br>0.633
0.166
0.4
56
WP
1
0.299
因此层次总排序:组合权向量为:
W
WP
2
0.246
WP
0.45
6
3
24
故最终决
策为
P
3
首选,
P
1
次之,
P
2
最后。
组合一致性检验:
m
j
CI
j
由
CR<
br>
a
j
1
m
可知:组合一致性检验结果为——层次总排序的一致性检验:
a
j
RI<
br>j
j
1
5
a
j
CI
3
C
R
j
1
5
a
j
RI
j
j
1
0.262
0.003
5
0.474
0.001
0.055
0
0.099
0.005
0.102
0
0.262
0.58
0.474
0.58
0.055
0.58
0.099
0.58
0.102
0.58
0
.0009
0.0005
0
0.0005
0
0.262
0.474
0.055
0.099
0.102
0.58
0.00190.0019
0.992
0.58
0.575
0.0033
0.1
故 一致性检验通过。
最层次总排序为
最决策为:
P
3
首选,
P
1
次之,
P
3
最后。
%层次分析法的matlab程序
disp('请输入判断矩阵A(n阶)');%
在屏幕显示这句话
A=input('A=');% 从屏幕接收判断矩阵
[n,n]=size(A);% 计算A的维度,这里是方阵,这么写不太好
x=ones(n,100);% x为n行100列全1的矩阵
y=ones(n,100);% y同x
m=zeros(1,100);%
m为1行100列全0的向量
m(1)=max(x(:,1));%
x第一列中最大的值赋给m的第一个分量
y(:,1)=x(:,1);%
x的第一列赋予y的第一列
x(:,2)=A*y(:,1);%
x的第二列为矩阵A*y(:,1)
m(2)=max(x(:,2));%
x第二列中最大的值赋给m的第二个分量
y(:,2)=x(:,2)m(2);%
x的第二列除以m(2)后赋给y的第二列
p=0.0001;i=2;k=abs(m(2)-m(1));%
初始化p,i,k为m(2)-m(1)的绝对值
while k>p% 当k>p是执行循环体
i=i+1;% i自加1
x(:,i)=A*y(:,i-1);%
x的第i列等于A*y的第i-1列
m(i)=max(x(:,i));%
m的第i个分量等于x第i列中最大的值
y(:,i)=x(:,i)m(i);%
y的第i列等于x的第i列除以m的第i个分量
25
k=abs(m(i)-m(i-1));% k等于m(i)-m(i-1)的绝对值
end
a=sum(y(:,i));% y的第i列的和赋予a
w=y(:,i)a;%
y的第i列除以a
t=m(i);% m的第i个分量赋给t
disp('权向量');disp(w);% 显示权向量w
disp('最大特征值');disp(t);% 显示最大特征值t
%以下是一致性检验
CI=(t-n)(n-1);%
t-维度再除以维度-1的值赋给CI
RI=[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];%
计算的标准
CR=CIRI(n);% 计算一致性
if CR<0.10
disp('此矩阵的一致性可以接受!');
disp('CI=');disp(CI);
disp('CR=');disp(CR);
else
disp('此矩阵的一致性不可以接受!');
end
26
层次分析法建模
层次分析法(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)机会多(如新单位或单位发展有后劲)等。
问题:现在有多个用人单位可供他选择,因此,他面临多种选择和决策,问题是他将如何
作出决
策和选择?——或者说他将用什么方法将可供选择的工作单位排序?
工作选择
贡献收入 发展 声誉 工作环境 生活环境
1
B.假期旅游地点选择 <
br>暑假有3个旅游胜地可供选择。例如:
P
1
:苏州杭州,
P
2
北戴河,
P
3
桂林,到底到哪个
地方去旅游最好?要作出决策和选择
。为此,要把三个旅游地的特点,例如:①景色;②费用;
③居住;④环境;⑤旅途条件等作一些比较—
—建立一个决策的准则,最后综合评判确定出一
个可选择的最优方案。
目标层
选择旅游地
景费居饮旅
准则层
色 用 住 食 途
方案层
P
1
P
2
P
3
C.资源开发的综合判断
7种金属可供开发,开发后对国家贡献可以通过两两比较
得到,决定对哪种资源先开发,
效用最用。
对经济发展、贡献U
经济价值开採费 风险费 要求量 战略重要性 交通条件
铁In 铜Co 磷酸盐 钿Ur 铝Al
金Go
可供选择的单位P
1’
P
2 ‘
-----
P
n
二、问题分析:
例如旅游地选择问题:一般说来,此决策问题可按如下步骤进行:
2
(S1)将决策解分解为三个层次,即:
目标层:(选择旅游地)
准则层:(景色、费用、居住、饮食、旅途等5个准则)
方案层:(有
P
1
,
P
2
,
P
3
三个选择地点)
并用直线连接各层次。
(S2)互相比较各准则对目标的权重,各方案对每一个准则的权重。
这些权限重在人的思维过
程中常是定性的。
例如:经济好,身体好的人:会将景色好作为第一选择;
中老年人:会将居住、饮食好作为第一选择;
经济不好的人:会把费用低作为第一选择。
而层次分析方法则应给出确定权重的定量分析方法。
(S3)将方案后对准则层的权重,及准则后对目标层的权重进行综合。
(S4)最终得出方案层对目标层的权重,从而作出决策。
以上步骤和方法即是AHP的决策分析方法。
三、确定各层次互相比较的方法——成对比较
矩阵和权向量
在确定各层次各因素之间
的权重时,如果只是定性的结果,则常常不容易被别人接受,因
而Santy等人提出:一致矩阵法
.....
即:1. 不把所有因素放在一起比较,而是两两相互比较
2.
对此时採用相对尺度,以尽可能减少性质不同的诸因素相互比较的困难,提高准确度。
因素比较方法 —— 成对比较矩阵法:
目的是,要比较某一层
n
个因素<
br>C
1
,C
2
,
, C
n
对上一
层因素O的影响(例如:旅游决策解
中,比较景色等5个准则在选择旅游地这个目标中的重要性)。 <
br>採用的方法是:每次取两个因素
C
i
和
C
j
比较其对
目标因素O的影响,并用
a
ij
表示,全部
比较的结果用成对比较矩阵表示,
即:
A(a
ij
)
nxn
,
a
ij
0, a
ji
1
(或a
ij
a
ij
1)
(1)
a
ij
1
a
ji
由于上述成对比较矩阵有特点:
A(a
ij
)
, a
ij
0,
a
ij
故可称
A
为正互反矩阵:显然,由
a
ij
1
,即:
a
ij
a
ji
1
,故有:
a
ji
1
a
ji
3
例如:在旅游决策问题中:
(O的重要性为1
C
C(景色)
1
景色)对目标
a
12
1
=
1
表示:
2
C(费用)
O的重要性为2
2
2
C(费用)对目标
1,费用重要性为2)
故:
a
12
1
(即景色重要性为
2
O的重要性为4
C<
br>1
(景色)对目标
C(景色)
1
4
a
13
4
=
表示:
1
C(居
住条件)
C(居住条件)对目标O的重要性为1
3
3
即:景色为4
,居住为1。
O的重要性为7
C
2
(费用)对目标
C(
费用)
2
a
23
7
7
=
表示:
1
C(居住条件)
O的重要性为1
3
3
C(居住条件)对目标
即:费用重要性为7,居住重要性为1。
1
2
1
因此有成对比较矩阵:
A
4
1
3
1
3
1
2
1
1
7
1
5
1
5
4
7
1
2
3
33
55
11
23
11
11
??问题:稍加分析就发现上述成对比较矩阵的问题:
①
即存在有各元素的不一致性,例如:
既然:
a
12
C
1
1
C
411
a
21
2;
a
13
1
a
31
C
2
2C
3
1a
13
4
所以应该有:
a23
C
2
a
21
C
1
2
8
8
1
1
C
3
a31
C
3
4
C
1
C
2
而不应为矩阵<
br>A
中的
a
23
7
1
②成对比较矩阵比较的次数要求太 ,因:
n
个元素比较次数为:
C
n
2
n(n
1)
次,
2!因此,问题是:如何改造成对比较矩阵,使由其能确定诸因素
C
1
,
, C
n
对上层因素O的
权重?
对此Saoty提出了:在成对比较出现不一致情况下,计算各因素
C
1
,
,
C
n
对因素(上层
因素)O的权重方法,并确定了这种不一致的容许误差范围。
为此,先看成对比较矩阵的完全一致性——成对比较完全一致性
4
四:一致性矩阵
Def
:
设有正互反成对比较矩阵:
a
W
1
W
11
<
br>W
1 a
12
1
1 ,
, a
W
1
1n
1
W
2
W
n
a
W
2W
2
W
2
21
a
22
1
,
A
W
1
W
,
a
2n
2
W
n
a
W
(4)
i
ij
W
j
a
W
n
W
n
W
n
n1
W
a
n2
a
nn
1
1
W
2W
n
除满足:(i)正互反性:即
a
1
ij
0 a
ij
a
( 或 a
ij
a
ji
1)
ji
而且还满足:(ii)一致性:即
a
a
i
ij
a
a
ah
i
k
a
kj
i
i, j1, 2, n
有点点
错误
j
a
j
h
则称满足上述条件的正互反对称矩阵A为一致性矩阵,
简称一致阵。
一致性矩阵(一致阵)性质:
性质1:
A
的秩
Rank(A)=1显然
A
的唯一非0的特征根为n
性质2:
A
的任一列(行)向量都是对应特征根
n
的特征向量:
即有(特征向量、特征值):
W
1
W
1
WW
W
1
W
<
br>W
1
W
2
W
n
A
22
2
W
W
1
1
W
W
W
2
,则向量
W
W<
br>n
W
2
n
W
n
<
br>
n
W<
br>1
W
2
W
n
W
3<
br>
5
W
1
W
1
W
1
W
2
满足:<
br>AW
W
n
W
n
WW
2
1
即:
(AnI)W0
W
1
W
nW
1
1
W
n
W
nW
<
br>
2
2
nW
W
n
W
nW
W
n
n
n
我的理解:通过A(变换A与W中的
元素
有关)变换将一致W矩阵变成权向量
W(特征向量),如果正互反矩阵W’接近一致矩
阵,同样
的道理变换A可以将W’变成权向量
(这里的权向量与W’稍有不同)
W
1
W
2
W
启发与思
考:既然一致矩阵有以上性质,即
n
个元素
W
1,
W
2,
W
3 ,
…
W
n
构成的向量
W
n
是一致矩阵A的特征向量,则可以把向量W归一化后的向量
,看成是诸元素
W
1,
W
2,
W
3 ,
…
W
n
目标的权向量,因此,可以用求
A
的特征根和特征
向量的办法,求出元素
W
1,
W
2,
W
3 ,
…
W
n
相对于目标O的劝向量。
解释:
一致矩阵即:
n
件物体
M
1
, M
2
,
,M
n
,它们重量分别为
W
1
,
W
2
,
,W
n
,将他们两两比较
W
1
W
2
重量,其比值构成一
致矩阵,若用重量向量
W
右乘
A
,则
W
n
6
A的特征根为n,
W
1
W
2
:重量向量
W=
,则归一化后的特征向量
以n为特征根的特征向量为
W
:
n
W
1
W=
W
i
=1
,就表示诸因素C
1
,C
2
,
,C
n
对上层因素O的权重,即为
W
权向量,此种用特征向量求权向量的方法 称特征根法,
分析:
W
1
W
2
若重量
向量
W
未知时,则可由决策者对物体
M
1
,
M
2
,
,M
n
之间两两相比关系,
W
n
主观作出比值的判
断,或用Delphi(调查法)来确定这些比值,使
A
矩阵(不一定有一致性)
为已
知的,并记此主观判断作出的矩阵为(主观)判断矩阵
A
,并且此
A
(不一致
)在不一致
的容许范围内,再依据:
A
ij
的特征根或和特征向量
W
连续地依赖于
ij
矩阵的元素
a
,即当
a
离一致性
的要求不太远时,
A
的特
征根
i
和特征值(向量)
W
与一致矩阵
A
的特征根
和特征
向量
W
也相差不
大的道理:由特征向量
W
求权向量
W
的方
法即为特征向量法,并由此
引出一致性检查的方法。
问题:Remark
以上讨论的用求特征根来求权向量
W
的方法和思路,在理论上应解决以下问题:
1. 一致阵的性质1是说:一致阵的最大特征根为
n
(即必要条件),但用特征根来
求特征向量
时,应回答充分条件:即正互反矩阵是否存在正的最大特征根和正的特征向量?且如果正互<
br>
7
反矩阵
A
的最大特征根
ma
x
n
时,
A
是否为一致阵?
2. 用主观判断矩阵
A<
br>的特征根
和特征向量
W
连续逼近一致阵
A
的特征根
和特征向量
W
时,即:
由
lim
k
k
k
得到:
limW
k
W
k
即:
limA
k
A
k
是否在理论上有依据。
3.一般情况下,主观判断矩阵
A<
br>在逼近于一致阵
A
的过程中,用与
A
接近的
A*
来代
替
A
,
即有
A*A
,这种近似的替代一致性矩阵
A
的作法,就导致了产生的偏差估计问题,即一致
性检验问题,即要确定一种一致性检验判断指标,由此
指标来确定在什么样的允许范围内,
主观判断矩阵是可以接受的,否则,要
两两比较构造主观判断矩阵。此问题即一致性检
验问题的内容。
以上三个问题:前两个问题由
数学严格比较可获得(见教材P325,定理1、定理2)。第
3个问题:Satty给出一致性指标(
TH1,TH2介绍如下:)
附:
Th1:(教材P326,perronTh 比隆
1970 )对于正矩阵
A
(
A
的所有元素为正数)
(1)
A
的最大特征根是正单根
;
(2)
<
br>对应正特征向量
W
(
W
的所有分量为正数)
1<
br>
k
Ae
1
(3)
lim
Tk
W
其中:
e
为半径向量
,
W
是对应
的归一化特征向量
k
eAe<
br>
1
证明:(3)可以通过将
A
化为标准形证明
Th2:
n
阶正互反阵A的最大特征根
n
;
当
n
时,
A
是一致阵
8
五、一致性检验——一致性指标:
1.一致性检验指标的定义和确定——
CI
(平均值)的定义:
当人们对
复杂事件的各因素,采用两两比较时,所得到的主观判断矩阵
A
,一般不可直接
保证正
互反矩阵
A
就是一致正互反矩阵
A
,因而存在误差(及误差估计问题)。这种
误差,必
然导致特征值和特征向量之间的误差
(
)及
W-W
。此时就导致问题
AW=
max
W
与问
题
AWnW
之间的差别。(上述问题中
max
是主观判断矩阵A
的特征值,
W
是带有偏差的相
对权向量)。这是由判断矩阵不一致性所
引起的。
因此,为了避免误差太大,就要衡量主观判断矩阵
A
的一致性。
因为:
①当主观判断矩阵
A
为一致阵
A
时就有:
=
k
k
1n
1
nn
k
a
k
k
1
n
A
为一致阵时有:
a
ii
1
(a[ii]为对角线k
1k
1
nn
上的值,按照一致性矩阵的理解,它
应该为1)
此时存在唯一的非
max
n
(由一致阵性质1:R
ark(4)=1,
A
有唯一非O最大特征根且
max
n
)
②当主观判断矩阵
A
不是一致矩阵时,此时一般有:
max
n
(Th2)
此时,应有:
max
即:
max
n
k
a
h
max
k
ii
n
(不大理解)
k
max
所以,可以取其平均值作为检验主观判断矩阵的准则,一致性的指标,
即:
C
I
显然:
(1)
当
max
n
时,有:
CI0
,
A
为完全一致性
max
n
n
1
k
max
k
n
1
9
(2)
CI值越大,主观判断矩阵
A
的完全一致性越差,即:
A
偏离
A越远(用特征向量
作为权向量引起的误差越大)
(3) 一般
CI01<
br>,认为主观判断矩阵
A
的一致性可以接受,否则应重新进行两两比
较,构造主观
判断矩阵。
2.随机一致性检验指标——
RI
问题:实际操作时发现:
主观判断矩阵
A
的维数越大,判断的一致性越差,故应放宽对高维矩
阵的一致性要求。
于是引入修正值
RI
来校正一致性检验指标:即定义
RI
的修正值表为:
A
的维数
1 2 3
4 5 6 7 8 9
0.00 0.00
0.58 0.96 1.12 1.24 1.32 1.41 1.45
RI
并定义新的一致性检验指标为:
C
R
C
I
R
I
随机一致性检验指标——
RI
的解释:
为确定
A
的不一致程度的容许范围,需要确定衡量
A
的一致性指示
C
I
的标准。于是Satty
又引入所谓随机一致性指标
RI
,其定义和计
算过程为:
(i
j)
从1~9和1~
1
中随机取值,①
对固定的
n
,随机构造正互反阵
A
,其元素
a
i
j
9
与
a
ji
的互反性,即:
aij
1
且满足
a
ij
a
ji
1
. ,且
a
ii
② 然后再计算
A<
br>
的一致性指标
CI
,因此
A
是非常不一致的,
此时,
CI
值相当大.
③ 如此构造相当多的
A
,再
用它们的
CI
平均值作为随机一致性指标。
④ Satty对于不同的
n
(n1
~11),用100~500个样本
A
计算出上表所列出的随机一
致性
指标
RI
作为修正值表。
3.
一致性检验指标的定义——一致性比率
CR
。
由随机性检验指标
CR
可知:
当
n1,
2
时,
RI0
,这是因为1, 2阶正互反阵总是一致阵。
对于
n3
的成对比较阵
A
,将它的一致性指标
CI
与同阶(指n
相同)的随机一致性指
标
RI
之比称为一致性比率——简称一致性指
标,
即有: 一致性检验指标的定义——一致性比率
10
定义:
C
R
C
IC
I
:
C
R
R
IR
I
C
I
当:
C
R
01
时,认为主观判断矩阵
A
的不一致程度在容许范围之内,R
I
可用其特征向量作为权向量。否则,对主观判断矩阵
A
重
新进行成对比较,构重新的主观
判断矩阵
A
。
注:上式
C
R
C
I
01
的选取是带有一定主观信度的。
R
I
六、标度——比较尺度解:
在构造正互反矩阵时,当比较两
个可能是有不同性质的因素
C
i
和
C
j
对于上层因素O的影响时,採用什么样的相对刻度较好,即
a
ij
的元素的值在(1~9)或(1
~
1
)或更多
9
的数字,Satty提出用1~9尺度最好,即
a<
br>ij
取值为1~9或其互反数1~
1
,心理学家也
9
提出:人
们区分信息等级的极限解能力为
7
±2。可见对
nn
阶矩阵,只需作出判断值即可
标度
a
ij
1
3
5
7
9
2,4,6,8,
倒数1,
定
义
n(n1)
个
2
因素
i
与因素
j
相
同重要
因素
i
比因素
j
稍重要
因素
i
比因素
j
较重要
因素
i
比因素
j
非常重要
因素
i
比因素
j
绝对重要
因素
i
与因素
j
的重要性的比较值介于上述两
个相邻等级之间
因素
j
与
因素
i
比较得到判断值为
a
ij
的互反
数,
111
11111
, , , , , , ,
23456789
a
ji
1
a
ii
1
a
ij
注:以上比较的标度Satt
y曾用过多种标度比较层,得到的结论认为:1~9尺度不仅在较简单
的尺度中最好,而且比较的结果并
不劣于较为复杂的尺度。Satty曾用的比较尺度为:
① 1~3, 1~5, 1~6,…,
1~11,以及
②
(d0.1)
~
(d0.9)
,其中
d1, 2, 3,
4
11
③
1
p
~
9
,其中
P2, 3, 4,
5
…
等共27种比较尺度,对放在不同距离处的光源亮度进行比较判断,并构造出成对比较矩
阵,计
算出权向量。同时把计算出来的这些权向量与按照物理学中光强度定律和其他物理知识得到的实际权向量进行对比。结果也发现1~9的比较标度不仅简单,而效果也较好(至少不比其他更
复杂
的尺度差)
因而用1~9的标度来构造成对比较矩阵的元素较合适。
P
七、组合权向量的计算——层次总排序的权向
量的计算
层次分析法的基本思想:
(1)
计算出下一层每个元素对上一层每个元素的权向量
W
def:层次总排序,计算同一层次所有元素对最高层相对重要性的排序权值。
当然要先:①构造下一层每个元素对上一次每个元素的成对比较矩阵
②计算出成对比较矩阵的特征向量(和法,根法,幂法)
③由特征向量求出最大特征根
max
(由和法,根法,幂法求得)
④用最大特征根
max
用方式
C
I
max
n
n
1
及
C
R
C
R
对成对比较矩
R
I
阵进行一致性检,并通过。
(2) 并把下层每个元素对上层每个元素的权向量按列排成以下
表格形式:例,假定:上
层
A
有
m
个元素,
A
1<
br>, A
2
,
,
A
m
,且其层次总排序权向量为
a
1
, a
2
,
, a
m
,
下层
B
有
n
个元
素
B
1
, B
2
,
,
B
n
,则按
B
j
对
A
i
个元素的单排
序权向量的列向
量为
b
ij
,即有:
层次
A
1
A
1
…
A
1
m
(权向
B
层总是排序权重
量、列向量)
12
a
1
a
2
…
a
m
B
1
B
2
B
n
b
11
b
12
b
1m
W1
a
j
b
1j
W
2
<
br>
a
j
b
2j
j
1
m
j
1
m
m
b
12
b
22
b
2m
b
n1
b
n2
b
nm<
br>
W
n
a
j
b<
br>nj
j
1
max
CI
计算出最大特根(方法:和法、根法、幂法)
一致性检验
C
I
一致性检验比率
max
n
n
1
CI
C
I
j
C
R
R
I
m
a
jCI
3
a
j
m
j
检验
CR01
否?
RI
j
注:①若下层元
素
B
k
与上层元素
A
j
无关系时,取
b
k
j
0
②总排序权向量各分量的计算公式:
W
i
<
br>
ab
j
1
m
jij
(i
1,
,n)
(3)
对层次总排序进行一致性检验:从高层到低层逐层进行,如果
如果
B
层次某些元素对
A
j
单的排序的一致性指标为
CI
j
,相应的平均随机一致
性指标为
RI
j
,则
B
层总排序随机一致性比率为:
C
R
aCI
j
1
m
j
m
j
a
j
1
j
R
I
j
当
CR01
时,认为层次总排序里有满意的一致性,否则应重新调整
判断矩阵的元素取值。
八、层次分析法的基本步骤:
(S1)建立层次结构模型
将有关因素按照属性自上而下地分解成若干层次:
同一层各因素从属于上一层因素,或对上层
因素有影响,同时又支配下一层的因素或受
到下层因素的影响。
13
最上层为目标层(一般只有一个因素),最下层为方案层或对象层决策层,中间可以有
1个 或几个层次,通常为准则层或指标层。
当准则层元素过多(例如多于9个)时,应进一步分解出子准则层。
(S2)构造成对比较矩 阵,以层次结构模型的第2层开始,对于从属于(或影响及)上一
层每个因素的同一层诸因素,用成对比 较法和1~9比较尺度构造成对比较矩阵,直到最下
层。
(S3)计算(每个成对比较矩阵的)权向量并作一致性检验
① 对每一个成对比较矩阵计算 最大特征根
max
及对应的特征向量(和法、根法、幂法等)
W
1
W
W
n
② 利用一致性指标
CI
,随机一致性指标
CR
和一致性比率作一致性检验
CR
C
I
R
I
W
1
③ 若通过检验(即
CR0.1
,或
CI0.1
)则将上层出权向量
W
归一化之后
W
n
作为(
B
j
到
A
j
)的权向量(即单排 序权向量)
④ 若
CR0.1
不成立,则需重新构造成对比较矩阵
(S4)计算组合权向量并作组合一致性检验——即层次总排序
W
1
① 利用单层权向量的权值
W
j
j
1,
, m
构组合权向量表:并计算出特征根,组
W
n
合特征向量,一致性
14
上
单 层
层 重
权
量
向
下层 量
层次
A
1
A
1
…
A
1
m
W
1
W
计算组合权向量
W
n
其中
W<
br>i
m
a
1
a
2
…
a
m
aW
j
1
j
m
ij
B
1
B
2
B
n
最大特征根
max
(i)W
11
W
12
W
1m
W
1
a
j
b
1j
W
2
a
j
b
2j
j
1
m
j
1
m
W
12
W
22
W
2m
W
n1
W
n2
W
nm
W
n
a
j
b<
br>nj
j
1
和法、根法、幂法
一致性检验
CI
CI
j
(j)max
n
n
1
CI0.1
?
一致性随机检验
RI
RI
j
对照表
一致性比率
CR
CI
CR
jRI
m
a
j
CI
j
CR01
?
a
j
m
j
RI
2j
W
1
W
1
② 若通过
一致性检验,则可按照组合权向量
W
的表示结果进
行决策(
W
W
W
n
n
中
W
i
中最大者的最优),即:
W*
maxW:W
i
W
1
,
,W
n
T
③ 若未能通过检验,则需重新考虑模型或重新构造那些一致性比率,
CR
较大
的成对比较
矩阵
九、特征根的近似求法(实用算法)
层次分析法的基本思路是计算上层每个元素对下一层次各元素的权向量(即最大特征根
15
W
1
m
ax
对应的特征向量
W
),以及组合权向
量及一致性检验问题。
W
n
计算判断矩
阵最大特征根和对应阵向量,并不需要追求较高的精确度,这是因为判断矩阵
本身有相当的误差范围。而
且优先排序的数值也是定性概念的表达,故从应用性来考虑也希望
使用较为简单的近似算法。常用的有以
下求特征根的近似求法:“和法”、“根法”、“幂法”,具
体如下:
1.“和法”求最大特征根和对应特征向量(近似解)
a
ij
~
(
S1)将矩阵
A(a
ij
)
nxm
的每一列向量的归一化得:W
ij
n
(利用数据验证即为:每个
a
ij
i
1
位置的数除以该列的合计)
n
~~
~
(S2)对
W
ij
按行求和得:
W
i
W
ij
j
1
(S3)将
W
i
归一化,即有:
W
i
~
~
W
i
W
1
~
,则有特征向量:
W
n
~
W
W
i
n
i
1
W
1
1
n
(AW)
i
(S4)计算与
特征向量
W
对应的最大特征根
max
的近似值:
max
nW
i
1
i
W
n
此
方法:实际上是将
A
的列向量归一化后取平均值作为
A
的特征向量。
解释:
当
A
为一致矩阵时,它的每一列向量都是特征向量
W
可以在
A
的不一致性不严重时,取
A
的列向量(归
一化后)的平均值作为近似特征向
量是合理的(有依据的)。
2.“根法”求最大特征根特征向量近似值:
步骤与“和法”相同,只是在(S2)时:对归
一化后的列向量按行“求和”改为按行“求
~
~
积”
再取
n
次方根,即:
W
i
W
ij
。
j
1
n
1
n
即有具体步骤:
16
(S1)将矩阵
A(a
ij
)
min
的每一列向量归一化得:
W
ij<
br>
~
a
ij
a
i
1
n
ij
(S2)对归一化以后的列向量各元素:
W
ij
~
a
ij
a
i
1
n
ij
n
~
~
按行“求和”并开
n次方根得:
W
i
W
ij
j
1
1
n
<
br>n
~
W
ij
~
W
i
j
1
~
(S3)再将<
br>W
i
归一化得:
W
i
n
1~
n
n
W
i
n
~
W
ij
i
1
i
1
j
1
W
1
W
2
W
得到特征向量近似值:
W
n
(S4)计算最大特征根:
max
<
br>1
n
1
(AW)
i
作为最大特征根的近似值。
nW
i
注:“根法”是将“和法”中求列向量的算术平均值改为求几何平均值。
3.“幂法”求最大特征根:
(S1)任取
n
维归一化初始向量
W
(S2)计算
W
(0)
~
(k
1)
AW
(k)
,
k
0, 1, 2,
~
(k
1)
~
(
k
1)
(k
1)
W
(S
3)
W
归一化,即令:
W
~
(k
1)
W
i
i
1
n
(S4)对预先给定的
,当
W
i
量;否则返回(S2)
(S5)计算最大特征根,
max
(k
1)
W
i
(k)
(i
1, 2,
, n)
时,
W
(k1)
即为所求的特征向
~
1
n
W
i
(k
1)
<
br>(k)
n
i
1
W
i
17
以上用幂法求最大特征根
max
对应特征向量的
迭代方法,其收敛性由TH1(教材P325)中
1
Ae<
br>的3)
lim
Tk
W
,其中
e
,
W
是对应
max
的归一化向量
特征。
k
eAe
1
k
(证明:可以将
A
化为标准形证明)保证。
W
(0)
任意选取,也可以取由“根法”、“和法”得
W
1
到的
W
W
2
W
n
注:在以上求特征根和特向量的方法中“和法”最简单。
例:在旅游问题中,求目标层到准则层的成对比较矩阵为
A
的特征向量和最大特征根:
选择旅游地
准则层:
景费居饮旅
色 用 住 食 途
方案层:
P
1
P
2
P
3
1
1
2
433
<
br>21755
11
10.5433
11
21755
A
4
1
23
10.50.333
<
br>
11
7
11
=
0.250.143
2
3
0.
3330.2211
0311
1
5
1
35
311
.3330.
2
18
1.312<
br>
2.37
AW
~~
i
0.273
W
0.493
0.511
W1
利用“和法”求
A
的特征向量
W
和特征根
max
W
n
(S1)将
A
Wij
nxn
的元素按列归一化得:
0.265
0.2450.2350.2860.29
0.4890.4110.4760.484
A
W
~
ij
<
br>0.510
nxn
0.0640.0700.0590.0480
.032
0.0850.0980.1180.0950.
097
0.0850.0980.1760.0950.097
2
1
2
0.25
0.333
0.333
3.917
2
0.5
1
0.143
0.2
0.2
2.043
3
4
7
1
2
3
17
各列归一化的分母
4
3
5
0.5
1
1
10.5
5
3
5
0.333
1
1
10.333
(S2)将
A
W
~<
br>~
~
n
~
ij
nxn
中元素
W<
br>ij
按行求和得各行元素之和:
W
i
W
ij
j
1
(S3)再将上述矩阵向量归一化得到特征向量近似值,
1.312
0.262
2.37
0.474
W
W
~
i
1
n
0.273
0.
W
4.999
.493
05
5
0.099
特征向量
i
0
i
1
0.511
0.102
5
其中
W
~
i
(1.312
2.37
0
.273
0.493
0.511)
4.999
1
(S4)计算与特征向量相对应最大特征根(的近似值)
max<
br>
1
n
AW
i
n
i
1
W
i
n
nnn
a
1
iaW
2j
W
i
a
3
j
W
i
a
4j
W
i
1j
ii
j
1i
j
1i
j
1
a
5j
W
i
5
W
i
j
1
1
W
2
W
3
W
4
W
5
19
0.262
0.262
0.262
0.4740.4740.4
74
1 0.5 4
33
0.055
21 7 5
5
0.055
0.25 0.1431
0.5 0.333
0.055
0.099
0.099
0.099
0.102
0.102
0.10
2
1
5
0.262
0.474
0.055
0.262
0.474
0.262
0.474
0.337 0.2
2 11
0.055
0.333 0.2
3 11
0.055
0.09
9
0.099
0.102
0.102
0.099
0.102<
br>
<
br>1
0.263
0.237
0.22
0.297
0.3060.524
0.474
0.385
0.495
0.
5
0.262
5
0.474
0.066
0.068
0.055
0.0495
0
.0340.087
0.095
0
0.055
.11
0.099
0.102
0.099
0.087
0.095
0.165
0.09
9
0.102
0.102
1<
br>
1.3232.3880.2730.4930
5
0.
262
0.474
0.055
0.099
.548
0.102
1
5
5.05
5.038
4.960
4.9
8
5.373
1
5
25.401
5.0802
0.262
0.
474
故有最大特征根
max
5.080
2 , W
0.055
0.099
0.102
20
对
A
一致性检验指标:
CI
max
n
n
1
5.0802<
br>
50.0802
0.02
44
RI
1.12
CR
故通过检验。
0.02
0.018
0.1
1.
12
十、应用实例
对前面旅游问题进行决策
目标层:
A
选择旅游地点
0.262 0.474 0.099 0.102
0.055
景费居饮
色 用 住 食
准则层:
B
1
B
2
B
3
B
4
0.595 0.129 0.129
0.277
决策层:
P
1
P
2
P
3
已知:①目标
A
对准则
B
i
i1, 2, 3,
4, 5
的权重向量为:
T
,并已通过一致性检
W
0.262 0.474
0.055 0.099 0.102
(由前面已算出)
旅
途
B
5
验。
②准则
B
1
,
B
2
, B
3
, B
4
,
B
5
相对于
P
1
, P
2
,
P
3
的成对比较矩阵为
B
1
对
P
1
,
P
2
, P
3
作用的成对比较矩阵为:
b
1
1
b
12
B
1
b
21b
22
b
31
b
32
b
13
12
b
23
<
br>
Y
2
1
1
b
33
Y
5
2
5
2
1
同样
B
2
对
P
1
, P
2
, P
3
作用的成对比较矩阵为:
21
1
B
1
<
br>
3
8
1
3
1
3<
br>1
1
8
1
B
1
3
3
Y
1
3
1
1
Y
3
3
3
1
134<
br>
11
B
4
1
11
B
5
1
1
3
1
44
11
4
解:
1
4
1
4
1
对以上每个比较矩阵都可计算出最大特征根
max
及
对象的特征向量
W
(即权重向量),并
进行一致性检验:
CIRI
CR
以
B
1
为例用“和法”求出
B
1
的
特征根
max
及对立的特征向量
W
1
25
1
B
1
0.512
0.20.51
0.5880.57
10.625
~
(S1)对
B
1
按列归一化
得:
B
1
W
ij
0.2940.2860.2
5
0.1180.1430.125
1.784
~
(S2)对按列归一化反向量再
按行求和:
W
W
ij
0.83<
br>
j
1
0.386
<
br>n
~
W
i
~
(S3)对
W
按行归一化得到特
征向量
W
W
n
~
W<
br>i
i
1
1.784
0.5
95
1.784
0.83
0.386
0.83
W
0.277
1.784
0.
83
0.386
0.386
0.129
1.784<
br>
0.83
0.386
1
(S4)计算特征根
max
(B)
max
1
BW
1
B
0.5
i
n
i
1
Wi
1
5
12
0.20.51
2
22
B)
(
max
1
0.595
0.595
0.595
1 2 50.2770.5 1 20.2770.2
0.5 10.277
0.129
0.129
0.129
1
3
0.5950.2770.12
9
1
0.595
0.554
0.645
0.2
98
0.277
0.258
0.119
0.139
0.129
<
br>
3
0.5950.2770.129
1
1.7940.8330.387
3
0.5950.2770.129
11
3.0
15
3.007
3
9.0223.00733
一致性检验:
CI
max
mn
1
RI
0.58
CI0.0035
CR
0.006
0.1
RI0.58
3.0
07
30.007
0.0035
0.1
3
12
故通过检验,既成对矩阵
B
1
可以接受。
同样步骤对
B
2
, B
3
,
B
4
, B
5
,对
P
1
,
P
2
, P
3
, P
4
, P
5
的影响
用特征向量
W
B
2
,
W
B
B
3
,
W
B
4
,
W
B
5
表示
35
24
,
max
,
max
,
max
最大特征根用:
max
表示
B
B
B
并分别计算一致性检
验指标:
CI
RI
B
2
3
CI
B
3
CI
B
4
CI
B
5
0.58 0.58
0.58
CR
列表如下:
权
准则层
值
决策层
(B
2
)
CI
CR
(B
3
)
CR
(B
4
)
CR
(B
5
)
RI
组合权向量
B1
B2 B3 B4 B5
0.262 0.474 0.055
0.099 0.102
W
i
a
j
b
ij
j
1
n
n
P
1
0.595
0.082 0.429 0.633 0.166
W
1
a
j
b
1j
0.299
j
1
23
P
2
0.277 0.236 0.429 0.193 0.166
W
2
a
j
b
2j
0.246
j
1
5
5
P
3
0.129
0.682 0.142 0.175 0.668
W
3
a
j
b
3j
0.456
j
1
max
CI
RI
CR
3.007 3.002 3 3.009 3
0.0035 0.001 0 0.005 0
0.58
0.58 0.58 0.58 0.58
0.006
n
其中
W
1
, W
2
,
W
3
的计算公式为:
W
i
ab
(i
1,
,n)
j
1
j
ij
0.595
0.082
5
W
1
a
j
b
1j
0.262, 0.474, 0.055, 0.099, 0.102
<
br>0.429
j
1
0.633
0.166
0.262
0.595
0.474
0.082.
0.055
0.429
0.099
0.633
0.1
02
0.166
0.156
0.039
<
br>0.024
0.063
0.017
0.299
0.595
0.082
<
br>
0.277 0.236 0.429 0.193 0.166<
br>
0.429
0.633
0.166
W
2
a
j
b
2j
j
1
5
<
br>0.246
0.595
0.082
5
W
3
a
j
b
3j
<
br>
0.129 0.682 0.142 0.175 0.668<
br>
0.429
j
1
<
br>0.633
0.166
0.4
56
WP
1
0.299
因此层次总排序:组合权向量为:
W
WP
2
0.246
WP
0.45
6
3
24
故最终决
策为
P
3
首选,
P
1
次之,
P
2
最后。
组合一致性检验:
m
j
CI
j
由
CR<
br>
a
j
1
m
可知:组合一致性检验结果为——层次总排序的一致性检验:
a
j
RI<
br>j
j
1
5
a
j
CI
3
C
R
j
1
5
a
j
RI
j
j
1
0.262
0.003
5
0.474
0.001
0.055
0
0.099
0.005
0.102
0
0.262
0.58
0.474
0.58
0.055
0.58
0.099
0.58
0.102
0.58
0
.0009
0.0005
0
0.0005
0
0.262
0.474
0.055
0.099
0.102
0.58
0.00190.0019
0.992
0.58
0.575
0.0033
0.1
故 一致性检验通过。
最层次总排序为
最决策为:
P
3
首选,
P
1
次之,
P
3
最后。
%层次分析法的matlab程序
disp('请输入判断矩阵A(n阶)');%
在屏幕显示这句话
A=input('A=');% 从屏幕接收判断矩阵
[n,n]=size(A);% 计算A的维度,这里是方阵,这么写不太好
x=ones(n,100);% x为n行100列全1的矩阵
y=ones(n,100);% y同x
m=zeros(1,100);%
m为1行100列全0的向量
m(1)=max(x(:,1));%
x第一列中最大的值赋给m的第一个分量
y(:,1)=x(:,1);%
x的第一列赋予y的第一列
x(:,2)=A*y(:,1);%
x的第二列为矩阵A*y(:,1)
m(2)=max(x(:,2));%
x第二列中最大的值赋给m的第二个分量
y(:,2)=x(:,2)m(2);%
x的第二列除以m(2)后赋给y的第二列
p=0.0001;i=2;k=abs(m(2)-m(1));%
初始化p,i,k为m(2)-m(1)的绝对值
while k>p% 当k>p是执行循环体
i=i+1;% i自加1
x(:,i)=A*y(:,i-1);%
x的第i列等于A*y的第i-1列
m(i)=max(x(:,i));%
m的第i个分量等于x第i列中最大的值
y(:,i)=x(:,i)m(i);%
y的第i列等于x的第i列除以m的第i个分量
25
k=abs(m(i)-m(i-1));% k等于m(i)-m(i-1)的绝对值
end
a=sum(y(:,i));% y的第i列的和赋予a
w=y(:,i)a;%
y的第i列除以a
t=m(i);% m的第i个分量赋给t
disp('权向量');disp(w);% 显示权向量w
disp('最大特征值');disp(t);% 显示最大特征值t
%以下是一致性检验
CI=(t-n)(n-1);%
t-维度再除以维度-1的值赋给CI
RI=[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];%
计算的标准
CR=CIRI(n);% 计算一致性
if CR<0.10
disp('此矩阵的一致性可以接受!');
disp('CI=');disp(CI);
disp('CR=');disp(CR);
else
disp('此矩阵的一致性不可以接受!');
end
26