对于R语言，optimize函数的一点疑问

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

今天看到optimize函数是用来求极值的，就找了个函数试了试，用R语言写了下，总共也就那么几行，
f=function(x) x*sin(10*pi*x)+1
curve(f,xlim=c(-2,1))
optimize(f,c(-2,1),tol=0.0001,maximum=T)
是求（-2，1）的极大值，随便画了张图，结果问题来了
这是函数给的极大值：
$maximum[1] -1.050968$objective[1] 2.050482
这是画出的图，很明显不在-1.05处啊， 请问，这是为什么呢，做极小值也是一样，也是不对的。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：OPTIMIZE Optim R语言 OPT IMI function

在R中输入?optimize，好好看看其用法和结果
maximum是指xlim的值
objective是函数的值，也就是你所求的值

hubifeng？ 发表于 2015-12-2 22:33
在R中输入?optimize，好好看看其用法和结果
maximum是指xlim的值
objective是函数的值，也就是你所求的值

这个我知道啊， 问题是他给出的值，你从图像上看的话，明显不是极大值啊

The first evaluation of f is always at x_1 = a + (1-φ)(b-a) where (a,b) = (lower, upper) and phi = (sqrt(5) - 1)/2 = 0.61803.. is the golden section ratio. Almost always, the second evaluation is at x_2 = a + phi(b-a). Note that a local minimum inside [x_1,x_2] will be found as solution, even when f is constant in there, see the last example.

上面这段出自optimize函数的帮助。所以是optimize函数的算法问题。

楼主找到解决方法了吗

suimong 发表于 2015-12-4 02:57
上面这段出自optimize函数的帮助。所以是optimize函数的算法问题。

有什么方法可以避免这个问题吗

那么多局部极值点，optimize怎么能一下子就找出全局最优点？你试个二次函数，肯定行。两种方法：1. 多试一些初始值。2. 选用其它方法。