我确实使用Gams解决了mcp优化问题。 我确实设置了所有方程式,但无法设置初始值。显然,第三个方程是试图对资源提取时的物流增长进行建模,但是由于我没有初始存货,所以它不起作用。
Equations
Objective_R(t) Objective Function
Movement_Stock(t) stock movement
Movement_ShadowPrice(t) shadow price
;
Objective_R(t)..
l(t)+c =g= (a - R(t))/b
;
Movement_Stock(t)..
l(t) + z*l(t) =g= l(t+1)
;
Movement_ShadowPrice(t)..
S(t-1)*g + S(t-1) -R(t) =e= S(t)
;
t0(t) = yes$(ord(t) eq 1);
S.L(t) = S0$t0(t)