层次分析法

%使用方法
%（1）构造判断矩阵A
%（2）将下文代码复制粘贴到Matlab中即可
%例如：A=[1 3 5;0.33 1 3;0.2 0.33,1]
%复制在命令窗口

disp('请输入准则层判断矩阵A(n阶)');
A=input('A=');
[n,n]=size(A);%得到矩阵的阶数
[V,D]=eig(A);%求得特征向量和特征值
            %求出最大特征值和它所对应的特征向量
            %V为特征向量，D为特征值，按对角线元素排序
tempNum=D(1,1); %特征值初值
pos=1;  %标记的第一个数
for h=1:n
    if D(h,h)>tempNum
        tempNum=D(h,h);
        pos=h;
    end
end
%以上得到最大特征值极其位置行列数pos  
w=abs(V(:,pos));%找出最大特征值对应的特征向量
w=w/sum(w);%归一化处理
t=D(pos,pos);%t为最大特征值
disp('准则层特征向量w=');disp(w);disp('准则层最大特征根t=');disp(t);
         %以下是一致性检验
CI=(t-n)/(n-1);RI=[0 0 0.52 0.89 1.12 1.26 1.36 1.41 1.46 1.49 1.52 1.54 1.56 1.58 1.59 1.60 1.61 1.615 1.62 1.63];
CR=CI/RI(n);
if CR<0.10
    disp('此矩阵的一致性可以接受!');
    disp('CI=');disp(CI);
    disp('CR=');disp(CR);
else disp('此矩阵的一致性验证失败，请重新进行评分!');
end