隨機(jī)一致性指標(biāo)求解

東南大學(xué)《數(shù)學(xué)實(shí)驗(yàn)》報(bào)告學(xué)號(hào)09008406姓名陳怡東 成績(jī)實(shí)驗(yàn)內(nèi)容：隨機(jī)一致性指標(biāo)求解一實(shí)驗(yàn)?zāi)康挠肕atlab求解隨機(jī)一致性指標(biāo)的方法加深對(duì)隨機(jī)一致性指標(biāo)概念的理解二預(yù)備知識(shí)熟悉層次分析法的含義了解RI的求解方法三實(shí)驗(yàn)內(nèi)容與要求用MATLAB編寫(xiě)程序分別計(jì)算n=2?30是的n的隨機(jī)一致性講演指標(biāo)的值RI，為保證隨機(jī)性，要求每階創(chuàng)建5000個(gè)矩陣四實(shí)驗(yàn)原理及分析層次分析法建模問(wèn)題中，需要用到對(duì)矩陣A的一致性檢驗(yàn)，然后對(duì)于一般的問(wèn)題，尤其當(dāng)考慮實(shí)際因素比較多時(shí)，很難保證判斷A為一致矩陣，因此在計(jì)算矩陣A的最大特征值之時(shí)，需要檢驗(yàn)矩陣A的一致程度。令：丹一I稱CI為一致性指標(biāo)，顯然CI=0是矩陣A為一致矩陣的必要條件，可以看出CI值越大，A的不一致程度越嚴(yán)重。RI是按照下面的方式選取的：對(duì)于固定的n，隨機(jī)構(gòu)造正互反矩陣A，他的元素是從1?9及其倒數(shù)中隨機(jī)選取的，因此A的一致性一般是很差的，取充分大的子樣得到A的最大特征值的平均值k，定義：CR稱為隨機(jī)一致性比率，RI稱為隨機(jī)一致性指標(biāo)。當(dāng)CR<0.1是樣板認(rèn)為矩陣A的不一致程度在容許范圍之內(nèi)，可以用其特征向量作為權(quán)向量代碼作為M文件結(jié)果functionri=RI(n)n=2RI=0n=ceil(n);n=3RI=0.5303sca=[1234567891/91/8n=4RI=0.88451/71/61/51/41/31/2];n=5RI=1.1143fora=2:n;n=6RI=1.2525stat=eye(a);n=7RI=1.3396maxr(5000)=0;n=8RI=1.4022forb=1:5000n=9RI=1.4490n=10RI=1.4841biaoji=ceil(17*rand(a));n=11RI=1.5079biaoji=triu(biaoji);n=12RI=1.5376fori=1:an=13RI=1.5502forj=1:i-1n=14RI=1.5701n=15RI=1.5819stat(j,i)=sca(biaoji(j,i));n=16RI=1.5923n=17RI=1.6084stat(i,j)=1/stat(j,i);n=18RI=1.6132endn=19RI=1.6220endr=eig(stat);maxr(b)=max(r);stat;endk=(sum(maxr))/5000;count=a-1;ri(count)=(k-a)/(a-1);endn=20 RI=1.6306n=21 RI=1.6347n=22 RI=1.6401n=23 RI=1.6482n=24 RI=1.6511n=25 RI=1.6550n=26RI=1.6586n=27 RI=1.6627n=28 RI=1.6662n=29 RI=1.6683n=30RI=1.6732代碼中忘記加入注釋了。六實(shí)驗(yàn)總結(jié)本實(shí)驗(yàn)通過(guò)時(shí)間學(xué)習(xí)了Matlab編程的