查看: 19801|回复: 5

【求助】29号错误！

字体大小: 正常放大

1 主题	0 听众	4 积分

升级 80%

该用户从未签到

自我介绍: aa

电梯直达

1^#

发表于 2012-12-25 16:32 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

model:
sets:
supply/p1..p12/:p;
demand/d1..d12/:d;
hub/h1..h5/:g,f;
tr/tr1..tr5/:z;
link1(supply,hub):c,w;
link2(hub,demand):h,x;
endsets
data:
p=75178 42603 36707 28838 28158 26031 21785 17587 17185 16143 13025 12976;
d=73948 40975 40578 31298 26713 25805 20192 18010 17068 15802 13206 12621;
g=3 2 3 2 3;
f=500000 400000 400000 300000 300000;
c=153 161 406 255 296
111 266 244 108 93
192 192 142 90 46
266 299 113 191 150
274 274 98 166 125
109 87 270 71 120
160 92 190 72 121;
h=153 111 192 266 274 109 160
161 266 192 299 274 87 92
406 244 142 113 98 270 190
255 108 90 191 166 71 72
296 93 46 150 125 120 121;
enddata
min=@sum(link1(k,i):c(k,i)*w(k,i))+@sum(link2(i,j):h(i,j)*x(i,j))+@sum(link1(k,i):g(i)*w(k,i))+@sum(hub(i):f(i)*z(i);
@for(supply(k):
@sum(link1(k,i):w(k.i))<=p(k));
@for(demand(j):
@sum(link2(i,j):x(i,j)>=d(j));
@for(hub(i):
@sum(link2(i,j):x(i,j)=@sum(link1(k,i):w(k,i)));
@sum(tr(i):z(i))<=3;
@sum(tr(i):@bin(z));
end

运行：
error code 29
invalid number of initialization values in a DATA/INIT/CALC statement.
the required number of value is : 60 21] 160 92 190 72 121;

哪里有问题，求大神帮忙···

zan

demand, supply, link1

转播0 淘帖0 分享0 收藏0 支持0 反对0 微信

相关帖子

使用道具举报

madio

3万主题	1310 听众	5万积分

TA的每日心情

	奋斗 2024-7-1 22:21

签到天数: 2014 天

[LV.Master]伴坛终老

自我介绍: 数学中国站长

群组: 数学建模培训课堂1

群组: 数学中国美赛辅助报名

群组: Matlab讨论组

群组: 2013认证赛A题讨论群组

群组: 2013认证赛C题讨论群组

2^#

发表于 2012-12-25 18:43 |只看该作者 |招呼Ta 关注Ta |邮箱已经成功绑定

你定义的时候定义了12个分量，给数据的时候只给了7个数据，当然会报错！

收起(2)

shsqtl :能不能帮我也看一下呀
model:
title 方案;
DATA:
r =0.9387, 0.8969, 0.9178, 0.8607, 0.9234, 0.9206, 0.9749, 1.0000
  0.5257,0.6711,0.6209,0.4336,0.5360,0.5936,0.7157,0.6702
  0.8496,0.6713,0.7680,0.7934,0.8889,0.8205,1.0000,0.9505
  0.0976,0.0081,0.2439,-0.2114,0.6992,0.5528,0.5528,1.0000
  0.2311,0.3023,1.0000,0.4101,0.6222,0.5236,0.3903,0.2481
  0.8187,0.7641,0.7376,0.5360,0.7198,0.8142,0.7842,1.0000
  0.0100,0.0096,0.0098,0.0098,0.0099,0.0103,0.0099,0.0102
  0.9227,1.0000,0.8390,0.8599,0.8535,0.7472,0.7697,0.7053
  0.8438,0.8484,0.8498,0.8945,0.9126,0.9459,0.9480,0.9309
  1.0000,0.4594,0.6762,0.6477,0.5762,0.5408,0.4781,0.5364
  0.6666,0.7404,0.6807,0.8865,0.8597,0.6902,0.8533,0.8740
  0.7571,0.7066,0.7508,0.7382,0.8013,0.7792,0.7003,0.6845
  0.8457,0.9298,0.8689,0.8154,1.0000,0.8196,0.8238,0.1023
  0.4407,0.7992,0.8538,0.7919,0.7472,0.7771,0.7836,0.6077
  0.3972,0.9519,0.5564,0.5084,0.5208,0.8167,1.0000,0.6363;
ENDDATA
sets:
canshu/1..15/:cityr,cityzongr,H,w;
city/1..9/;
link(sanshu,city):f;
endsets
min = (DI-1)^2+(CI-1)^2
k=1/@log(9);
DI=F;
CI=1-S/fmean;
@for(canshu(i):cityzongr(i)=@sum(city(j)|j#lt#9:r(i,j))+cityr(i));
@for(link(i,j)|j#lt#9:f(i,j)=r(i,j)/cityzongr(i));
@for(canshu(i):f(i,9)=cityr(i)/cityzongr(i));
@for(canshu(i):H(i)=-k*@sum(city(j):f(i,j)*@log(f(i,j))));
@for(canshu(i):w(i)=(1-H(i))/@sum(canshu(j)1-H(j))));
Fec=w(1)*f(1,9)+w(2)*f(2,9);
Fen=w(3)*f(3,9)+w(4)*f(4,9)+w(5)*f(5,9);
Fso=w(6)*f(6,9)+w(7)*f(7,9)+w(8)*f(8,9);
Fpr=w(9)*f(9,9)+w(10)*f(10,9)+w(11)*f(11,9)+w(12)*f(12,9)+w(13)*f(13,9)+w(14)*f(14,9)+w(15)*f(15,9);
DI=Fec+Fen+Fso+Fpr;
fmean=(Fec+Fen+Fso+Fpr)/4;
S=@sqrt(((Fec-fmean)^2+(Fso-fmean)^2+(Fen-fmean)^2+(Fpr-fmean)^2)/4);
end

2017-01-22 23:18 回复

shsqtl :model:
title 方案;
DATA:
r =0.9387, 0.8969, 0.9178, 0.8607, 0.9234, 0.9206, 0.9749, 1.0000
  0.5257,0.6711,0.6209,0.4336,0.5360,0.5936,0.7157,0.6702
  0.8496,0.6713,0.7680,0.7934,0.8889,0.8205,1.0000,0.9505
  0.0976,0.0081,0.2439,-0.2114,0.6992,0.5528,0.5528,1.0000
  0.2311,0.3023,1.0000,0.4101,0.6222,0.5236,0.3903,0.2481
  0.8187,0.7641,0.7376,0.5360,0.7198,0.8142,0.7842,1.0000
  0.0100,0.0096,0.0098,0.0098,0.0099,0.0103,0.0099,0.0102
  0.9227,1.0000,0.8390,0.8599,0.8535,0.7472,0.7697,0.7053
  0.8438,0.8484,0.8498,0.8945,0.9126,0.9459,0.9480,0.9309
  1.0000,0.4594,0.6762,0.6477,0.5762,0.5408,0.4781,0.5364
  0.6666,0.7404,0.6807,0.8865,0.8597,0.6902,0.8533,0.8740
  0.7571,0.7066,0.7508,0.7382,0.8013,0.7792,0.7003,0.6845
  0.8457,0.9298,0.8689,0.8154,1.0000,0.8196,0.8238,0.1023
  0.4407,0.7992,0.8538,0.7919,0.7472,0.7771,0.7836,0.6077
  0.3972,0.9519,0.5564,0.5084,0.5208,0.8167,1.0000,0.6363;
ENDDATA
sets:
canshu/1..15/:cityr,cityzongr,H,w;
city/1..9/;
link(sanshu,city):f;
endsets
min = (DI-1)^2+(CI-1)^2
k=1/@log(9);
DI=F;
CI=1-S/fmean;
@for(canshu(i):cityzongr(i)=@sum(city(j)|j#lt#9:r(i,j))+cityr(i));
@for(link(i,j)|j#lt#9:f(i,j)=r(i,j)/cityzongr(i));
@for(canshu(i):f(i,9)=cityr(i)/cityzongr(i));
@for(canshu(i):H(i)=-k*@sum(city(j):f(i,j)*@log(f(i,j))));
@for(canshu(i):w(i)=(1-H(i))/@sum(canshu(j)1-H(j))));
Fec=w(1)*f(1,9)+w(2)*f(2,9);
Fen=w(3)*f(3,9)+w(4)*f(4,9)+w(5)*f(5,9);
Fso=w(6)*f(6,9)+w(7)*f(7,9)+w(8)*f(8,9);
Fpr=w(9)*f(9,9)+w(10)*f(10,9)+w(11)*f(11,9)+w(12)*f(12,9)+w(13)*f(13,9)+w(14)*f(14,9)+w(15)*f(15,9);
DI=Fec+Fen+Fso+Fpr;
fmean=(Fec+Fen+Fso+Fpr)/4;
S=@sqrt(((Fec-fmean)^2+(Fso-fmean)^2+(Fen-fmean)^2+(Fpr-fmean)^2)/4);
end

2017-01-22 23:19 回复

我也说一句

发表

使用道具举报

丄課覺覺.

1 主题	0 听众	4 积分

升级 80%

该用户从未签到

自我介绍: aa

3^#

发表于 2012-12-25 19:31 |只看该作者 |招呼Ta 关注Ta

madio 发表于 2012-12-25 18:43
你定义的时候定义了12个分量，给数据的时候只给了7个数据，当然会报错！

谢谢····我把矩阵改过之后，又有新问题了，

error code 50
improper use of @for( ) function 34] @for(supply(k):

这个程序是我模仿其他人论文上写的，其实不懂lingo，所以很多不懂·····

点评

madio model: sets: supply/p1..p7/:p; demand/d1..d7/:d; hub/h1..h5/:g,f; tr/tr1..tr5/:z; link1(supply,hub):c,w; link2(hub,demand):h,x; endsets data: p=75178 42603 36707 28838 28158 26031 21785; d= 详情回复发表于 2012-12-25 21:05

我也说一句

发表

使用道具举报

madio

3万主题	1310 听众	5万积分

TA的每日心情

	奋斗 2024-7-1 22:21

签到天数: 2014 天

[LV.Master]伴坛终老

自我介绍: 数学中国站长

群组: 数学建模培训课堂1

群组: 数学中国美赛辅助报名

群组: Matlab讨论组

群组: 2013认证赛A题讨论群组

群组: 2013认证赛C题讨论群组

4^#

发表于 2012-12-25 21:05 |只看该作者 |招呼Ta 关注Ta |邮箱已经成功绑定

丄課覺覺. 发表于 2012-12-25 19:31
谢谢····我把矩阵改过之后，又有新问题了，

error code 50

model:
sets:
supply/p1..p7/:p;
demand/d1..d7/:d;
hub/h1..h5/:g,f;
tr/tr1..tr5/:z;
link1(supply,hub):c,w;
link2(hub,demand):h,x;
endsets
data:
p=75178 42603 36707 28838 28158 26031 21785;
d=73948 40975 40578 31298 26713 25805 20192 ;
g=3 2 3 2 3;
f=500000 400000 400000 300000 300000;
c=153 161 406 255 296
111 266 244 108 93
192 192 142 90 46
266 299 113 191 150
274 274 98 166 125
109 87 270 71 120
160 92 190 72 121;
h=153 111 192 266 274 109 160
161 266 192 299 274 87 92
406 244 142 113 98 270 190
255 108 90 191 166 71 72
296 93 46 150 125 120 121;
enddata
min=@sum(link1(k,i):c(k,i)*w(k,i))+@sum(link2(i,j):h(i,j)*x(i,j))+@sum(link1(k,i):g(i)*w(k,i))+@sum(hub(i):f(i)*z(i));
@for(supply(k):@sum(link1(k,i):w(k,i))<=p(k));
@for(demand(j):@sum(link2(i,j):x(i,j))>=d(j));
@for(hub(i):@sum(link2(i,j):x(i,j))=@sum(link1(k,i):w(k,i)));
@sum(tr(i):z(i))<=3;
@for(tr(i):@bin(z(i)));
end

错误很多，到处丢括号，我都改了，现在语法没有错误了，就是数据上我为了保证一致，删除了一些数据，最后的结果是没有可行解，我估计和数据或者模型有关，你自己对一下！