快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

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

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

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

7

主题

5

听众

28

积分

升级  24.21%

该用户从未签到
跳转到指定楼层
1#
发表于 2009-9-5 07:49 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
model:; U$ R. ^; d  b- X2 w$ V. P
sets:
8 r* H' B1 M2 h. xshengshi/1..13/;
# |, R2 w4 u& }$ n) Lzhibiao/1..6/:a;
  K2 Z* W8 p, M5 i7 y3 O0 h! Qlinks(shengshi,zhibiao):x;
$ b9 u4 k( f: Uendsets
0 E, C6 O  c/ b" i& _# n% U! W4 c! Qdata:
$ ]$ u6 l( E2 {+ ~x=
  \4 K+ \/ ^1 k/ {5 G  {2 @$ Z" P% e( _0.099599  0.104095  0.021186  0.489849  0.031419  0.399601 7 c  O6 F# e) p& [) z3 ?5 T8 g
0.012222  0.006320  0.035075  0.558043  0.253776  0.465170 : ^) ^) L! b% |
0.016260  0.004213  0.084423  0.631442  0.375174  0.535579 5 f+ T8 k7 q+ g/ O7 K" M' e" J
0.004587  0.000000  0.086154  0.557522  0.210591  0.464721
) L; x: D$ b- B1 m0.039880  0.015333  0.077921  1.000000  0.488760  0.888523
2 z9 i5 f' t* B6 e+ }4 |' ~0.003290  0.002399  0.023441  0.299844  1.000000  0.217814
0 j0 j8 H" `8 E: p+ z0.021750  0.023296  0.055727  0.635086  0.657561  0.538723
; m, d6 J7 S. g1.000000  1.000000  1.000000  0.000000  0.000000  0.430589
( i5 m/ w. A+ x- C4 o0.070963  0.277539  0.385947  0.153566  0.014757  0.000000 + o* Q9 A9 @! F' \4 C' m2 ]* m
0.038164  0.038373  0.024990  0.489849  0.031159  1.000000
+ ]( h# l& e* l- R2 g; A0.000000  0.008735  0.015259  0.288912  0.037009  0.040918 # X9 v8 C& [; O! R
0.057181  0.038508  0.254109  0.184279  0.093160  0.052246
* g+ G. g! T9 l* w$ x; E0.048211  0.009641  0.000000  0.387819  0.135834  0.118513;
1 p/ u& [& @& D+ |+ S+ a4 menddata/ C( M& b9 @$ U2 r6 B2 ?( ^, ]! I
max=% `3 C' a8 V- p/ T) l7 D" m
(@sum(shengshi(i):
/ H1 ]$ l) J$ B. q+ t$ U, x/ K; ?, v(@sum(links(i,j):a(j)*x(i,j))-@sum(shengshi(i):
0 M3 I8 B3 d$ g" X1 r" \@sum(links(i,j):a(j)*x(i,j)))/13)^2)/12)^(1/2);; j) |, Q, k+ P7 D- d6 O* U
& N3 Z- a7 v. p8 m1 M8 K: s( P
@sum(zhibiao(j):a(j)^2)=1;$ s; m' o% R3 d" ~! D; }8 y5 J
@for(zhibiao(j):a(j)>-1);
9 P8 W! U* c4 `) i" \% q2 b: Z@for(zhibiao(j):a(j)<1);3 m8 e+ O" x- w6 V
@for(zhibiao(j)free(a(j)));. g9 y% ~( G; R
end ) v, c9 \! A* Q; R! @7 I% r" }
求不出全局最优解，请大侠帮帮忙！！！
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.
    9 [* ^0 {* o/ |' X7 B. e, P( [1 x* L$ |  Objective value:                             0.7272518
    7 `. D# f* w5 K, q6 b  Objective bound:                             0.7272527
    , _- x6 k2 m0 P0 G2 A' P  Infeasibilities:                              0.000000
    % a3 K1 T8 {0 ~5 v8 s$ M: G  Extended solver steps:                           38843
    & e* ~8 h6 ]; X1 p: S  Total solver iterations:                       75613473 K3 t8 H6 k6 O
    1 W. f( R0 Q4 \9 w3 W. r- ~2 K2 i
    & X- T) j3 P0 c/ t& B/ p9 p# ]
                           Variable           Value        Reduced Cost
    - g# Z- ]% D# x7 z* u                          A( 1)       0.1248020            0.000000
    / |" q3 M6 J+ _. ]7 o* v                          A( 2)       0.1216772            0.000000
      l# o! d" p  ]) H6 K                          A( 3)       0.1621721            0.000000$ K% L5 q* N+ A2 D7 q
                              A( 4)       0.6318675            0.000000: R  u! z  S+ ~
                              A( 5)       0.3868056            0.000000) X. m% u/ N- n* ~
                              A( 6)       0.6280479            0.000000
    6 e4 L! w1 [( Y                       X( 1, 1)       0.9959900E-01        0.000000
    $ X0 `8 V9 P6 r                       X( 1, 2)       0.1040950            0.000000
    7 u0 U  W3 r- m5 ]8 X' V                       X( 1, 3)       0.2118600E-01        0.000000
    ' j/ d3 L- t! ^# i1 f1 }. L. }                       X( 1, 4)       0.4898490            0.000000
    ! I/ J% G3 M+ ^; g8 s2 D                       X( 1, 5)       0.3141900E-01        0.000000
    + q' C; q4 H& c- @  G% t                       X( 1, 6)       0.3996010            0.000000
    5 X+ i. w* |, _; i9 N                       X( 2, 1)       0.1222200E-01        0.000000! o. k) |* p) `) s- }/ F
                           X( 2, 2)       0.6320000E-02        0.000000
    : |, W4 L6 \0 x+ B/ O5 t                       X( 2, 3)       0.3507500E-01        0.000000$ N6 T) }# }6 Q) ]% m
                           X( 2, 4)       0.5580430            0.000000
    , n7 @, n/ s8 x2 }+ G                       X( 2, 5)       0.2537760            0.000000
    * }2 Y& f7 }7 u2 v. N                       X( 2, 6)       0.4651700            0.000000
    9 Q/ ?" e( i$ U- A* Q                       X( 3, 1)       0.1626000E-01        0.000000
    ! ~5 G/ v2 Z* R$ j) G& |1 L                       X( 3, 2)       0.4213000E-02        0.000000
    ' ~) |) Y* t6 @) x# Z' L6 l                       X( 3, 3)       0.8442300E-01        0.000000
    % T2 i6 x5 ]3 F" r3 c                       X( 3, 4)       0.6314420            0.000000
    ; D& G4 ]7 A. ~; W/ u) @" Y. E" F                       X( 3, 5)       0.3751740            0.000000, K! X) m5 v: @, Y- }7 a
                           X( 3, 6)       0.5355790            0.000000% y) t6 x0 g( v3 A+ ~! U
                           X( 4, 1)       0.4587000E-02        0.000000
    1 {5 h( P# n$ J6 x$ y                       X( 4, 2)        0.000000            0.000000% c* g3 o" {: f; c! I: F2 c
                           X( 4, 3)       0.8615400E-01        0.000000
    - n# M  m" u8 u( w                       X( 4, 4)       0.5575220            0.000000
    % D4 e5 ?! _% @8 o0 }* d8 |; j4 R% @                       X( 4, 5)       0.2105910            0.0000002 Z8 b$ f( e( C$ ]2 A' C" C- i
                           X( 4, 6)       0.4647210            0.000000
    % F4 L+ ~0 w& c% B0 N( b# }                       X( 5, 1)       0.3988000E-01        0.000000
    0 N0 [2 Q* i" `2 U$ l                       X( 5, 2)       0.1533300E-01        0.0000003 W! t% O5 X6 b+ C" [' y
                           X( 5, 3)       0.7792100E-01        0.000000
    4 V% d: g  p$ {$ n                       X( 5, 4)        1.000000            0.000000
    ( }2 d- z5 t) v: B+ _* p1 n                       X( 5, 5)       0.4887600            0.000000
    / k+ u. c- H# e& O1 h2 d3 u( k                       X( 5, 6)       0.8885230            0.000000
    1 Y. F( J$ v4 P% G                       X( 6, 1)       0.3290000E-02        0.000000" ~9 r1 q. M6 N2 |, P$ q9 r4 \, X
                           X( 6, 2)       0.2399000E-02        0.000000$ Z+ o8 i7 {' d( D+ g! E! N1 z
                           X( 6, 3)       0.2344100E-01        0.0000005 p2 ~2 M0 G. o( j1 f" y6 X
                           X( 6, 4)       0.2998440            0.000000
    ( T, ~+ u, J( G4 N+ F5 v1 u1 b6 W                       X( 6, 5)        1.000000            0.0000004 W, c3 O# w9 C1 L3 ?7 N8 w
                           X( 6, 6)       0.2178140            0.000000
    * @: o+ s$ S& c+ g+ g5 R7 q                       X( 7, 1)       0.2175000E-01        0.000000
    ' D8 ~1 T9 A  p- ?                       X( 7, 2)       0.2329600E-01        0.000000% z! Z' \, B9 W
                           X( 7, 3)       0.5572700E-01        0.000000
    2 T  |+ j, ~& |                       X( 7, 4)       0.6350860            0.000000
    3 o7 K9 K' w. Q0 |/ D7 S                       X( 7, 5)       0.6575610            0.000000
    6 {" k$ X; a9 j5 h9 {/ l7 W! ^                       X( 7, 6)       0.5387230            0.000000/ o2 F) k, e& l7 _* C
                           X( 8, 1)        1.000000            0.000000* M, R' t. I, s
                           X( 8, 2)        1.000000            0.0000004 V% c1 i" W  q3 X7 o
                           X( 8, 3)        1.000000            0.000000
    " \! q9 ]; o. Q0 f                       X( 8, 4)        0.000000            0.000000
    6 G2 w# V7 x; ~                       X( 8, 5)        0.000000            0.000000
    1 U/ [& r* P; V                       X( 8, 6)       0.4305890            0.000000( g7 g: Y5 `  r3 P" i
                           X( 9, 1)       0.7096300E-01        0.000000
    8 g, n5 q4 o. V3 p                       X( 9, 2)       0.2775390            0.000000
    5 k* \/ R# k5 R  E* W& j                       X( 9, 3)       0.3859470            0.000000
    ) D4 z( r  ]; ~                       X( 9, 4)       0.1535660            0.0000009 m# M7 x+ S5 I8 f, ^4 ~
                           X( 9, 5)       0.1475700E-01        0.0000005 |: l. H6 E& n; K- }. x. W
                           X( 9, 6)        0.000000            0.000000
    9 v# w! m- W, H. @' T! _$ {. D                      X( 10, 1)       0.3816400E-01        0.0000006 a/ {& ?/ ?0 a4 T% p7 K
                          X( 10, 2)       0.3837300E-01        0.000000
    ) o2 A/ J0 ^1 b- ^2 _) o! }                      X( 10, 3)       0.2499000E-01        0.000000
    : `; X: f/ e& t8 L+ A( h                      X( 10, 4)       0.4898490            0.000000
    $ i# ?0 u/ p2 ]9 ~- d3 |                      X( 10, 5)       0.3115900E-01        0.000000/ ~9 n6 z  s6 I# u, f
                          X( 10, 6)        1.000000            0.000000
    4 @' q- S0 {2 t* p1 ]                      X( 11, 1)        0.000000            0.0000001 J( F4 ^3 h% H  W) W' U  R
                          X( 11, 2)       0.8735000E-02        0.0000004 \7 y( j. |* u2 M" P% E$ i
                          X( 11, 3)       0.1525900E-01        0.000000
    4 M* k, W# L- U: o" ?                      X( 11, 4)       0.2889120            0.0000002 e3 g5 |' M( @; Q: @( x) {
                          X( 11, 5)       0.3700900E-01        0.0000001 u3 X  M$ ~2 x5 t' b: d
                          X( 11, 6)       0.4091800E-01        0.000000  B2 F7 e1 O, p% z, K5 c! A/ h
                          X( 12, 1)       0.5718100E-01        0.000000- N0 r* B; r+ c
                          X( 12, 2)       0.3850800E-01        0.000000
    4 ~+ g% C# Y+ T  Z9 |, m8 h                      X( 12, 3)       0.2541090            0.0000007 f" ~5 j  ^/ @/ T
                          X( 12, 4)       0.1842790            0.000000
    : U: @* L1 [3 d8 f7 J( L! H                      X( 12, 5)       0.9316000E-01        0.000000
    7 ~/ F: D/ z9 ~2 ]6 i9 p* v                      X( 12, 6)       0.5224600E-01        0.000000% u6 @& @5 p6 l  V) L9 Z+ N
                          X( 13, 1)       0.4821100E-01        0.000000
    * W! K6 O' ]0 |( U9 r  ^9 B; x                      X( 13, 2)       0.9641000E-02        0.000000
    " ?+ v( p( P: i/ t                      X( 13, 3)        0.000000            0.000000* M2 S5 B2 C( d1 i
                          X( 13, 4)       0.3878190            0.000000
    5 Q: [5 c; b0 x& F$ y' O! _                      X( 13, 5)       0.1358340            0.000000- \. d7 i4 N) I: M! f
                          X( 13, 6)       0.1185130            0.0000003 k+ c* i( q0 {1 ?! r8 H, }/ H* M
    7 c* O/ N8 F3 m. Q2 ?% b
                                Row    Slack or Surplus      Dual Price
    2 I2 ?$ U; p. S7 C$ D' N% ~' c                              1       0.7272518            1.000000' [6 C3 j) Q5 N0 G7 P) x$ T8 Z
                                  2        0.000000           0.3636259
    ; E$ S8 s! g8 p+ K4 B                              3        1.124802            0.0000004 q6 ~+ @3 L9 ]7 k6 C
                                  4        1.121677            0.000000
    ! E7 @: B9 p! ~7 `                              5        1.162172            0.000000) j/ z/ K6 S9 E6 W$ {
                                  6        1.631867            0.000000
    7 ?$ ~. g8 g+ A/ J8 h" E! T- n                              7        1.386806            0.000000  s4 Q& y5 k* |4 G% P+ ]3 [5 T
                                  8        1.628048            0.000000
    - Q2 f' J) g( J                              9       0.8751980            0.000000( _7 T3 h( Y. e# a
                                 10       0.8783228            0.000000! F' F! R, [# _, R3 A6 H! f' g
                                 11       0.8378279            0.000000# G+ h! ~2 H2 y
                                 12       0.3681325            0.000000
    . l) [' U7 E  R9 C; d0 j7 y                             13       0.6131944            0.000000
    ! x. L- Y' I. x  N                             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
    但是 我设置全局求解时：两三秒就停止了& A  t7 s" b1 a* m4 [
    出现错误提示9 }: L6 }3 O+ A
    还有 上面的 全局的做大值还没有局部做大的大  为什么！！！
    回复

    使用道具 举报

    返回列表 发新帖
    高级模式
    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 23:37 , Processed in 0.632193 second(s), 102 queries .

    回顶部