问题摘要:
在
matlab中计算含有第一类修正的贝塞尔函数的积分时,出现了在部分值能计算出结果,部分值无法计算出结果的问题。
一、变量定义:
syms x
y=1.0+1i*x; %积分路径
xitap、E、B、R、me、e、yipesiu0、c、lamuda
gamma~E
yita=e*B*R*xitap/(me*c*gamma);
aa=4*pi*R/(3*lamuda*gamma^3*sqrt(1+yita^2));
a=aa*yita/(1+yita^2);
二、积分计算
xx=int((1-2*y^2)*besseli(0,a*y^3)*exp(-1.5*aa*(y-y^3/3))*diff(y)/y,x,-inf,inf);
所定义的变量中,B、R、me、e、yipesiu0、c都取定值,E=30,xitap=0.19时,lamuda=0.5能计算出积分,但lamuda=0.9时却计算不出来,请问可能会是哪些原因导致的?