这里的平均随机一致性指标(RI) 是由大量随机试验得出的平均数据
代码部分
function RI=getRI(n)
times=10000;
CIS=zeros(1,times);
for t=1:times
CI=getCI(n);
CIS(t)=CI;
end
RI=sum(CIS)/times;
function CI=getCI(n)
A=eye(n);
mat=[9 8 7 6 5 4 3 2 1 1/2 1/3 1/4 1/5 1/6 1/7 1/8 1/9];
for i=1:n-1
for j=i+1:n
randnum=randi(17);
A(i,j)=mat(randnum);
A(j,i)=mat(18-randnum);
end
end
Asum=sum(A,1);
Aprogress=A./(ones(n,1)*Asum);
W=sum(Aprogress,2)./n;
w=A*W;
lam=sum(w./W)/n;
CI=(lam-n)/(n-1);
end
end
经过计算得到矩阵大小为(1-33)时RI大小 可以看出与标准值还是由一定的出入的,可能是由于matlab随机数生成的机制导致的