D = input('请输入判断矩阵D:');
% 一致性检验
a_max = max(eig(D));% a_max为最大特征值
disp(['λmax = ',num2str(a_max)])
[~,n] = size(D); % 使用不到的数可以替换为~
CI = (a_max-n)/(n-1);
disp('CI = ');disp(CI)
% RI为平均随机一致性指标表
RI = [0 0 0.52 0.89 1.12 1.26 1.26 1.36 1.41 1.49 1.52 1.54 1.56 1.58 1.59] ;
CR = CI/RI(n);
disp('CR = '); disp(CR)
if CR<0.1
disp('**********CR<0.1,判断矩阵通过一致检验!**********')
else
error('一致检验不通过,请修改判断矩阵!')
end
输入判断矩阵即可