数学建模社区-数学中国

标题: 一个模型的编程。我自己是在是编不出来了。（初学者） [打印本页]

作者: `Ｍr.Z.. 时间: 2012-10-13 18:20
标题: 一个模型的编程。我自己是在是编不出来了。（初学者）
model:
SETS:
ZD/D1..D8/:Pi;
ZL/L1..L8/:Qj;
links1(ZD,ZL)|&1#eq#&2:x;
TD/TD1..TD8/:Ai;
TL/TL1..TL8/:Bj;
links2(TD,TL)|&1#eq#&2:y;
tx/xt1..xt8/:ri;
xs/xs1..xs8/:si;
ENDSETS
DATA:
Pi=60,55,51,43,41,52,50,51;
Qj=35,37,22,32,41,32,43,38;
Bj=38,37,23,38,44,33,45,38;
k=0.1;
ENDDATA
max=k*(@sum(links1(i,j):Pi*Qj)-@sum(links2(i,j):Ai*Bj));
@for(xs(i):si=1+((Bi-Qi)/Qi));
@for(tx(i):ri<1.1);
@for(tx(i):ri>0.9);
@for(TD(i):Ai=ri*Pi);
@for(TL(i):Bi=yi*Qi);
@sum(links1(i,j):Pi*Qj)=@sum(links2(i,j):Ai*Bj);
求不出来解。看看我哪写错了。谢谢。
上面的衍生集合。不会用啊。

作者: madio 时间: 2012-10-13 18:28
@sum(links1(i,j):Pi*Qj)=@sum(links2(i,j):Ai*Bj);
最后一句后面少了一个分号，已经错误提示了，很好找的。

作者: `Ｍr.Z.. 时间: 2012-10-13 19:20
好像没有少。是你少复制了吧

作者: `Ｍr.Z.. 时间: 2012-10-13 19:25

madio 发表于 2012-10-13 18:28
@sum(links1(i,j)i*Qj)=@sum(links2(i,j):Ai*Bj);
最后一句后面少了一个分号，已经错误提示了，很好找的 ...

你用这个程序，能算出数来？

作者: madio 时间: 2012-10-13 19:40

`Ｍr.Z.. 发表于 2012-10-13 19:25
你用这个程序，能算出数来？

能找到局部最优解，说明程序没有问题！

作者: `Ｍr.Z.. 时间: 2012-10-13 20:02
本帖最后由 `Ｍr.Z.. 于 2012-10-13 20:03 编辑

madio 发表于 2012-10-13 19:40
能找到局部最优解，说明程序没有问题！

我这出来是这个

作者: madio 时间: 2012-10-13 20:03

`Ｍr.Z.. 发表于 2012-10-13 20:02
我这出来是这个。

看不到图片，你得将图片当附件上传后再点击放在贴子中，要不看不见！

作者: `Ｍr.Z.. 时间: 2012-10-13 20:04

madio 发表于 2012-10-13 20:03
看不到图片，你得将图片当附件上传后再点击放在贴子中，要不看不见！

好了。帮我看一下，衍生集合，定义后，在程序里应用时正确吗？

作者: Vir 时间: 2012-10-13 21:04
感觉应该配套原有公式，否则不知道你的lingo到底是否对应上正确的公式没

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)