如何在Pyomo中建模流量平衡方程？

Question

我需要在我的Pyomo模型中指定以下方程式：

w[t]=w[t-1]+y[t]+z[t]-v[t]

t是時間步長， w ， y ， z和v是隨時間步長變化的變量。 我應該將其建模為約束嗎？ 我贊賞任何類似的方程式（也許不完全相同）建模的例子。

Answer 1

是的，您應該將此建模為約束。 這是一個幫助您入門的小示例：

m.t = Set(initialize=[1,2,3,4])
m.w = Var(m.t)
m.y = Var(m.t)

def _flow_rule(m, t):
   if t == 1:
      return Constraint.Skip
   return m.w[t] == m.w[t-1] + m.y[t]
m.flow = Constraint(m.t, rule=_flow_rule)