我有一个GAMS模型,我有
Set i / 1*6 /;
Alias (i,ip,il) ;
Variables
x(i,ip) ;
我想生成运算于x中所有向量的标量积的方程,不包括相同向量的乘积。类似的东西:
scalarProduct(i)..
sum(ip,x(i,ip)x(i,il)) =e= someConstant;
然而,这并不排除相同载体的产物。如何添加?我能以某种方式用美元声明吗?无论如何,这个陈述中可能有一些错误,我没有尝试过,因为我认为我想要的排除缺失。
答案 0 :(得分:1)
所以我想做的是:
Sets
i / 1*13 /
ii(i,i) diagonal elements / #i:#i /
ij(i,i) all elements / #i.#i /
ij_wo_ii(i,i);
获取没有对角线元素的所有组合:
ij_wo_ii(i,j) = ij(i,j) - ii(i,j);
然后我在这样的等式中使用它:
equation(j,k)..
sum(i,x(i,j)*x(i,k)$ij_wo_ii(j,k)) =l= 1;
这类似于正交性,除了矩阵中向量的乘积必须小于某个值,而不一定是0.不知道是否有一个术语。希望它对其他人也有用。