数模实验报告
计算n阶矩阵的随机一致性指标RI
实验内容:
用matlab或C++编写程序分别计算n=3-30时的n阶矩阵的随机一致性检验指标的值RI。
>> RI=zeros(1,30);
for n=3:30
%n=20;
times=10000;
enum=[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];
L=zeros(1,times);
A=ones(n,n);
for num=1:times
for i=1:n
for j=i+1:n
A(i,j)=enum(ceil(17*rand(1)));
A(j,i)=1/A(i,j);
A(i,i)=1;
end
正则化一个5 5随机矩阵 end
V=eig(A);
L(num)=max(V);
end
k=sum(L)/times;
RI(n)=(k-n)/(n-1);
end
RI
RI =
Columns 1 through 15
0 0 0.5104 0.8845 1.1108 1.2493 1.3416 1.4046 1.4530 1.4851 1.5152 1.5360 1.5542 1.5700 1.5851
Columns 16 through 30
1.5947 1.6068 1.6135 1.6212 1.6283 1.6346 1.6402 1.6471 1.6508 1.6552 1.6601 1.6630 1.6661 1.6699 1.6729
>>
版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198,我们将在24小时内删除。
发表评论