e01sef.m

matlab经典算法的程序-数据拟合

function [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