【源代码】AHP法算权重

%此子程序来源参考：《层次分析法在确定企业家评价指标权重中的应用》
function [W,CR]=AHP(P)
%计算指标权重
[m,n]=size(P);
for i=1:m 
    M_t(i)=prod(P(i,:));%计算每一行元素的积
    W_t(i)=M_t(i)^(1/n);
end
W=W_t/sum(W_t);
%计算判断矩阵的最大特征值
PW=P*W';
PW_t=PW'./W/n;
lamb=sum(PW_t);

%一致性检验
%R1为平均一致性指标表：第一行为阶数
R1=[1 2 3 4 5 6 7 8 9;
    0 0 0.58 0.90 1.12 1.24 1.32 1.41 1.45];
C1=(lamb-n)/(n-1);
CR=C1/R1(2,n);