快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

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

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

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

7

主题

5

听众

28

积分

升级  24.21%

该用户从未签到
跳转到指定楼层
1#
发表于 2009-9-5 07:49 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
model:
: [: h) l) u" e" |; V& ~sets:
9 P, A* Q  q5 R" pshengshi/1..13/;
" r& P0 d; q# x& n% ]% \zhibiao/1..6/:a;
+ P( \) d) v& a" U& l% G" S4 |& Alinks(shengshi,zhibiao):x;
$ S6 r. m  \  o* W% @endsets& P3 i; m; c0 P* Z
data:
! E. i" ]( T7 W; tx=
1 a) s# u1 I% k0 i+ R0.099599  0.104095  0.021186  0.489849  0.031419  0.399601 8 B9 j/ ~3 S# b( i6 Z  R
0.012222  0.006320  0.035075  0.558043  0.253776  0.465170 # }% ?5 [9 v1 x2 T+ E2 q
0.016260  0.004213  0.084423  0.631442  0.375174  0.535579
% Q6 X8 }7 c$ T0.004587  0.000000  0.086154  0.557522  0.210591  0.464721 , P6 J) h  t' y6 L# P7 ?4 @
0.039880  0.015333  0.077921  1.000000  0.488760  0.888523
0 c+ V6 @! _: Y) I. }- k0 Q0.003290  0.002399  0.023441  0.299844  1.000000  0.217814 5 `: m  N, V( ?# x* l+ u
0.021750  0.023296  0.055727  0.635086  0.657561  0.538723
$ Q9 P- ^3 e/ p4 ^0 M1.000000  1.000000  1.000000  0.000000  0.000000  0.430589 ' ~( }; q% w; y: m3 m
0.070963  0.277539  0.385947  0.153566  0.014757  0.000000 6 x; [# }( w6 }8 @7 Z
0.038164  0.038373  0.024990  0.489849  0.031159  1.000000 . I, @% W1 [6 A0 P$ x
0.000000  0.008735  0.015259  0.288912  0.037009  0.040918 * q3 q) J) m" x% o3 {
0.057181  0.038508  0.254109  0.184279  0.093160  0.052246 * r! }/ S6 l4 V0 `% ~5 l
0.048211  0.009641  0.000000  0.387819  0.135834  0.118513;
$ {$ H: J! v! R$ h+ E, penddata
2 I1 @( T' |1 p  S, v* x7 G* umax=" e2 A5 j4 O1 D4 a
(@sum(shengshi(i):& e& ]+ _* g2 o
(@sum(links(i,j):a(j)*x(i,j))-@sum(shengshi(i):# c2 k5 V  R1 `
@sum(links(i,j):a(j)*x(i,j)))/13)^2)/12)^(1/2);) [: q0 `) k6 \- H! L# e2 t4 b

" D" H' G' W9 z+ C7 u1 f8 O@sum(zhibiao(j):a(j)^2)=1;6 e# }4 I7 _9 }0 o& N/ j. n
@for(zhibiao(j):a(j)>-1);
" i' c/ O. ~2 ^. I# U- g. g0 [@for(zhibiao(j):a(j)<1);
9 a4 b, {# K7 |' O4 \! H# e@for(zhibiao(j)free(a(j)));4 k& p( c  P% X3 y; f* A% ^
end
( h0 j3 O: I$ B- W7 m求不出全局最优解，请大侠帮帮忙！！！
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.
    6 {4 y$ `! S" V8 |4 Z  Objective value:                             0.7272518
    ; A/ M+ W: [# G( u# N- Y0 x  Objective bound:                             0.7272527
    ) F' H3 Y7 g# h! z  Infeasibilities:                              0.000000
    ! {) M5 p: |  ]+ Z  Extended solver steps:                           38843
      E. i' p4 \( S9 o0 j! Z  Total solver iterations:                       7561347  q, ?' d! f, b* Y$ C2 Y9 C7 @
    3 ]/ W  A9 e% P

    * |( W0 u& [1 V9 b                       Variable           Value        Reduced Cost
    9 n# r0 e- {  ]3 w' k                          A( 1)       0.1248020            0.000000
    $ B. u, ?  C. S$ W3 k) ]' \* r& [* p                          A( 2)       0.1216772            0.000000
    . T) b  P% ^/ L$ |                          A( 3)       0.1621721            0.000000
    : k9 D4 D" n8 j% q& [3 y                          A( 4)       0.6318675            0.000000
    $ O/ r( o1 M5 d                          A( 5)       0.3868056            0.000000, @8 t, g! i# X2 k) C9 B+ j6 b! K8 R
                              A( 6)       0.6280479            0.0000007 }/ U7 `2 J( [' x
                           X( 1, 1)       0.9959900E-01        0.000000
      j- o' Y) T& S. w" X. k$ F* A                       X( 1, 2)       0.1040950            0.000000
    . T6 ]. a( ?  d# C  H: N$ [  c                       X( 1, 3)       0.2118600E-01        0.000000# a0 X# \9 o" J% e
                           X( 1, 4)       0.4898490            0.000000. N3 l7 G) v" T5 l2 i# o4 J9 Q
                           X( 1, 5)       0.3141900E-01        0.000000" _5 i* k" P# M- ]& Y  R! A$ f; e
                           X( 1, 6)       0.3996010            0.000000
    5 T7 ]6 M! ^2 @# j; }4 V8 O                       X( 2, 1)       0.1222200E-01        0.000000; ?- j6 R4 G8 `1 q2 k
                           X( 2, 2)       0.6320000E-02        0.000000
    7 Q. s7 e. V" |. H2 L                       X( 2, 3)       0.3507500E-01        0.000000
    6 X7 h& K2 |; f/ e" S" ?: s* T                       X( 2, 4)       0.5580430            0.000000
    7 @$ c4 A% t+ d+ X2 Q2 |                       X( 2, 5)       0.2537760            0.000000
    ( T+ P3 \3 O' S* u1 t' ?                       X( 2, 6)       0.4651700            0.000000) c) A8 c7 S1 a4 p0 J
                           X( 3, 1)       0.1626000E-01        0.000000
    4 J. z) v) a0 A/ e9 _& N" J" A& M; X                       X( 3, 2)       0.4213000E-02        0.000000  K! k6 v9 \- p! h) I3 g
                           X( 3, 3)       0.8442300E-01        0.000000
    & @# [! G' D/ o# j3 t                       X( 3, 4)       0.6314420            0.000000% H# o& C7 L, Q$ Q* t: g
                           X( 3, 5)       0.3751740            0.000000* n$ L, E% g; k# _$ g4 O
                           X( 3, 6)       0.5355790            0.000000% Z4 u9 P4 B$ U: f: r4 W3 x% U
                           X( 4, 1)       0.4587000E-02        0.000000% Z  _, }# O2 B, C1 |! F# q
                           X( 4, 2)        0.000000            0.000000
    ! n! `6 {( i6 O4 r, o( p                       X( 4, 3)       0.8615400E-01        0.000000- O- N' h5 j& O% m7 \. a+ B
                           X( 4, 4)       0.5575220            0.000000
      ]$ A, g0 u, j: V                       X( 4, 5)       0.2105910            0.0000002 q0 e7 S+ U/ [  m! N" i# s
                           X( 4, 6)       0.4647210            0.000000
    " y% c& Z: z0 o2 U1 u                       X( 5, 1)       0.3988000E-01        0.000000
    5 z  n: V2 j  S% Y* S: Y                       X( 5, 2)       0.1533300E-01        0.000000
    6 @) }0 E* c  @( [                       X( 5, 3)       0.7792100E-01        0.000000
    / R# u7 C2 V  g3 R# U                       X( 5, 4)        1.000000            0.000000, a+ o* v+ w% ]0 u$ r
                           X( 5, 5)       0.4887600            0.000000( w; S* ~% A1 e9 u
                           X( 5, 6)       0.8885230            0.000000
    ) m9 i% M( a* q. j( ?                       X( 6, 1)       0.3290000E-02        0.000000
    ; K( d8 L8 {% ~                       X( 6, 2)       0.2399000E-02        0.000000
    ! K. I5 \5 `6 Y! J) p7 G                       X( 6, 3)       0.2344100E-01        0.000000
    ' A2 ^! Q- Q$ e2 ^# l  Z5 P                       X( 6, 4)       0.2998440            0.000000
    . R1 B. D' \9 s& {  z+ J* v1 \                       X( 6, 5)        1.000000            0.000000. ~' _$ N* c$ C- _3 T
                           X( 6, 6)       0.2178140            0.000000
    + A5 {. f* l/ V5 f: T9 q: E                       X( 7, 1)       0.2175000E-01        0.000000
    , G& u& {, L: ?( U% }                       X( 7, 2)       0.2329600E-01        0.000000
    # P6 R$ R2 h" t8 z                       X( 7, 3)       0.5572700E-01        0.000000
    : z8 V+ K# ^. F, w# r1 N( X# ^                       X( 7, 4)       0.6350860            0.000000
    ! n6 D# e5 _6 S- p" }: `2 r2 h                       X( 7, 5)       0.6575610            0.000000
    ; @" x1 y3 |* L  Y4 K7 q" b                       X( 7, 6)       0.5387230            0.000000+ ~+ t" i  E6 @% z
                           X( 8, 1)        1.000000            0.000000
    & |* W' K+ M8 q+ g( f* F2 x$ P  M- l                       X( 8, 2)        1.000000            0.000000" C  X  s0 U1 W+ e3 E
                           X( 8, 3)        1.000000            0.000000
    . d) @6 ]9 ^3 _+ G                       X( 8, 4)        0.000000            0.000000
    " R9 ?# w4 x$ S4 b                       X( 8, 5)        0.000000            0.0000006 e0 h5 l, I$ t) Y- s/ J6 K, z
                           X( 8, 6)       0.4305890            0.000000
    ( D* ^  n7 \# y9 l" Q' }1 o                       X( 9, 1)       0.7096300E-01        0.0000002 A6 H3 B3 C6 l" d9 \; }1 ^
                           X( 9, 2)       0.2775390            0.000000
    7 m7 K% z7 a7 ~" A( {9 o                       X( 9, 3)       0.3859470            0.0000006 n/ f3 B! c; \7 l1 C
                           X( 9, 4)       0.1535660            0.000000' @4 A5 w- O- q! x
                           X( 9, 5)       0.1475700E-01        0.000000  p8 k' J  d. d/ M1 [7 K
                           X( 9, 6)        0.000000            0.0000001 A' V8 F6 j( n. ]' E. a) J6 O
                          X( 10, 1)       0.3816400E-01        0.000000
    . G. n" b9 v. \                      X( 10, 2)       0.3837300E-01        0.000000
    $ B  _! n8 `- i9 ?0 ~% E6 h% r4 {6 O$ H                      X( 10, 3)       0.2499000E-01        0.0000003 ^( D8 C2 s9 K+ e
                          X( 10, 4)       0.4898490            0.000000
    / [5 z, P4 w9 z* X                      X( 10, 5)       0.3115900E-01        0.000000# o) c* P  Q! w2 w9 K; {
                          X( 10, 6)        1.000000            0.000000" Z: Z9 Z8 v& m4 K% C
                          X( 11, 1)        0.000000            0.000000' I. L$ ]0 A5 X
                          X( 11, 2)       0.8735000E-02        0.000000
    - H( _# j1 w; j% _                      X( 11, 3)       0.1525900E-01        0.000000
    ( A1 z7 |. h) _                      X( 11, 4)       0.2889120            0.000000: S& b5 e9 ?1 x) }1 y# Q+ D
                          X( 11, 5)       0.3700900E-01        0.000000
    7 x* L+ [- g  }1 k' ?. K                      X( 11, 6)       0.4091800E-01        0.000000
    " A* C  F5 b0 n% g                      X( 12, 1)       0.5718100E-01        0.000000, A/ b( q8 D5 }
                          X( 12, 2)       0.3850800E-01        0.000000/ O, {8 s1 C  g2 v
                          X( 12, 3)       0.2541090            0.000000" U5 p4 L! l+ D& V# ]
                          X( 12, 4)       0.1842790            0.0000005 C+ H/ D2 M" j  N+ T' }# a% R+ G: g
                          X( 12, 5)       0.9316000E-01        0.000000* D! [) C& w# t
                          X( 12, 6)       0.5224600E-01        0.000000  V0 k4 }# g) t7 e3 S+ K# A
                          X( 13, 1)       0.4821100E-01        0.000000! x  c2 L2 O7 s2 I  Q
                          X( 13, 2)       0.9641000E-02        0.000000' T3 \1 p) r  X" H4 ?) W
                          X( 13, 3)        0.000000            0.000000# V# [! b0 d8 x. e
                          X( 13, 4)       0.3878190            0.0000001 ~+ T4 ~; l7 \
                          X( 13, 5)       0.1358340            0.000000& |  a% b7 r- H9 q5 g: }- m
                          X( 13, 6)       0.1185130            0.000000& Y# `, H2 L4 g2 x  p0 n* W; {1 A6 t
    & y, _2 Y' e3 ?+ t) N/ S- X
                                Row    Slack or Surplus      Dual Price
    : J/ p+ Y1 y9 D! B4 C6 z                              1       0.7272518            1.000000
    # S, I2 Q/ Y- a3 v* _( c                              2        0.000000           0.3636259: k4 _8 k- s3 _  _& w( R
                                  3        1.124802            0.0000001 j0 h2 s; ?) I8 c, h5 Z! r
                                  4        1.121677            0.0000002 g2 t0 W% [0 ?; e- C  Z
                                  5        1.162172            0.0000006 h9 n8 E  `: ?1 L- x
                                  6        1.631867            0.000000/ o$ p1 v/ D- n) ?4 l
                                  7        1.386806            0.000000
    6 r9 a4 a' ]. \/ z                              8        1.628048            0.000000
    + G$ Q$ _8 T- T                              9       0.8751980            0.000000
    6 r& Y8 \* X, s2 {2 }& K% g1 m                             10       0.8783228            0.0000001 y( c; K9 _0 S, M
                                 11       0.8378279            0.000000
    8 L: {3 x, }7 u                             12       0.3681325            0.000000- O* L7 O  E& W
                                 13       0.6131944            0.000000: \( T0 Y! |2 }1 W  D) B% n* c
                                 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
    但是 我设置全局求解时：两三秒就停止了
    7 I7 W. M. S2 l) h4 }, _0 k; k+ x出现错误提示4 y/ ?8 S! E0 H0 M
    还有 上面的 全局的做大值还没有局部做大的大  为什么！！！
    回复

    使用道具 举报

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

    回顶部