优化超过3维分段常数函数

Question

我正在研究一个具有三维分段常数函数的模拟项目，并且我试图找到最大化输出的输入。使用optim()或Nelder-Mead算法在 R 中使用SANN似乎最好（他们不要求函数可以区分），但我＆＃39 ; m发现optim()最终确切地返回了我的起始值。这个起始值是使用网格搜索获得的，所以它可能相当不错，但如果它是确切的峰值，我会感到惊讶。

我怀疑optim()没有测试点远离初始猜测，导致所有测试点都给出相同的输出。

1. 这是一个合理的担忧吗？
2. 如何调整optim()在搜索时正在测试的值的广度？