层次分析法步骤解析—根法、和法、幂法
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
层次分析法(AHP)
AHP(Analytic Hierarchy Process)方法,是由20世纪70年代由美国著名运筹学学家T.L.Satty提出的。
它是指将决策问题的有关元素分解成目标、准则、方案等层次,在此基础上进行定性分析和定量分析的一种决策方法。
这一方法的特点,是在对复杂决策问题的本质、影响因素及其内在关系等进行深入分析之后,构建一个层次结构模型,然后利用较少的定量信息,把决策的思维过程数学化,从而为求解多准则或无结构特性的复杂决策问题提供了一种简便的决策方法。
AHP十分适用于具有定性的,或定性定量兼有的决策分析。
这是一种十分有效的系统分析和科学决策方法,现在已广泛地应用在企业信用评级、经济管理规划、能源开发利用与资源分析、城市产业规划、企业管理、人才预测、科研管理、交通运输、水资源分析利用等方面。
一、递阶层次结构的建立
一般来说,可以将层次分为三种类型:
(1)最高层:只包含一个元素,表示决策分析的总目标,因此也称为总目标层。
(2)中间层:包含若干层元素,表示实现总目标所涉及的各子目标,包含各种准则、约束、策略等,因此也称为目标层。
(3)最低层:表示实现各决策目标的可行方案、措施等,也称为方案层。
典型的递阶层次结构如下:
总目标
m
一个好的递阶层次结构对解决问题极为重要,因此在建立递阶层次结构时,应注意到:
(1)从上到下顺序地存在支配关系,用直线段(作用线)表示上一层次因素与下一层次因素之间的关系,同一层次及不相邻元素之间不存在支配关系。
(2)整个结构不受层次限制。
(3)最高层只有一个因素,每个因素所支配元素一般不超过9个,元素过多可进一步分层。
(4)对某些具有子层次结构可引入虚元素,使之成为典型递阶层次结构。
二、构造比较判断矩阵
设有m个目标(方案或元素),根据某一准则,将这m个目标两两进行比较,把第i个目标(i=1,2,…,m)对第j个目标的相对重要性记为a ij,(j=1,2,…,m),这样构造的m阶矩阵用于求解各个目标关于某准则的优先权重,成为权重解析判断矩阵,简称判断矩阵,记作A=(a ij)m×m。
Satty 于1980年根据一般人的认知习惯和判断能力给出了属性间相对重要性等级表(见表4-4)。
利用该表取a ij 的值,
称为1-9标度方法。
表4-4 目标重要性判断矩阵A 中元素的取值
若决策者能够准确估计a ij (i,j,k=1,2,…,m),则有: a ij =1/a ji a ij= a ik ·a kj a ii =1
定义4-1 设A=(a ij )m×m ,A>0,(即a ij >0;i,j=1,2,…,m),如果满足条件(1)a ii =1(i =1,2,…,m);(2)a ij =1/a ji
(i,j =1,2,…,m),则称矩阵A 为互反正矩阵。
定义4-2 设A=(a ij )m×m ,A>0,如果满足条件a ij= a ik ·a kj (i,j,k=1,2,…,m)则称矩阵A 为一致性矩阵。
定理4-1 对于任何一个m 阶互反正矩阵A ,均有
≥m,其中
是矩阵A 的最大特征值。
定理4-2 m 阶互反正矩阵A 为一致性矩阵的充分必要条件是A 的最大特征根为m 。
三、单准则下的排序
层次分析法的信息基础是比较判断矩阵。
由于每个准则都支配下一层若干因素,这样对于每一个准则及它所支配的因素都可以得到一个比较判断矩阵。
因此根据比较判断矩阵如何求得各因素w 1,w 2, …,w m 对于准则A 的相对排序权重的过程称为单准则下的排序。
这里设A=(a ij )m×m ,A>0。
(一)本征向量法 利用AW=
W 求出所有的值,其中为的最大值,求出对应的特征向量W *
,然后把特征向量W *
规一化为
向量W ,则W=[w 1,w 2, …w m ]T
为各个目标的权重。
求需要解m 次方程,当m≥3时,计算比较麻烦,可以利用matlab 来求解。
(二)判断矩阵的近似解法
判断矩阵是决策者主观判断的定量描述,求解判断矩阵不要求过高的精度。
这里,介绍三种近似计算方法:根法、和法及幂法。
幂法适于在计算机上运算。
1、根法
(1)A 中每行元素连乘并开m 次方,得到向量
其中,
(2)对W *作归一化处理,得到权重向量W=(w 1,w 2, …w m )T
,其中
(3)对A 中每列元素求和,得到向量S=(s 1,s 2, …s m ),其中s j =
max
λmax
λλλmax λλmax λλT m
w w w W )
,...,,(**2*1*
=m m
j ij
i
a w ∏==1
*∑==m
i i i
i w w w 1
*
*/∑=m
i ij
a
1
(4)计算的值,
=
2、和法
(1)将A 的元素按列作归一化处理,得矩阵Q=(q ij )m×m 。
其中,
(2)将Q 的元素按行相加,得向量。
其中,
(3)对向量
作归一化处理,得权重向量W=(w 1
,w 2
, …w m
)T
,其中
(4)求出最大特征值
3、幂法
幂法是一种逐步迭代的方法,经过若干次迭代计算,按照规定的精度,求出判断矩阵A 的最大特征值及其对应的特征向
量。
定理 3 设矩阵A=(a ij )m×m ,A>0,则,其中,W 是A 的最大特征值对应的的特征向量,C 为常
数,向量e=(1,1,…,1)T。
幂法的计算步骤是:
①任取初始正向量X (0)=(x 1(0), x 2(0), …, x m (0))T
,计算
②迭代计算,对于k=0,1,2, …计算
③精度检查。
当
时,转入步骤④;否则,令k=k+1,转入步骤②。
④求最大特征值和对应的特征向量,将Y (k+1)归一化,即
例 判断矩阵 1 2 5 A = 1/2 1 7 1/5 1/7 1
用幂法计算A 的最大特征值
及其对应额特征向量。
精度
=0.0001。
解:取初始向量X (0)
=(1,1,1)T
,迭代过程见下表
max
λSW w s i m
i i ==∑=1max
λ∑=m i i i
w AW m 1)(1∑==m
k kj
ij ij a a q 1
/T
m ),...,,(21αααα=∑==m
j ij
i q 1
αα∑==m
k k
i i w 1
/αα∑==m i i
i w AW m 1max
)(1λCW e A e e
A k
T k k =∞
→lim
)0()0()0()
0(0/},{max m X Y x X m i i
===∞
1
)1()1()1()
1(1)()1(/},{,max ++++∞
+++====k k k k i i
k k k k m X Y x X m AY X ε
<-+k k m m 11
max 1
)1()
1(,/+=++==∑k m
i k i k m y Y
W λmax
λε
由上表看出,当k=7时,|m 8-m 7|=|3.1189-3.1189|=0<0.0001,迭代终止,得到=3.1189,W=
(0.5415,0.3816,0.0769)T
四、单准则下的一致性检验
由于客观事物的复杂性,会使我们的判断带有主观性和片面性,完全要求每次比较判断的思维标准一致是不太可能的。
因此在我们构造比较判断矩阵时,我们并不要求n(n-1)/2次比较全部一致。
但这可能出现甲与乙相比明显重要,乙与丙相比极端重要,丙与甲相比明显重要,这种比较判断会出现严重不一致的情况。
我们虽然不要求判断具有一致性,但一个混乱的,经不起推敲的比较判断矩阵有可能导致决策的失误,所以我们希望在判断时应大体一致。
而上述计算权重的方法,当判断矩阵过于偏离一致性时,其可靠程度也就值得怀疑了。
因此,对于每一层次作单准则排序时,均需要作一致性的检验。
一致性指标(Consistency Index,CI ):
随机指标(Random Index,RI )
一致性比率(Consistency Rate,CR ):CR=CI/RI 当CR 取0.1时,最大特征值=CI·(m -1)+m=0.1·RI·(m -1)+m
表中当n=1,2时,RI=0,这是因为1,2阶判断矩阵总是一致的。
当n≥3时,若CR<0.1即
<
,认为比较判断矩阵的一致性可以接受,否则应对判断矩阵作适当的修正,直
到
小于
通过一致性检验时,求得的W 才有效。
五、层次总排序
计算同一层次中所有元素对最高层(总目标)的相对重要性标度(又称权重向量)称为层次总排序。
1、层次总排序的步骤为:
(1)计算同一层次所有因素对最高层相对重要性的权重向量,这一过程是自上而下逐层进行; (2)设已计算出第k-1层上有n k-1个元素相对总目标的权重向量为w
(k-1)
=(w 1
(k-1)
, w 2
(k-1)
,…, w n(k-1)
(k-1))T
(3)第k 层有个n k 个元素,他们对于上一层次(第k-1层)的某个元素j 的单准则权重向量为p j (k)
=(w 1j (k)
, w 2j (k)
,…, w nkj)(k))T
(对于与k-1层第j 个元素无支配关系的对应w ij 取值为0); (4)第k 层相对总目标的权重向量为w k
= (p 1(k)
, p 2(k)
,…p k-1(k)
,)w
(k-1)
max
λ1
max --=
m m
CI λ'max λ'
max
λ
max λ'max
λmax
λ'max
λ
2、层次总排序的一致性检验
人们在对各层元素作比较时,尽管每一层中所用的比较尺度基本一致,但各层之间仍可能有所差异,而这种差异将随着层次总排序的逐渐计算而累加起来,因此需要从模型的总体上来检验这种差异尺度的累积是否显著,检验的过程称为层次总排序的一致性检验。
第k层的一致性检验指标CIk=(CI1(k-1), CI2(k-1),…, CIn K(k-1))w(k-1)
RI k=(RI1(k-1), RI2(k-1),…, RIn K(k-1))w(k-1)
CR k=CR k-1+CI k/RI k(3≤k≤n)
当CR k <0.1,可认为评价模型在第k层水平上整个达到局部满意一致性。
六、递阶层次结构权重解析过程
1、树状结构目标体系
目标可分为多个层次,每个下层目标都隶属于一个而且只隶属一个上层目标,下层目标是对上层目标的具体说明。
对于树状结构的目标体系,需由上而下逐步确定权重,即由树干向树梢,求树杈各枝相对于树杈的权重。
2、网状结构目标体系
网状结构的目标也分为多个层次,每个下层目标隶属于某几个上层目标(至少有一个下层目标隶属于不止一个上层目标)。
七、AHP方法的基本步骤
层次分析法大体分为以下六个步骤:
1、明确问题;建立层次结构;
2、两两比较,建立判断矩阵;2、层次单排序及其一致性检验;
3、层次总排序及其一
致性检验;4、根据分析计算结果,考虑相应的决策。
例4-4
(1)建立方案评价的递阶层次结构模型
设以A
为比较准则,B、
C、D、E
的两两比较判断矩阵为A,类似地构造矩阵B
、C、D
、E。
对于上述各比较判断矩阵,用Matlab数学软件求出其最大特征值及其对应的特征向量,将特征向量经归一化后,即可得相应的层次单排序的相对重要性权重向量,以及一致性指标CI和一致性比率CR,列表如下:
层次单排序及一致性计算结果
max max
由此可见,所有五个层次单排序的的值均小于,CR的值均小于0.1,符合一致性要求。
(4)层次总排序
已知第二层(B、C、D、E)相对于总目标层A的排序向量为W(2)= (0.4694,0.2788,0.1491,0.1027)T,而第三层(X、Y、Z)以第二层各个因素为准则的排序向量分别为:
P b(3)= (0.3429,0.2286,0.4286)T
P c(3)= (0.2849,0.2174,0.4977)T
P d(3)= (0.5472,0.1897,0.2631)T
P e(3)= (0.2727,0.1818,0.5455)T
则第三层(X、Y、Z)相对于总目标的排序向量为
W(3)=( P b(3), P c(3), P d(3), P e(3))·W(2)=(0.3500,0.2149,0.4352)T
(5)层次总排序的一致性检验
由于CI(2)=(CIb,CIc,CId,CIe)=(O,0.0091,0.0008,0.0000)
RI(2)=(RIb,RIc,RId,RIe)=(0.5800,0.5800,0.5800,0.5800)
CR2=CRa=0.0061
因此,CI3= CI(2)·W(2)=0.0027
RI3= RI(2)·W(2)=0.5800
CR3= CR2+ CI3/ RI3=0.0107<0.1
当CR3 <0.1,可以认为评价模型在第3层水平上整个达到局部满意一致性
(6)决策
所考虑的三种方案排序的相对优先排序为Z>X>Y,应该选择企业Z作为结盟伙伴。
利用MATLAB求解例4-4
(1)第二层对目标层的权重与一致性检验
A=[1 2 3 4;1/2 1 2 3;1/3 1/2 1 3/2;1/4 1/3 2/3 1] %输入比较判断矩阵A
a=eig(A) %求出A的所有的特征值
[X,D]=eig(A) %求出A的所有特征向量及对角矩阵
a1=a(1,:) %在A的所有特征值中取出最大的特征值(第1个)a2=X(:,1) %求最大特征值所对应的特征向量
a3=ones(1,4) %构造一个其中元素全为1的1×4矩阵
a4=a3*a2 %求出a2中所有元素的和
W2=1/a4*a2 %求出矩阵A的权重向量
CIa=(a1-4)/3 %求出A的一致性指标
CRa=CIa/0.90 %求出A的一致性比率
(2)求出第三层对第二层的各个因素的权重向量及一致性检验
程序基本和(1)中相似,求出B、C、D、E的最大特征值、权重向量和一致性比率
B b1 Pb CIb CRb
C c1 Pc CIc CRc
D d1 Pd CId CRd
E e1 Pe CIe CRe
(3)总排序
P=[Pb Pc Pd Pe] %构造矩阵W3
W3=P*W2 %求出总的排序向量
(4)总排序的一致性检验
CI2=[CIb CIc CId CIe]
RI2=[0.58 0.58 0.58 0.58]
CR2=CRa
CI3=CI2*W2 %第3层的一致性指标
RI3=RI2*W2 %第3层随机指标
CR3=CR2+CI3/RI3 %第3层的一致性比率
参考文献:
1.张强:《决策理论与方法》,大连,东北财经大学出版社,2009。
2.李柏年:《模糊数学及其应用》,合肥,合肥工业大学出版社,2007。
3.赵新泉、彭勇行:《管理决策分析》,北京,科学出版社,2008。
4.杜栋、庞庆华:《现代综合评价方法与案例精选》,北京,清华大学出版社,2005。