数学建模社区-数学中国

标题: lingo模型一直报错，求解不出来，跪求想请大佬帮帮忙！！！ [打印本页]

作者: 沈宁 时间: 2022-5-18 21:05
标题: lingo模型一直报错，求解不出来，跪求想请大佬帮帮忙！！！
这是lingo中输入的代码：

sets:
period/1..3/: ;    !t;
reservoir/1..3/: ; !i;
level/1..5/: ;    !h;
S1(period,reservoir):BM, SM, BI, SI, BA, SA, BH, SH, SDM, SDI, SDA, LDM, DM, LDI, DI,
LDA, DA, LDH, DH, LFM, FM, LFI, FI, LFA, FA;  !t,i;
S2(period,level):  R_1, R_2, R_3,Q_1, Q_2, Q_3, Sr_1, Sr_2, Sr_3, Ee_1, Ee_2, Ee_3; !h,t;
S3(level,period,reservoir):Va, p, W, G, RR;  !h,t,i;

endsets

data:
BM, BI, BA, BH, SM, SI, SA, SH, LDM, LDI, LDA, LDH, LFM, LFI, LFA, p, Q_1, Q_2, Q_3, Sr_1, Sr_2, Sr_3,
Ee_1, Ee_2, Ee_3, W, G=@ole('C:mydata.xlsx');
@ole('C:mydata.xlsx')=DM, DI, DA, DH, FM, FI, FA;
enddata

max=F1+F2-F3-F4-F5-F6;
F1=@sum(S1(t,i):BM(t,j)*SM(t,i)+BI(t,i)*SI(t,i));
F2=@sum(S1(t,i):BA(t,i)*SA(t,i))+@sum(S1(t,i):BH(t,i)*0.6*SH(t,i)*8.5);
F3=@sum(S3(h,t,i):p(h,t,i))*@sum(S1(t,i)

DM(t,i)*DM(t,i)+LDI(t,i)*DI(t,i));
F4=@sum(S3(h,t,i):p(h,t,i))*@sum(S1(t,i)

DA(t,i)*DA(t,i));
F5=@sum(S3(h,t,i):p(h,t,i))*@sum(S1(t,i)

DH(t,i)*0.6*DH(t,i)*8.5);
F6=@sum(S3(h,t,i):p(h,t,i))*@sum(S1(t,i)

FM(t,i)*FM(t,i)+LFI(t,i)*FI(t,i)+LFA(t,i)*FA(t,i));

@sum(S3(h,t,i):RR(h,t,i)-W(h,t,i)-G(h,t,i))=@sum(S3(h,t,i):Va(h,t,i));

@sum(S1(t,i):BM(t,i)+BI(t,i)+BA(t,i))
+@sum(S1(t,i):0.6*BH(t,i)*8.5)
-@sum(S3(h,t,i):p(h,t,i))*@sum(S1(t,i)

M(t,i)+DI(t,i)+DA(t,i))
-@sum(S3(h,t,i):p(h,t,i))*@sum(S1(t,i):0.6*DH(t,i)*8.5)<=Va(h,t,i);

@sum(S1(t,i):BM(t,i)+BI(t,i)+BA(t,i))+@sum(S1(t,i):0.6*BH(t,i)*8.5)
+@sum(S3(h,t,i):p(h,t,i))*@sum(S1(t,i):FM(t,i)+FI(t,i)+FA(t,i))>=Va(h,t,i);

@sum(link2(h,t)*@sum(reservoir(i)|i#eq#1

_1(h,t)+Sr_1(h,t)-Ee_1(h,t))=@sum(link2(h,t)

sum(reservoir(i)|i#eq#1:R_1(h,t));

@sum(link2(h,t)*@sum(reservoir(i)|i#eq#2

_2(h,t)+Sr_2(h,t)-Ee_2(h,t))=@sum(link2(h,t)

sum(reservoir(i)|i#eq#2:R_2(h,t));

@sum(link2(h,t)

sum(reservoir(i)|i#eq#3

_3(h,t)+Sr_3(h,t)-Ee_3(h,t))=@sum(link2(h,t)

sum(reservoir(i)|i#eq#3:R_3(h,t));

end

模型数学公式如下：

作者: 沈宁 时间: 2022-5-18 21:06
希望有懂得的老哥帮帮忙！小弟必有重谢！！！

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