東南大學(xué)數(shù)模實驗報告隨機一致性指標求解一、實驗?zāi)康?）掌握用matlab求解隨機一致性指標的方法2）加深對隨機一致性指標概念的理解二、實驗內(nèi)容用matlab或C++編寫程序分別計算n=3-30時的n階矩陣的隨機一致性檢驗指標的值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)三、實驗結(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ù)據(jù)與書上給的前11個誤差不大，由于選用了較為笨拙的循環(huán)算法，使計算高階的矩陣時耗時很長。曲線插值一、實驗?zāi)康?）熟悉一般的曲線插值的方法2）熟悉“\”、polyfit、polyval、interp1、spline、cscvn等Matlab命令3）學(xué)會用常見插值函數(shù)的求解及應(yīng)用二、實驗內(nèi)容(1)已知某平原地區(qū)的一條公路經(jīng)過如下坐標所表示的點，請用樣條插值繪出這條公路（不考慮公路的寬度）。X(m)03050708090120148170180Y(m)80644742486680120121138X(m)202212230248268271280290300312Y(m)160182200208212210200196188186X(m)320340360372382390416430478440Y(m)200184188200202240246280296308X(m)420380360340320314280240200Y(m)334328334346356360392390400(2)對于上表給出的數(shù)據(jù)，估計公路長度程序如下：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