LINGO8.0 for windows软件及应用（编译）

水葫芦娃 · 发表于 2005-9-5 01:43:02

<

>下载了lingo,就看到了你的热帖,麻烦你也发到我的有箱一份.非常感激!

水葫芦娃 · 发表于 2005-9-5 01:46:21

<

><a href="mailthuluwa1961@sina.com" target="_blank" >huluwa1961@sina.com</A>

lishawn · 发表于 2005-9-5 02:11:25

<

>请教一个问题，lingo程序如下：
<

>model: sets: ai/1..10/:a; bi/1..5/:b; endsets init: b=10; endinit @for(ai(i):a(i)=i); @for(bi(i)

bnd(1,b(i),10)); @for(bi(i)

gin(b(i))); <a href="mailtmin=@sum(bi(i)

sum(ai(j):@if(j" target="_blank" >min=@sum(bi(i):@sum(ai(j):@if(j</A> #eq# b(i),a(j),0))); end
<

>本想结果是b(i)都等于1，min=5。但运行后min=0，何解？
请指教。

wanbaocheng · 发表于 2005-9-5 02:20:46

水葫芦娃:已发送，请查收！

lishawn · 发表于 2005-9-5 03:16:26

刚才的问题解决了，谢谢～

wanbaocheng · 发表于 2005-9-5 03:19:01

<

>model: 用这段代码也表达了同样的意思
<

>sets: bi/1..5/:b; endsets
<

>@for(bi(i)

bnd(1,b(i),10)); @for(bi(i)

gin(b(i)));
<a href="mailtmin=@sum(bi:b" target="_blank" >min=@sum(bi:b</A>); end
当然，你那样做也许有别的目的。我试了，确实有问题。问题就出在@if那里，这里我通过引入0－1指示变量，就实现了@if的功能，运行结果是正确的。对于新加入的代码已经做了注释。
model: data: M1=10; !b(i)-j的上界; M2=-10;!b(i)-j的下界; enddata sets: ai/1..10/:a; bi/1..5/:b; abi(bi,ai):del1,del2,del; !0-1指示变量,共计150个; endsets @for(ai(i):a(i)=i); @for(bi(i)

bnd(1,b(i),10)); @for(bi(i):@gin(b(i))); @for(abi: @bin(del1);@bin(del2);@bin(del););
!为了引入指示变量而加入的约束; @for(abi(i,j): b(i)-j<=M1*del1(i,j)+0.5*(del1(i,j)-1); b(i)-j>=M2*del2(i,j)-0.5*(del2(i,j)-1); del1(i,j)+del2(i,j)-del(i,j)<=1; b(i)-j<=M1*(1-del(i,j)); b(i)-j>=M2*(1-del(i,j)); );
<a href="mailt!min=@sum(bi(i):@sum(ai(j):@if(j" target="_blank" >!min=@sum(bi(i):@sum(ai(j):@if(j</A> #eq# b(i),a(j),0))); !引入指示变量后的目标函数,实际上代替了@if的功能; <a href="mailtmin=@sum(bi(i" target="_blank" >min=@sum(bi(i</A>): @sum(ai(j): del(i,j)*a(j)) ); end

wanbaocheng · 发表于 2005-9-5 03:21:33

你是如何解决的，可以给大家说一说，交流交流嘛。

lishawn · 发表于 2005-9-5 03:50:18

<

>我发觉是@gin(x)这个函数的问题。lingo对整数有一个精度限制，存在一个绝对误差限和相对误差限。例如绝对误差限默认值为10的－6次方，则lingo会把1.999999当成整数2，所以程序中b(i)未必就取真正意义上的整数。于是，虽然有了
<

>@for(bi(i)

bnd(1,b(i),10)); @for(bi(i)

gin(b(i))); 但对于<a href="mailti@if" target="_blank" >@if</A>函数里面的判断条件j #eq# b(i),就算j取遍1到10，条件仍可能全都判假。(只要b(i)不是真正的整数）
<

>于是改写一下程序就可以了：
model: sets: ai/1..10/:a; bi/1..5/:b; endsets <a href="mailtmin=@sum(bi(i)

sum(ai(j):@if(j-b(i" target="_blank" >min=@sum(bi(i):@sum(ai(j):@if(j-b(i</A>) #le# 0.1 #or# j-b(i) #le# -0.1,a(j),0))); @for(ai(i):a(i)=i); @for(bi(i):@bnd(1,b(i),10)); @for(bi(i):@gin(b(i))); end

lishawn · 发表于 2005-9-5 04:17:59

<

>sorry,上面的代码好似还有问题~
<

>不过下面的例子可以说明一点问题：
<

>sets: ai/1..10/:a; endsets @for(ai(i):a(i)=i); @gin(c); @bnd(1,c,10); <a href="mailtmin=@sum(ai(i)

if(i" target="_blank" >min=@sum(ai(i)

if(i</A> #eq# c,a(i),0)); <a href="mailt!min=@sum(ai(i)

if(i-c" target="_blank" >!min=@sum(ai(i):@if(i-c</A> #le# 0.1 #and# i-c #ge# -0.1,a(i),0));
本来两个min的结果应该一致的，答案是1。
但实验发现上面的min取值是0,用下面的min可以得到正确结果。

lishawn · 发表于 2005-9-5 05:51:11

<

>第429楼的程序应该如下：
<

>model: sets: ai/1..10/:a; bi/1..5/:b; endsets <a href="http://www.shumo.com/bbs/mailtmin=@sum(bi(i)

sum(ai(j)

if((j-b(i" target="_blank" >min=@sum(bi(i)

sum(ai(j):@if((j-b(i</A>) #le# 0.1) #and# (j-b(i) #ge# -0.1),a(j),0))); @for(ai(i):a(i)=i); @for(bi(i):@bnd(1,b(i),10);); @for(bi(i):@gin(b(i))); end
<

>同时还需要修改options里面的整数的绝对误差限，把-6改为-8次方。
我晕了~

[此贴子已经被作者于2005-9-4 22:34:39编辑过]

		自动登录	找回密码
密码			注-册-帐-号