向dblquad中带入参数出错

liujian

我在书上看到dblquad的用法如下
q = dblquad(fun,xmin,xmax,ymin,ymax,tol,method,p1,p2,…)   %将可选参数p1,p2,..等传递给函数fun(x,y,p1,p2,…)。若tol=[]，method=[]，则使用缺省精度和算法quad

我需要求得A  B 分别取0-100时的积分值，其中在A=0，B=0的时候
A=0；B=0；t=dblquad(inline('(x.^2+y.^2<=1).*exp(i*(-0.346181.*(A.*x+B.*y)-0.1730905.*(x.^2+y.^2)))'),-1,1,-1,1,0.001,@quad,A,B);
result=((abs(t)).^2)/(pi.^2)
得到的结果是
result =
   1.62113893827740

如果我把A=0  B=0直接带入的话也就是下面的式子
t=dblquad(inline('(x.^2+y.^2<=1).*exp(i*(-0.1730905.*(x.^2+y.^2))+0)'),-1,1,-1,1,0.001,@quad);
result=((abs(t)).^2)/(pi.^2)
得到的结果是
result =
   0.98525945483360

怎么结果不一样呢 ，第二个就是我要的结果，请问哪位高手知道是怎么回事吗？第一个式子哪里不对了吗？实在是不知道哪里不对了？已经苦恼了一个星期了，拜托各位高手了

sigma665

第一个'(x.^2+y.^2<=1).*exp(i*(-0.346181.*(A.*x+B.*y)-0.1730905.*(x.^2+y.^2)))'，用到了'   '，里面的A B不是你原来的A B,也就是说，不等于0。
如果非要用这种表示，可以参考帖子：关于带参数积分，自行搜索


第2个是正确的