主页 / coding / 问题 / 79285727
Accepted
garej
Asked: 2024-12-17 02:29:24 +0800 CST

如何正确循环具有方程动态的 GAMS 模型？

  • 772

我正在尝试了解在 GAMS 中创建递归动态模型的正确方法。为此，我制作了一个由 2 个线性方程组成的玩具示例（又称蛛网模型）。现在它从随机点开始P(t=1)=40并移动到平衡点（系统解）。

t在循环内运行带有索引的简单代码似乎有些棘手。GAMS 总是抱怨 (!)：

错误 50：循环控制指数出现在模型方程中

我一直在查看其他代码片段（例如），但找不到完整的示例来重现。玩弄代码时，我刚刚将另一个z和 循环放入 中z。我工作正常。

到目前为止一切顺利，但请注意 howt仍然在循环内。我有点困惑。（至少我不知道这种编码的副作用。）

Set t /t1*t10/
 * HERE IS THE TRICK
    z /z1*z10/;

Parameters a /100/, b /2/, c /20/, d /1/;

Variables P(t), Qd(t), Qs(t), obj;

Positive Variables P, Qd, Qs;

Equations Demand(t), Supply(t), Equilibrium(t), Objective; 
    Demand(t).. Qd(t) =E= a - b * P(t);
    Supply(t).. Qs(t) =E= c + d * P(t-1);
    Equilibrium(t).. Qd(t) =E= Qs(t);
    Objective.. obj =E= 1; 

P.l('t1') = 40;

Model cobwebModel /all/;

Loop(z,
    Solve cobwebModel using LP maximizing obj;
    P.l(t+1) = (a - c + d * P.l(t)) / (b + d);
);

问题是：如何重构上述代码以使其符合适当的 GAMS 风格并准备扩展？

loops

1 个回答

  1. Best Answer

    我不太明白你最后想解决什么问题，但如果我理解正确的话，你t一次只想解决 中的一个元素？在你的代码中，这不会发生，相反，你会对 中的所有 10 个元素求解 10 次t。但你可以使用动态子集，以不同的方式执行此操作，如下所示：

    Set t /t1*t10/
    tt(t);
    
Alias (t,ta);

Parameters a /100/, b /2/, c /20/, d /1/;

Variables P(t), Qd(t), Qs(t), obj;

Positive Variables P, Qd, Qs;

Equations Demand(t), Supply(t), Equilibrium(t), Objective; 
    Demand(tt(t)).. Qd(t) =E= a - b * P(t);
    Supply(tt(t)).. Qs(t) =E= c + d * P(t-1);
    Equilibrium(tt(t)).. Qd(t) =E= Qs(t);
    Objective.. obj =E= 1; 

P.l('t1') = 40;

Model cobwebModel /all/;

Loop(ta,
    tt(ta)=yes;
    Solve cobwebModel using LP maximizing obj;
    P.l(ta+1) = (a - c + d * P.l(ta)) / (b + d);
    tt(ta)=no;
);
    • 1

相关问题