計(jì)算隨機(jī)一致性指標(biāo)RI作業(yè)

1、東南大學(xué)數(shù)學(xué)實(shí)驗(yàn)報告學(xué)號學(xué)010314姓名梅杰成績實(shí)驗(yàn)內(nèi)容：計(jì)算隨機(jī)一致性指標(biāo)RI一實(shí)驗(yàn)?zāi)康挠?jì)算n=330時的n階矩陣的隨機(jī)一致性指標(biāo)RI二預(yù)備知識(1)熟悉隨機(jī)一致性指標(biāo)的含義及計(jì)算方法(2)熟悉eig、rand等Matlab命令三實(shí)驗(yàn)內(nèi)容與要求用MATLAB編制程序，(要求采用和法計(jì)算最大特征值)，分別計(jì)算n=330時的n階矩陣的隨機(jī)一致性指標(biāo)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;%定義結(jié)果輸出格式并初始化，RI(1)直接賦值為0%循環(huán)計(jì)算階數(shù)3到30的隨機(jī)正互反矩陣的RI%起初以20階矩陣為例測試%10000個子樣，應(yīng)該夠多了吧enum=9

4、876543211/21/31/41/51/61/71/81/9;%矩陣元素從enum中取得lamda=zeros(1,times);A=ones(n,n);%最大特征值向量初始化%初始化相應(yīng)階數(shù)的矩陣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%循環(huán)%把矩陣A賦值為正互反矩陣%矩陣的上半部分從enum中隨機(jī)取值%矩陣的下半部分與上半部分成倒數(shù)%矩陣對角線為1%求得A的特征向量%以最大特征值給lamda向量賦值k=sum(lamda)/times;%最大特征值的平均值endRIRI(n)=(k-n)/(n-1);%得出對應(yīng)的RI(n)%最后輸出RI向量，即130階矩陣的平均隨機(jī)一致性指標(biāo)四實(shí)驗(yàn)心得一開始對matlab命令不熟悉，感覺無從下手，通過查閱matlab相關(guān)資料，進(jìn)一步了解matlab后，找到門路和方法，由于對程序和函