?? sea.m
字號:
clear;
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,'+');
z=[-4 -8 -6 -8 -6 -8 -8 -9 -9 -8 -8 -9 -4 -9];
[fnodes,minnq,rnw,rnq,ifail]=e01sef(x,y,z);%求插值函數的參數
nx=100;
px=linspace(75,200,nx);
ny=200;
py=linspace(-50,150,ny);
for i=1:ny
for j=1:nx
[pf(i,j),ifail]=e01sff(x,y,z,rnw,fnodes,px(j),py(i));%求插值值
end
end
figure(2)
meshz(px,py,pf+5)%作海底地貌圖
figure(3)
contour(px,py,pf,[-5 -5]) %作深度為5的海底等值線圖
grid
[i1,j1]=find(pf<-5);
for k=1:length(1)
pf(i1(k),j1(k))=-5;
end
figure(4)
meshc(px,py,pf) %作水深低于5英尺的部分海底曲面圖
rotate3d
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -