在R中迭代处理数值积分发散问题：逐步调下限至收敛

针对你遇到的积分下限为0时发散的问题，可以通过循环测试逐步增大的下限值，结合R的错误处理机制来找到能让积分收敛的最小下限。以下是具体实现代码：

代码实现 #

# 定义被积函数（将H和a作为参数传入，避免全局变量依赖）
integrand <- function(u, H, a) {
  log(u) * u^(2 * (H - 1)) * (exp(a * u) - exp(-a * u))
}

# 设置参数值
H <- 0.1
a <- 0.1
upper_limit <- 100

# 生成候选下限序列：从1e-1到1e-13，每次步长为10倍缩小（对应一位小数级调整）
lower_candidates <- 10^(-1:-13)

# 初始化结果变量
converged_result <- NULL
successful_lower <- NULL

# 遍历测试每个候选下限
for (lower in lower_candidates) {
  # 用tryCatch捕获积分的错误与警告
  integration_attempt <- tryCatch(
    {
      # 尝试执行积分
      integrate(f = integrand, lower = lower, upper = upper_limit,
                rel.tol = .Machine$double.eps^0.01, subdivisions = 1e6,
                H = H, a = a)
    },
    error = function(e) {
      # 积分发散或报错时，打印信息并返回NULL
      message(sprintf("下限 %.1e 积分失败：%s", lower, e$message))
      return(NULL)
    },
    warning = function(w) {
      # 出现警告时（比如精度不足），打印信息并返回积分结果
      message(sprintf("下限 %.1e 积分出现警告：%s", lower, w$message))
      return(integrate(f = integrand, lower = lower, upper = upper_limit,
                       rel.tol = .Machine$double.eps^0.01, subdivisions = 1e6,
                       H = H, a = a))
    }
  )
  
  # 检查是否成功得到有效结果
  if (!is.null(integration_attempt)) {
    converged_result <- integration_attempt
    successful_lower <- lower
    message(sprintf("找到收敛下限：%.1e", lower))
    break  # 找到第一个收敛下限后停止循环
  }
}

# 输出最终结果
if (!is.null(converged_result)) {
  cat("\n成功收敛的积分结果：\n")
  print(converged_result)
  cat("使用的下限：", successful_lower, "\n")
} else {
  cat("\n所有测试的下限均无法使积分收敛\n")
}

关键说明 #

• 错误处理：使用tryCatch捕获integrate函数抛出的错误（如下限导致积分发散），避免程序直接崩溃，同时可以继续测试下一个下限。
• 候选下限序列：生成1e-1到1e-13的序列，每次下限缩小10倍，符合你“一位小数步长调整”的需求。
• 参数传递：将H和a作为参数传入被积函数，避免依赖全局变量，让代码更健壮。
• 终止逻辑：找到第一个收敛的下限就停止循环，若需要测试所有可能的下限，可移除break语句并记录所有收敛结果。

内容的提问来源于stack exchange，提问作者Mr Frog