将矩阵与 GAMS 中的向量相乘

问题描述

在 GAMS 中考虑下表：

Set i /i1*i4 /;
Table a(i,j) 'original matrix'
        i1         i2     i3           i4
   i1    2          0      0            0
   i2    0.272727   2      0.727273     0
   i3    0          0.8    2            0.2
   i4    0          0      0            2;

我想将 a 的倒矩阵与以下向量相乘：

parameters v(i) /
i1 0
i2 -0.03413
i3 -0.01174
i4 0 /;

所以我最终会得到以下向量：

parameters B(i) /
i1 0
i2 -0.01747
i3 0.001116
i4 0

我尝试了以下方法：

Alias (i,j);
Parameter inva(i,j) 'inverse of a';

execute_unload 'a.gdx',i,a;
execute '=invert.exe a.gdx i a b.gdx inva';
execute_load 'b.gdx',inva;

parameter B(i);

loop(j,loop(i,B(i) = inva(i,j)*v(i)
    )
);

我的逆矩阵看起来不错，但是矩阵/向量乘法出了点问题。

解决方法

只是简单地查看了您的代码，在我看来，您似乎希望将最后一个循环设为：

loop(j,B(i) = inva(i,j)*v(j)
     );

因为您只想对矩阵的列求和。