请大家帮忙看哈该lingo程序错在哪里

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

由于急于做个东西，刚开始学lingo，程序总是编不对，请大家帮忙看看错在哪里，谢谢！
model:
sets:
dmu/a1..a39/:ee;
output/o1..o4/;
do(dmu,output):f;
wdo(dmu,output):w;
endsets
data:
n=4;
f=41.374 0.76 0.01214 10.7516   
7.5057 0.692 0.00089 0.05207   
64.977 2.85 0.01989 54.8142   
160.75 1.252 0.02476 53.4763   
39.785 5.563 0.02646 5.57088   
408.5 1.298 0.09064 17.1816   
391.58 19.82 0.01033 0.20892   
188.83 29.45 0.00955 0.76588   
250.24 4.361 0.01199 0.05094   
5.3832 0.058 0.00084 0.30294   
115.85 5.526 0.01165 0.0156   
8.173 0.345 0.00093 0.04326   
92.117 9.302 0.00211 0.20876   
83.279 5.697 0.01185 0.03052   
5.7716 0.076 0.00156 4.93978   
2018.5 39.17 0.04799 0.05951   
6.2366 0.179 0.00086 0.00539   
2.8194 0.132 0.00042 5.47648   
42.478 10.07 0.03094 5.20861   
168.74 65.45 0.03044 0.09468   
164.89 4.455 0.00601 0.25325   
322.08 9.937 0.03895 0.24721   
18.543 1.284 0.01 0.08281   
3.064 0.26 0.00113 2.90088   
21.045 0.93 0.08005 0.73361   
43.218 4.642 0.02656 1.09635   
35.178 10.17 0.0581 0.0869   
11.856 0.793 0.00219 0.16704   
7.3769 0.264 0.00224 0.14108   
18.473 4.966 0.00324 0.33516   
17.068 2.088 0.00218 0.01928   
5.4787 0.196 0.00081 0.02805   
3.0761 0.435 0.00047 0.01578   
19.818 0.757 0.00199 0.01822   
5.1734 0.343 0.00085 0.02332   
109.35 11.87 0.01363 0.8445   
32.229 0.812 0.21538 4.7541   
59.355 33.63 0.03417 1.57874   
133.09 3.28 0.0013 0.02362;
enddata
MIN=@SUM(DMU:EE);
@for(dmu(i):
ee(i)=@sum(f(j)|j #le# n:f(i,j)*w(i,j));
@sum(w(i,j))=1;
@for(dmu(k):
  @sum(f(j)|j #le# n:f(k,j)*w(k,j))>=1
)
);
@for(wdo(i,j):@bnd(0.000001,x,1));
end

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：lingo ING Lin NGO output 程序 lingo 帮忙 哈该

没人回应啊，自己顶起。
我将程序改了哈，但还是不对（目标函数不对）：
model:
sets:
dmu/a1..a39/;
out/o1..o4/;
do(dmu,out):f,w;
endsets
data:
f=41.374 0.76 0.01214 10.7516   
7.5057 0.692 0.00089 0.05207   
64.977 2.85 0.01989 54.8142   
160.75 1.252 0.02476 53.4763   
39.785 5.563 0.02646 5.57088   
408.5  1.298 0.09064 17.1816   
391.58 19.82 0.01033 0.20892   
188.83 29.45 0.00955 0.76588   
250.24 4.361 0.01199 0.05094   
5.3832 0.058 0.00084 0.30294   
115.85 5.526 0.01165 0.0156   
8.173  0.345 0.00093 0.04326   
92.117 9.302 0.00211 0.20876   
83.279 5.697 0.01185 0.03052   
5.7716 0.076 0.00156 4.93978   
2018.5 39.17 0.04799 0.05951   
6.2366 0.179 0.00086 0.00539   
2.8194 0.132 0.00042 5.47648   
42.478 10.07 0.03094 5.20861   
168.74 65.45 0.03044 0.09468   
164.89 4.455 0.00601 0.25325   
322.08 9.937 0.03895 0.24721   
18.543 1.284 0.01   0.08281   
3.064  0.26 0.00113 2.90088   
21.045 0.93 0.08005 0.73361   
43.218 4.642 0.02656 1.09635   
35.178 10.17 0.0581 0.0869   
11.856 0.793 0.00219 0.16704   
7.3769 0.264 0.00224 0.14108   
18.473 4.966 0.00324 0.33516   
17.068 2.088 0.00218 0.01928   
5.4787 0.196 0.00081 0.02805   
3.0761 0.435 0.00047 0.01578   
19.818 0.757 0.00199 0.01822   
5.1734 0.343 0.00085 0.02332   
109.35 11.87 0.01363 0.8445   
32.229 0.812 0.21538 4.7541   
59.355 33.63 0.03417 1.57874   
133.09 3.28 0.0013 0.02362;   
enddata
@for(dmu(m)|m #le# 39:
max=1/(@sum(out(i):f(m,i)*w(m,i)));
@sum(out(i):w(m,i))=1;
@for(dmu(k):
        @sum(out(j):f(k,j)*w(m,j))>=1;
  )
);
end

部分扩展模型如下，从中清楚地看出目标函数有问题：
MODEL:
   [_1559] MAX= 1 / ( ( 133.09 * W_A39_O1 + 3.28 * W_A39_O2 +
      0.0013 * W_A39_O3 + 0.02362 * W_A39_O4 ) ) ;
   [_2] W_A1_O1 + W_A1_O2 + W_A1_O3 + W_A1_O4 = 1 ;
   [_3] 41.374 * W_A1_O1 + 0.76 * W_A1_O2 + 0.01214 * W_A1_O3 + 10.7516 *
   W_A1_O4 >= 1 ;
   [_4] 7.5057 * W_A1_O1 + 0.692 * W_A1_O2 + 0.00089 * W_A1_O3 + 0.05207 *
   W_A1_O4 >= 1 ;
   [_5] 64.977 * W_A1_O1 + 2.85 * W_A1_O2 + 0.01989 * W_A1_O3 + 54.8142 *
   W_A1_O4 >= 1 ;
   [_6] 160.75 * W_A1_O1 + 1.252 * W_A1_O2 + 0.02476 * W_A1_O3 + 53.4763 *
   W_A1_O4 >= 1 ;
   [_7] 39.785 * W_A1_O1 + 5.563 * W_A1_O2 + 0.02646 * W_A1_O3 + 5.57088 *
   W_A1_O4 >= 1 ;
   [_8] 408.5 * W_A1_O1 + 1.298 * W_A1_O2 + 0.09064 * W_A1_O3 + 17.1816 *
   W_A1_O4 >= 1 ;
   [_9] 391.58 * W_A1_O1 + 19.82 * W_A1_O2 + 0.01033 * W_A1_O3 + 0.20892 *
   W_A1_O4 >= 1 ;
   [_10] 188.83 * W_A1_O1 + 29.45 * W_A1_O2 + 0.00955 * W_A1_O3 + 0.76588 *
   W_A1_O4 >= 1 ;
   [_11] 250.24 * W_A1_O1 + 4.361 * W_A1_O2 + 0.01199 * W_A1_O3 + 0.05094 *
   W_A1_O4 >= 1 ;
   [_12] 5.3832 * W_A1_O1 + 0.058 * W_A1_O2 + 0.00084 * W_A1_O3 + 0.30294 *
   W_A1_O4 >= 1 ;
   [_13] 115.85 * W_A1_O1 + 5.526 * W_A1_O2 + 0.01165 * W_A1_O3 + 0.0156 *
   W_A1_O4 >= 1 ;
   [_14] 8.173 * W_A1_O1 + 0.345 * W_A1_O2 + 0.0009300000000000001 *
   W_A1_O3 + 0.04326 * W_A1_O4 >= 1 ;
   [_15] 92.117 * W_A1_O1 + 9.302 * W_A1_O2 + 0.00211 * W_A1_O3 + 0.20876 *
   W_A1_O4 >= 1 ;
   [_16] 83.279 * W_A1_O1 + 5.697 * W_A1_O2 + 0.01185 * W_A1_O3 + 0.03052 *
   W_A1_O4 >= 1 ;
   [_17] 5.7716 * W_A1_O1 + 0.076 * W_A1_O2 + 0.00156 * W_A1_O3 + 4.93978 *
   W_A1_O4 >= 1 ;
   [_18] 2018.5 * W_A1_O1 + 39.17 * W_A1_O2 + 0.04799 * W_A1_O3 + 0.05951 *
   W_A1_O4 >= 1 ;
   [_19] 6.2366 * W_A1_O1 + 0.179 * W_A1_O2 + 0.00086 * W_A1_O3 + 0.00539 *
   W_A1_O4 >= 1 ;
   [_20] 2.8194 * W_A1_O1 + 0.132 * W_A1_O2 + 0.00042 * W_A1_O3 + 5.47648 *
   W_A1_O4 >= 1 ;
   [_21] 42.478 * W_A1_O1 + 10.07 * W_A1_O2 + 0.03094 * W_A1_O3 + 5.20861 *
   W_A1_O4 >= 1 ;
   [_22] 168.74 * W_A1_O1 + 65.45 * W_A1_O2 + 0.03044 * W_A1_O3 + 0.09468 *
   W_A1_O4 >= 1 ;
   [_23] 164.89 * W_A1_O1 + 4.455 * W_A1_O2 + 0.00601 * W_A1_O3 + 0.25325 *
   W_A1_O4 >= 1 ;
   [_24] 322.08 * W_A1_O1 + 9.936999999999999 * W_A1_O2 + 0.03895 * W_A1_O3
   + 0.24721 * W_A1_O4 >= 1 ;
   [_25] 18.543 * W_A1_O1 + 1.284 * W_A1_O2 + 0.01 * W_A1_O3 + 0.08281 *
   W_A1_O4 >= 1 ;
   [_26] 3.064 * W_A1_O1 + 0.26 * W_A1_O2 + 0.00113 * W_A1_O3 + 2.90088 *
   W_A1_O4 >= 1 ;
   [_27] 21.045 * W_A1_O1 + 0.9300000000000001 * W_A1_O2 + 0.08005 *
   W_A1_O3 + 0.73361 * W_A1_O4 >= 1 ;
   [_28] 43.218 * W_A1_O1 + 4.642 * W_A1_O2 + 0.02656 * W_A1_O3 + 1.09635 *
   W_A1_O4 >= 1 ;
   [_29] 35.178 * W_A1_O1 + 10.17 * W_A1_O2 + 0.0581 * W_A1_O3 +
   0.08690000000000001 * W_A1_O4 >= 1 ;
   [_30] 11.856 * W_A1_O1 + 0.793 * W_A1_O2 + 0.00219 * W_A1_O3 + 0.16704 *
   W_A1_O4 >= 1 ;
   [_31] 7.3769 * W_A1_O1 + 0.264 * W_A1_O2 + 0.00224 * W_A1_O3 + 0.14108 *
   W_A1_O4 >= 1 ;
   [_32] 18.473 * W_A1_O1 + 4.966 * W_A1_O2 + 0.00324 * W_A1_O3 + 0.33516 *
   W_A1_O4 >= 1 ;
   [_33] 17.068 * W_A1_O1 + 2.088 * W_A1_O2 + 0.00218 * W_A1_O3 + 0.01928 *
   W_A1_O4 >= 1 ;
   [_34] 5.4787 * W_A1_O1 + 0.196 * W_A1_O2 + 0.00081 * W_A1_O3 + 0.02805 *
   W_A1_O4 >= 1 ;
   [_35] 3.0761 * W_A1_O1 + 0.435 * W_A1_O2 + 0.00047 * W_A1_O3 + 0.01578 *
   W_A1_O4 >= 1 ;
   [_36] 19.818 * W_A1_O1 + 0.757 * W_A1_O2 + 0.00199 * W_A1_O3 + 0.01822 *
   W_A1_O4 >= 1 ;
   [_37] 5.1734 * W_A1_O1 + 0.343 * W_A1_O2 + 0.00085 * W_A1_O3 + 0.02332 *
   W_A1_O4 >= 1 ;
   [_38] 109.35 * W_A1_O1 + 11.87 * W_A1_O2 + 0.01363 * W_A1_O3 + 0.8445 *
   W_A1_O4 >= 1 ;
   [_39] 32.229 * W_A1_O1 + 0.8120000000000001 * W_A1_O2 + 0.21538 *
   W_A1_O3 + 4.7541 * W_A1_O4 >= 1 ;
   [_40] 59.355 * W_A1_O1 + 33.63 * W_A1_O2 + 0.03417 * W_A1_O3 + 1.57874 *
   W_A1_O4 >= 1 ;
   [_41] 133.09 * W_A1_O1 + 3.28 * W_A1_O2 + 0.0013 * W_A1_O3 + 0.02362 *
   W_A1_O4 >= 1 ;
   [_42] 1 / ( ( 7.5057 * W_A2_O1 + 0.692 * W_A2_O2 + 0.00089 *
      W_A2_O3 + 0.05207 * W_A2_O4 ) ) ;
   [_43] W_A2_O1 + W_A2_O2 + W_A2_O3 + W_A2_O4 = 1 ;
   [_44] 41.374 * W_A2_O1 + 0.76 * W_A2_O2 + 0.01214 * W_A2_O3 + 10.7516 *
   W_A2_O4 >= 1 ;

请高手指正，在下感激不尽！