disp('请输入判断矩阵A(n阶)');
A=input('A=');
%矩阵改良
B=log10(A);
C=zeros(5);
for j=1:5
for i=1:5
s=0;
for k=1:5
a=B(j,k)-B(i,k);
s=s+a;
end
C(j,i)=1/5*s;
end
end
A1=10.^C;
%改良后的矩阵为A1,以下为求特征值和特征向量.
[n,n]=size(A1);
x=ones(n,100);
y=ones(n,100);
m=zeros(1,100);
m(1)=max(x(:,1));
y(:,1)=x(:,1);
x(:,2)=A1*y(:,1);
m(2)=max(x(:,2));
y(:,2)=x(:,2)/m(2);
p=0.0001;
i1=2;
k1=abs(m(2)-m(1));
while k1>p
i1=i1+1;
x(:,i1)=A1*y(:,i1-1);
m(i1)=max(x(:,i1));
y(:,i1)=x(:,i1)/m(i1);
k1=abs(m(i1)-m(i1-1));
end
a1=sum(y(:,i1));
w1=y(:,i1)/a1;
t1=m(i1);
disp(w1);disp(t1);
%以下是一致性检验
CI=(t1-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 disp('此矩阵的一致性可以接受!');
disp('CI=');disp(CI);
disp('CR=');disp(CR);
end
