Pyomo多重积分

Question

是否存在与Pyomo的求和函数的连续（积分）模拟，可以一次进行多次求和

sum(model.x[i,j] for i in model.index_is, for j in model.index_js)

dae示例中的Heat_Conduction.py有一个对嵌套集成的引用，但它被注释掉并抛出一个ValueError：每当我将它注释回来时，必须用连续S索引积分。

# def _obj(m):
#     return Integral(Integral(expr=m.intExp,wrt=m.x,bounds=(0,10)), wrt=m.y,bounds=(0,5))
# m.obj = Objective(rule=_obj)

我到处搜索，但我找不到一个有效的例子，我无法理解文档。

提前致谢

Answer 1

Integral组件尚未完全开发，我不建议使用它。使用Pyomo.DAE组件制定连续积分的更好方法是区分积分方程并将其转换为微分方程。请参阅此处的讨论：How to define an Integral as an objective function in pyomo?