最优化问题

nntx

f<-function(x)x^2

optimize(f,lower=-1,upper=1)

为什么得不到最小值呢，请问有些什么好方法求global extreme value呢，谢谢

Ihavenothing

回复第1楼的 nntx：

你不是已经得到它了吗？

nntx

回复第2楼的 Ihavenothing：这个得出来的不对啊，很显然这个最小值应该是x=0,y=0,可是这个得出来不对。麻烦您看看，谢谢

suckbunny

回复第3楼的 nntx：

<br />
f<-function(x)x^2<br />
optimize(f,lower=-1,upper=1)<br />
$minimum<br />
[1] 0</p>
<p>$objective<br />
[1] 0<br />

</p>

Ihavenothing

回复第3楼的 nntx：

suckbunny的结果就是啊，要不你把你的程序和输出结果贴出来。

suckbunny

回复第5楼的 Ihavenothing：其实站长你可以称呼我4楼的[s:18]

youngtf

是浮点数的问题么。。。

nan.xiao

<br />
f <- function(x) x^2<br />
optimize(f, lower = -1, upper = 1)<br />
$minimum<br />
[1] -2.775558e-17</p>
<p>$objective<br />
[1] 7.70372e-34<br />

</p>

R 3.0.0 x86_64 linux-gnu

nntx

回复第5楼的 Ihavenothing：麻烦您看8楼，我得到的结果也是8楼。谢谢

nntx

回复第5楼的 Ihavenothing：请问您是如何设置的精度呢？我发现可能是精度的问题，但我试了试还是不对。谢谢了

helend

回复第10楼的 nntx：是浮点数的精度问题。如果结果的绝对值小于.Machine$double.eps就可以认为是等于0了。

<br />
> f <- function(x) x^2<br />
> a <- optimize(f, lower = -1, upper = 1)<br />
><br />
> abs(unlist(a)) < .Machine$double.eps<br />
  minimum objective<br />
     TRUE      TRUE<br />
> unlist(a) == 0<br />
  minimum objective<br />
    FALSE     FALSE<br />

</p>

如果你一定想让结果等于0，可以把tol设定的更小一些，例如

<br />
> a <- optimize(f, lower = -1, upper = 1,<br />
+ tol = .Machine$double.eps ^ 3)<br />
><br />
> unlist(a) == 0<br />
  minimum objective<br />
     TRUE      TRUE<br />

</p>