数学建模 海底测量.docVIP

例9.5（海底测量）下表给出水面直角坐标（x,y）处水深z,这是在低潮时测得的。如果船的吃水深度为5米，试问在矩形域75x200,-50y150中行船应避免进入那些区域？ x(m) 129 140 108 88 185 195 105 y(m) 7 141 28 147 22 137 85 z(m) 4 8 6 8 6 8 8 x(m) 157 107 77 145 162 162 117 y(m) -6 -81 3 45 -66 84 -38 z(m) 9 9 8 8 9 4 9 解 我们先看看测量点的位置。 clear;close; x=[129 140 108 88 185 195 105 157 107 77 145 162 162 117]; y=[7 141 28 147 22 137 85 -6 -81 3 45 -66 84 -38]; plot(x,y,o); 这是一批不规则数据。由于没有先验函数，我们使用插值法。为了使结果更为直观, 我们考虑将z的数据转化为相对于海面的高度。下面绘出所考虑区域海底地形图,可以清晰看见在（129,7.5）和（162,84）附近各有一块暗礁（图9.11）。 z=[4 8 6 8 6 8 8 9 9 8 8 9 4 9]; h=-z; xi=75:5:200;yi=[-50:10:150]; HI=griddata(x,y,h,xi,yi,cubic); mesh(xi,yi,HI);view(-60,30);