Блог

Я хочу иметь возможность определять интеграл pyomo как часть целевой функции.

Я не могу понять, какое выражение необходимо для интеграла.
Вот мое лучшее предположение:

 model = ConcreteModel()
model.t = ContinuousSet(bounds = (0,1))
model.y = Var(model.t)
model.dydt = DerivativeVar(model.y, wrt=(model.t))

def myintegral(model,i):
    return model.dydt[i]

model.n = Integral(model.t, wrt=model.t, rule=myintegral)  # this line is the trouble

def myobjective(model):
     return model.n

model.obj = Objective(rule=myobjective)
  

Ошибка: TypeError: A callable type that is not a Pyomo expression can not be used to initialize an Expression object. Use 'rule' to initalize with function types.

Но я не понимаю, почему выражение внутри интеграла является проблемой, поскольку эти переменные, по-видимому, полностью индексируются индексом model.t :

 # but, this is totally fine:  
print model.dydt[0]
print model.dydt[1]
  

Я что-то недопонимаю в этом?

Вот некоторые ресурсы, с которыми я консультировался до сих пор:

https://groups.google.com/forum /#!тема/pyomo-форум/6RhEXEMDTPc
https://software.sandia.gov/downloads/pub/pyomo/PyomoOnlineDocs.html#_parameters
https://projects.coin-or.org/Coopr/browser/pyomo/trunk/examples/dae/Heat_Conduction.py?rev=9315

Я открыт для предложений / ссылок о других ресурсах pyomo .

Гейб прав, это действительно ошибка в классе Integral, и она была исправлена в репозитории github. Еще одна ошибка в вашем примере модели — это спецификация целевого компонента. Вы должны использовать ключевое слово ‘rule’ вместо ‘expr’

 def myobjective(model):
     return model.n
model.obj = Objective(rule=myobjective)
  

Кроме того, я хочу повторить кое-что, упомянутое в онлайн-документации для pyomo.dae. Компонент Integral является прототипом и не полностью разработан. Я не рекомендую использовать его для сложных интегралов или моделей, требующих решений с высокой точностью. Класс Integral использует правило трапеции для численного интегрирования. Я бы рекомендовал преобразовать любые интегралы в вашей задаче в дифференциальные уравнения и решить их, используя предоставленные автоматические преобразования дискретизации.

Это похоже на ошибку. Вы должны открыть заявку здесь: https://github.com/Pyomo/pyomo/issues