我目前正在努力解决这个问题。我需要最大化这家公司的利润。
这就是我目前拥有的代码:
from pyomo.environ import *
from pyomo.opt import *
opt = solvers.SolverFactory("ipopt")
model = ConcreteModel()
model.x1 = Var(within=NonNegativeIntegers)
model.x2 = Var(within=NonNegativeIntegers)
model.y1 = Var(within=NonNegativeIntegers)
model.y2 = Var(within=NonNegativeIntegers)
model.b1 = Var(within=Boolean)
model.b2 = Var(within=Boolean)
model.c1 = Constraint(expr = model.x1 + model.x2 + model.y1 + model.y2 <= 7000)
model.c2 = Constraint(expr = 2*model.x1 + 2*model.x2 + model.y1 + model.y2 <= 10000)
model.c3 = Constraint(expr = model.x1 <= 2000)
model.c4 = Constraint(expr = model.x2 <= 1000)
model.c5 = Constraint(expr = model.y1 <= 2000)
model.c6 = Constraint(expr = model.y2 <= 3000)
model.z = Objective(expr= (150*model.x1 + 180*model.x2*model.b1 + 100*model.y1 + 110*model.y2*model.b2), sense=maximize)
results = opt.solve(model)
这是我尝试为约束编写的代码,只要它不超过 2000 个产品,它就只使用第一个斜率:
def ObjRule(model):
if model.x1 >= 2000:
return model.b1==1
if model.x2 >= 2000:
return model.b2 == 1`
如果有人有提示,我该如何继续那就太好了。
先感谢您,
帕特里克
在 Pyomo 中,规则不是发送到求解器的回调。对每个索引调用一次它们以获得一组静态表达式。这组表达式被发送到求解器。您在规则内部使用的任何 if 逻辑不应涉及变量的值(除非它基于变量的初始值,在这种情况下,无论您在主表达式之外使用变量,您都可以将变量包装在 value() 中返回)。
如果您想对分段函数进行建模,则需要应用某种建模技巧来实现。在某些情况下,这涉及引入离散变量(参见examples https://github.com/Pyomo/pyomo/tree/master/examples/pyomo/piecewise对于分段分量),在其他情况下则不然(例如,当最大化可表示为有限数量仿射函数的最小值的分段函数时)。
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)