资源描述：
《基于MATLAB的AHP层次分析法资料的.ppt》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、基于MATLAB的AHP（层次分析法）的实现案例：整治路口交通秩序方案的评估决策AHP（层次分析法）介绍及步骤步骤:1、构造层次结构模型2、建立判断矩阵3、检验判断矩阵的一致性4、计算权重，选择决策层次分析法是一种定性与定量分析相结合的多因素决策分析方法。这种方法将决策者的经验判断给于数量化，在目标因素结构复杂且缺乏必要数据的情况下使用更为方便，因而在实践中得到广泛应用。关于MATLAB的相关介绍Matlab是矩阵实验室（MatrixLaboratory）的简称，是美国MathWorks公司出品的数学软件，用于算法开
2、发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境。Matlab可以进行矩阵运算、绘制函数和数据图像、设计算法、创建用户界面、连接用其他编程语言编写的程序等。Matlab以矩阵为计算单位，采用M语言作为程序语言，与C语言有诸多相似之处，并可方便地与C/C++、MicrosoftExcel等工具和软件进行结合并进行代码共享和数据交换，可以方便地进行数值分析、图像处理等功能，配合功能强大的统计和金融工具箱，Matlab已经可以在概率统计、经济管理等方面发挥强大的作用。案例简介某市一十字路口常常因行人过街拥
3、挤，存在安全隐患，市政部门欲对该路口进行改造，现提出了3套改造方案：方案1（S1）：建地下通道；方案2（S2）：建人行天桥；方案3（S3）：拆除周围的旧建筑，拓宽街面。市政部门认为，该改造工程需考虑如下几个方面的指标：指标1（P1）：通车能力的大小；指标2（P2）：交通安全系数的高低；指标3（P3）：建筑费用的高低；指标4（P4）：群众出行方便度的大小；指标5（P5）：市容整洁程度的高低。现在需要就以上问题进行决策，需决定在三套方案（S1~S3）中选用最优方案。构建结构层次模型整治路口交通秩序通车能力P1交通安全系数
4、P2建筑费用P3群众出行方便度P4市容整洁程度P5地下通道S1拆除建筑S3过街天桥S2目标层准则层措施层形成判断矩阵标度aij定义（比较因素i与j）1因素i与j同样重要3因素i与j稍微重要5因素i与j较强重要7因素i与j强烈重要9因素i与j绝对重要2,4,6,8上述两相邻判断的中间值1~9的倒数若因素j与因素i比较，得到的判断值为aji=1/aij，aii=1Saaty九级标度法及其含义形成判断矩阵目标P1P2P3P4P5P113456P21/31333P3¼1/311/33P41/51/3313P51/61/31/
5、31/31指标间重要程度判断矩阵（表一）P1S1S2S3S1113S2113S31/31/31通车能力（P1）判断矩阵（表二）形成判断矩阵P1S1S2S3S1121S21/211/3S3131交通安全系数（P2）判断矩阵（表3）P1S1S2S3S111/41/7S2411/4S3741建筑费用（P3）判断矩阵（表4）P1S1S2S3S111/31/5S231½S3521群众方便度（P4）判断矩阵（表5）P1S1S2S3S1123S2½11S31/311市容整洁度（P5）判断矩阵（表6）计算-MATLABMatlab程
6、序设计思路一致性不接受开始输入判断矩阵计算权重向量和最大特征值输出权重向量和最大特征值计算CRCR<0.1Y一致性接受输出CI、CRN结束程序的设计目标首先是计算判断矩阵的权重向量和最大特征值，然后进行CI和CR的计算，进而判断矩阵的一次性。程序流程图如右图所示程序代码disp('请输入判断矩阵A(n阶)');A=input('A=');[n,n]=size(A);x=ones(n,100);y=ones(n,100);m=zeros(1,100);m(1)=max(x(:,1));y(:,1)=x(:,1);x(:
7、,2)=A*y(:,1);m(2)=max(x(:,2));y(:,2)=x(:,2)/m(2);p=0.0001;i=2;k=abs(m(2)-m(1));whilek>pi=i+1;x(:,i)=A*y(:,i-1);m(i)=max(x(:,i));y(:,i)=x(:,i)/m(i);k=abs(m(i)-m(i-1));enda=sum(y(:,i));w=y(:,i)/a;t=m(i);disp('权向量');disp(w);disp('最大特征值');disp(t);%以下是一致性检验CI=(t-n)/
8、(n-1);RI=[000.520.891.121.261.361.411.461.491.521.541.561.581.59];CR=CI/RI(n);ifCR<0.10disp('此矩阵的一致性可以接受!');disp('CI=');disp(CI);disp('CR=');disp(CR);elsedisp('此矩阵的一致性不可以接受!
