数学建模社区-数学中国

标题: 请教一段lingo程序，多谢各位大侠帮忙啦~~帮顶也好~~~ [打印本页]

---

作者: yansichong    时间: 2010-8-5 11:27
标题: 请教一段lingo程序，多谢各位大侠帮忙啦~~帮顶也好~~~
sets:
supply/1..7/:p,s;
need/1..15/:y,t,l;
link/supply,need/:x,c;
endsets
data:
s=800,800,1000,2000,2000,2000,3000;
p=160,155,155,160,155,150,160;
c=264.3 160.3 140.2 98.6 38.0 20.5 3.1 21.2 64.2 92.0 96.0 106.0 121.2 128.0 142.0
  309.3 205.3 190.2 171.6 111.0 95.5 86.0 71.2 114.2 142.0 146.0 156.0 171.2 178.0 192.0
  324.3 220.3 200.2 181.6 121.0 105.5 96.0 86.2 48.2 82.0 86.0 96.0 111.2 118.0 132.0
  354.3 250.3 235.2 216.6 156.0 140.5 131.0 116.2 84.2 62.0 51.0 61.0 76.2 83.0 97.0
  349.3 245.3 225.2 206.6 146.0 130.5 121.0 111.2 79.2 57.0 33.0 51.0 71.2 73.0 87.0
  359.3 255.3 235.2 216.6 156.0 140.5 131.0 121.2 84.2 62.0 51.0 45.0 26.2 11.0 28.0
  369.3 265.3 245.2 226.6 166.0 150.5 141.0 131.2 99.2 77.0 66.0 56.0 38.2 26.0 2.0;
enddata
init:
endinit
min=@sum(links(i,j):(p(i)+c(i,j))*x(i,j))+0.05*@sum(need(i):y*(y+1)+(l-y)*(l-y-1))
@for(supply(i):@sum(need(j):x(i,j))>=500*t(i));
     @for(supply(i):@sum(need(j):x(i,j))<=s(i)*t(i));
     @for(supply(i):@bin(t(i)));
@for(need(j)|j#NE#1 #AND# j#NE#15:@sum(supply(i):x(i,j))=y(j)+(l(j-1)-y(j-1)));
@sum(link(i,j)|j#EQ#15:x(i,j))=l(14)-y(14);
y(1)=0;
@for(need(i):y(i)<=l(i));
@for(need(i):y(i)>=0);
@for(link(i,j):x(i,j)>=0);
@sum(link(i,j):x(i,j))=5171;
end
总是对c 数据的第一行的140.2报错……说是不可用的没有初始化的数据…………

---

作者: yansichong    时间: 2010-8-5 11:28
自己顶个~~~~~~~~~

---

作者: yansichong    时间: 2010-8-5 11:28
呃……怎么有这么多表情………………

---

作者: linmatsas    时间: 2010-8-5 14:55
对呀……下次直接顺手传个附件吧……要不这么多表情怎么看呀……而且而且……能不能把方程一起发上来~~~~~

---

作者: lixiaoyong    时间: 2010-8-5 19:32
不好意思，学得不精，不能帮上你

---

作者: tiancai1234    时间: 2010-8-10 18:40
将问题也一起发出来吧

---

作者: tiancai1234    时间: 2010-8-11 12:53
能将题目，问题一并拿出来吗

---

作者: wuguishouling    时间: 2010-8-11 16:01
你把lg4文件发上来

---

作者: wuguishouling    时间: 2010-8-11 17:18
晕···link(supply,need)!

---

作者: yansichong    时间: 2010-8-16 16:00
这两个一个是论文，另一个是根据论文的第一个模型编的程序~~编译通过了，但是结果不对……不知道应该怎么改了……

钢管订购和运输问题优化——安冉.pdf

2010-8-16 15:53 上传

点击文件名下载附件

下载积分: 体力 -2 点

143.07 KB, 下载次数: 0, 下载积分: 体力 -2 点

LINGO4.lg4

2010-8-16 15:58 上传

点击文件名下载附件

下载积分: 体力 -2 点

8 KB, 下载次数: 0, 下载积分: 体力 -2 点

---

作者: yansichong    时间: 2010-8-16 16:21
求解释……

---

作者: zhao6884269    时间: 2010-8-26 11:33
试试运气啦~~~~~~~~~~~

---

作者: lilijie_shangha    时间: 2010-8-26 11:34
来报道!!!!!!!!!!!

---

作者: asdffdsa    时间: 2010-8-26 11:35
我要把这个帖子一直往上顶，往上顶！

---

作者: mintw    时间: 2010-8-26 11:36
强烈支持。楼主万岁

---

作者: z123456789    时间: 2010-8-26 11:38
强烈支持。楼主万岁

---

作者: 547896321    时间: 2010-8-26 11:46
我来了~~~~~~~~~ 闪人~~~~~~~~~~~~~~~~  

---

作者: mybaby    时间: 2010-8-26 13:52
哦～～

---

作者: 4tm    时间: 2010-8-26 18:59
留个脚印```````

---

作者: zhyg    时间: 2010-8-26 19:17
强烈支持。楼主万岁

---

作者: shbjava    时间: 2010-8-27 00:00
鉴定完毕！  

---

作者: sigh    时间: 2010-8-27 12:00
强烈支持。楼主万岁

---

作者: arthuryuan    时间: 2010-8-27 15:00
楼主，你写得实在是太好了。我惟一能做的，就只有把这个帖子顶上去这件事了

---

作者: abottle    时间: 2010-8-27 20:00
我来了~~~~~~~~~ 闪人~~~~~~~~~~~~~~~~  

---

作者: move170    时间: 2010-8-28 12:00
鉴定完毕！  

---

作者: gctygl    时间: 2010-8-28 15:00
我来了~~~~~~~~~ 闪人~~~~~~~~~~~~~~~~  

---

作者: wushenlijessie    时间: 2010-8-28 20:00
强烈支持。楼主万岁

---

作者: zhaojie_1999    时间: 2010-8-28 23:59
哦～～

---

作者: saimons    时间: 2010-8-29 08:00
来报道!!!!!!!!!!!

---

作者: hcq0_0    时间: 2010-8-29 12:00
试试运气啦~~~~~~~~~~~

---

作者: smsmsm    时间: 2010-8-29 15:00
楼主的帖子实在是写得太好了。可是我立刻想到，这么好的帖子，倘若别人看不到，那么不是浪费楼主的心血吗？经过痛苦的思想斗争，我终于下定决心，牺牲小我，奉献大我。我要拿出这帖子奉献给世人赏阅，我要把这个帖子一直往上顶，往上顶！顶到所有人都看到为止！  

---

作者: luohan    时间: 2010-8-29 20:00
我来了~~~~~~~~~ 闪人~~~~~~~~~~~~~~~~  

---

作者: why8255    时间: 2010-8-30 12:00
声明一下：本人看贴和回贴的规则，好贴必看，精华贴必回。

---

作者: 天涯客    时间: 2010-8-30 15:00
试试运气啦~~~~~~~~~~~

---

作者: surfw    时间: 2010-8-30 20:00
鉴定完毕！  

---

作者: wz    时间: 2010-8-31 08:00
我回不抢呢 考虑再三 还是不抢了吧 ^_^

---

作者: hayate    时间: 2010-8-31 12:00
鉴定完毕！  

---

作者: myflying    时间: 2010-8-31 15:00
顶顶更健康，越顶吃的越香。

---

作者: dargal    时间: 2010-8-31 20:02
强烈支持。楼主万岁

---

作者: fbxyxy    时间: 2010-9-1 12:00
呵呵 大家好奇嘛 来观看下~~~~  

---

作者: zhaoxilin    时间: 2010-9-1 12:00
不错不错,我喜欢看  

---

作者: equan4321    时间: 2010-9-1 15:00
强烈支持。楼主万岁

---

作者: ywycmm    时间: 2010-9-1 20:00
楼主，你写得实在是太好了。我惟一能做的，就只有把这个帖子顶上去这件事了

---

作者: UserQi    时间: 2010-9-2 12:00
哦～～

---

作者: cugbhedayi    时间: 2010-9-2 15:00
鉴定完毕！  

---

作者: julpee    时间: 2010-9-2 20:00
顶顶更健康，越顶吃的越香。

---

作者: huzh369    时间: 2010-9-3 12:00
试试运气啦~~~~~~~~~~~

---

作者: viviny    时间: 2010-9-3 15:00
强人，佩服死了。呵呵，不错啊

---

作者: zhengyx    时间: 2010-9-3 20:00
哦～～

---

作者: alair002    时间: 2012-1-13 20:46
没有体力啦，资料能发给我一份吗？我的邮箱是18633525948圈163邮箱，谢啦

---

作者: hyltt7636    时间: 2012-2-1 00:34
试试这个吧
sets:
supply/1..7/:p,s;
need/1..15/:y,t,l;
link(supply,need):x,c;
endsets
data:
s=800,800,1000,2000,2000,2000,3000;
p=160,155,155,160,155,150,160;
c=264.3 160.3 140.2 98.6 38.0 20.5 3.1 21.2 64.2 92.0 96.0 106.0 121.2 128.0 142.0
  309.3 205.3 190.2 171.6 111.0 95.5 86.0 71.2 114.2 142.0 146.0 156.0 171.2 178.0 192.0
  324.3 220.3 200.2 181.6 121.0 105.5 96.0 86.2 48.2 82.0 86.0 96.0 111.2 118.0 132.0
  354.3 250.3 235.2 216.6 156.0 140.5 131.0 116.2 84.2 62.0 51.0 61.0 76.2 83.0 97.0
  349.3 245.3 225.2 206.6 146.0 130.5 121.0 111.2 79.2 57.0 33.0 51.0 71.2 73.0 87.0
  359.3 255.3 235.2 216.6 156.0 140.5 131.0 121.2 84.2 62.0 51.0 45.0 26.2 11.0 28.0
  369.3 265.3 245.2 226.6 166.0 150.5 141.0 131.2 99.2 77.0 66.0 56.0 38.2 26.0 2.0;
enddata
init:
endinit
min=@sum(link(i,j)p(i)+c(i,j))*x(i,j))+0.05*@sum(need(i):y*(y+1)+(l-y)*(l-y-1));
@for(supply(i)sum(need(j):x(i,j))>=500*t(i));
      @for(supply(i)sum(need(j):x(i,j))<=s(i)*t(i));
      @for(supply(i)bin(t(i)));
@for(need(j)|j#NE#1 #AND# j#NE#15sum(supply(i):x(i,j))=y(j)+(l(j-1)-y(j-1)));
@sum(link(i,j)|j#EQ#15:x(i,j))=l(14)-y(14);
y(1)=0;
@for(need(i):y(i)<=l(i));
@for(need(i):y(i)>=0);
@for(link(i,j):x(i,j)>=0);
@sum(link(i,j):x(i,j))=5171;
end

---

作者: hyltt7636    时间: 2012-2-1 00:42
试试这个
sets:
supply/1..7/:p,s;
need/1..15/:y,t,l;
link(supply,need):x,c;
endsets
data:
s=800,800,1000,2000,2000,2000,3000;
p=160,155,155,160,155,150,160;
c=264.3 160.3 140.2 98.6 38.0 20.5 3.1 21.2 64.2 92.0 96.0 106.0 121.2 128.0 142.0
  309.3 205.3 190.2 171.6 111.0 95.5 86.0 71.2 114.2 142.0 146.0 156.0 171.2 178.0 192.0
  324.3 220.3 200.2 181.6 121.0 105.5 96.0 86.2 48.2 82.0 86.0 96.0 111.2 118.0 132.0
  354.3 250.3 235.2 216.6 156.0 140.5 131.0 116.2 84.2 62.0 51.0 61.0 76.2 83.0 97.0
  349.3 245.3 225.2 206.6 146.0 130.5 121.0 111.2 79.2 57.0 33.0 51.0 71.2 73.0 87.0
  359.3 255.3 235.2 216.6 156.0 140.5 131.0 121.2 84.2 62.0 51.0 45.0 26.2 11.0 28.0
  369.3 265.3 245.2 226.6 166.0 150.5 141.0 131.2 99.2 77.0 66.0 56.0 38.2 26.0 2.0;
enddata
init:
endinit
min=@sum(links(i,j)p(i)+c(i,j))*x(i,j))+0.05*@sum(need(i):y*(y+1)+(l-y)*(l-y-1))
@for(supply(i)sum(need(j):x(i,j))>=500*t(i));
      @for(supply(i)sum(need(j):x(i,j))<=s(i)*t(i));
      @for(supply(i)bin(t(i)));
@for(need(j)|j#NE#1 #AND# j#NE#15sum(supply(i):x(i,j))=y(j)+(l(j-1)-y(j-1)));
@sum(link(i,j)|j#EQ#15:x(i,j))=l(14)-y(14);
y(1)=0;
@for(need(i):y(i)<=l(i));
@for(need(i):y(i)>=0);
@for(link(i,j):x(i,j)>=0);
@sum(link(i,j):x(i,j))=5171;
end

---

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

Powered by Discuz! X2.5