层次分析法程序代码

层次分析法1）建立层次结构模型：（2）构造判断矩阵判断矩阵ijAa应为正互反矩阵，而且ija的判断如下（1~9尺度法）：标度aij含义1的影响与CCji相同3的影响比CCji稍强5的影响比CCji强7的影响比CCji明显的强9的影响比CCji绝对的强2，4，6，8的影响与CCji之比在上述两个相邻的等级之间91,......21的影响与CCji之比为上面aij的复反数（3）单层排序及一致性检验1、单层排序求解判断矩阵A的最大特征值max，再由最大特征值求出对应的特征向量决策目标P1准则C1准则C2准则C3准则C4P2P3P4P5P6maxA，并将标准化，即为同一层相对于上一层某一因素的权重，根据此权重的大小，便可确定该层因素的排序。2、一致性检验取一致性指标max1nCIn，（n为A的阶数）取随机性指标RI如下：n123456789……RI000.580.91.121.241.321.411.45……令CICRRI，若0.1CR，则认为A具有一致性。否则，需要对A进行调整，直到具有满意的一致性为止。（4）层次总排序及一致性检验假定准则层12,,,nCCC排序完成，其权重分别为12,,,naaa，方案层P包含m个方案：12,,,mPPP。其相对于上一层的1,2,,jCjn对方案层P中的m个方案进行单层排序，其排序权重记为12,,,jjmjbbb1,2,,jn，则方案层P中第i个方案Pi的总排序权重为1njijjab，见下表：层次C层次P1C2CnCP层总排序权重1a2ana1P11b12b1nb11njjjab2P21b22b2nb21njjjabnP1mb2mbmnb1njmjjab从而确定P层的排序。例：纯文本文件txt3.txt中的数据格式如下：111411/2112411/211/21531/21/41/41/511/31/3111/331122233111/41/241321/3111/41/5411/2521131/31/311/737111/353171/51/711171171/71/711791/7111/911matlab程序：fid=fopen('txt3.txt','r');n1=6;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:n1[x,y]=eig(eval(char(['b',int2str(i)])));lamda=max(diag(y));num=find(diag(y)==lamda);w1(:,i)=x(:,num)/sum(x(:,num));cr1(i)=(lamda-n2)/(n2-1)/ri(n2);endcr1,ts=w1*w0,cr=cr1*w0