二维插值实验——船在该海域会搁浅吗？

在某海域测得一些点(x,y)处的水深z由下表给出，船的吃水深度为5英尺，在矩形区域（75，200）*（-50，150）里的哪些地方船要避免进入。

 .作出测量点的分布图； 

clear;clc
x=[129 140 103.5 88 185.5 195 105.5 157.5 107.5 77 81 162 162 117.5];
y=[7.5 141.5 23 147 22.5 137.5 85.5 -6.5 -81 3 56.5 -66.5 84 -33.5];
z=[-4 -8 -6 -8 -6 -8 -8 -9 -9 -8 -8 -9 -4 -9];
plot3(x,y,z,'r.','markersize',25);
grid on

 

clear;clc
x=[129 140 103.5 88 185.5 195 105.5 157.5 107.5 77 81 162 162 117.5];
y=[7.5 141.5 23 147 22.5 137.5 85.5 -6.5 -81 3 56.5 -66.5 84 -33.5];
plot(x,y,'+');

 .求出矩形区域（75，200）*（-50，150）的细分网格节点之横、纵坐标向量；

close all
clear;clc;
x0=[129 140 103.5 88 185.5 195 105 157.5 107.5 77 81 162 162 117.5]
y0=[7.5 141.5 23 147 22.5 137.5 85.5 -6.5 -81 3 56.5 -66.5 84 -33.5]
z0=-[4 8 6 8 6 8 8 9 9 8 8 9 4 9]
x=75:1:200
y=-50:1:150
z=griddata(x0,y0',z0,x,y','v4')
contour3(x,y,z,12,'k')

作出海底曲面图形和等高线图

clear;clc;
x0=[129 140 103.5 88 185.5 195 105 157.5 107.5 77 81 162 162 117.5]
y0=[7.5 141.5 23 147 22.5 137.5 85.5 -6.5 -81 3 56.5 -66.5 84 -33.5]
z0=-[4 8 6 8 6 8 8 9 9 8 8 9 4 9]
x=77:1:195
y=-81:1:147
z=griddata(x0,y0',z0,x,y','v4')
meshc(x,y,z)

.作出水深小于5的海域范围。

clear;clc;
x=[129 140 103.5 88 185.5 195 105 157.5 107.5 77 81 162 162 117.5];
y=[7.5 141.5 23 147 22.5 137.5 85.5 -6.5 -81 3 56.5 -66.5 84 -33.5];
z=-[4 8 6 8 6 8 8 9 9 8 8 9 4 9];
x1=75:200;
y1=-50:150;
z1=griddata(x,y,z,x1,y1','v4')
mesh(x1,y1,z1)
k=find(z1<-5);
z1(k)=nan;
mesh(x1,y1,z1)

 

clear;clc;
x=[129 140 103.5 88 185.5 195 105 157.5 107.5 77 81 162 162 117.5];
y=[7.5 141.5 23 147 22.5 137.5 85.5 -6.5 -81 3 56.5 -66.5 84 -33.5];
z=-[4 8 6 8 6 8 8 9 9 8 8 9 4 9];
x1=75:200;
y1=-50:150;
z1=griddata(x,y,z,x1,y1','v4')
mesh(x1,y1,z1)
k=find(z1<-5);
z1(k)=-5;
mesh(x1,y1,z1)

close all
clear;clc;
x0=[129 140 103.5 88 185.5 195 105 157.5 107.5 77 81 162 162 117.5]
y0=[7.5 141.5 23 147 22.5 137.5 85.5 -6.5 -81 3 56.5 -66.5 84 -33.5]
z0=-[4 8 6 8 6 8 8 9 9 8 8 9 4 9]
x=75:1:200
y=-50:1:150
z=griddata(x0,y0',z0,x,y','v4')
contour3(x,y,z,12,'k')