查看单个帖子
旧 2008-07-18, 20:37   #1
rain_fall
初级会员
 
注册日期: 2008-04-02
年龄: 43
帖子: 5
声望力: 0
rain_fall 正向着好的方向发展
默认 [求助]dblquad带入数值出错

我在书上看到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
怎么结果不一致呢,第二个就是我要的结果,请问哪位高手知道是怎么回事吗?第一个式子哪里不对了呢?已经苦恼了一个星期了,拜托各位高手了
rain_fall 当前离线   回复时引用此帖