如何解决鲍威尔方法中的边界Scipy
让功能最小化
f =lambda x: (x+1)**2
在scipy中使用鲍威尔方法
如果我们使用
scipy.optimize.minimize(f,1,method='Powell',bounds=None)
返回是
direc: array([[1.]])
fun: array(0.)
message: 'Optimization terminated successfully.'
nfev: 20
nit: 2
status: 0
success: True
x: array(-1.)
即最小值应为-1。如果我们提供界限
scipy.optimize.minimize(f,bounds=[(0,2)])
再来一次
direc: array([[1.]])
fun: array(0.)
message: 'Optimization terminated successfully.'
nfev: 20
nit: 2
status: 0
success: True
x: array(-1.)
现在错误了!正确答案应该为0。这就像没有将边界考虑在内一样。我正在使用scipy'1.4.1'和python 3.7.6。 有人有线索吗?
解决方法
使用scipy 1.4.x,鲍威尔方法无法处理约束或界限,如您所见here。更新到scipy 1.5.x,它可以处理边界,请参见here:
In [11]: scipy.optimize.minimize(f,x0=1.0,method='Powell',bounds=[(0.0,2.0)])
Out[11]:
direc: array([[1.64428414e-08]])
fun: array(1.)
message: 'Optimization terminated successfully.'
nfev: 103
nit: 2
status: 0
success: True
x: array([2.44756652e-12])
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。