资源描述：
《最新层次分析法matlab源程序.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、__________________________________________________层次分析法matlab源程序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(:,2)=Ay(:,1);m(2)=max(x(:,2));y(:,2)=x(:,2)/m(2);p=0.0001;i=2;k=abs(m(2)-m(1));while k>p i=i+1; x(:,i)
2、=Ay(:,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(w);disp(t); %以下是一致性检验CI=(t-n)/(n-1);RI=[000.520.891.121.261.361.411.461.491.521.541.561.581.59];CR=CI/RI(n)
3、;ifCR<0.10 disp('此矩阵的一致性可以接受!'); disp('CI=');disp(CI); disp('CR=');disp(CR);endfunctionAHPInit1(x,y)%层次分析的初始化%默认只有两层x为准则数，y为方案数%CToT为准则对目标生成的比较阵%EigOfCri为准则层的特征向量%EigOfOpt为选项层的特征向量EigOfCri=zeros(x,1);%准则层的特征向量EigOfOpt=zeros(y,x);dim=x;%维度RI=[000.580.901.121.241.321.411.451.491.51];%RI
4、标准%生成成对比较阵fori=1:dim CToT(i,:)=input('请输入数据:');endCToT收集于网络，如有侵权请联系管理员删除__________________________________________________%输出pause,tempmatrix=zeros(x+1);tempmatrix=AHP1(dim,CToT);EigOfCri=tempmatrix(1:x);ci1=tempmatrix(1+x);EigOfCrici1pause,matrix=cell(x);%元胞数组ci=zeros(1,x);dim=y;fork
5、=1:x matrix{k}=zeros(dim,dim);%生成成对比较阵fori=1:dim 收集于网络，如有侵权请联系管理员删除__________________________________________________matrix{k}(i,:)=input('请输入数据:');end%判断该比较阵是不是一致阵tempmatrix=zeros(y+1);tempmatrix=AHP1(dim,matrix{k});EigOfOpt(:,k)=tempmatrix(1:y);ci(k)=tempmatrix(y+1);EigOfOpt(:,k)c
6、i(k)pause,end%下面进行组合一致性检查RI=[000.580.901.121.241.321.411.451.491.51];CR=ci1/RI(x)+ciEigOfCri/RI(y);CRifCR>0.1 disp('组合一致性不通过，请重新评分') returnend%下面根据比较阵的结果进行组合result=EigOfOptEigOfCri;resultfunctionf=AHP1(dim,CmpMatrix)RI=[000.580.901.121.241.321.411.451.491.51];%判断该比较阵是不是一致阵%判断该比较阵是不是一致阵[
7、V,D]=eig(CmpMatrix);%求得特征向量和特征值%求出最大特征值和它所对应的特征向量tempNum=D(1,1);pos=1;forh=1:dim ifD(h,h)>tempNum 收集于网络，如有侵权请联系管理员删除__________________________________________________tempNum=D(h,h); pos=h; endend eigVector=V(:,pos);maxeig=D(pos,pos);max
