层次分析法简单程序.docVIP

层次分析法简单程序 毕业后第一次做的项目是关于飞机首飞的风险管理。老实说，工厂对飞机试飞的工作并不太重视，这牵扯到多方面的原因了，不想在此发牢骚。此项目走走停停，到这个月也几近尾声了，剩下的一步是风险监控的计划与实施。其间深切感受到了要做好一件工作是多么的不易。下面是在做风险调查问卷资料分析用的一个小程度，没有优化过算法。数据都用excel储存，这程序可写成一ABH（Analytical Hierarchy Process,层次分析法）函数,在matlab主函数中直接将要分析的判断矩阵从excel导入到matlab中，再调用此函数分析即可。另外在A—B—C—D等多层子因素的情况下还要对C，D层次排序做一致性检查，具体的原理方法就查阅相关的书籍吧。程序如下： %层次分析法简单程序format shortclearclcdisp(Please Input The Matrix:);A=input(A=);[m,n]=size(A);if m~=n??? return;endWpj=ones(m,1); %用于存储判断矩阵每一行的几何平均值Wgy=ones(m,1); %用于存储Wpj归一后值，意义指每一行在判断矩阵中所占的比例（即判断矩阵中各影响因素所占的比重，以行表示影响因素，列表示判断分值）Aw=ones(m,1);? %判断矩阵与Wgy相乘后的特征向量，其用途是之后用于求最大特征值a=1;for i=1:m??? for j=1:n? %其实判断矩阵m,n是相等的,也可用m??????? a=a*A(i,j);??? end??? Wpj(i)=power(a,1/n);??? a=1;ends=sum(Wpj);for i=1:m??? Wgy(i)=Wpj(i)/s;endAw=A*Wgy;LANDAmax=0;%用于存储最大特征值for i=1:m??? LANDAmax=LANDAmax+Aw(i)/(n*Wgy(i));%求得最大特征值endCI=(LANDAmax-n)/(n-1);RIn=[0 0 0.58 0.90 1.12 1.24 1.32 1.41 1.45];RI=RIn(n);Judge=CI/RI;disp(Wgy=);disp(Wgy);disp(CI=);disp(CI);disp(RI=);disp(RI);disp(Judge=);disp(Judge);if Judge0.1??? disp(判断矩阵满足一致性检验要求);else??? disp(判断矩阵错误，不满足一致性检验要求);end %%%%% ? 未做将数据导回excel的工作，工作中用的是ctrl+C和ctrl+V，呵呵，也蛮快。 例子（在Command Window的显示）： Please Input The Matrix:A=[1 5 3;1/5 1 1/2;1/3 2 1]Wgy=??? 0.6483??? 0.1220??? 0.2297 CI=??? 0.0018 RI=??? 0.5800 Judge=??? 0.0032 判断矩阵满足一致性检验要求