matlab层次分析法的组合一致性检验

functions(x,y)%层次分析的初始化%默认只有两层x为准则数，y为方案数%A为准则层对目标层生成的比较阵%B为方案层对准则层生成的比较阵%EigA为准则层的特征向量%EigB为方案层的特征向量x=5;y=3;EigA=zeros(x,1);EigB=zeros(y,x);dim=x;RI=[000.580.901.121.241.321.411.451.491.51];%RI标准%准则层对目标层的成对比较阵A=[11/2433217551/41/711/21/31/31/52111/31/5311]tempB=zeros(x+1);tempB=jianyan(dim,A);EigA=tempB(1:x);ci1=tempB(1+x);'准则层的特征向量:'EigA'准则层的ci值:'ci1ci=zeros(1,x);dim=y;%方案层对准则层各成对比较阵B{1}=[1251/2121/51/21];B{2}=[11/31/8311/3831];B{3}=[1131131/31/31];B{4}=[1341/3111/411];B{5}=[111/4111/4441];'方案层对准则层各成对比较阵:'B{:}%判断该比较阵是不是一致阵fork=1:xtempB=zeros(y+1);tempB=jianyan(dim,B{k});EigB(:,k)=tempB(1:y);ci(k)=tempB(y+1);['方案层对准则层的第',num2str(k),'因素的特征向量和CI值:']EigB(:,k)ci(k)end%下面进行组合一致性检查RI=[000.580.901.121.241.321.411.451.491.51];CR=ci1/RI(x)+ci*EigA/RI(y);'组合一致性比率CR值'CRifCR0.1disp('组合一致性不通过，请重新评分')returnend%下面根据比较阵的结果进行组合result=EigB*EigA;'方案层对目标层的组合权向量为：'resultfunctionf=jianyan(dim,CmpMatrix)RI=[000.580.901.121.241.321.411.451.491.51];%判断该比较阵是不是一致阵[V,D]=eig(CmpMatrix);%求得特征向量和特征值%求出最大特征值和它所对应的特征向量tempNum=D(1,1);pos=1;forh=1:dimifD(h,h)tempNumtempNum=D(h,h);pos=h;endendeigVector=V(:,pos);maxeig=D(pos,pos);CI=(maxeig-dim)/(dim-1);CR=CI/RI(dim);ifCR0.1disp('准则对目标影响度评分生成的矩阵不是一致阵，请重新评分')returnend%归一化sum=0;forh=1:dimsum=sum+eigVector(h);endforh=1:dimeigVector(h)=eigVector(h)/sum;endf=[eigVector;CI];结果：A=1.00000.50004.00003.00003.00002.00001.00007.00005.00005.00000.25000.14291.00000.50000.33330.33330.20002.00001.00001.00000.33330.20003.00001.00001.0000ans=准则层的特征向量:EigA=0.26360.47580.05380.09810.1087ans=准则层的ci值:ci1=0.0180ans=方案层对准则层各成对比较阵:ans=1.00002.00005.00000.50001.00002.00000.20000.50001.0000ans=1.00000.33330.12503.00001.00000.33338.00003.00001.0000ans=1.00001.00003.00001.00001.00003.00000.33330.33331.0000ans=1.00003.00004.00000.33331.00001.00000.25001.00001.0000ans=1.00001.00000.25001.00001.00000.25004.00004.00001.0000ans=方案层对准则层的第1因素的特征向量和CI值:ans=0.59540.27640.1283ans=0.0028ans=方案层对准则层的第2因素的特征向量和CI值:ans=0.08190.23630.6817ans=7.7081e-004ans=方案层对准则层的第3因素的特征向量和CI值:ans=0.42860.42860.1429ans=-4.4409e-016ans=方案层对准则层的第4因素的特征向量和CI值:ans=0.63370.19190.1744ans=0.0046ans=方案层对准则层的第5因素的特征向量和CI值:ans=0.16670.16670.6667ans=-4.4409e-016ans=组合一致性比率CR值CR=0.0188ans=方案层对目标层的组合权向量为：result=0.29930.24530.4554