数学建模社区-数学中国

标题: lingo高手帮忙调试下！！！万分感谢 [打印本页]
作者: jiiing    时间: 2009-9-4 14:18
标题: lingo高手帮忙调试下！！！万分感谢
model:
sets:
shengshi/1..13/;
zhibiao/1..6/:a;
links(shengshi,zhibiao):x;
endsets
data:
x=
0.099599  0.104095  0.021186  0.489849  0.031419  0.399601
0.012222  0.006320  0.035075  0.558043  0.253776  0.465170
0.016260  0.004213  0.084423  0.631442  0.375174  0.535579
0.004587  0.000000  0.086154  0.557522  0.210591  0.464721
0.039880  0.015333  0.077921  1.000000  0.488760  0.888523
0.003290  0.002399  0.023441  0.299844  1.000000  0.217814
0.021750  0.023296  0.055727  0.635086  0.657561  0.538723
1.000000  1.000000  1.000000  0.000000  0.000000  0.430589
0.070963  0.277539  0.385947  0.153566  0.014757  0.000000
0.038164  0.038373  0.024990  0.489849  0.031159  1.000000
0.000000  0.008735  0.015259  0.288912  0.037009  0.040918
0.057181  0.038508  0.254109  0.184279  0.093160  0.052246
0.048211  0.009641  0.000000  0.387819  0.135834  0.118513;
enddata
[obj]max=
@sum(shengshi@sum(links(i,j):a(j)*x(i,j))-@sum(shengshi(i)sum(links(i,j):a(j)*x(i,j)))/13)^2/12)^(1/2);
@sum(zhibiao(j):a(j)^2)=1;
@for(zhibiao(j):-1<a(j)<1);
end
作者: jiiing    时间: 2009-9-4 14:20
表示：      （
   表示   ：   @
作者: shumo_bin    时间: 2009-9-4 15:12
model:
sets:
shengshi/1..13/;
zhibiao/1..6/:a;
links(shengshi,zhibiao):x;
endsets
data:
x=
0.099599  0.104095  0.021186  0.489849  0.031419  0.399601
0.012222  0.006320  0.035075  0.558043  0.253776  0.465170
0.016260  0.004213  0.084423  0.631442  0.375174  0.535579
0.004587  0.000000  0.086154  0.557522  0.210591  0.464721
0.039880  0.015333  0.077921  1.000000  0.488760  0.888523
0.003290  0.002399  0.023441  0.299844  1.000000  0.217814
0.021750  0.023296  0.055727  0.635086  0.657561  0.538723
1.000000  1.000000  1.000000  0.000000  0.000000  0.430589
0.070963  0.277539  0.385947  0.153566  0.014757  0.000000
0.038164  0.038373  0.024990  0.489849  0.031159  1.000000
0.000000  0.008735  0.015259  0.288912  0.037009  0.040918
0.057181  0.038508  0.254109  0.184279  0.093160  0.052246
0.048211  0.009641  0.000000  0.387819  0.135834  0.118513;
enddata
max=@sum(shengshisum(links(i,j):a(j)*x(i,j))-@sum(shengshi(i)sum(links(i,j):a(j)*x(i,j))/13)^2/12)^(1/2);
@sum(zhibiao(j):a(j)^2)=1;
@for(zhibiao(j):a(j)>-1);
@for(zhibiao(j):a(j)<1);
end
作者: shumo_bin    时间: 2009-9-4 15:15
Local optimal solution found.
  Objective value:                              10.70675
  Infeasibilities:                              0.000000
  Total solver iterations:                            19


                       Variable           Value        Reduced Cost
                          A( 1)       0.1594391          -0.7405255E-08
                          A( 2)       0.1725755          -0.1392315E-07
                          A( 3)       0.2330697          -0.3587752E-07
                          A( 4)       0.6408935          -0.6208556E-07
                          A( 5)       0.3758956          -0.1745871E-07
                          A( 6)       0.5817503            0.000000
                       X( 1, 1)       0.9959900E-01        0.000000
                       X( 1, 2)       0.1040950            0.000000
                       X( 1, 3)       0.2118600E-01        0.000000
                       X( 1, 4)       0.4898490            0.000000
                       X( 1, 5)       0.3141900E-01        0.000000
                       X( 1, 6)       0.3996010            0.000000
                       X( 2, 1)       0.1222200E-01        0.000000
                       X( 2, 2)       0.6320000E-02        0.000000
                       X( 2, 3)       0.3507500E-01        0.000000
                       X( 2, 4)       0.5580430            0.000000
                       X( 2, 5)       0.2537760            0.000000
                       X( 2, 6)       0.4651700            0.000000
                       X( 3, 1)       0.1626000E-01        0.000000
                       X( 3, 2)       0.4213000E-02        0.000000
                       X( 3, 3)       0.8442300E-01        0.000000
                       X( 3, 4)       0.6314420            0.000000
                       X( 3, 5)       0.3751740            0.000000
                       X( 3, 6)       0.5355790            0.000000
                       X( 4, 1)       0.4587000E-02        0.000000
                       X( 4, 2)        0.000000            0.000000
                       X( 4, 3)       0.8615400E-01        0.000000
                       X( 4, 4)       0.5575220            0.000000
                       X( 4, 5)       0.2105910            0.000000
                       X( 4, 6)       0.4647210            0.000000
                       X( 5, 1)       0.3988000E-01        0.000000
                       X( 5, 2)       0.1533300E-01        0.000000
                       X( 5, 3)       0.7792100E-01        0.000000
                       X( 5, 4)        1.000000            0.000000
                       X( 5, 5)       0.4887600            0.000000
                       X( 5, 6)       0.8885230            0.000000
                       X( 6, 1)       0.3290000E-02        0.000000
                       X( 6, 2)       0.2399000E-02        0.000000
                       X( 6, 3)       0.2344100E-01        0.000000
                       X( 6, 4)       0.2998440            0.000000
                       X( 6, 5)        1.000000            0.000000
                       X( 6, 6)       0.2178140            0.000000
                       X( 7, 1)       0.2175000E-01        0.000000
                       X( 7, 2)       0.2329600E-01        0.000000
                       X( 7, 3)       0.5572700E-01        0.000000
                       X( 7, 4)       0.6350860            0.000000
                       X( 7, 5)       0.6575610            0.000000
                       X( 7, 6)       0.5387230            0.000000
                       X( 8, 1)        1.000000            0.000000
                       X( 8, 2)        1.000000            0.000000
                       X( 8, 3)        1.000000            0.000000
                       X( 8, 4)        0.000000            0.000000
                       X( 8, 5)        0.000000            0.000000
                       X( 8, 6)       0.4305890            0.000000
                       X( 9, 1)       0.7096300E-01        0.000000
                       X( 9, 2)       0.2775390            0.000000
                       X( 9, 3)       0.3859470            0.000000
                       X( 9, 4)       0.1535660            0.000000
                       X( 9, 5)       0.1475700E-01        0.000000
                       X( 9, 6)        0.000000            0.000000
                      X( 10, 1)       0.3816400E-01        0.000000
                      X( 10, 2)       0.3837300E-01        0.000000
                      X( 10, 3)       0.2499000E-01        0.000000
                      X( 10, 4)       0.4898490            0.000000
                      X( 10, 5)       0.3115900E-01        0.000000
                      X( 10, 6)        1.000000            0.000000
                      X( 11, 1)        0.000000            0.000000
                      X( 11, 2)       0.8735000E-02        0.000000
                      X( 11, 3)       0.1525900E-01        0.000000
                      X( 11, 4)       0.2889120            0.000000
                      X( 11, 5)       0.3700900E-01        0.000000
                      X( 11, 6)       0.4091800E-01        0.000000
                      X( 12, 1)       0.5718100E-01        0.000000
                      X( 12, 2)       0.3850800E-01        0.000000
                      X( 12, 3)       0.2541090            0.000000
                      X( 12, 4)       0.1842790            0.000000
                      X( 12, 5)       0.9316000E-01        0.000000
                      X( 12, 6)       0.5224600E-01        0.000000
                      X( 13, 1)       0.4821100E-01        0.000000
                      X( 13, 2)       0.9641000E-02        0.000000
                      X( 13, 3)        0.000000            0.000000
                      X( 13, 4)       0.3878190            0.000000
                      X( 13, 5)       0.1358340            0.000000
                      X( 13, 6)       0.1185130            0.000000
作者: shumo_bin    时间: 2009-9-4 15:17
不过不是全局最优
作者: jiiing    时间: 2009-9-4 15:24
谢谢高手帮忙！！！
作者: jiiing    时间: 2009-9-4 15:42
怎样才能得到全局最优解
作者: tiger659    时间: 2009-9-4 15:52
model:
sets:
shengshi/1..13/;
zhibiao/1..6/:a;
links(shengshi,zhibiao):x;
endsets
data:
x=
0.099599  0.104095  0.021186  0.489849  0.031419  0.399601
0.012222  0.006320  0.035075  0.558043  0.253776  0.465170
0.016260  0.004213  0.084423  0.631442  0.375174  0.535579
0.004587  0.000000  0.086154  0.557522  0.210591  0.464721
0.039880  0.015333  0.077921  1.000000  0.488760  0.888523
0.003290  0.002399  0.023441  0.299844  1.000000  0.217814
0.021750  0.023296  0.055727  0.635086  0.657561  0.538723
1.000000  1.000000  1.000000  0.000000  0.000000  0.430589
0.070963  0.277539  0.385947  0.153566  0.014757  0.000000
0.038164  0.038373  0.024990  0.489849  0.031159  1.000000
0.000000  0.008735  0.015259  0.288912  0.037009  0.040918
0.057181  0.038508  0.254109  0.184279  0.093160  0.052246
0.048211  0.009641  0.000000  0.387819  0.135834  0.118513;
enddata
[obj]max=
(@sum(shengshi(i)sum(links(i,j):a(j)*x(i,j))-@sum(shengshi(i)sum(links(i,j):a(j)*x(i,j)))/13)^2/12)^(1/2);
@sum(zhibiao(j):a(j)^2)=1;
@for(zhibiao(j)bnd(-1,a(j),1));
end
作者: tiger659    时间: 2009-9-4 15:53
设置global就可以
作者: jiiing    时间: 2009-9-5 07:25
设置global还是求不出。。。
大侠帮忙
作者: shumo_bin    时间: 2009-9-5 20:16
非线性规划 难求全局最优 lingo难以判断=====个人见解
作者: lirui0081    时间: 2009-9-5 20:56
这个问题 可以尝试用matlab 的遗传 蚁群算法之类的尝试下



欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5