(完整版)层次分析法实例讲解学习

1、层次分析法实例讲解学习生活实际例题：旅游实例，有三个旅游地点供游客们选择，连云港，常州，徐州。影响游客们决策的因素主要有以下五项：景色、费用、居住、饮食、旅途。请根据个人偏好选择最佳旅游地点。分析：旅游点是方案层，将它们分别用B,B,B表示，影响旅游决策的因素为准123则层A,A,A,A,A；目标层为选择旅游地，即可以建立以下模型：12345费用饮食常州徐州选择旅游地建立判断矩阵：准则层判断矩阵（即各种因素在旅客偏好选择中所占有的不同比重）11/243321755A=1/41/711/21/31/31/52111/31/5311方案层判断矩阵建立（针对每一个影响因素来对方案层建立）：-125j

2、1/31/8-113_B=1/212B=311/3B=11311/51/21183111/31/31134111/4B=1/311B=111/411/4111441求准则层判断矩阵A的特征值：Matlab运行程序：a,b=eig（A）b'矩阵的对角线为准则层判断矩阵A的特征值:00000.03100000.03100000.00500000.0055.0730000=0.031,九4=0.005,九=0.0055即九=5.073,九=0.031,九123选出最大特征值:九=max(九,九,九,九,九)=X123451最大特征值的特征向量即为准则层的影响因素所占的权重，X所对应的特征向量

3、为:w=1-0.4658-0.8409-0.0951-0.1733-0.19201归一化(最简matlab程序为w=w1./sum(w1)w=b.26360.47590.05380.09810.1087一致性指标的检验:由maxX是否等于5来检验判断矩阵A是否为一致矩阵。由于特征根连续地依赖于矩阵A中的值，故maxX比5大得越多，A的非一致性程度也就越严重，maxX对应的标准化特征向量也就越不能真实地反映出对因素A(i=1,5)的影i响中所占的比重。计算一致性指标CI:CI=Xnn1此题的一致性指标为CI=Xmaxn=5.073-5=0.018n15-1平均随机一致性指标RI相对固定，如下表:

4、RI随机一致性指标n1234567891011RI000.580.901.121.241.321.411.451.491.51计算一致性比例CR:CR=馬当CR时，认为判断矩阵的一致性是可以接受的，否则应对判断矩阵作适当修正。本题:CR=RI吧=0-016<0-1可行。按照如上方式处理矩阵B,B,B,B,B得:12345B：1九=3.005maxCI=0.003RIw=O.595,O.277,O.1291=O.58CR二0.005V0.1B：2九=3.002maxCI=0.001RIw=O.O82,O.236,O.6822=O.58CR二0.002V0.1B：3九=3maxCI=ORIw

5、=O.429,O.429,O.1423=O.58CR二OvO.1B4：九=3.009maxCI=0.005RIw=O.633,O.193,O.1754=O.58CR二0.009V0.1B：5九=3maxCI=0RIw=O.166,O.166,O.6685=O.58CR二OvO.1计算方案层j(j=1,n)对总目标的权值：Q(j)=工w(i)w(j)ii=1即连云港对总排序的权值为：Q(1)=工w(i)w(1)=ii=10.263x0.595+0.475x0.082+0.055x0.429+0.099x0.633+0.110x0.116=0.3同理可计算出常州，徐州对总排序的权值：Q(2)=工w

6、(i)w(2)=0.246Q(3)=工w(i)w=0.456i=1iii=1可建立如下表格清晰数据的排列将准则层矩阵的最大特征值的特征向量代表的权值放于准则层权值的右方，准则层每个因素所对应的方案层矩阵的最大特征值的特征向量放于所对应的影响因素的下面。层次总排序合成表准侧层权值风景费用居住饮食旅游总排序权值0.26360.47590.05380.09810.1087Ji=r方案层连云港0.5950.0820.4290.6330.1660.300单排序常州0.2770.2360.4290.1930.1660.246权值徐州0.1290.6820.1420.1750.6680.456对于整个算法编

7、写matlab程序如下：fid=fopen('txt2.txt','r');n1=5;n2=3;a=;fori=1:n1tmp=str2num(fgetl(fid);a=a;tmp;%读准则层判断矩阵endfori=1:n1str1=char('b',int2str(i),'=;');str2=char('b',int2str(i),'=b',int2str(i),'tmp;');eval(str1);forj=1:n2tmp=str2num(fgetl(fid);eval(str2);%读方案层的判断矩阵endendri=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)fori=1:n1x,y=eig(eval(char('b',int2str(i);lamda=m