计算随机一致性指标实验报告.doc

PPT
东南大学《数学实验》报告
学号 姓名 成绩
实验内容：计算随机一致性指标RI
一 实验目的
计算n=2~30时的n阶矩阵的随机一致性指标RI
二 预备知识
(1)熟悉随机一致性指标的PPT
东南大学《数学实验》报告
学号 姓名 成绩
实验内容：计算随机一致性指标RI
一 实验目的
计算n=2~30时的n阶矩阵的随机一致性指标RI
二 预备知识
(1)熟悉随机一致性指标的含义及计算方法
(2)熟悉eig、rand等Matlab命令
三 实验内容与要求
用MATLAB或C++编制程序，（要求采用和法计算最大特征值），分别计算n=2~30时的n阶矩阵的随机一致性指标RI
命令：
clear
for k=1:50
for n=1:10000
b=[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];
randommatrix=ceil(17*rand(k));
for i=1:k
for j=i:k
matrix(i,j)=b(randommatrix(i,j));
matrix(j,i)=1/matrix(i,j);
matrix(i,i)=1;
end
end
matrixone=ones(1,k);
for i=1:k
for j=1:k
matrixone(i)=matrixone(i)*matrix(i,j);
end
end
matrixone=matrixone.^(1./k);
matrixone=matrixone./sum(matrixone);
L=sum((matrix*matrixone')./(k*matrixone'));
CI(n)=(L-k)/(k-1);
end
RI(k)=sum(CI)./10000;
end
结果：
2
0
3

4

5
**********
6

7
**********
8

9