没人回应啊，自己顶起。
我将程序改了哈，但还是不对（目标函数不对）：
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 ;

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