新手求大侠帮忙看看代码！谢谢！！！

terryniu01

看了篇文献想用matlab来显现它里面的一个重要公式，计算出最后结果。附件是要是实现的公式。里面含有贝塞尔函数，要对它进行0到无穷大的积分。以下是我写的代码：

1. function T=T_3D(L,a,f)
   
2. syms r lamda;
   
3. P=1;
   
4. k=[0.026,65,51.9];
   
5. alpha=[2.08e-5,17.14e-6,13.6e-6];
   
6. omiga=2*pi*f;
   
7. det=(1+i).*sqrt(omiga./(2*alpha));
   
8. sigma=(lamda^2+det.^2).^(1/2);
   
9. b10=k(2)*sigma(2)/(k(1)*sigma(1));
   
10. b21=k(3)*sigma(3)/(k(2)*sigma(2));
    
11. gama10=(b10-1)/(b10+1);
    
12. gama21=(b21-1)/(b21+1);
    
13. part1=1/(k(1)*sigma(1)*(1+b10));
    
14. part2=(1-gama21*exp(-2*sigma(2)*L))/(1+gama10*gama21*exp(-2*sigma(2)*L));
    
15. F=part1*part2*exp(-lamda^2*a^2/4)*besselj(0,lamda*r)*lamda;
    
16. T=P*int(F,lamda,0,inf)
    
17. 
    

复制代码

运行的结果是matlab不报错也不给出结果，比较郁闷。想请大家帮忙看一看是不是不能用int函数来实现，这个公式里面带有贝塞尔0阶函数，以前没有接触过它的积分，还是0到无穷大限，望赐教，万分感激！