如何在R中计算含变量上限的二重积分？

针对内层积分上限包含变量的二重积分，有两种实用的计算方法：

方法一：嵌套使用基础包的integrate()函数 #

R基础包的integrate()可处理单变量积分，通过嵌套调用，先计算每个x对应的内层y积分，再对x完成外层积分。

假设你的积分形式为：

∫(x=0到0.5) [∫(y=0到x) f(x,y) dy] dx

其中f(x,y)为被积函数，对应目标结果的被积函数为y/2，代码示例如下：

# 定义内层积分：给定x，计算y的积分结果
inner_integral <- function(x) {
  integrate(function(y) y/2, lower = 0, upper = x)$value
}

# 计算外层积分：对x积分得到最终结果
result <- integrate(inner_integral, lower = 0, upper = 0.5)
print(result$value)

运行后将得到预期的0.0104166666666667。

方法二：使用cubature包的自适应积分函数 #

针对更复杂的多重积分，cubature包的adaptIntegrate()支持直接处理多维积分，且能灵活设置变量约束。

首先安装并加载包：

install.packages("cubature")
library(cubature)

接着定义被积函数（输入为二维向量，x[1]对应原变量x，x[2]对应原变量y），通过condition参数指定y ≤ x的约束：

# 定义被积函数
f <- function(x) {
  x[2]/2
}

# 计算二重积分
result <- adaptIntegrate(f, lower = c(0, 0), upper = c(0.5, 0.5), 
                         condition = function(x) x[2] <= x[1])
print(result$integral)

运行后同样会输出目标结果。

注意事项 #

• 若被积函数存在奇点或不连续，可调整adaptIntegrate()的tol参数提升计算精度。
• 嵌套integrate()适合简单二重积分场景，cubature包则更适配高维或复杂约束的积分需求。

内容的提问来源于stack exchange，提问作者Homer Jay Simpson