求助matlab探讨：quad和quadl为何积分结果相差极大？

soberprogress

1. clear all;
   
2. syms n r lamda z rho alpha w;
   
3. k=2*pi/lamda;
   
4. sgb=exp(-(r/w)^n)*besselj(0,alpha*r);%SGB
   
5. intgrnd=k/z*exp(i*k*rho^2/(2*z))*sgb*exp(-i*k*r^2/(2*z))*besselj(0,k*rho*r/z)*r;
   
6. 
   
7. N=8;% the order of SG
   
8. Z=100:2:106;
   
9. Rho=linspace(-0.01,0.01);
   
10. intgrnd=subs(intgrnd,{n lamda z alpha w},{N,1550*10^-9,Z,600,0.05});%关于Z的向量
    
11. 
    
12. intValue=zeros(length(Z),length(Rho));%预分配空间
    
13. for iValue=1:length(Z)
    
14. disp(iValue);
    
15. intgrndRho=subs(intgrnd(iValue),rho,Rho);
    
16. for jValue=1:length(Rho)
    
17. disp(jValue);
    
18. intf=inline(vectorize(char(intgrndRho(jValue))));
    
19. intValue(iValue,jValue)=abs(quad(intf,0,1));
    
20. end;%内循环结束
    
21. plot(Rho,intValue.^2);
    
22. hold on
    
23. end %外循环结束

复制代码

如果将上面程序中的quad 换为quadl ，则得到很好的图像，但是数量级不正确，请问咋回事？

sigma665

应该差不了多少
quadl更精确些
你的积分方程是不是奇异了

soberprogress

都可以给出结果，就是数量级相差极大！
试试看，就明白了！

一个图形漂亮，但是数量级不对，一个数量级正确但是图形不对！

[ 本帖最后由 soberprogress 于 2008-10-15 10:44 编辑 ]

sigma665

quad(intf,0,1,1e-15)时，会提示奇异
所以应该2个算的都不对

我倾向于相信quadl

soberprogress

我用的是matlab6.0 没有提示任何错误啊！:@( 谢谢您！
是不是因为你设置了精度，才会出现这样的结果，我也倾向于相信quadl，因为它精度高。

sigma665

我设置了精度

soberprogress

我给 quad命令设置了 精度后，很好！