快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

留学>>
返回列表 发新帖
查看: 10726|回复: 9
打印 上一主题 下一主题

[问题求助] lingo编程求不出全局最优解，请大侠帮帮忙！！！

[复制链接]
字体大小: 正常 放大
jiiing        

7

主题

5

听众

28

积分

升级  24.21%

该用户从未签到
跳转到指定楼层
1#
发表于 2009-9-5 07:49 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
model:/ k! r. w# z& Q3 [
sets:
* S5 Z5 i: I% M' lshengshi/1..13/;
& D* f- Z9 J9 T2 L$ l8 u% Nzhibiao/1..6/:a;
' I- x* J: Y: h, R5 v6 H1 m3 N: k, Mlinks(shengshi,zhibiao):x;: u! P$ S+ W  E3 T& c6 y  E
endsets
( c% q( h% v9 d" i8 V' J6 Gdata:  _5 t; c, e! [" v5 v
x=
, [" X- k4 R7 U3 a+ o0.099599  0.104095  0.021186  0.489849  0.031419  0.399601 ' y  I" L  E9 f; ~
0.012222  0.006320  0.035075  0.558043  0.253776  0.465170
, j1 Z/ O5 l0 k3 S3 ^( Q5 k0.016260  0.004213  0.084423  0.631442  0.375174  0.535579
" A9 g8 w$ P$ x6 T0.004587  0.000000  0.086154  0.557522  0.210591  0.464721   v6 D) v# T* y; s: B7 h# }
0.039880  0.015333  0.077921  1.000000  0.488760  0.888523 ' ^5 H" {) w% p* g7 w- b8 `& i, K
0.003290  0.002399  0.023441  0.299844  1.000000  0.217814 , v" z8 @" H7 G( ~5 E: v
0.021750  0.023296  0.055727  0.635086  0.657561  0.538723 - v/ t( k3 i$ }: }, d* B' i
1.000000  1.000000  1.000000  0.000000  0.000000  0.430589 ' f+ F4 r3 I. Q% p
0.070963  0.277539  0.385947  0.153566  0.014757  0.000000 0 Q  `, M" m& K' x# ^4 R( S
0.038164  0.038373  0.024990  0.489849  0.031159  1.000000
+ `8 n! l2 d8 x' q' I( E0.000000  0.008735  0.015259  0.288912  0.037009  0.040918 ; h, y) p% d- S8 O! C) U* T
0.057181  0.038508  0.254109  0.184279  0.093160  0.052246 0 Z: S$ r9 [! H( Q* L/ h
0.048211  0.009641  0.000000  0.387819  0.135834  0.118513;3 V! P8 I) ?/ k  d
enddata
* A, v# y0 Y0 A, r5 r4 A4 cmax=
' N8 f" _* K5 n; f. @(@sum(shengshi(i):, |+ V2 a$ m$ h$ d! b; [
(@sum(links(i,j):a(j)*x(i,j))-@sum(shengshi(i):/ }# |- `5 Y; j
@sum(links(i,j):a(j)*x(i,j)))/13)^2)/12)^(1/2);- L' D$ N" i( Q9 H2 L! ]" O" f
6 {- t8 C3 Z" M& w$ K+ l
@sum(zhibiao(j):a(j)^2)=1;3 r, J2 C6 J" c% h. d( z
@for(zhibiao(j):a(j)>-1);) I( x% t. ?* K+ y
@for(zhibiao(j):a(j)<1);
+ s( a. O/ U7 M  I$ ]; X@for(zhibiao(j)free(a(j)));
( a( {/ u2 n( _' _7 D* U  W; ~end & b: A" ~+ ?( j9 T9 Y$ W/ r5 I
求不出全局最优解，请大侠帮帮忙！！！
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
回复

使用道具 举报

kongfanqing 实名认证       

45

主题

5

听众

2009

积分

  • TA的每日心情
    郁闷
    2020-3-29 17:53

    • 签到天数: 71 天

    [LV.6]常住居民II

    新人进步奖

    群组: 数学建模

    群组: 数学趣味、游戏、IQ等

    群组: LINGO

    群组: 学术交流D
    2#
    发表于 2009-9-5 08:32 |只看该作者 |招呼Ta 关注Ta
    运算时间太长，你能否增加一个初始值就好了
    回复

    使用道具 举报

    kongfanqing 实名认证       

    45

    主题

    5

    听众

    2009

    积分

  • TA的每日心情
    郁闷
    2020-3-29 17:53

    • 签到天数: 71 天

    [LV.6]常住居民II

    新人进步奖

    群组: 数学建模

    群组: 数学趣味、游戏、IQ等

    群组: LINGO

    群组: 学术交流D
    3#
    发表于 2009-9-5 08:47 |只看该作者 |招呼Ta 关注Ta
    Global optimal solution found.# W( ~* F5 _, Z: G$ h
      Objective value:                             0.7272518
    + W7 R  s3 k8 D+ |8 N* g4 j& E  Objective bound:                             0.7272527' H, f( V4 s) x1 g# M9 G8 R* p
      Infeasibilities:                              0.0000000 h5 F7 p  z( [; e6 x* P7 m6 S$ a. C2 `
      Extended solver steps:                           388431 r* u6 ^7 F7 s; j" z. {
      Total solver iterations:                       7561347
    9 g% I, L4 h! T
    : \; e8 e: U1 B( W. h# Q* W: O$ ^+ i' h* O
                           Variable           Value        Reduced Cost8 x7 W8 s6 s( m9 S
                              A( 1)       0.1248020            0.0000006 x! r1 r$ T/ u
                              A( 2)       0.1216772            0.000000
    4 h. t/ K+ ?/ k; n! X4 J+ w                          A( 3)       0.1621721            0.000000
    1 a$ s8 `$ [4 y( f- w5 L/ p$ S                          A( 4)       0.6318675            0.0000004 }3 |4 X' b0 ~! q( H
                              A( 5)       0.3868056            0.0000002 Z) i) g* J! l- U$ B
                              A( 6)       0.6280479            0.000000: T1 a" [; X: U# U$ {6 ?; e
                           X( 1, 1)       0.9959900E-01        0.000000
    8 v( d) a2 D3 r5 x# p; c                       X( 1, 2)       0.1040950            0.000000
    ( }, q2 H/ W7 M7 }8 w4 b# h" B                       X( 1, 3)       0.2118600E-01        0.000000' k, y5 T/ v% a( [% W5 r6 |- r
                           X( 1, 4)       0.4898490            0.000000
    4 G" D* J) c5 U  k1 r                       X( 1, 5)       0.3141900E-01        0.000000) U5 p9 M4 y& G; h/ P$ n
                           X( 1, 6)       0.3996010            0.000000
    0 v6 E# X  a1 S* b                       X( 2, 1)       0.1222200E-01        0.000000( K* g2 x/ B: A& x, ]% \9 m- G( a+ P
                           X( 2, 2)       0.6320000E-02        0.000000
    ! e" |: J( y  v" \/ c                       X( 2, 3)       0.3507500E-01        0.000000) X% z9 `, g, Y; v; S: k! W
                           X( 2, 4)       0.5580430            0.000000
    7 o0 Z7 D( j/ H/ j  ]                       X( 2, 5)       0.2537760            0.000000! o, f# e/ ?  E# _$ C& L* b
                           X( 2, 6)       0.4651700            0.000000
    ' O' \! J( G* L8 i) x1 |( @                       X( 3, 1)       0.1626000E-01        0.000000( @6 i, W8 B8 o
                           X( 3, 2)       0.4213000E-02        0.000000
    % w7 |% o4 @: O( Q5 }+ N& h( M                       X( 3, 3)       0.8442300E-01        0.000000
    . k: j9 T- `9 W* b* m& T                       X( 3, 4)       0.6314420            0.000000$ m1 x* ^7 _# @( P/ p
                           X( 3, 5)       0.3751740            0.000000& T& I3 z! O4 |  A6 Z) M
                           X( 3, 6)       0.5355790            0.000000
    8 Z9 B3 k) l. X# x, S                       X( 4, 1)       0.4587000E-02        0.000000
    : \/ E, Q9 |2 Y3 X' n                       X( 4, 2)        0.000000            0.000000. D, P$ W% O: G& M9 I
                           X( 4, 3)       0.8615400E-01        0.000000  N/ g1 ^5 d) c; y" R
                           X( 4, 4)       0.5575220            0.000000
    * v# u% Z) G! H$ \: a) s8 k7 c                       X( 4, 5)       0.2105910            0.000000" W* Y% R" J) ^$ {& U% b3 v4 f. G
                           X( 4, 6)       0.4647210            0.000000* k: t+ k2 ?  Z* D; H2 I, ?! ?# n2 D
                           X( 5, 1)       0.3988000E-01        0.000000* Y( j( y* V8 Z( z' K" y9 b
                           X( 5, 2)       0.1533300E-01        0.0000001 y, W4 G" s, ~& n6 j* S
                           X( 5, 3)       0.7792100E-01        0.0000002 M4 t3 T% @: W: A  N0 R/ l
                           X( 5, 4)        1.000000            0.000000; u  P! Q' ?; ~' @% U2 Z% E- P/ @
                           X( 5, 5)       0.4887600            0.000000* W2 r$ u  T4 n, O9 |4 r% w
                           X( 5, 6)       0.8885230            0.000000
    3 H3 i* @9 i/ B' T8 a                       X( 6, 1)       0.3290000E-02        0.000000
    # @# J  H4 q8 q+ H                       X( 6, 2)       0.2399000E-02        0.000000
    + i) s! ]( x" H% r/ k                       X( 6, 3)       0.2344100E-01        0.000000
    2 A( j) k( P9 o2 \1 g$ c. R                       X( 6, 4)       0.2998440            0.000000
    4 G$ X, N' z0 L1 Q                       X( 6, 5)        1.000000            0.0000003 _1 m6 W, H1 P
                           X( 6, 6)       0.2178140            0.000000/ ~" L8 ~+ `! ?8 R. a" k% g. @
                           X( 7, 1)       0.2175000E-01        0.000000
    , B1 w+ [8 W" D) p& C, a                       X( 7, 2)       0.2329600E-01        0.0000006 J4 k2 I; |$ X4 P2 G5 V
                           X( 7, 3)       0.5572700E-01        0.000000
    $ N1 G+ J- n/ v$ s1 k: ]1 O                       X( 7, 4)       0.6350860            0.000000
    + M% _% @' C3 g$ c, y( t) o                       X( 7, 5)       0.6575610            0.000000
    8 m+ r6 D% h( b  i                       X( 7, 6)       0.5387230            0.000000. S: g( u8 v; E
                           X( 8, 1)        1.000000            0.000000* o* z4 e4 l3 x4 M8 ?9 D
                           X( 8, 2)        1.000000            0.000000
    ) _7 w, n, V# o+ Z0 ?; i# q8 @                       X( 8, 3)        1.000000            0.000000! _: g! j! M8 [& H2 T
                           X( 8, 4)        0.000000            0.0000004 V+ m. z0 ]' z5 j1 I( T
                           X( 8, 5)        0.000000            0.000000
    ( W7 Y. F$ S* N, L0 j3 o* W4 m# N                       X( 8, 6)       0.4305890            0.000000
    - V. h. N" X7 Y5 s0 }4 m' }                       X( 9, 1)       0.7096300E-01        0.000000
    6 j, _7 n" j  v6 S, T& @                       X( 9, 2)       0.2775390            0.000000, t9 ?2 L% p5 I7 b  C" H8 L4 e
                           X( 9, 3)       0.3859470            0.000000+ b) a* y9 {' [, k' Y% _
                           X( 9, 4)       0.1535660            0.000000, D1 u/ q! n# o9 t
                           X( 9, 5)       0.1475700E-01        0.0000009 a% t3 }5 j% {5 q  r/ ^
                           X( 9, 6)        0.000000            0.000000
    . r  F- f! Y- Z& f% M) _                      X( 10, 1)       0.3816400E-01        0.0000000 q0 F1 J2 O1 X0 F7 Y3 n$ q# V
                          X( 10, 2)       0.3837300E-01        0.000000
    " _. o! F$ ^1 T2 R. N; X% s4 D0 m! ?                      X( 10, 3)       0.2499000E-01        0.000000
    ! F% m# m: l( _7 A9 u. e) I0 m7 U( P- r                      X( 10, 4)       0.4898490            0.000000
      g: w5 Y2 l" ~# V                      X( 10, 5)       0.3115900E-01        0.000000* F$ M  t- j5 k7 L( ~1 W% {) \8 n7 x+ c! [
                          X( 10, 6)        1.000000            0.0000008 |3 {, t. W1 W+ y) M6 Y9 C  O& L1 h
                          X( 11, 1)        0.000000            0.000000
    1 L, {5 J/ }: ?# ^9 X                      X( 11, 2)       0.8735000E-02        0.000000
    # O, ^4 s! p6 K8 W$ F8 W3 m                      X( 11, 3)       0.1525900E-01        0.000000/ M3 \9 Z/ l! O
                          X( 11, 4)       0.2889120            0.000000
    ; F% p4 i. C! @5 a                      X( 11, 5)       0.3700900E-01        0.000000
    ) x& m) U# y% O4 }                      X( 11, 6)       0.4091800E-01        0.000000, a3 a# o3 p. f! x
                          X( 12, 1)       0.5718100E-01        0.000000
    4 H, q% |/ N" C9 ~                      X( 12, 2)       0.3850800E-01        0.000000
    7 R7 X  [9 X/ v  f                      X( 12, 3)       0.2541090            0.0000001 S/ O1 E6 j1 N
                          X( 12, 4)       0.1842790            0.0000001 \' \4 [; d6 Q/ J6 m$ Y& e
                          X( 12, 5)       0.9316000E-01        0.000000
    3 z% {+ s  J) D+ K0 e                      X( 12, 6)       0.5224600E-01        0.000000# |4 K) \( Y7 y  {/ A, [
                          X( 13, 1)       0.4821100E-01        0.000000
    / n- q& j! z2 Y0 G  P1 J                      X( 13, 2)       0.9641000E-02        0.000000
    ( i( G: d0 L" N% C6 C1 e. x                      X( 13, 3)        0.000000            0.000000
    0 a+ N! k9 B+ Z                      X( 13, 4)       0.3878190            0.000000
    - l5 @* h4 {) S- Q) D                      X( 13, 5)       0.1358340            0.000000
    8 E" o1 Y- m: P* V7 T                      X( 13, 6)       0.1185130            0.000000; N3 v. J0 j- z5 U* I

    - f. U- h/ y! @3 d6 \                            Row    Slack or Surplus      Dual Price
    ( f$ |* z2 c7 ]  {+ P, v# L                              1       0.7272518            1.000000
    3 o: Y2 c& t$ B( N0 k                              2        0.000000           0.3636259
    ) g' [* ~( z" z1 R) v                              3        1.124802            0.000000
    0 n% V, |) Y3 O& M2 g8 o                              4        1.121677            0.000000  Q1 |- K3 O2 x  M
                                  5        1.162172            0.0000007 o- u1 @2 Q( M
                                  6        1.631867            0.000000
    $ b& b8 [" x8 @  I8 |* ~# m                              7        1.386806            0.0000003 u( q; {; p2 _. n
                                  8        1.628048            0.000000! L: E  Q5 ]! b6 D
                                  9       0.8751980            0.000000
    4 \5 Y6 ]3 S6 x* o                             10       0.8783228            0.0000002 R8 V3 [- v# y- L5 }8 H- `
                                 11       0.8378279            0.000000
    ! I. m- J8 s& ?1 p6 ^* J" `                             12       0.3681325            0.000000
      j+ E" x* y  D6 X3 P8 E                             13       0.6131944            0.000000
    5 Z4 G/ c7 D! @7 z" {, Y# B: G! V% S! G                             14       0.3719521            0.000000
    回复

    使用道具 举报

    kongfanqing 实名认证       

    45

    主题

    5

    听众

    2009

    积分

  • TA的每日心情
    郁闷
    2020-3-29 17:53

    • 签到天数: 71 天

    [LV.6]常住居民II

    新人进步奖

    群组: 数学建模

    群组: 数学趣味、游戏、IQ等

    群组: LINGO

    群组: 学术交流D
    4#
    发表于 2009-9-5 08:48 |只看该作者 |招呼Ta 关注Ta
    运算了12分钟才计算出来，给他一个初始值就好了！
    回复

    使用道具 举报

    ychexingwuji        

    4

    主题

    4

    听众

    70

    积分

    升级  68.42%

    该用户从未签到
    5#
    发表于 2009-9-5 08:49 |只看该作者 |招呼Ta 关注Ta
    学习了 呵呵
    回复

    使用道具 举报

    formulaone9275        

    1

    主题

    4

    听众

    188

    积分

    升级  44%

    该用户从未签到
    6#
    发表于 2009-9-5 08:54 |只看该作者 |招呼Ta 关注Ta
    有很多时候都是求不出最优解的，局部最优就可以用了 ，我看了你的程序，其中max里的sum有对i求了两次和，你看看能不能再改改，求求试试
    回复

    使用道具 举报

    yuan07150215        

    0

    主题

    4

    听众

    168

    积分

    升级  34%

    该用户从未签到
    7#
    发表于 2009-9-5 08:58 |只看该作者 |招呼Ta 关注Ta
    先在限制条件下得出最低值，之后运行lingo程序，当运行到20-30分钟时，人为的终止求解（按下Interrupt Solver），此时运行的解与目标函数的最优解很相近！可认为是最优解！……
    回复

    使用道具 举报

    xuyin        

    0

    主题

    4

    听众

    198

    积分

    升级  49%

  • TA的每日心情
    开心
    2016-3-23 08:32

    • 签到天数: 26 天

    [LV.4]偶尔看看III
    8#
    发表于 2009-9-5 09:06 |只看该作者 |招呼Ta 关注Ta
    哇，都是高手啊
    回复

    使用道具 举报

    liulongfei 实名认证       

    2

    主题

    2

    听众

    100

    积分

    升级  0%

    该用户从未签到

    群组: 数学趣味、游戏、IQ等

    群组: 数学建摸协会

    群组: 华中科技大学
    9#
    发表于 2009-9-5 09:35 |只看该作者 |招呼Ta 关注Ta
    选择全局最优解，然后就是时间的问题了
    回复

    使用道具 举报

    jiiing        

    7

    主题

    5

    听众

    28

    积分

    升级  24.21%

    该用户从未签到
    10#
    发表于 2009-9-5 12:16 |只看该作者 |招呼Ta 关注Ta
    但是 我设置全局求解时：两三秒就停止了: M5 ?1 f6 V) l' q* b2 h0 o
    出现错误提示
    ' ^8 ^- A0 `! r2 s  G! w还有 上面的 全局的做大值还没有局部做大的大  为什么！！！
    回复

    使用道具 举报

    返回列表 发新帖
    高级模式
    B Color Image Link Quote Code Smilies
    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩

    • 电话咨询

    • 04714969085

    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问：王兆丰

    GMT+8, 2026-10-7 10:31 , Processed in 3.617367 second(s), 103 queries .

    回顶部