数学建模社区-数学中国

标题: lingo运行问题 [打印本页]

作者: tiant0928 时间: 2011-1-19 14:15
标题: lingo运行问题
同样的一个程序，为什么第一次运行可以，再连续进行第二次运行就出现1017的错误～～
请教高手解答！

作者: madio 时间: 2011-1-19 14:57
把程序发上来看看

作者: tiant0928 时间: 2011-1-19 16:40
model:
Title: 单周期蒸汽动力系统运行优化;

!B1约束;
95.8*YB1+0.8488*DB1O=8.89*COAL1*YB1;
160*YB1<DB1O;
DB1O<320*YB1;
@BIN(YB1);
!B2约束;
95.8*YB2+0.8488*DB2O=8.89*COAL2*YB2;
160*YB2<DB2O;
DB2O*YB2<320;
DB1O=DB2O;
@BIN(YB2);

Init:
DB1O=200;
DB2O=200;
pp=0;
sp=0;
DL4i=0;
endinit

!B3约束;
19.47*YB3+0.931*DB3O=13.5*FOB3+13.5*GB3;
80*YB3<DB3O;
140*YB3>DB3O;
@BIN(YB3);
!B4约束;
5.818*YB4+1.021*DB4O=13.5*FOB4+13.5*GB4;
DB4O>YB4*35;
DB4O<YB4*75;
@BIN(YB4);
!B5约束;
5.159*YB5+1.032*DB5O=13.5*FOB5+13.5*GB5;
DB5O>YB5*50;
DB5O<YB5*65;
@BIN(YB5);
!瓦斯约束;
YB3*GB3+YB4*GB4+YB5*GB5<12;
!C1约束;
PCC1=252*DC1I-235*DC1C1-102*DC1C2-179*YC1;
YC1*100<DC1I;
YC1*380>DC1I;
YC1*10<DC1C1;
YC1*190>DC1C1;
YC1*10<DC1C2;
YC1*180>DC1C2;
DC1C1<5*DC1C2;
DC1C2<5*DC1C1;
DC1I>DC1C1+DC1C2+10;
PCC1<50000*YC1;
PCC1>20000*YC1;
@BIN(YC1);
!C2约束;
PCC2=252*DC2I-235*DC2C1-102*DC2C2-179*YC2;
YC2*100<DC2I;
YC2*380>DC2I;
YC2*10<DC2C1;
YC2*190>DC2C1;
YC2*10<DC2C2;
YC2*180>DC2C2;
YC2*50000>PCC2;
DC2I>DC2C1+DC2C2+10;
PCC2<50000*YC2;
PCC2>20000*YC2;
@BIN(YC2);
DC2I=DC1I;
DC1C1=DC2C1;
DC1C2=DC2C2;
!T1约束;
PT1=86.2*DT1I-1459*YT1;
!YT1*50<DT1I;
!YT1*86.5>DT1I;
PT1<6500*YT1;
PT1>1500*YT1;
@BIN(YT1);
!T2约束;
PT2=88.6*DT2I-1650*YT2;
!YT2*50<DT2I;
!YT2*100>DT2I;
PT2<6500*YT2;
PT2>1500*YT2;
@BIN(YT2);
!T3约束;
PT3=73.37*DT3I-23.3*DT3C-116.74*YT3;
YT3*60<DT3I;
YT3*100>DT3I;
!YT3*10<DT3C;
!YT3*30>DT3C;
PT3<6000*YT3;
PT3>3500*YT3;
DT3C<DT3I;
@BIN(YT3);
!减温器约束;
!DL1I=0.952*DL1O;
!DL2I=0.952*DL2O;
DL1I=0.989*DL1O;
DL2I=0.989*DL2O;

DL3I=0.938*DL3O;
DL4I=0.987*DL4O;
YL4*20<DL4I;
YL4*80>DL4I;
@Bin(YL4);
!HP约束;
DB1O+DB2O>DL1I+DL2I+DC1I+DC2I;
HPI=DB1O+DB2O;
HPO=DL1I+DL2I+DC1I+DC2I;
!MP约束，购入PMP，外输SMP，工艺产生PGMP＝57，工艺需求PDMP＝336.8;
DL1O+DL2O+DC1C1+DC2C1+YB3*DB3O+YB4*DB4O+YB5*DB5O+57+PMP>SMP+236.8+DL3I+DT1I+DT2I;
MPI=DL1O+DL2O+DC1C1+DC2C1+YB3*DB3O+YB4*DB4O+YB5*DB5O+57+PMP;
MPO=SMP+236.8+DL3I+DT1I+DT2I;
!LP约束，购入PLP，外输SLP，工艺产生PGLP＝76.8，工艺需求PDLP＝204.7;
DL3O+DC1C2+DC2C2+DT1I+DT2I+76.8+PLP>204.7+DL4I+DT3I+SLP;
SLP<50;
LPI=DL3O+DC1C2+DC2C2+DT1I+DT2I+76.8+PLP;
LPO=204.7+DL4I+DT3I+SLP;
!LLP约束，购入PLLP，外输SLLP，工艺产生PGLLP＝34.9，工艺需求PDLLP＝121.3;
DL4O+DT3C+34.9+PLLP>121.3+SLLP;
SLLP<50;
LLPI=DL4O+DT3C+34.9+PLLP;
LLPO=121.3+SLLP;

!电力约束,购入PP,外输SP，工艺需求P**＝75476;
PCC1+PCC2+PT1+PT2+PT3+PP>75476+SP;
pp<50000;
PI=PCC1+PCC2+PT1+PT2+PT3+PP;
PO=75476+SP;
DATA:
!Power_Price_Purchase=?;
ENDDATA
!目标函数;
PzichanMP=(3500*(FOB3+FOB4+FOB5)+1200*(GB3+GB4+GB5))/(DB3O+DB4O+DB5O);
!       购入蒸汽费用                                        外输蒸汽费用    ;
min=180*PMP+135*PLP+80*PLLP-(110*SMP+85*SLP+45*SLLP)
!    燃煤费用                         燃油费用                   燃用瓦斯费用 ;
  +1000*(COAL1+COAL2)+2500*(FOB3+FOB4+FOB5)+1200*(GB3+GB4+GB5)
!    购电费用外输电费用                                  折旧费;
            +0.45*PP-0.35*SP
!                锅炉折旧费                ;
+320*YB1+320*YB2+200*YB3+185*YB4+170*YB5+200*YC1+200*YC2+160*YT1+140*YT2+110*YT3;

end

作者: chenhan 时间: 2011-1-22 00:23
太多了凤飞飞

作者: chenhan 时间: 2011-1-22 00:31
哈哈哈哈哈哈哈哈哈哈哈

作者: chenhan 时间: 2011-1-22 01:02
哈哈哈哈哈哈哈哈哈

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