數(shù)學(xué)建模實驗報告之曲線插值

PAGEPAGE3東南大學(xué)《數(shù)學(xué)實驗》報告學(xué)號09008226姓名畢斌成績實驗內(nèi)容：曲線插值一實驗?zāi)康某Ｒ姴逯岛瘮?shù)的求解及應(yīng)用二預(yù)備知識(1)熟悉一般的曲線插值的方法(2)熟悉“\”、polyfit、polyval、interp1、spline、cscvn等Matlab命令三實驗內(nèi)容與要求(1)已知某平原地區(qū)的一條公路經(jīng)過如下坐標(biāo)所表示的點，請用樣條插值繪出這條公路（不考慮公路的寬度）。X(m)03050708090120148170180Y(m)80644742486680120121138X(m)202212230248268271280290300312Y(m)160182200208212210200196188186X(m)320340360372382390416430478440Y(m)200184188200202240246280296308X(m)420380360340320314280240200Y(m)334328334346356360392390400完整代碼：functioninterpolation%路線x軸上單調(diào)遞增到最大值road_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];%插值0到478x1=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);end%做圖并保持到下一個，一起繪制圖形plot(road_x1,road_y1,'.',x1,y1);holdon;%x軸上單調(diào)遞減，第二段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];%插值200到478x2=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);end%繪圖plot(road_x2,road_y2,'.',x2,y2);%plot(road_x2,road_y2