解决R中的数值不稳定积分

问题描述

我有一个功能想集成到x和-Inf之间的Inf。我在R中使用函数integrate。但是，我确实收到一条错误消息，说Non-finite function value。

test_func<-function(x,t,A,sigma,y){
  
  n1 = exp(-((y - x)^2)/(4 * sigma * t)) + exp(-((y + x)^2)/(4 * sigma * t))
  d1 = sqrt(4 * pi * sigma * t)
  a2 = 2 * A /(sigma) * exp((4 * A *(y + x +(4 * A * t)))/(sigma))
  b2 = pracma::erfc(((y + x)+(8 * A * t))/(2 * sqrt(sigma * t)))
  res = (n1/d1) - (a2 * b2)
  return(res)
}

t<-1000
A<-0.01
sigma<-1
y<-1

integrate(test_func,lower = -Inf,upper = Inf,y)$value

我确实收到以下错误：

Error in integrate(test_func,: 
  non-finite function value

我知道有一些方法可以使这种类型的积分数值稳定，但是我不知道，在这种情况下哪种方法可以提供帮助。

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）