数学建模竞赛-海底测量

例（海底测量）下表给出水面直角坐标（x,y）处水深z,这是在低潮时测得的。如果船的吃水深度为5米，试问在矩形域75

解我们先看看测量点的位置。

>>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);