我使用pyomo.dae解决微分方程优化问题我定义了一个集合
m.e = ContinuousSet(bounds=(e0, ef))
但是我希望ef作为变量自由。这是一个灵活的端点微分问题。我该怎么办?
ef
您更改了微分方程
y'(x) = f(x,y(x))
在弹性间隔[e0, ef]到超过标准间隔[0,1]的版本,通过
[e0, ef]
[0,1]
u'(s) = T*f(e0+T*s, u(s))
现在e0和T=ef-e0可以像其他任何参数一样对待。
e0
T=ef-e0