東南大學(xué)數(shù)模實(shí)驗(yàn)報(bào)告

東南大學(xué)數(shù)模實(shí)驗(yàn)報(bào)告隨機(jī)一致性指標(biāo)求解一、實(shí)驗(yàn)?zāi)康?）掌握用matlab求解隨機(jī)一致性指標(biāo)的方法2）加深對(duì)隨機(jī)一致性指標(biāo)概念的理解二、實(shí)驗(yàn)內(nèi)容用matlab或C++編寫程序分別計(jì)算n=3-30時(shí)的n階矩陣的隨機(jī)一致性檢驗(yàn)指標(biāo)的值RI。程序如下：B=[1/9,1/8,1/7,1/6,1/5,1/4,1/3,1/2,1,1,2,3,4,5,6,7,8,9]k=0n=3forp=1:1000A=zeros(n)fori=1:nforj=i:nb=randi(18,1)A(i,j)=B(1,b)A(j,i)=1/A(i,j)A(i,i)=1endendk=k+max(eig(A))endk=k/pr=(k-n)/(n-1)三、實(shí)驗(yàn)結(jié)果改變程序中n的值，求的結(jié)果為：n12345678910RI000.530.881.101.251.341.401.451.49n11121314151617181920RI1.511.541.551.571.581.601.611.621.631.63n21222324252627282930RI1.641.641.651.661.671.671.681.681.681.68四、實(shí)驗(yàn)分析實(shí)驗(yàn)所得數(shù)據(jù)與書上給的前11個(gè)誤差不大，由于選用了較為笨拙的循環(huán)算法，使計(jì)算高階的矩陣時(shí)耗時(shí)很長(zhǎng)。曲線插值一、實(shí)驗(yàn)?zāi)康?）熟悉一般的曲線插值的方法2）熟悉“\”、polyfit、polyval、interp1、spline、cscvn等Matlab命令3）學(xué)會(huì)用常見插值函數(shù)的求解及應(yīng)用二、實(shí)驗(yàn)內(nèi)容(1)已知某平原地區(qū)的一條公路經(jīng)過如下坐標(biāo)所表示的點(diǎn)，請(qǐng)用樣條插值繪出這條公路（不考慮公路的寬度）。X(m)03050708090120148170180Y(m)80644742486680120121138X(m)202212230248268271280290300312Y(m)160182200208212210200196188186X(m)320340360372382390416430478440Y(m)200184188200202240246280296308X(m)420380360340320314280240200Y(m)334328334346356360392390400(2)對(duì)于上表給出的數(shù)據(jù)，估計(jì)公路長(zhǎng)度程序如下：functioninterpolationroad_x1=[0,30,50,70,80,90,120,148,170,180,202,212,230,248,268,271,280,290,300,312,320,340,360,372,382,390,416,430,478];road_y1=[80,64,47,42,48,66,80,120,121,138,160,182,200,208,212,210,200,196,188,186,200,184,188,200,202,240,246,280,296];x1=0:478;y1=interp1(road_x1,road_y1,x1,'spline');length1=0;fori=0:477y_i=abs(interp1(road_x1,road_y1,i+1)-interp1(road_x1,road_y1,i));length1=length1+sqrt(1+(y_i)^2);endplot(road_x1,road_y1,'.',x1,y1);holdon;road_x2=[478,440,420,380,360,340,320,314,280,240,200];road_y2=[296,308,334,328,334,346,356,360,392,390,400];x2=200:478;y2=interp1(road_x2,road_y2,x2,'spline');length2=0;forj=200:477y_j=abs(interp1(road_x2,road_y2,j+1)-interp1(road_x2,road_y2,j));length2=length2+sqrt(1+(y_j)^2);endplot(road_x2,road_y2,'.',x2,y2);h