用GAMS如何编写如下约束条件，在线等。

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

p(j)+q(i,j)=l=0.05*x(i)*a(i,j)这样写为什么会报错？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：gams 约束条件 GAM AMS 在线等 在线 如何

你定义的约束方程是不是忘了写维度？

EQUATIONS
        CONS(i,j)
;
CONS(i,j)..    z(i) + p(i,j) =G= a(i,j) * x(j);

qianhaoqi 发表于 2014-9-8 23:03
你定义的约束方程是不是忘了写维度？

EQUATIONS

嗯，非常感谢啊，但还是会出现domain violation for set的报错，还有，当j=1的时候需要单独列出，如何写？

qianhaoqi 发表于 2014-9-8 23:03
你定义的约束方程是不是忘了写维度？

EQUATIONS

face(1)..z- sum(m,u(m)*do1(m))+2*p(1)+sum(i,q(i,1)) =l= 0
或 scalar b /1/;
face(b)..z- sum(m,u(m)*do1(m))+2*p(b)+sum(i,q(i,b)) =l= 0都会出现domain violation for set，如何改写？

【炊烟】 发表于 2014-9-9 10:26
嗯，非常感谢啊，但还是会出现domain violation for set的报错，还有，当j=1的时候需要单独列出，如何写？ ...

需要单独列出的，然后定义一个子集的映射，用自己去conditional控制equations
举例来说：

SETS
        i       / 1*2 /
        j       / 1*5 /
        ij(i,j) / 1*2.1 /
;

EQUATIONS
        CONS1(i,j)
        CONS2(i,j)
;
CONS1(i,j)$ij(i,j)..    z(i) + p(i,j) =G= a(i,j) * x(j);
CONS2(i,j)$(not ij(i,j))..    z(i) + p(i,j) =G= a(i,j) * x(j);

上面定义的ij(i,j)是一个映射集合，这个集合定义的是（i=1,j=1）和（i=2,j=1），所以对应的CONS1方程就仅仅包含j=1的所有情况，而CONS2方程对应的是剩下的所有情况~

【炊烟】 发表于 2014-9-9 10:34
face(1)..z- sum(m,u(m)*do1(m))+2*p(1)+sum(i,q(i,1)) =l= 0
或 scalar b /1/;
face(b)..z- sum(m,u ...

conditional equation要用conditional的写法，GAMS里用的是dollar符号$，见上楼的例子可参考~或者查看GAMS的users guide里也有说明~

qianhaoqi 发表于 2014-9-10 13:42
conditional equation要用conditional的写法，GAMS里用的是dollar符号$，见上楼的例子可参考~或者查看GAM ...

请问$符号有什么用？

例如等式里 z=e=sum((i,j) $sam(i,j),(Q(i,j)-sam(i,j))**2);

bkjg 发表于 2015-6-1 20:02
请问$符号有什么用？

例如等式里 z=e=sum((i,j) $sam(i,j),(Q(i,j)-sam(i,j))**2);

$符号是表示行内条件判断，$sam(i,j)这个条件就是说只有sam这个二维变量中，第i行第j列的数据不为0时，才会对(i,j)这组维度进行加总，否则就跳过这一“sam单元格”。