关于程序员:一致性检验方法

45次阅读

共计 797 个字符，预计需要花费 2 分钟才能阅读完成。

一致性测验 是指对由不同样本计算的各平均值或方差进行测验。

disp('请输出判断矩阵 A(n 阶)');
A=input('A=');
[n,n]=size(A);
x=ones(n,100);
y=ones(n,100);
m=zeros(1,100);
m(1)=max(x(:,1));
y(:,1)=x(:,1);
x(:,2)=A*y(:,1);
m(2)=max(x(:,2));
y(:,2)=x(:,2)/m(2);
p=0.0001;i=2;k=abs(m(2)-m(1));
while  k>p
  i=i+1;
  x(:,i)=A*y(:,i-1);
  m(i)=max(x(:,i));
  y(:,i)=x(:,i)/m(i);
  k=abs(m(i)-m(i-1));
end
a=sum(y(:,i));
w=y(:,i)/a;
t=m(i);
disp(w);
         % 以下是一致性测验
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];
CR=CI/RI(n);
if CR<0.10
    disp('此矩阵的一致性能够承受!');
    disp('λmax=');disp(t);
    disp('CI=');disp(CI);
    disp('CR=');disp(CR);
end


A=[1,2,3,5,7,9,8,5,6;
1/2,1,2,3,6,8,7,4,5;
1/3,1/2,1,2,5,7,6,3,4;
1/5,1/3,1/2,1,3,6,5,2,3;
1/7,1/6,1/5,1/3,1,3,3,1/2,1/2;
1/9,1/8,1/7,1/6,1/3,1,1/2,1/5,1/3;
1/8,1/7,1/6,1/5,1/3,2,1,1/4,1/2;
1/5,1/4,1/3,1/2,2,5,4,1,2;
1/6,1/5,1/4,1/3,2,3,2,1/2,1;]

正文完