在某海域测得一些点(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')