1、山区地貌：在某山区测得一些地点的高程如下表：平面区域1200x40001200y3600，试作出该山区的地貌图和等高线图，并对几种插值方法进行比较。
3600320028002400200016001200Yx
14801500150015001390132011301200
15001550151014301300120098015501600155016001600160015501200110015501600155013801070120011001350145012001150101015001500140090011001060950145014201400130070090085012501280123010409005007001600200024002800320036004000
方法一：利用插值的方法，绘制山区的地貌图和等高线，采用了5种插值方法，分别是最邻近插值、线性插值、三次样条插值、立方插值、分段线性插值，得到如图15所示的图像：
图1最邻近插值地貌图（左），等高线（右）
图2线性插值地貌图（左），等高线（右）
1
f图3三次样条插值地貌图（左），等高线（右）
图4立方插值地貌图（左），等高线（右）
图5分段线性插值地貌图（左），等高线（右）比较由以上五种插值方法得到的地貌图和等高线图，可以看出，由于两个高度之间直线为最短距离，因此利用最邻近插值得到的地貌图和等高线为直
2
f线，描述的山地地貌为陡崖，对于一般山区的地貌是不符合的；分段线性插值得到的图像随着分段数目的增多，而更加平缓，棱角更加不明显；利用线性插值、三次样条插值和立方插值所得到的图像，较为平滑，更加适合描述该区山地的地貌。
图像绘制程序：x12004004000y12004003600z11301250128012301040900500700132014501420140013007009008501390150015001400900110010609501500120011001350145012001150101015001200110015501600155013801070150015501600155016001600160015501480150015501510143013001200980
figure1meshzxyzxlabelXylabelYzlabelZtitle网格面
xi1200404000yi1200403600
figure2z1ii
terp2xyzxiyi
earest最邻近插值subplot121surfcxiyiz1ixlabelXylabelYzlabelZtitle最邻近插值subplot122co
tourxiyiz1i10r
figure3z2ii
terp2xyzxiyisubplot121surfcxiyiz2ixlabelXylabelYzlabelZ分段线性插值title分段线性插值subplot122co
tourxiyiz2i10r
figure4z3ii
terp2xyzxiyicubic
3
fsurfcxiyiz3ixlabelXylabelYzlabelZ立方插值title立方插值
figure5z4ii
terp2xyzxiyispli
esurfcxiyiz4ixlabelXylabelYzlabelZ三次样条插值title三次样条插值
figure6z5ii
terp2xyzxiyili
earsurfcxiyiz4ixlabelXylabelYzlabelZ线性插值title线性插值
figure7subplot321co
tourxiyiz1i10rsubplot322co
tourxiyiz2i10r