a=[1,1,1,4,1,1/2
1,1,2,4,1,1/2
1,1/2,1,5,3,1/2
1/4,1/4,1/5,1,1/3,1/3
1,1,1/3,3,1,1
2,2,2,3,3,1];
[x,y]=eig(a); %求a矩阵的所有特征值构成对角阵y,求a的特征向量构成列向量x
eigenvalue=diag(y); %返回y主对角线上的元素构成矩阵eigenvalue
lamda=eigenvalue(1); %lamda=矩阵eigenvalue的第一个元素
ci1=(lamda-6)/5; %这是层次分析法里计算平均一致性检验系数
cr1=ci1/1.24 %这块应该是层次分析法里面计算随机一致性检验系数
w1=x(:,1)/sum(x(:,1)) %对x进行归一化处理
b1=[1,1/4,1/2;4,1,3;2,1/3,1]; %从这里往后和前面一样,只是比较矩阵变了
[x,y]=eig(b1);
eigenvalue=diag(y);
lamda=eigenvalue(1);
ci21=(lamda-3)/2;
cr21=ci21/0.58
w21=x(:,1)/sum(x(:,1));
如果你这就是层次分析法,建议你在百度上看看层次分析法的原理,经典运筹学方法。一看就懂。看完这程序你自己就会编