如何解决解决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 举报,一经查实,本站将立刻删除。