【发布时间】:2022-01-05 07:09:56
【问题描述】:
我正在考虑 Nelder Mead 优化技术,我需要返回第一列中具有最小值的矩阵。
我的代码:
fun <- function(x){
3*(sin(0.5+0.25*x[2]*x[1]))*cos(x[1])
}
out <- matrix(NA, nrow=100, ncol=3)
for (i in 1:100) {
x <- runif(1, -7, 7)
y <- runif(1, -7, 7)
x0 <- c(x,y)
res <- optim(x0,fun,method="Nelder-Mead")
out[i,] <- round(c(res$value,res$par) , digits = 5)
}
out
【问题讨论】: