2

我想将值分配给 GAMS 中的 3-D 表。但它似乎不像在 Matlab 中那样工作......运气好吗?代码如下,问题出在最后几行:

    Sets
         n           nodes                       / Sto , Lon , Par , Ber , War , Mad , Rom /
         i           scenarios                   / 1 * 4 /
         k           capacity level              / L, N, H /  ;

alias(n,m);

Table balance(n,i) traffic balance for different nodes
         1       2       3       4
Sto      50      50      -50     -50
Lon      -40     40      -40     40
Par      0       0       0       0
Ber      0       0       0       0
War      40      -40     40      -40
Mad      0       0       0       0
Rom      -50     -50     50      50  ;

Scalar r  fluctuation rate of the capacity level
/0.15/;

Parameter p(k) probability of each level
/ L   0.25
  N   0.5
  H   0.25 / ;

Table nor_cap(n,m) Normal capacity level from n to m
          Sto    Lon     Par     Ber     War     Mad     Rom
Sto       0      11      14      25      30      0       0
Lon       11     0       21      0       0       14      0
Par       14     21      0       22      0       31      19
Ber       25     0       22      0       26      0       18
War       30     0       0       26      0       18      22
Mad       0      14      31      0       18      0       15
Rom       0      0       19      18      22      15      0  ;

Table  max_cap(n,m,k) capacity level under each k
max_cap(n,m,'N')=nor_cap(n,m)
max_cap(n,m,'L')=nor_cap(n,m)*(1-r)
max_cap(n,m,'H')=nor_cap(n,m)*(1+r);
4

1 回答 1

0

对 3-D 矩阵的最终分配应使用 PARAMETER 而不是 TABLE 完成。一般来说,我还会注意到 TABLE 非常严格(二维,代码内的文本输入)。您可能需要考虑使用 $GDXIN(或 EXECUTE_LOAD)和一些 GAMS 实用程序来加载 xls 或 csv 文件。

作为 MATLAB 和 GAMS 的用户,我会注意到 GAMS 依赖于每个数组的“索引”,但除此之外它们可能非常相似。在您的情况下, max_cap(n,m,k) 类似于每个容量级别方案下 from_city 和 to_city 之间的最大容量。您的矩阵需要声明为 PARAMETER,它可以是任何 n 维(索引)矩阵,甚至包括 SCALAR。

另外,如果您真的需要快速回答,请尝试 GAMS 邮件列表,全球精通 GAMS 用户的数量不能超过几千,因此在 StackOverflow 上可能很难找到快速答案 - 太棒了更常见的语言。

于 2015-04-05T20:28:44.190 回答