查看: 2140|回复: 2

求高人相助LINGO 编程

字体大小: 正常放大

3 主题	4 听众	18 积分

升级 13.68%

TA的每日心情

	开心 2012-2-8 20:27

签到天数: 3 天

[LV.2]偶尔看看I

电梯直达

1^#

发表于 2012-8-2 14:49 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

5体力

model:
sets:
product /1..7/:value，capacity;
need /1..14/:demand;
link(product,need):cost,number;
endsets
data:
value=160 155 155 160 155 150 160;
cost= 85.3 85.2 140 38 20.5 3.1 21.2 64.2 92 96 96 116.2 123 132
135.3 135.2 185 111 95.5 86 71.2 114.2 142 146 146 166.2 173 182
145.3 145.2 200 121 105.5 96 86.2 48.2 82 86 86 106.2 108 122
180.3 180.2 230 156 140.5 131 116.2 84.2 62 51 51 71.2 73 87
175.3 120.2 225 141 125.5 116 106.2 74.2 51 33 78 61.2 63 77
175.3 170.2 225 146 130.5 121 111.2 79.2 57 45 45 26.2 11 28
190.3 185.2 240 151 135.5 136 126.5 94.2 67 61 45 38.2 26 2;
capacity=800 800 1000 2000 2000 2000 3000;
demand=500 920 1130 1350 1650 2130 2810 3011 3216 3410 4016 4766 5067 5171;
enddata
!目标函数;
min=@sum(link(i,j):cost(i,j)*number(i,j))+@sum(product(i):value(i)*@sum(need(j):number(i,j)));
!number各行和要求
@for(product(i)|i#eq#1:@sum(need(j):number(i,j))>=500;@sum(need(j):number(i,j))<=capacity(i));
@for(product(i)|i#eq#2:@sum(need(j):number(i,j))>=500;@sum(need(j):number(i,j))<=capacity(i));
@for(product(i)|i#eq#3:@sum(need(j):number(i,j))>=500;@sum(need(j):number(i,j))<=capacity(i));
@for(product(i)|i#eq#4:@sum(need(j):number(i,j))>=500;@sum(need(j):number(i,j))<=capacity(i));
@for(product(i)|i#eq#5:@sum(need(j):number(i,j))>=500;@sum(need(j):number(i,j))<=capacity(i));
@for(product(i)|i#eq#6:@sum(need(j):number(i,j))>=500;@sum(need(j):number(i,j))<=capacity(i));
@for(product(i)|i#eq#7:@sum(need(j):number(i,j))>=500;@sum(need(j):number(i,j))<=capacity(i));
!number各列和要求
@for(need(j)|j#eq#14:@sum(product(i):number(i,j))>=0; @sum(product(i):number(i,j))<=500);
@for(need(j)|j#ge#13 #and# j#le#14:@sum(product(i):number(i,j))>=500; @sum(product(i):number(i,j))<=920);
@for(need(j)|j#ge#12 #and# j#le#14:@sum(product(i):number(i,j))>=920; @sum(product(i):number(i,j))<=1130);
@for(need(j)|j#ge#11 #and# j#le#14:@sum(product(i):number(i,j))>=1130; @sum(product(i):number(i,j))<=1350);
@for(need(j)|j#ge#10 #and# j#le#14:@sum(product(i):number(i,j))>=1350; @sum(product(i):number(i,j))<=1650);
@for(need(j)|j#ge#9 #and# j#le#14:@sum(product(i):number(i,j))>=1650; @sum(product(i):number(i,j))<=2130);
@for(need(j)|j#ge#8 #and# j#le#14:@sum(product(i):number(i,j))>=2130; @sum(product(i):number(i,j))<=2810);
@for(need(j)|j#ge#7 #and# j#le#14:@sum(product(i):number(i,j))>=2810; @sum(product(i):number(i,j))<=3011);
@for(need(j)|j#ge#6 #and# j#le#14:@sum(product(i):number(i,j))>=3011; @sum(product(i):number(i,j))<=3216);
@for(need(j)|j#ge#5 #and# j#le#14:@sum(product(i):number(i,j))>=3216; @sum(product(i):number(i,j))<=3410);
@for(need(j)|j#ge#4 #and# j#le#14:@sum(product(i):number(i,j))>=3410; @sum(product(i):number(i,j))<=4016);
@for(need(j)|j#ge#3 #and# j#le#14:@sum(product(i):number(i,j))>=4016; @sum(product(i):number(i,j))<=4766);
@for(need(j)|j#ge#2 #and# j#le#14:@sum(product(i):number(i,j))>=4766; @sum(product(i):number(i,j))<=5067);
@for(need(j)|j#ge#1 #and# j#le#14:@sum(product(i):number(i,j))>=5067; @sum(product(i):number(i,j))<=5171);
!number中所有元素和要求
@for(link(i,j):@sum(link(i,j):number(i,j))=5171);
!number中的整数线性规划
@for(link (i,j): @gin(number (i,j)));
end

复制代码

程序结果是找不到可行解，是编程出错了吧~~~求高人指教，高人啊高人
是不是这里三句有错：（1）@for(link(i,j):@sum(link(i,j):number(i,j))=5171);本意是求所有元素和
（2）@sum(product(i):value(i)*@sum(need(j):number(i,j))；本意是求（value(i)*（（i）行的和））的（i）从1到7的和
（3）写各列和的意思是（第3列列和+第4列列和+……第14列列和）的总和大于小于某值

zan

160, 155, product, number, demand