数模论坛

 找回密码
 注-册-帐-号
搜索
热搜: 活动 交友 discuz
查看: 2707|回复: 4

LINGO的问题,大家帮我看看,谢谢

[复制链接]
发表于 2004-9-22 23:38:22 | 显示全部楼层 |阅读模式
有一个新建油库,需要购买油罐来储油备用,总共有12种大小的油罐可供选择,油罐总储油量应该大于该油库三天的供应量Qd,而且规定油库油罐的种类不能超过两种,而且油罐总数不能少于两个,建立模型,使油库购买油罐的费用最小。
已知,油库每天的供应量为Qd=14863公斤/天;各种油罐的容积vpst(立方米)、价格mpst(千元);储存天数T=3;油的密度ry=534.89公斤/立方米;油罐的充装系数为0.9;xpst(I)为0-1变量,表示是否选用第i种大小的油罐;Npst(I)为整数,表示油库第i种油罐的数量。我建立了下面的模型,但是求出来的结果明显不正确,大家帮我看看问题出在哪里?
 楼主| 发表于 2004-9-22 23:38:44 | 显示全部楼层
model:
sets:
LPGST/1..12/:mpst,vpst,xpst,Apst,Npst,Npst_1;
endsets
data:
vpst=5,10,20,25,32,50,60,80,100,150,200,400;
mpst=45,58,68,75,85,110,120,150,175,250,380,700;
enddata
data:
Qd=14863;ry=534.89;
T=3;
enddata
min=(0.10+(1+0.30)/@fpa(0.15,20))*@sum(LPGST:mpst*xpst*Npst);
<a href="mailtry*@sum(LPGST:vpst*0.90*xpst*Npst) target="_blank" >=T*Qd">ry*@sum(LPGST:vpst*0.90*xpst*Npst)&gt;=T*Qd</A>;
<a href="mailtry*@sum(LPGST:vpst*0.90*xpst*Npst)<=T*Qd+ry*@min(LPGST:vpst*0.90*Apst" target="_blank" >ry*@sum(LPGST:vpst*0.90*xpst*Npst)&lt;=T*Qd+ry*@min(LPGST:vpst*0.90*Apst</A>);
@sum(LPGST:xpst)&gt;=1;
@sum(LPGST:xpst)&lt;=2;
@for(LPGST(I):Apst(I)=@if(xpst(I)#eq#1,1,1000000));
@for(LPGST(I):Npst(I)=@if(xpst(I)#eq#0,0,1+Npst_1(I)));
@sum(LPGST:xpst*Npst)&gt;=2;
@for(LPGST:Npst_1&gt;=0);
@for(LPGSTgin(Npst));
@for(LPGSTgin(Npst_1));
@for(LPGSTbin(xpst));
end
注:<a href="mailtry*@sum(LPGST:vpst*0.90*xpst*Npst)<=T*Qd+ry*@min(LPGST:vpst*0.90*Apst" target="_blank" >ry*@sum(LPGST:vpst*0.90*xpst*Npst)&lt;=T*Qd+ry*@min(LPGST:vpst*0.90*Apst</A>)这是我加的储油上限约束,<a href="mailtry*@min(LPGST:vpst*0.90*Apst" target="_blank" >ry*@min(LPGST:vpst*0.90*Apst</A>)这是所选储罐中单罐容积最小的储罐全储满油的油量,这个上限约束可能有点问题,但是不给上限约束,求出来的结果更加离谱。
@for(LPGST(I):Npst(I)=@if(xpst(I)#eq#0,0,1+Npst_1(I)))表示如果xpst=0,则Npst=0,否则,Npst&gt;=1。
还有一个疑问,用Lingo8.0求解时,我调换约束的前后位置后,以前能得出局部最优解的,竟然得不出结果了,为什么?
发表于 2004-9-23 15:27:06 | 显示全部楼层
请首先把语法问题解决吧
 楼主| 发表于 2004-9-23 17:12:12 | 显示全部楼层
不太明白,楼上能否讲清楚一点?
发表于 2004-9-24 03:09:12 | 显示全部楼层
<>你把网页上得程序拷入lingo,运行一次,看看</P>
您需要登录后才可以回帖 登录 | 注-册-帐-号

本版积分规则

小黑屋|手机版|Archiver|数学建模网 ( 湘ICP备11011602号 )

GMT+8, 2025-5-9 05:22 , Processed in 0.050367 second(s), 18 queries .

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表