我的脚本是:
from __future__ import division
import numpy
import scipy
from pyomo.environ import *
from pyomo.dae import *
from pyomo.opt import SolverFactory
m=ConcreteModel()
m.x3=Var(within=NonNegativeReals)
m.u=Var(within=NonNegativeReals)
def _con(m):
return m.x3 >=3
m.con=Constraint(rule=_con)
def _con2(m):
return 4 >= m.u >=1
m.con2=Constraint(rule=_con2)
m.obj=Objective(expr=m.x3*m.u)
opt = SolverFactory("Ipopt", executable = "/Ipopt-3.12.6/bin/ipopt")
results = opt.solve(m)
results.write()
虽然这是一个非常简单的问题,尽管程序说明了这一点 它找到了最佳解决方案, 解决方案的数量为0,并且没有显示任何解决方案。
任何想法??
非常感谢。
答案 0 :(得分:5)
在Pyomo的最新版本中,默认情况下将解决方案加载到模型中。结果对象应该用于检查状态信息。如果您想询问解决方案,可以直接访问模型组件并检查其值,或者您可以执行以下操作:
m.solutions.store_to(results)
results.write()