二维离散数据插值后数值积分出现问题

susu0526

程序如下，不知道为什么就是出错，一维的可以算出来，二维不知道是为什么出错了，调了好久了，还是出不来
下面是一维的：
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;
[t11,t22]=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);


怎么办呢？
请牛人指导，谢谢

zhouyang664

dblquad不是这么用的，看看help

ChaChing

水平有限, 但感觉lz是否交代原工程问题, 比较好些!?