LINGO 下标越界

再回首0000

model:
sets:
HOU/1..13/;
DIS/1..5/:Fd,Dd,yd;
FAC/1..3/:Ff,Df,yf;
DEA/1,2/;
CUS/1/;
CAR/1..3/:g1;
PART1/1..2/:g2;
PART2/1/:g3;
VOL(CAR,HOU):Ami;
BS1(CAR,DIS):Bmj,Smj;
BS2(PART1,FAC):Bnk,Snk;
BS3(PART2,FAC):Bok,Sok;
COST1(HOU,DIS):m1;
COST2(DIS,FAC):m2;
COST3(DIS,DEA):m3;
COST4(FAC,CUS):m4;
LONG1(HOU,DIS):d1;
LONG2(DIS,FAC):d2;
LONG3(DIS,DEA):d3;
LONG4(FAC,CUS):d4;
TIME1(HOU,DIS):h1;
TIME2(DIS,FAC):h2;
TIME3(DIS,DEA):h3;
TIME4(FAC,CUS):h4;
LINK1(CAR,HOU,DIS):c1,x1;
LINK21(PART1,DIS,FAC):c21,x2;
LINK22(PART2,DIS,FAC):c22,u2;
LINK31(PART1,DIS,DEA):c31,x3;
LINK32(PART2,DIS,DEA):c32,u3;
LINK41(PART1,FAC,CUS):c41,x4;
LINK42(PART2,FAC,CUS):c42,u4;
TLINK1(CAR,HOU,DIS):T1,v1;
TLINK21(PART1,DIS,FAC):T21,v2;
TLINK22(PART2,DIS,FAC):T22,w2;
TLINK31(PART1,DIS,DEA):T31,v3;
TLINK32(PART2,DIS,DEA):T32,w3;
TLINK41(PART1,FAC,CUS):T41,v4;
TLINK42(PART2,FAC,CUS):T42,w4;
LLINK1(CAR,HOU,DIS):a1,f1;
LLINK21(PART1,DIS,FAC):a21,f2;
LLINK22(PART2,DIS,FAC):a22,j2;
LLINK31(PART1,DIS,DEA):a31,f3;
LLINK32(PART2,DIS,DEA):a32,j3;
LLINK41(PART1,FAC,CUS):a41,f4;
LLINK42(PART2,FAC,CUS):a42,j4;


endsets
data:
t=0.8;   r=0.7;  b=0.5;  tm=0.1;  tn=0.08; to=0.06; u=200;
g1=3 10 5;  g2=0.3,0.5;  g3=4.2;   e=45;

Ami=
45000     62000     33000      22000     9500    26000     13000      13400     12000     8200     17400     8000     8000
3000       3000        2100         1700       800      1700       2000         1500       1300       1100      2600      1300      1500
243          216          163            93           74        175          153            304         134          192        475        657        475
;
Smj=
0.37    0.38    0.33    0.35   0.34
0.39    0.43    0.35    0.38   0.36
0.36    0.34    0.35    0.32   0.31
;
Snk=
0.68    0.64    0.66
0.58    0.54    0.56
;
Sok=
0.43    0.39     0.41
;

Bmj=
15  10  6  5  5
2   0.5 0.5  0.2   0.2
0.3  0.3  0.5 0.2  0.2
;
Bnk=
20  10    5
10  20    2
;
Bok=150 150 15;
Fd,Dd=4200 380  4200 380  3700 250  3900 280  3800 270;  
Ff,Df=6000 380    4500 330    5000 360;
m1=
5     86.8  139.2      110   40.4
86.8  5     228    52.5  98.5
74.8  16.8  198    62.5  84
52    36.8  182    73.5  56.5
30    66    154.8  85    13.5
110   52.5  196.8  5     68.8
62.4  82    152.4  53.2  16.4
134   125   128    75.6  83.5
40.4  98.5  182.5  68.8  5
76    165   97.5   118.4 86.5
139.2 228   5      246   182.5
124.5 210   65.5   193   131.5
162.5 217.5 86.4   181.5 139
;
m2=
5       134     52
86.8    125     36.8
139.2   128     182
110     75.6    73.5
40.4    83.5    56.5
;
m3=
124.5    95
210      165
65.5     97.5
193      148
131.5    86.5
;
m4=
5   134    52
;
h1=
0.125   1.8    2.9    2.3    0.8
1.8     0.125  4.75   1.3    2.5
1.6     0.35   4.125  1.5625 2.1
1.1     0.8    3.8    1.8375 1.4125
0.625   1.375  3.225  2.125  0.3375
2.3     1.3125 4.1    0.125  1.4
1.3     2.05   3.175  1.1    0.3
3.35    3.125  2.7    1.575  2.0875
0.8     2.4625 4.5625 1.4    0.125
1.6     4.125  2.4375 2.5    2.1625
2.9     4.75   0.125  6.15   4.5625
3.1125  5.25   1.6375 4.825  3.275
4.0625  5.4375 1.8    4.5375 3.475
;
h2=
0.125    3.35   1.1  
1.8      3.125  0.8
2.9      2.7    3.8
2.3      1.575  1.8375
0.8      2.0875 1.4125
;
h3=
3.1125  2.375
5.25    4.125
1.6375  2.4375
4.825   3.7
3.2875  2.1625
;
h4=
0.125  3.35  1.1
;
d1=
10   217  348  275  101
217  10   570  105  197
187  42   495  125  168
130  92   455  147  113
75   165  387  170  27
275  105  492  10   172
156  164  381  133  41
268  250  320  189  167
101  197  365  172  10
190  330  195  296  173
348  570  10   492  365
249  420  131  386  263
325  435  216  363  278
;

d2=
10  268  130
217 250  92
348 320  455
275 189  147
101 167  113
;

d3=
249  190
420  330
131  195
386  296
263  173
;
d4=
10 268 130;

Nd=3;
Nf=2;
enddata

MIN=@SUM(DIS(j):Fd(j)*yd(j))+@sum(FAC(k):Ff(k)*yf(k)) +@SUM(DIS(j):Dd(j)*yd(j))+@sum(FAC(k):Df(k)*yf(k)) +@SUM(DIS(j):@SUM(CAR(m):(Smj(m,j)*@SUM(HOU(i):x1(m,i,j)))))
+@SUM(FAC(k):@SUM(PART1(n):(Snk(n,k)*@SUM(DIS(j):x2(n,j,k)))))+@SUM(FAC(k):@SUM(PART2(o):(Sok(o,k)*@SUM(DIS(j):u2(o,j,k)))))
+(@SUM(CAR(m):@SUM(HOU(i):@SUM(DIS(j):c1*x1)))+@SUM(PART1(n):@SUM(DIS(j):@SUM(FAC(k):c21*x2)))+@SUM(PART2(o):@SUM(DIS(j):@SUM(FAC(k):c22*u2)))
+@SUM(PART1(n):@SUM(DIS(j):@SUM(DEA(l):c31*x3)))+@SUM(PART2(o):@SUM(DIS(j):@SUM(DEA(l):c32*u3)))+@SUM(PART1(n):@SUM(FAC(k):@SUM(CUS(p):c41*x4)))
+@SUM(PART2(o):@SUM(FAC(k):@SUM(CUS(p):c42*u4))))/10000
+(@SUM(CAR(m):@SUM(HOU(i):@SUM(DIS(j):T1*v1)))+@SUM(PART1(n):@SUM(DIS(j):@SUM(FAC(k):T21*v2)))+@SUM(PART2(o):@SUM(DIS(j):@SUM(FAC(k):T22*w2)))
+@SUM(PART1(n):@SUM(DIS(j):@SUM(DEA(l):T31*v3)))+@SUM(PART2(o):@SUM(DIS(j):@SUM(DEA(l):T32*w3)))+@SUM(PART1(n):@SUM(FAC(k):@SUM(CUS(p):T41*v4)))
+@SUM(PART2(o):@SUM(FAC(k):@SUM(CUS(p):T42*w4))))/10000
+u*(tm*(@SUM(DIS(j):@SUM(CAR(m):@SUM(HOU(i):a1*f1))))
+tn*(@SUM(FAC(k):@SUM(PART1(n):@SUM(DIS(j):a21*f2)))+@SUM(DEA(l):@SUM(PART1(n):@SUM(DIS(j):a31*f3))))
+@SUM(CUS(p):@SUM(PART1(n):@SUM(FAC(k):a41*f4))))
+to*(@SUM(FAC(k):@SUM(PART2(o):@SUM(DIS(j):a22*j2)))+@SUM(DEA(l):@SUM(PART2(o):@SUM(DIS(j):a32*j3)))
+@SUM(CUS(p):@SUM(PART2(o):(@SUM(FAC(k):a42*j4)))))
+e*(@SUM(DIS(j):@SUM(CAR(m):@SUM(HOU(i):f1/a1)))
+@SUM(FAC(k):@SUM(PART1(n):@SUM(DIS(j):f2/a21))))+@SUM(DEA(l):@SUM(PART1(n):@SUM(DIS(j):f3/a31)))
+@SUM(CUS(p):@SUM(PART1(n):@SUM(FAC(k):f4/a41)))
+@SUM(FAC(k):@SUM(PART2(o):@SUM(DIS(j):j2/a22)))+@SUM(DEA(l):@SUM(PART2(o):@SUM(DIS(j):j3/a32)))
+@SUM(CUS(p):@SUM(PART2(o):@SUM(FAC(k):j4/a42)))
;


@FOR(CAR(m):@For(HOU(i):@FOR(DIS(j):c1(m,i,j)=g1(m)*m1(i,j))));
@FOR(PART1(n):@For(DIS(j):@FOR(FAC(k):c21(n,j,k)=g2(n)*m2(j,k))));
@FOR(PART1(n):@For(DIS(j):@FOR(DEA(l):c31(n,j,l)=g2(n)*m3(j,l))));
@FOR(PART1(n):@For(FAC(k):@FOR(CUS(p):c41(n,k,p)=g2(n)*m4(k,p))));
@FOR(PART2(o):@For(DIS(j):@FOR(FAC(k):c22(o,j,k)=g3(o)*m2(j,k))));
@FOR(PART2(o):@For(DIS(j):@FOR(DEA(l):c32(o,j,l)=g3(o)*m3(j,l))));
@FOR(PART2(o):@For(FAC(k):@FOR(CUS(p):c42(o,k,p)=g3(o)*m4(k,p))));
@FOR(CAR(m):@For(HOU(i):@FOR(DIS(j):T1(m,i,j)=b*g1(m)*h1(i,j))));
@FOR(PART1(n):@For(DIS(j):@FOR(FAC(k):T21(n,j,k)=b*g2(n)*h2(j,k))));
@FOR(PART1(n):@For(DIS(j):@FOR(DEA(l):T31(n,j,l)=b*g2(n)*h3(j,l))));
@FOR(PART1(n):@For(FAC(k):@FOR(CUS(p):T41(n,k,p)=b*g2(n)*h4(k,p))));
@FOR(PART2(o):@For(DIS(j):@FOR(FAC(k):T22(o,j,k)=b*g3(o)*h2(j,k))));
@FOR(PART2(o):@For(DIS(j):@FOR(DEA(l):T32(o,j,l)=b*g3(o)*h3(j,l))));
@FOR(PART2(o):@For(FAC(k):@FOR(CUS(p):T42(o,k,p)=b*g3(o)*h4(k,p))));
@FOR(CAR(m):@For(HOU(i):@FOR(DIS(j):a1=d1(i,j))));
@FOR(PART1(n):@For(DIS(j):@FOR(FAC(k):a21=d2(j,k))));
@FOR(PART1(n):@For(DIS(j):@FOR(DEA(l):a31=d3(j,l))));
@FOR(PART1(n):@For(FAC(k):@FOR(CUS(p):a41=d4(k,p))));
@FOR(PART2(o):@For(DIS(j):@FOR(FAC(k):a22=d2(j,k))));
@FOR(PART2(o):@For(DIS(j):@FOR(DEA(l):a32=d3(j,l))));
@FOR(PART2(o):@For(FAC(k):@FOR(CUS(p):a42=d4(k,p))));


@FOR(CAR(m):@FOR(HOU(i):@SUM(DIS(j):x1(m,i,j))=Ami(m,i)));
@FOR(DIS(j):@SUM(CAR(m):@SUM(HOU(i):x1(m,i,j)))=@SUM(FAC(k):x2(1,j,k))+@SUM(DEA(l):x3(1,j,l)));
@FOR(DIS(j):@SUM(CAR(m):@SUM(HOU(i):x1(m,i,j)))=@SUM(FAC(k):x2(2,j,k))+@SUM(DEA(l):x3(2,j,l)));
@for(dis(j):r*@sum(car(m):g1(m)*@sum(hou(i):x1(m,i,j)))=@SUM(FAC(k):u2(1,j,k))+@SUM(DEA(l):u3(1,j,l)));
@FOR(DIS(j):t*@SUM(CAR(m):@SUM(HOU(i):x1(m,i,j)))<=@SUM(FAC(k):x2(1,j,k)));
@FOR(DIS(j):t*@SUM(CAR(m):@SUM(HOU(i):x1(m,i,j)))<=@SUM(FAC(k):x2(2,j,k)));
@for(dis(j):t*r*@sum(car(m):g1(m)*@sum(hou(i):x1(m,i,j)))<=@SUM(FAC(k):u2(1,j,k)));
@for(fac(k):@sum(dis(j):x2(1,j,k))=@sum(cus(p):x4(1,k,p)));
@for(fac(k):@sum(dis(j):x2(2,j,k))=@sum(cus(p):x4(2,k,p)));
@for(fac(k):@sum(dis(j):u2(1,j,k))=@sum(cus(p):u4(1,k,p)));
@for(dis(j):@for(car(m):@sum(hou(i):x1(m,i,j))<=10000*bmj(m,j)*yd(j)));
@for(fac(k):@for(part1(n):@sum(dis(j):x2(n,j,k))<=10000*bnk(n,k)*yf(k)));
@for(fac(k):@for(part2(o):@sum(dis(j):u2(o,j,k))<=10000*bok(o,k)*yf(k)));

@FOR(CAR(m):@FOR(HOU(i):@SUM(DIS(j):v1(m,i,j))=Ami(m,i)));
@FOR(DIS(j):@SUM(CAR(m):@SUM(HOU(i):v1(m,i,j)))=@SUM(FAC(k):v2(1,j,k))+@SUM(DEA(l):v3(1,j,l)));
@FOR(DIS(j):@SUM(CAR(m):@SUM(HOU(i):v1(m,i,j)))=@SUM(FAC(k):v2(2,j,k))+@SUM(DEA(l):v3(2,j,l)));
@for(dis(j):r*@sum(car(m):g1(m)*@sum(hou(i):v1(m,i,j)))=@SUM(FAC(k):w2(1,j,k))+@SUM(DEA(l):w3(1,j,l)));
@FOR(DIS(j):t*@SUM(CAR(m):@SUM(HOU(i):v1(m,i,j)))<=@SUM(FAC(k):v2(1,j,k)));
@FOR(DIS(j):t*@SUM(CAR(m):@SUM(HOU(i):v1(m,i,j)))<=@SUM(FAC(k):v2(2,j,k)));
@for(dis(j):t*r*@sum(car(m):g1(m)*@sum(hou(i):v1(m,i,j)))<=@SUM(FAC(k):w2(1,j,k)));
@for(fac(k):@sum(dis(j):v2(1,j,k))=@sum(cus(p):v4(1,k,p)));
@for(fac(k):@sum(dis(j):v2(2,j,k))=@sum(cus(p):v4(2,k,p)));
@for(fac(k):@sum(dis(j):w2(1,j,k))=@sum(cus(p):w4(1,k,p)));
@for(dis(j):@for(car(m):@sum(hou(i):v1(m,i,j))<=10000*bmj(m,j)*yd(j)));
@for(fac(k):@for(part1(n):@sum(dis(j):v2(n,j,k))<=10000*bnk(n,k)*yf(k)));
@for(fac(k):@for(part2(o):@sum(dis(j):w2(o,j,k))<=10000*bok(o,k)*yf(k)));

@FOR(CAR(m):@FOR(HOU(i):@SUM(DIS(j):f1(m,i,j))=Ami(m,i)));
@FOR(DIS(j):@SUM(CAR(m):@SUM(HOU(i):f1(m,i,j)))=@SUM(FAC(k):f2(1,j,k))+@SUM(DEA(l):f3(1,j,l)));
@FOR(DIS(j):@SUM(CAR(m):@SUM(HOU(i):f1(m,i,j)))=@SUM(FAC(k):f2(2,j,k))+@SUM(DEA(l):f3(2,j,l)));
@for(dis(j):r*@sum(car(m):g1(m)*@sum(hou(i):f1(m,i,j)))=@SUM(FAC(k):j2(1,j,k))+@SUM(DEA(l):j3(1,j,l)));
@FOR(DIS(j):t*@SUM(CAR(m):@SUM(HOU(i):f1(m,i,j)))<=@SUM(FAC(k):f2(1,j,k)));
@FOR(DIS(j):t*@SUM(CAR(m):@SUM(HOU(i):f1(m,i,j)))<=@SUM(FAC(k):f2(2,j,k)));
@for(dis(j):t*r*@sum(car(m):g1(m)*@sum(hou(i):f1(m,i,j)))<=@SUM(FAC(k):j2(1,j,k)));
@for(fac(k):@sum(dis(j):f2(1,j,k))=@sum(cus(p):f4(1,k,p)));
@for(fac(k):@sum(dis(j):f2(2,j,k))=@sum(cus(p):f4(2,k,p)));
@for(fac(k):@sum(dis(j):j2(1,j,k))=@sum(cus(p):j4(1,k,p)));
@for(dis(j):@for(car(m):@sum(hou(i):f1(m,i,j))<=10000*bmj(m,j)*yd(j)));
@for(fac(k):@for(part1(n):@sum(dis(j):f2(n,j,k))<=10000*bnk(n,k)*yf(k)));
@for(fac(k):@for(part2(o):@sum(dis(j):j2(o,j,k))<=10000*bok(o,k)*yf(k)));



@for(dis(j):@bin(yd(j)));
@for(fac(k):@bin(yf(k)));
@for(car(m):@for(hou(i):@for(dis(j):@gin(x1(m,i,j)))));
@for(part1(n)|n#lt#3:@for(dis(j):@for(fac(k):@gin(x2(n,j,k)))));
@for(part1(n)|n#lt#3:@for(dis(j):@for(dea(l):@gin(x3(n,j,l)))));
@for(part1(n)|n#lt#3:@for(fac(k):@for(cus(p):@gin(x4(n,k,p)))));
@for(car(m):@for(hou(i):@for(dis(j):@gin(v1(m,i,j)))));
@for(part1(n)|n#lt#3:@for(dis(j):@for(fac(k):@gin(v2(n,j,k)))));
@for(part1(n)|n#lt#3:@for(dis(j):@for(dea(l):@gin(v3(n,j,l)))));
@for(part1(n)|n#lt#3:@for(fac(k):@for(cus(p):@gin(v4(n,k,p)))));
@for(car(m):@for(hou(i):@for(dis(j):@gin(f1(m,i,j)))));
@for(part1(n)|n#lt#3:@for(dis(j):@for(fac(k):@gin(f2(n,j,k)))));
@for(part1(n)|n#lt#3:@for(dis(j):@for(dea(l):@gin(f3(n,j,l)))));
@for(part1(n)|n#lt#3:@for(fac(k):@for(cus(p):@gin(f4(n,k,p)))));

end

madio

问题出在目标函数上！写成下面的就可以运行，找到最优解，所以可能是一部分有问题。


MIN=@SUM(DIS(j):Fd(j)*yd(j))
+@sum(FAC(k):Ff(k)*yf(k))
+@SUM(DIS(j):Dd(j)*yd(j))
+@sum(FAC(k):Df(k)*yf(k))
+@SUM(DIS(j):@SUM(CAR(m):(Smj(m,j)*@SUM(HOU(i):x1(m,i,j)))))
+@SUM(FAC(k):@SUM(PART1(n):(Snk(n,k)*@SUM(DIS(j):x2(n,j,k)))))
+@SUM(FAC(k):@SUM(PART2(o):(Sok(o,k)*@SUM(DIS(j):u2(o,j,k)))))
+@SUM(CAR(m):@SUM(HOU(i):@SUM(DIS(j):c1*x1)))
+@SUM(PART1(n):@SUM(DIS(j):@SUM(FAC(k):c21*x2)))
+@SUM(PART2(o):@SUM(DIS(j):@SUM(FAC(k):c22*u2)))
+@SUM(PART1(n):@SUM(DIS(j):@SUM(DEA(l):c31*x3)))
+@SUM(PART2(o):@SUM(DIS(j):@SUM(DEA(l):c32*u3)))
+@SUM(PART1(n):@SUM(FAC(k):@SUM(CUS(p):c41*x4)))
+@SUM(PART2(o):@SUM(FAC(k):@SUM(CUS(p):c42*u4)))/10000
+@SUM(CAR(m):@SUM(HOU(i):@SUM(DIS(j):T1*v1)))
+@SUM(PART1(n):@SUM(DIS(j):@SUM(FAC(k):T21*v2)))
+@SUM(PART2(o):@SUM(DIS(j):@SUM(FAC(k):T22*w2)))
+@SUM(PART1(n):@SUM(DIS(j):@SUM(DEA(l):T31*v3)))
+@SUM(PART2(o):@SUM(DIS(j):@SUM(DEA(l):T32*w3)))
+@SUM(PART1(n):@SUM(FAC(k):@SUM(CUS(p):T41*v4)))
+@SUM(PART2(o):@SUM(FAC(k):@SUM(CUS(p):T42*w4)))/10000

;



出现错误的就在这一段目标函数的代码中：
+u*(tm*(@SUM(DIS(j):@SUM(CAR(m):@SUM(HOU(i):a1*f1)))))
+tn*(@SUM(FAC(k):@SUM(PART1(n):@SUM(DIS(j):a21*f2))))
+@SUM(DEA(l):@SUM(PART1(n):@SUM(DIS(j):a31*f3)))
+@SUM(CUS(p):@SUM(PART1(n):@SUM(FAC(k):a41*f4)))
+to*(@SUM(FAC(k):@SUM(PART2(o):@SUM(DIS(j):a22*j2))))
+@SUM(DEA(l):@SUM(PART2(o):@SUM(DIS(j):a32*j3)))
+@SUM(CUS(p):@SUM(PART2(o):@SUM(FAC(k):a42*j4)))
+e*(@SUM(DIS(j):@SUM(CAR(m):@SUM(HOU(i):f1/a1))))
+@SUM(FAC(k):@SUM(PART1(n):@SUM(DIS(j):f2/a21)))
+@SUM(DEA(l):@SUM(PART1(n):@SUM(DIS(j):f3/a31)))
+@SUM(CUS(p):@SUM(PART1(n):@SUM(FAC(k):f4/a41)))
+@SUM(FAC(k):@SUM(PART2(o):@SUM(DIS(j):j2/a22)))
+@SUM(DEA(l):@SUM(PART2(o):@SUM(DIS(j):j3/a32)))
+@SUM(CUS(p):@SUM(PART2(o):@SUM(FAC(k):j4/a42)))

再回首0000

madio 发表于 2017-3-2 22:16
问题出在目标函数上！写成下面的就可以运行，找到最优解，所以可能是一部分有问题。

确实是的  我在检查这段程序 谢谢！

再回首0000

madio 发表于 2017-3-2 22:16
问题出在目标函数上！写成下面的就可以运行，找到最优解，所以可能是一部分有问题。

谢谢！按照你的说法  我重新敲入代码  问题已经解决！！