数模论坛

 找回密码
 注-册-帐-号
搜索
热搜: 活动 交友 discuz
查看: 5463|回复: 3

帮帮忙呀,关于e01sff,e01sef函数

[复制链接]
发表于 2003-9-15 05:10:18 | 显示全部楼层 |阅读模式
各位高手帮帮忙呀,,小弟我写了一个matlab的函数,是求画出海底地貌的。但是总是提是我ndefined function or variable 'e01sef'.这是怎么回事呀,,怎么会没有这个函数??我可是找着参考书上一步一步打出来的呀。
程序如下:
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]);
grid,%做深度为5的海底等值线图
[i1,j1]=find(pf<-5);
for k=1;length(i1)
    pf(i1(k),j1(k))=-5;
end
figure(4),meshc(px,py,pf),rotate3d%做水深低于5英尺的部分海地区面图
发表于 2003-9-15 16:14:16 | 显示全部楼层
[fnodes,minnq,rnw,rnq,ifail]=e01sef(x,y,z);%求差值函数的参数
[pf(i,j),ifail]=e01sff(x,y,z,rnw,fnodes,px(j),py(i));%求插值值


  以上两个函数并非MatLab自带的函数,而是二次开发者自己编的,你在你那本书上找找,应有e01sef和e01sff函数的matlab源代码

发表于 2003-9-15 19:00:29 | 显示全部楼层
我的机子上也没有e01sef这个函数
发表于 2003-9-16 01:54:15 | 显示全部楼层
lear
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];
%a=linspace(75,200,100);
%b=linspace(-50,150,200);
[x,y]=meshgrid(75:.5:200,-50:.5:150);
z=griddata(X,Y,Z,x,y,'cubic');
figure(2);
meshz(x,y,z+5);
figure(3);
contour(x,y,z,[-5,5]);
figure(4);
meshz(x,y,z);
rotate3d
同样可以作出海底地貌图,只是效果不是太好(没有课件中的好)
您需要登录后才可以回帖 登录 | 注-册-帐-号

本版积分规则

小黑屋|手机版|Archiver|数学建模网 ( 湘ICP备11011602号 )

GMT+8, 2024-11-26 23:17 , Processed in 0.056385 second(s), 19 queries .

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表