层次分析法实例讲解学习

r
层次分析法实例讲解学习
生活实际例题：
旅游实例，有三个旅游地点供游客们选择，连云港，常州，徐州。

影响游客们决策的因素主要有以下五项：景色、费用、居住、饮食、旅途。

请根据个人偏好选择最佳旅游地点。

分析: 旅游点是方案层，将它们分别用表示，影响旅游决策的因素为准
321,,B B B 则层；目标层为选择旅游地，即可以建立以下模型：
54321,,,,A A A A A 建立判断矩阵：
准则层判断矩阵（即各种因素在旅客偏好选择中所占有的不同比重）：
选择旅游地
景色徐州
连云港
旅途
饮食
居住
费用
常州
⎥⎥⎥
⎥
⎥⎥⎦
⎤⎢⎢⎢⎢⎢⎢⎣⎡=113
5
/13
/11125/13/13/12/117/14
/15571
233
42/11A 方案层判断矩阵建立（针对每一个影响因素来对方案层建立）:
⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=12/15/1212/15211B ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=1381/3131/81/311B ⎥⎥
⎥⎦
⎤⎢⎢⎢⎣⎡=13/13/1311311
1B ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=114/1113/14311B ⎥⎥
⎥⎦
⎤
⎢⎢⎢⎣⎡=1441/4111/4111B 求准则层判断矩阵的特征值：A Matlab 运行程序： [a,b]=eig(A)
‘b’矩阵的对角线为准则层判断矩阵的特征值：
A 即005.0,005.0,031.0,031.0,073.554321-=-==-==λλλλλ选出最大特征值：
1
54321,,,,max λλλλλλλ==）（最大特征值的特征向量即为准则层的影响因素所占的权重，所对应的特征向
λ量为：
[]
0.1920-0.1733-0.0951-0.8409-0.4658-1=w 归一化（最简matlab 程序为 w=w1./sum(w1)）
[]
1087.00981.00538.04759.02636.0=w 一致性指标的检验：
由 是否等于5来检验判断矩阵 A 是否为一致矩阵。

由于特征根连续地依λmax 赖于矩阵A 中的值 ，故 比5大得越多， A 的非一致性程度也就越严重，
λmax 对应的标准化特征向量也就越不能真实地反映出对因素的影
λmax )5,,1( =i A i 响中所占的比重。

计算一致性指标：
CI 1
max --=
n n
CI λ此题的一致性指标为
0.0181
-55
-5.0731
max ==
--=
n n
CI λ平均随机一致性指标相对固定，如下表：
RI 随机一致性指标
RI n 12345
6
7891011RI
0.58
0.90
1.12
1.24
1.32
1.41
1.45
1.49
1.51
计算一致性比例：
CR RI
CI CR =
当时，认为判断矩阵的一致性是可以接受的，否则应对判断矩阵作适当修正。

CR
本题：
1.00.0161.12
0.018<===
RI CI CR 可行。

按照如上方式处理矩阵得：
54321,,,,B B B B B ：
1B
3.005 [0.595,0.277,0.129]
=max λ=1w 0.003 0.58
0.005<0.1
=CI =RI =CR ：
2B 3.002 [0.082,0.236,0.682]
=max λ=2w 0.001 0.58 0.002<0.1
=CI =RI =CR ：
3B 3 [0.429,0.429,0.142]
=max λ=3w 0 0.58 0<0.1
=CI =RI =CR ：
4B 3.009 [0.633,0.193,0.175]
=max λ=4w 0.005 0.58
0.009<0.1
=CI =RI =CR ：
5B 3
[0.166,0.166,0.668]
=max λ=5w 0 0.58
0<0.1
=CI =RI =CR 计算方案层对总目标的权值：
),,1(n j j =∑==n
i i j w i w j Q 1
)
()()(即连云港对总排序的权值为：
=
=∑=5
1)1()()1(i i w i w Q
0.3
0.1160.1100.6330.0990.4290.0550.0820.4750.5950.263=⨯+⨯+⨯+⨯+⨯同理可计算出常州，徐州对总排序的权值：
0.246)2()()2(5
1==∑=i i w i w Q 0.456
)3()()3(5
1
==∑=i i w i w Q 可建立如下表格清晰数据的排列：
将准则层矩阵的最大特征值的特征向量代表的权值放于准则层权值的右方，准则层每个因素所对应的方案层矩阵的最大特征值的特征向量放于所对应的影响因素的下面。

层次总排序合成表
风景费用居住
饮食
旅游准侧层权值0.26360.47590.05380.09810.1087
总排序权值连云港0.5950.0820.4290.6330.1660.300常州0.2770.2360.4290.1930.1660.246方案层单排序权值
徐州
0.129
0.682
0.142
0.175
0.668
0.456
对于整个算法编写matlab 程序如下：
fid=fopen('txt2.txt','r');
n1=5;n2=3;
a=[];
for i=1:n1
tmp=str2num(fgetl(fid));
a=[a;tmp]; %读准则层判断矩阵end
for i=1:n1
str1=char(['b',int2str(i),'=[];']);
str2=char(['b',int2str(i),'=[b',int2str(i),';tmp];']); eval(str1);
for j=1:n2
tmp=str2num(fgetl(fid));
eval(str2); %读方案层的判断矩阵
end
end
ri=[0,0,0.58,0.90,1.12,1.24,1.32,1.41,1.45]; %一致性指标[x,y]=eig(a);
lamda=max(diag(y));
num=find(diag(y)==lamda);
w0=x(:,num)/sum(x(:,num));
cr0=(lamda-n1)/(n1-1)/ri(n1)
for i=1:n1
[x,y]=eig(eval(char(['b',int2str(i)])));
lamda=max(diag(y));
num=find(diag(y)==lamda);
w1(:,i)=x(:,num)/sum(x(:,num));
cr1(i)=(lamda-n2)/(n2-1)/ri(n2);
end
cr1, ts=w1*w0, cr=cr1*w0。