含min()函数的积分结果与内部函数积分结果不符的技术问询

你遇到的问题核心在于**min()是标量函数，而integrate要求被积函数支持向量化输入**：

• R的integrate函数为提升计算效率，会向被积函数传入向量形式的采样x值，而非逐个传入标量x。
• min(a,b)的逻辑是返回两个输入的全局最小值——当你传入min(-x+10,10)时，-x+10是一个向量，10是标量，min会把整个向量和10比较，返回整个向量的全局最小值（在x∈[0,10]区间里，这个最小值是0，当x=10时取到）。
• 这就导致第一个积分的被积函数实际变成了x * 0（全局最小值），积分结果自然完全偏离预期。

验证与解决方法 #

• 验证：运行x <- c(0,5,10); min(-x+10,10)，返回的是0，而不是预期的c(10,5,0)；但用pmin(-x+10,10)则会返回c(10,5,0)，这才是逐个元素取最小的正确结果。
• 修复代码：将min替换为向量化的pmin函数，重新计算积分：

  integrate(\(x) {x * pmin(-x+10, 10)}, lower = 0, upper = 10)$value
  # 输出结果：[1] 166.6667
  

  此时结果就和第二个积分完全一致了。

内容的提问来源于stack exchange，提问作者Victor Nielsen