数学建模题目.docxVIP

PAGE PAGE # 10QD 1000（左），等高线（右）图1最邻近插值地貌图竝性插值1000 1000塢性更值1、山区地貌：在某山区测得一些地点的高程如下表：（平面区域 1200=xv=4000,1200=yv=3600）,试作出该山区的地貌图和等高线图， 并对几种 插值方法进行比较。 10QD 1000 （左），等高线（右） 图1最邻近插值地貌图 竝性插值 1000 1000 塢性更值 3600 1480 1500 1550 1510 1430 1300 1200 980 3200 1500 1550 1600 1550 1600 1600 1600 1550 2800 1500 1200 1100 1550 1600 1550 1380 1070 2400 1500 1200 1100 1350 1450 1200 1150 1010 2000 1390 1500 1500 1400 900 1100 1060 950 1600 1320 1450 1420 1400 1300 700 900 850 1200 1130 1250 1280 1230 1040 900 500 700 Y/x 1200 1600 2000 2400 2800 3200 3600 4000 方法一：利用插值的方法，绘制山区的地貌图和等高线，采用了 5种插值方 法，分别是最邻近插值、线性插值、三次样条插值、立方插值、分段线性插 值，得到如图1-5所示的图像： I I II I 图2线性插值地貌图（左），等高线（右） 三发择条死值 三发择条死值 图3三次样条插值地貌图（左），等高线（右） 3X015001COO500 3X0 1500 1COO 500 0 4i：i：0 立方播葫 图4立方插值地貌图（左），等高线（右） 图5分段线性插值地貌图（左），等高线（右） 比较由以上五种插值方法得到的地貌图和等高线图，可以看出，由于两个 高度之间直线为最短距离，因此利用最邻近插值得到的地貌图和等高线为直线，描述的山地地貌为陡崖 ，对于一般山区的地貌是不符合的；分段线性插 值得到的图像随着分段数目的增多，而更加平缓，棱角更加不明显；禾I」用线性 插值、三次样条插值和立方插值所得到的图像，较为平滑，更加适合描述该区 山地的地貌。 图像绘制程序: x=1200:400:4000; y=1200:400:3600; z=[1130 1250 1280 1230 1040 900 500 700; 1320 1450 1420 1400 1300 700 900 850; 1390 1500 1500 1400 900 1100 1500 1200 1100 1350 1450 1200 1500 1200 1100 1500 1550 1600 1480 1500 1550 1550 1600 1550 1550 1600 1600 1510 1430 1300 1060 950; 1150 1010; 1380 1070; 1600 1550; 1200 980]; 最邻近插值分段线性插值 最邻近插值 分段线性插值 figure(1); meshz(x,y,z) xlabel(X),ylabel(Y),zlabel(Z) title( 网格面) xi=1200:40:4000; yi=1200:40:3600; figure(2) z1i=i nterp2(x,y,z,xi,yi, nearest);% % subplot(1,2,1), surfc(xi,yi,z1i) xlabel(X),ylabel(Y),zlabel(Z) title( 最邻近插值) % subplot(1,2,2),co ntour(xi,yi,z1i,10,r); figure(3) z2i=i nterp2(x,y,z,xi,yi); % subplot(1,2,1), surfc(xi,yi,z2i) xlabel(X),ylabel(Y),zlabel(Z)% title( 分段线性插值) % subplot(1,2,2),co ntour(xi,yi,z2i,10,r); figure(4) z3i=i nterp2(x,y,z,xi,yi,cubic); surfc(xi,yi,z3i) xlabel(X),ylabel(Y),zlabel(Z)% title( 立方插值) 立方插值 figure(5) z4i=i nterp2(x,y,z,xi,yi,spli ne); surfc(xi,yi,z4i) xlabel(X),ylabel(Y),zlabel(Z)% title( 三次样条插值) 三次样条插值% figure(6) z5i=i nterp2