R具有多个参数的双积分

Question

我正在尝试在R中转换以下双整数。

n <- 1
InnerFunc <- function(x, y) {f(y - x)*g(x)}
InnerIntegral <- Vectorize(function(y) {integrate(InnerFunc, -Inf, Inf)$value})
integrate(InnerIntegral, n, Inf)

我试图从另一个question中寻找灵感，但是我对这些论点感到困惑。感谢您的帮助。

Answer 1

感谢G5W向我推荐这篇文章： Numerical Triple Integration in R

我想我现在可以回答自己的问题。

n <- 1
integrate(Vectorize(function(y) {
integrate(Vectorize(function(x) {f(y - x)*g(x)}), -Inf, Inf)$value }), n, Inf)

我用

测试了语法

f <- function(x) {x^2}
g <- function(x) {sin(x)}

integrate(Vectorize(function(y) {
integrate(Vectorize(function(x) {f(y-x)*g(y)}), 0, 1)$value }), 0, pi)

哪个给出{@ 3}}中的3.394678