我想模拟电池的放电。 我想在该时刻（t）上存储m.energy_remain [t]上剩余的能量。直到现在，我只能计算最终时间中剩余的能量，因为integrate（）函数会计算所有时间水平（m.t）。 我该如何改善它使其每次也具有“ i”？
def _energy(M,t):
return m.power[1,t]+m.power[2,t]
m.energy = Integral(m.t, wrt=m.t,rule=_energy)
def _energy_consumed(M,t):
return M.energy_remain[t] - 40000 + M.energy ==0
m.energy_remain_cons = Constraint(m.t,rule=_energy_consumed)