关于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
}