快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

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

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

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

7

主题

5

听众

28

积分

升级  24.21%

该用户从未签到
跳转到指定楼层
1#
发表于 2009-9-5 07:49 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
model:
: d& W) j: g# Z" f1 Xsets:" m0 P0 m) W6 q/ |/ e! \
shengshi/1..13/;
/ Z, @9 z7 B, V2 jzhibiao/1..6/:a;9 `! _* ^* ], h+ t& R  x$ d, N
links(shengshi,zhibiao):x;! j* `( o' V0 W, N) Z
endsets
2 h; H' S' r; ^5 Udata:
* n. `3 b  o% A5 b; S. \: bx=* O5 |/ V. \; N8 y' b+ k; \) {
0.099599  0.104095  0.021186  0.489849  0.031419  0.399601 ( d, M0 I; ~% H* y/ E
0.012222  0.006320  0.035075  0.558043  0.253776  0.465170 , h' }; {3 Y* E8 d" ]2 o* E4 \
0.016260  0.004213  0.084423  0.631442  0.375174  0.535579 " d7 A8 i2 i. D9 Q5 m' ]* \. [0 v
0.004587  0.000000  0.086154  0.557522  0.210591  0.464721 2 ^; U/ x! V2 O4 q4 ?1 r/ h; Y
0.039880  0.015333  0.077921  1.000000  0.488760  0.888523
6 l0 s. Q$ |5 e& Q# m/ ?0.003290  0.002399  0.023441  0.299844  1.000000  0.217814
: L" U: g* r2 e9 @" Z$ W1 X0.021750  0.023296  0.055727  0.635086  0.657561  0.538723
/ h: M8 n) C+ [/ |1 m! }, b' x1.000000  1.000000  1.000000  0.000000  0.000000  0.430589 # z, K0 I; q! H5 u& H
0.070963  0.277539  0.385947  0.153566  0.014757  0.000000 - L4 X& C3 {( b5 w* R7 Q/ I; I
0.038164  0.038373  0.024990  0.489849  0.031159  1.000000 0 H1 d" l$ ^3 `7 Q! L( u9 R2 L
0.000000  0.008735  0.015259  0.288912  0.037009  0.040918
7 E+ x1 c5 b% U! F- ?0.057181  0.038508  0.254109  0.184279  0.093160  0.052246
& R& L, l3 E. _0 j1 Y0.048211  0.009641  0.000000  0.387819  0.135834  0.118513;
0 m+ `- @, s2 b- i: b4 N# kenddata
2 }/ r  T1 A- smax=. z1 L, j  g8 `& C+ x
(@sum(shengshi(i):) e) n% i1 i8 ?$ a& z9 G2 |' C
(@sum(links(i,j):a(j)*x(i,j))-@sum(shengshi(i):1 y& m4 z' T9 o+ O/ b, Y
@sum(links(i,j):a(j)*x(i,j)))/13)^2)/12)^(1/2);+ W8 w8 v8 Q  b- i) G

9 Q8 _" G5 B" ?- P- j8 X4 f@sum(zhibiao(j):a(j)^2)=1;
0 _. N1 ^, H7 P- T@for(zhibiao(j):a(j)>-1);
! y0 \6 j3 S' g4 M1 ~3 v- H9 Q@for(zhibiao(j):a(j)<1);) B( ^& W# t& K/ q) U
@for(zhibiao(j)free(a(j)));
2 T% ?! H: U) @end / e5 O& z) @) ]/ S1 E/ G
求不出全局最优解，请大侠帮帮忙！！！
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.! l3 R5 r7 o+ }) T9 j/ R+ k/ t0 h
      Objective value:                             0.7272518) p+ V7 j/ |' t1 N
      Objective bound:                             0.7272527
    2 ?/ C& U5 H1 i6 b2 R! j: \  Infeasibilities:                              0.000000" X* B+ Q% P% y$ r6 O" p
      Extended solver steps:                           38843
    ) ~4 s- N4 ^2 z9 \5 n  Total solver iterations:                       7561347
    ) p4 {, {5 H& x6 C; ]/ [) P1 o
    # X6 A& c6 ]9 \2 f" x" `4 V7 M% a+ F9 S+ P7 P
                           Variable           Value        Reduced Cost- D# @4 X: J! _* e9 s/ p% b% _7 a
                              A( 1)       0.1248020            0.000000& g/ G5 ~5 q: U  F) m
                              A( 2)       0.1216772            0.000000- y6 i1 B1 b+ y& @6 z( g5 P
                              A( 3)       0.1621721            0.0000002 v  \4 Q  U7 E
                              A( 4)       0.6318675            0.000000" _# g6 j+ j1 M0 c% [5 ~9 J
                              A( 5)       0.3868056            0.000000
    4 {1 j3 U2 |% O0 M  W                          A( 6)       0.6280479            0.000000
    . ^$ c$ v9 D  f6 p# |1 q+ l1 f- u9 h                       X( 1, 1)       0.9959900E-01        0.000000
    & f- x+ _, Z/ @$ G3 @                       X( 1, 2)       0.1040950            0.000000! d$ g0 ~" H( R, h7 h% A8 u' h
                           X( 1, 3)       0.2118600E-01        0.000000+ j( x9 E& h; i8 @0 Y! b, f# g' I+ M
                           X( 1, 4)       0.4898490            0.000000
    + X$ I9 |; M3 }) f  O                       X( 1, 5)       0.3141900E-01        0.000000+ q& `, p4 e' N1 o" w) u0 U4 |
                           X( 1, 6)       0.3996010            0.000000
    5 A3 K; n+ p( ?$ G% q7 R9 y                       X( 2, 1)       0.1222200E-01        0.0000007 n7 L+ N* P3 ^
                           X( 2, 2)       0.6320000E-02        0.000000
    : f! X4 y3 {+ T0 k, L% {) j                       X( 2, 3)       0.3507500E-01        0.0000006 f8 X) [: i# x( R
                           X( 2, 4)       0.5580430            0.000000% x; H. I% }, g
                           X( 2, 5)       0.2537760            0.000000: e# M1 U* _5 D- W# Y, |" w
                           X( 2, 6)       0.4651700            0.000000  @3 {4 n1 ], I
                           X( 3, 1)       0.1626000E-01        0.000000+ X4 c8 J0 ?/ Z
                           X( 3, 2)       0.4213000E-02        0.000000
    % R. [) E# j1 ~" y                       X( 3, 3)       0.8442300E-01        0.000000
    , ~- P) G  Z% |  _- \. d                       X( 3, 4)       0.6314420            0.0000009 F6 d! t( |& v; i3 y2 j. t, Q
                           X( 3, 5)       0.3751740            0.000000
    ( Y" p1 V7 G% f; V2 @! O1 b% r) {                       X( 3, 6)       0.5355790            0.000000# T1 y3 ]2 D1 |: u3 C
                           X( 4, 1)       0.4587000E-02        0.000000
    ; R6 ~, X% V+ J3 R: x9 i0 T* v4 c! e                       X( 4, 2)        0.000000            0.000000
    0 |5 N; ^& N! S. b                       X( 4, 3)       0.8615400E-01        0.000000/ h: D+ L  [. O) G
                           X( 4, 4)       0.5575220            0.000000
    $ W0 F8 Z0 N' N9 T7 j# j8 a/ K( f                       X( 4, 5)       0.2105910            0.000000
    4 ^2 D# j1 \3 \* Q# {+ r5 r                       X( 4, 6)       0.4647210            0.000000% |8 }; m. v9 K5 ?3 m7 t! q6 V
                           X( 5, 1)       0.3988000E-01        0.0000000 k1 D: b" B& J5 B6 Q* z
                           X( 5, 2)       0.1533300E-01        0.000000
    " p) V$ s0 M, K                       X( 5, 3)       0.7792100E-01        0.000000
    $ U. l5 d3 X6 E! S& @/ Y                       X( 5, 4)        1.000000            0.000000
    * P; ]3 S% U% }. f% X                       X( 5, 5)       0.4887600            0.000000# Q4 S( x0 G+ \; P" K1 t  Y8 r
                           X( 5, 6)       0.8885230            0.000000
    & D4 Y' I1 ]5 W& D                       X( 6, 1)       0.3290000E-02        0.000000" P# L7 H0 O4 T- r/ c
                           X( 6, 2)       0.2399000E-02        0.000000
    & S* m" t$ ~2 c5 f                       X( 6, 3)       0.2344100E-01        0.000000$ O# y; `7 T5 h' Q
                           X( 6, 4)       0.2998440            0.0000000 |* m( s7 O$ ?! k4 }/ i4 P
                           X( 6, 5)        1.000000            0.000000
    . {& a$ ~. j8 Z+ b, G                       X( 6, 6)       0.2178140            0.000000
    8 P# N/ d, c8 X- Q                       X( 7, 1)       0.2175000E-01        0.000000
    . T( }4 g* n* W* I                       X( 7, 2)       0.2329600E-01        0.000000
    ' d9 g. |  F% f. h4 a7 Y                       X( 7, 3)       0.5572700E-01        0.000000+ g; t% f* M$ ^# C. X
                           X( 7, 4)       0.6350860            0.000000
    * W* R. `( ~+ h7 D                       X( 7, 5)       0.6575610            0.000000
    6 [, ]$ f7 `6 d. h% U                       X( 7, 6)       0.5387230            0.000000
    ' d! q/ ~6 R9 X6 x( _0 l$ F' I) U                       X( 8, 1)        1.000000            0.000000
    + n8 F" Y9 W8 ]5 h/ h2 g                       X( 8, 2)        1.000000            0.000000
    & K  m1 {4 ~. K4 `: D1 `! Q                       X( 8, 3)        1.000000            0.0000003 ?1 N! _  E. g: V4 i2 A
                           X( 8, 4)        0.000000            0.000000
    ! o% j  |8 K# d                       X( 8, 5)        0.000000            0.0000000 I8 k' ?# c  `4 U7 r
                           X( 8, 6)       0.4305890            0.000000
    3 f$ c' y  [. X8 J' ^                       X( 9, 1)       0.7096300E-01        0.000000! O% P* ~$ m" D6 x
                           X( 9, 2)       0.2775390            0.000000
    2 T/ S8 }* C1 R' X! \                       X( 9, 3)       0.3859470            0.000000, M* G5 E* e  m8 h+ ~* w9 Y
                           X( 9, 4)       0.1535660            0.000000
    ' |8 |3 M/ }5 S                       X( 9, 5)       0.1475700E-01        0.000000
    ( o& e$ H, a7 g5 H! n' D                       X( 9, 6)        0.000000            0.000000
    " d2 U. f$ V+ l, ~0 P; {+ M% b                      X( 10, 1)       0.3816400E-01        0.000000
    + A! |4 v7 g. \0 ~# a6 w3 [                      X( 10, 2)       0.3837300E-01        0.000000
    ! i( R6 N5 Z+ K: R                      X( 10, 3)       0.2499000E-01        0.000000
    5 ]* g7 v) C/ |; Y2 X                      X( 10, 4)       0.4898490            0.000000
    8 F3 G1 k! a6 G9 @5 k                      X( 10, 5)       0.3115900E-01        0.000000
    6 E9 P( g( {4 T! x9 x. t5 A2 b                      X( 10, 6)        1.000000            0.0000008 ^7 s) w- I4 D# ^' m- q
                          X( 11, 1)        0.000000            0.000000
    # d" l$ ]& ], f% |* N                      X( 11, 2)       0.8735000E-02        0.0000007 D: t0 E! {" `6 f7 L# T3 E- G# X
                          X( 11, 3)       0.1525900E-01        0.000000' B2 X2 G1 I1 b0 ?& D1 n% E4 q( ?
                          X( 11, 4)       0.2889120            0.000000
    % T* X: e! X" D0 k                      X( 11, 5)       0.3700900E-01        0.000000
    " l7 t9 X% L0 j( A: m' @                      X( 11, 6)       0.4091800E-01        0.000000
    # _& A  N: q- V( T( h0 j                      X( 12, 1)       0.5718100E-01        0.0000008 b# k# ]' t' |% o$ Y
                          X( 12, 2)       0.3850800E-01        0.000000
    9 d( M9 [1 s9 U  y                      X( 12, 3)       0.2541090            0.000000
    ' `4 d2 z3 \* J4 n8 B: ~                      X( 12, 4)       0.1842790            0.000000$ E/ t- \6 ]( h2 }! e
                          X( 12, 5)       0.9316000E-01        0.000000
    9 z. h5 ]( Z9 p1 e7 U                      X( 12, 6)       0.5224600E-01        0.000000! W* U8 N7 x0 s8 R3 |
                          X( 13, 1)       0.4821100E-01        0.000000) K! {& ^, a( c: Q: p& W  V
                          X( 13, 2)       0.9641000E-02        0.0000003 j. _, I8 a1 l3 X. j6 s6 v; ^
                          X( 13, 3)        0.000000            0.000000  Q4 c5 z" u7 H. |
                          X( 13, 4)       0.3878190            0.0000003 T. V( i: g! }& ~: F
                          X( 13, 5)       0.1358340            0.000000* R( Q2 B( r8 K& z! |
                          X( 13, 6)       0.1185130            0.000000
    ' \9 }1 a- f5 S# K
    6 b$ u8 ]/ ]8 |& f" j                            Row    Slack or Surplus      Dual Price
    ' s: s9 [1 t  z                              1       0.7272518            1.000000, }4 J, `2 y( O/ j7 _1 A$ D: F
                                  2        0.000000           0.3636259+ r  V! D: ?! o' Q# r3 ]- a
                                  3        1.124802            0.000000
    % c& g7 S6 e4 B! O& @                              4        1.121677            0.000000
    & f) U) [& I; N( M                              5        1.162172            0.000000
    % f& F9 Y, M: p1 ]& u7 E5 X/ {                              6        1.631867            0.000000
    % D- ~; Z3 p$ s% I7 u6 I: p                              7        1.386806            0.000000
    + D. V9 N8 s2 ~                              8        1.628048            0.000000
    - Z; {# x* u0 K/ A0 o# k: b9 p                              9       0.8751980            0.000000
    ) S3 h$ N0 X) I: O5 ^                             10       0.8783228            0.000000
    , a3 O, l7 Q6 b" T! W. K( B  _                             11       0.8378279            0.000000
    ! h4 {7 k' Y, M/ N7 f1 e: M7 o) J  J                             12       0.3681325            0.000000: _* t. B' G" u2 }& V
                                 13       0.6131944            0.000000
    " Y4 J& v) n) U) t2 I* p. I                             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
    但是 我设置全局求解时：两三秒就停止了
    9 z( s. B) _5 b1 g: C出现错误提示
    5 p/ p& s- c5 \1 x还有 上面的 全局的做大值还没有局部做大的大  为什么！！！
    回复

    使用道具 举报

    返回列表 发新帖
    高级模式
    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, 2025-9-17 08:47 , Processed in 0.758452 second(s), 102 queries .

    回顶部