R中集成功能的优化

Posted 2023-03-24

【中文标题】R中集成功能的优化

【英文标题】：Optimization of an integrated function in R

【发布时间】：2020-12-09 06:23:33

【问题描述】：

我想为下面的积分找到a 的最佳值，它是x 的函数

  integrand <- function(x,a) 
          D=a/((x+1)*sqrt(x+a))
          D
        

我可以对a 的固定值进行积分。例如，如果 a=5

> integrate(integrand,0,5,a=5)$value
[1] 3.490687

但我想使用 R 中的 optim() 函数或任何可用的内置优化函数找到最佳值。我试过下面的代码，不幸的是它不起作用，

optim(5,integrate(integrand, lower = 0, upper = 5))

感谢任何帮助。

【问题讨论】：

如果积分是某个函数下面的区域，为什么要优化一个数字？

【参考方案1】：

您可以使用“优化”功能，该功能非常适合单参数优化：

optimize(f = integrand, interval = c(0, 5), a = 5)

【讨论】：

我怎么知道a的值是5？我以它为例。我的目标是找到a的最佳值