[问题类型]:
R的optim套件在某些时候无法找到极值
举例: f(x1,x2) = x1^10000+(1-x1)^10000+x2^10000+(1-x2)^10000, 0<x1<1 0<x2<1
这东西直接手算微分很容易
但optim套件会因为f几乎都趋近于零导致无法迭代
请板上高手大大帮看
谢谢
[软件熟悉度]:
使用者(已经有用R 做过不少作品)
[问题叙述]:
如问题类型所说
[程式范例]:
f = function(x){
x1 = x[1]
x2 = x[2]
return(x1^10000+(1-x1)^10000+x2^10000+(1-x2)^10000)
}
optim(par = c(0.3,0.3), fn = f)
[环境叙述]:
output 结果如下(事实上手算很容易知道极值为0.5 0.5)
$par
[1] 0.3 0.3
$value
[1] 0
$counts
function gradient
9 NA
$convergence
[1] 0
$message
NULL