% 输入判断矩阵
disp('请输入矩阵A: ')
A = input(‘判断矩阵A=’)
% 求权重(V特征值-列向量,D特征向量)
[V,D] = eig(A)
% max(D)每一列最大值 两次max
Max_eig = max(max(D))
% 最大特征值所在位置
D == Max_eig
[r,c] = find(D == Max_eig , 1)
% 最大特征值所在列(权重)
V(:,c)
disp(‘权重结果为:’);
disp( V(:,c) ./ sum(V(:,c)) )
% 一致性检验
% 一致性比例CI
CI = (Max_eig - 4) / (4 - 1);
% 一致性指标RI
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];
% 如果 CR < 0.1 ,则可认为判断矩阵的一致性可以接受; 否则需要对判断矩阵进行修正
CR = CI / RI(4);
disp('一致性指标CI=: ');disp(CI);
disp(‘一致性比例CR=:’);disp(CR);
if CR < 0.10
disp(‘因为CR<0.1,所以该判断矩阵A的一致性可以接受!’);
else
disp('注意CR>=0.1,因此该判断矩阵A需要进行修改! ');
end