为了账号安全,请及时绑定邮箱和手机立即绑定

不满足约束条件时,成功进行Scipy Optimize.minimize退出

不满足约束条件时,成功进行Scipy Optimize.minimize退出

蛊毒传说 2021-03-29 16:19:10
我一直在使用scipy.optimize.minimize (docs)当我定义一个无法满足约束的问题时,发现一些奇怪的行为。这是一个例子:from scipy import optimize# minimize f(x) = x^2 - 4xdef f(x):    return x**2 - 4*xdef x_constraint(x, sign, value):    return sign*(x - value)# subject to x >= 5 and x<=0 (not possible)constraints = []constraints.append({'type': 'ineq', 'fun': x_constraint, 'args': [1, 5]})constraints.append({'type': 'ineq', 'fun': x_constraint, 'args': [-1, 0]})optimize.minimize(f, x0=3, constraints=constraints)结果输出:fun: -3.0     jac: array([ 2.]) message: 'Optimization terminated successfully.'    nfev: 3     nit: 5    njev: 1  status: 0 success: True       x: array([ 3.])没有满足约束条件的解决方案,但是,使用初始条件作为最佳解决方案,minimum()成功返回。这种行为是故意的吗?如果是这样,如果最佳解决方案不满足约束条件,是否有办法强制失败?
查看完整描述

1 回答

?
千巷猫影

TA贡献1829条经验 获得超7个赞

这似乎是一个错误。我在github上的问题上添加了带有您的示例变体的评论。


如果您使用其他方法(例如COBYLA),则该函数将无法正确找到解决方案:


In [10]: optimize.minimize(f, x0=3, constraints=constraints, method='COBYLA')

Out[10]: 

     fun: -3.75

   maxcv: 2.5

 message: 'Did not converge to a solution satisfying the constraints. See `maxcv` for magnitude of violation.'

    nfev: 7

  status: 4

 success: False

       x: array(2.5)


查看完整回答
反对 回复 2021-04-20
  • 1 回答
  • 0 关注
  • 829 浏览
慕课专栏
更多

添加回答

举报

0/150
提交
取消
微信客服

购课补贴
联系客服咨询优惠详情

帮助反馈 APP下载

慕课网APP
您的移动学习伙伴

公众号

扫描二维码
关注慕课网微信公众号