GAMS，矩阵直接赋值答案

【问题标题】：GAMS, matrix direct assignmentGAMS，矩阵直接赋值
【发布时间】：2011-09-21 22:49:55
【问题描述】：

我想为 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);

【问题讨论】：

标签： matrix variable-assignment gams-math

【解决方案1】：

应使用 PARAMETER 而不是 TABLE 对 3-D 矩阵进行最终分配。一般来说，我还会注意到 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 上可能很难找到快速答案 - 真棒它适用于更常见的语言。

【讨论】：