如何解决Gekko Solver 不满足约束给出错误的解决方案
我使用 Gekko 求解器来优化函数,但即使在不满足给定约束的简单问题中,它也给出了错误的解决方案。
from gekko import GEKKO
m = GEKKO(remote=False)
a = m.Var(value=0,integer=True)
b = m.Var(value=0,integer=True)
# constraints
m.Equation([a + b > 7,a**2 + b**2 < 40])
# Objective function
m.Maximize(a**3 + b**3)
m.solve(disp=False)
print(a.value[0])
print(b.value[0])
max_value = a.value[0]**3 + b.value[0]**3
print(max_value)
输出:
2.0
6.0
224.0
解决方法
在末尾添加检查表明 Gekko 在请求的容差范围内找到了正确的解,尽管默认求解器是 IPOPT,即使在请求 integer=True
时也能找到连续解。
a = a.value[0]; b = b.value[0]
print('a+b>7',a+b)
print('a^2+b^2<40',a**2+b**2)
# 0.71611780666
# 6.2838821838
# 248.50000026418317
# a+b>7 6.99999999046
# a^2+b^2<40 40.00000001289459
尝试使用 m.options.SOLVER=1
切换到 MINLP 求解器 APOT。不等式 <
和 <=
在 Gekko 中是等价的,因为它是一个数值解。
2.0
6.0
224.0
a+b>7 8.0
a^2+b^2<40 40.0
如果在数学意义上是 <
或 >
7 和 40 是不允许的,则在约束上向上或向下移动一个整数,例如 >=8
和 {{1 }}。
<=39
结果正确:
m.Equation([a + b >= 8,\
a**2 + b**2 <= 39])
这里还有什么遗漏吗?为什么声称没有可行的解决方案?
a=3.0 b=5.0 Objective: 152.0
a+b>=8 with a+b=8.0 (constraint satisfied)
a^2+b^2=<39 with a^2+b^2=34.0 (constraint satisfied)
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。