二维离散数据插值后数值积分出现问题
程序如下,不知道为什么就是出错,一维的可以算出来,二维不知道是为什么出错了,调了好久了,还是出不来下面是一维的:
t1=0:0.001:1;
y1=cos(2*t1);
ff=quadspln(t1,y1,0,2);
function y=quadspln(x0,y0,a,b)
f=inline('interp1(x0,y0,x,''spline'')','x','x0','y0');
y=quadl(f,a,b,1e-8,[],x0,y0);
下面是二维的:
t10=0:0.001:1;
t20=0:0.002:2;
=meshgrid(t10,t20);
y10=cos(2*t11)+cos(t22);
f=inline('interp2(t10,t20,y10(.*exp(-1i*t10-1i*t20),t1,t2,''spline'')','t1','t2');
ff=dblquad(f,0,1,1e-8,0,2,1e-8);
运行后
??? Error using ==> fcnchk at 108
FUN must be a function, a valid string expression, or an inline function object.
Error in ==> dblquad at 54
quadf = fcnchk(quadf);
Error in ==> Untitled2 at 9
ff=dblquad('f',0,1,1e-8,0,2,1e-8);
怎么办呢?
请牛人指导,谢谢 dblquad不是这么用的,看看help 水平有限, 但感觉lz是否交代原工程问题, 比较好些!?
页:
[1]