Warning: file_get_contents(/data/phpspider/zhask/data//catemap/1/wordpress/12.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181
关于R中的积分给出错误“;积分可能是发散的”;_R_Integration - Fatal编程技术网

关于R中的积分给出错误“;积分可能是发散的”;

关于R中的积分给出错误“;积分可能是发散的”;,r,integration,R,Integration,所讨论的积分是: integrand<-function(y){ exp(-sqrt(2*y + alpha^2)*abs(x))/ (pi^2 * y * ((besselJ(delta*sqrt(2*y), lambda))^2) + (besselY(delta*sqrt(2*y), lambda))^2) } integral<-function(x){integrate(integrand, lower=0, upper=Inf, subdivisions=2

所讨论的积分是:

integrand<-function(y){
   exp(-sqrt(2*y + alpha^2)*abs(x))/ (pi^2 * y * ((besselJ(delta*sqrt(2*y), lambda))^2)
   + (besselY(delta*sqrt(2*y), lambda))^2)
}

integral<-function(x){integrate(integrand, lower=0, upper=Inf, subdivisions=20000)$value}

被积函数好的,我自己解决了

贝塞尔函数的加法缺少一对括号;该函数应为:

integrand<-function(y){
   exp(-sqrt(2*y + alpha^2)*abs(x))/ (pi^2 * y * (((besselJ(delta*sqrt(2*y), lambda))^2)
   + (besselY(delta*sqrt(2*y), lambda))^2))
}

integral<-function(x){
   integrate(integrand, lower=0, upper=Inf, subdivisions=20000)$value
}
被积函数
integrand<-function(y){
   exp(-sqrt(2*y + alpha^2)*abs(x))/ (pi^2 * y * (((besselJ(delta*sqrt(2*y), lambda))^2)
   + (besselY(delta*sqrt(2*y), lambda))^2))
}

integral<-function(x){
   integrate(integrand, lower=0, upper=Inf, subdivisions=20000)$value
}