计算随机一致性指标RI作业

1、东南大学数学实验报告学号学010314姓名梅杰成绩实验内容：计算随机一致性指标RI一实验目的计算n=330时的n阶矩阵的随机一致性指标RI二预备知识(1)熟悉随机一致性指标的含义及计算方法(2)熟悉eig、rand等Matlab命令三实验内容与要求用MATLAB编制程序，(要求采用和法计算最大特征值)，分别计算n=330时的n阶矩阵的随机一致性指标RI。CommandWindow'f-NewtoMATLAB?WatchthisVidechseeDbedslorreadGettincaStarted»Rlsatrasd,30)；farn=3:30Kn=20;tiraes=tOO

2、OG;enunt=9876543211/21/31/4i/S1/6i/71/B1/9;lamda=ztEGS(1times),A=ones(iijn);farnuft=l:ti»esfciri=I:nfarj=i+l:nA(ijj)=enwt(ceil(17*ranii(l)：A(Li)=l/A(ivj):A(iji)=;endendV=eie(A):1amda(ninn1三max(¥);endk=snj»(1amda)/1xaes:RI(n)=(k-n)/(n-l);endRIRI=000.51S20.89421.1104L24301.33241.4028t.4

3、520E4875Calunns11through20515Q1,53751.55781.57011.53371.5955I.A0391.6144L6215K6308Coluhns21through301,63571.64141.6458R65031.65511.65991,66221.66631.66881.6723Columns1through10RI=zeros(1,30);forn=3:30%n=20;times=10000;%定义结果输出格式并初始化，RI(1)直接赋值为0%循环计算阶数3到30的随机正互反矩阵的RI%起初以20阶矩阵为例测试%10000个子样，应该够多了吧enum=9

4、876543211/21/31/41/51/61/71/81/9;%矩阵元素从enum中取得lamda=zeros(1,times);A=ones(n,n);%最大特征值向量初始化%初始化相应阶数的矩阵fornum=1:timesfori=1:nforj=i+1:nA(i,j)=enum(ceil(17*rand(1);A(j,i)=1/A(i,j);A(i,i)=1;endendV=eig(A);lamda(num)=max(V);end%循环%把矩阵A赋值为正互反矩阵%矩阵的上半部分从enum中随机取值%矩阵的下半部分与上半部分成倒数%矩阵对角线为1%求得A的特征向量%以最大特征值给lamda向量赋值k=sum(lamda)/times;%最大特征值的平均值endRIRI(n)=(k-n)/(n-1);%得出对应的RI(n)%最后输出RI向量，即130阶矩阵的平均随机一致性指标四实验心得一开始对matlab命令不熟悉，感觉无从下手，通过查阅matlab相关资料，进一步了解matlab后，找到门路和方法，由于对程序和函