快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

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

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

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

7

主题

5

听众

28

积分

升级  24.21%

该用户从未签到
跳转到指定楼层
1#
发表于 2009-9-5 07:49 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
model:
& X$ {. h6 t' C  {sets:& S! s+ \% b  p# p
shengshi/1..13/;5 [! l0 X/ h3 X2 B
zhibiao/1..6/:a;- Z/ D1 g) t! w
links(shengshi,zhibiao):x;
  V) x9 ?# i8 A9 B+ fendsets
5 Z. B3 V" l6 R1 mdata:4 X/ e* A/ Y* ?; ]3 G# @, U. C
x=
3 i: @( i9 R' {# D/ u0 F. [0.099599  0.104095  0.021186  0.489849  0.031419  0.399601 * L7 S+ e1 \5 g9 b# b6 {' q% v
0.012222  0.006320  0.035075  0.558043  0.253776  0.465170 0 u, [% F! l2 a6 c, _+ C
0.016260  0.004213  0.084423  0.631442  0.375174  0.535579 6 y# F7 g7 a/ e5 o
0.004587  0.000000  0.086154  0.557522  0.210591  0.464721 ! q  l8 o9 ^% W% U
0.039880  0.015333  0.077921  1.000000  0.488760  0.888523
9 Z. Y* d# {4 i" i0.003290  0.002399  0.023441  0.299844  1.000000  0.217814
* O0 Q; w- U, i( N0.021750  0.023296  0.055727  0.635086  0.657561  0.538723 ) M8 E3 K9 F1 k1 w
1.000000  1.000000  1.000000  0.000000  0.000000  0.430589
# E2 D- Z# e4 x8 L  \( O: t3 N! ?0.070963  0.277539  0.385947  0.153566  0.014757  0.000000
7 I7 A! S9 C( c- I* E6 Y0.038164  0.038373  0.024990  0.489849  0.031159  1.000000
! {4 W8 M3 `7 R0 d  {0.000000  0.008735  0.015259  0.288912  0.037009  0.040918
1 R! a- x1 N0 K: ]0.057181  0.038508  0.254109  0.184279  0.093160  0.052246 5 Y! ^" d5 s" o9 ]+ y$ k3 j9 p
0.048211  0.009641  0.000000  0.387819  0.135834  0.118513;) n% z1 l  u4 P3 d/ m7 T" `
enddata) r; W+ y) T. n5 L9 d* |
max=/ |) R# i2 V' q8 E; B0 e- i
(@sum(shengshi(i):" w" D$ J5 M5 Q- g/ D$ I
(@sum(links(i,j):a(j)*x(i,j))-@sum(shengshi(i):
; z( [3 l! p- }9 Q4 T@sum(links(i,j):a(j)*x(i,j)))/13)^2)/12)^(1/2);/ Z* E9 w' ~9 P, O. H* J' ?& L

3 [% D4 ?4 Z: O& Q0 \. _4 U@sum(zhibiao(j):a(j)^2)=1;/ D6 o" H& h) U2 J& @; v6 ^# V: h/ h
@for(zhibiao(j):a(j)>-1);
2 D: F$ Q8 C  S@for(zhibiao(j):a(j)<1);- e1 W+ K% d2 W! b3 I
@for(zhibiao(j)free(a(j)));; w) i  p9 |8 |
end ! D6 O6 i) [( C' z( G$ v
求不出全局最优解，请大侠帮帮忙！！！
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.
    5 e3 q0 ^- ?& E0 ~( \* @! z) c  Objective value:                             0.7272518
    . Y, L, U# |+ R  Objective bound:                             0.72725271 f7 |/ [1 g) A5 Z
      Infeasibilities:                              0.000000
      u+ @; G# {& q+ t5 h  Extended solver steps:                           38843" D# S' T- L# B. ^0 Q. @8 [' @% ^
      Total solver iterations:                       7561347: p5 i. j8 d8 E- o

    " n+ J8 L3 `+ M; |+ D7 Q9 i* E: e5 G
                           Variable           Value        Reduced Cost% j3 r% y1 D% O, Z4 T0 H
                              A( 1)       0.1248020            0.000000
    + q  l5 s% Z" Y9 \! h                          A( 2)       0.1216772            0.000000$ v9 U. V9 `( v7 Y6 z' x
                              A( 3)       0.1621721            0.000000
      H' W5 D" P: U1 X: q& r& j                          A( 4)       0.6318675            0.000000
    & y9 q; ~2 }7 A# N& U* d4 d                          A( 5)       0.3868056            0.000000( w3 Y, N2 w/ o8 b6 P1 a
                              A( 6)       0.6280479            0.0000004 f5 i  b9 p  N( M! k1 J6 D
                           X( 1, 1)       0.9959900E-01        0.000000
    , [7 ?4 P. v. B$ y  ?. P1 s                       X( 1, 2)       0.1040950            0.0000000 n& X6 f1 V3 Z5 u7 a
                           X( 1, 3)       0.2118600E-01        0.000000- ^6 R  Q9 T+ j6 Z- \
                           X( 1, 4)       0.4898490            0.0000009 t1 a' |) y2 i
                           X( 1, 5)       0.3141900E-01        0.000000$ Y+ ?& S- N. i, x6 S2 i* R
                           X( 1, 6)       0.3996010            0.000000% {9 s) N& A2 h
                           X( 2, 1)       0.1222200E-01        0.000000
    8 C7 U* B1 }+ q" }+ t; U4 P6 B  p                       X( 2, 2)       0.6320000E-02        0.000000$ U6 ~( U3 e3 E" _+ L+ E, E
                           X( 2, 3)       0.3507500E-01        0.0000002 m" L, }7 Z: M/ L. {! b
                           X( 2, 4)       0.5580430            0.000000
    2 ?0 k* O5 ^3 w! B* m                       X( 2, 5)       0.2537760            0.000000& ]7 }& B5 B( `5 t' d8 ^& s' [
                           X( 2, 6)       0.4651700            0.000000
    - F9 P8 E" q' |5 `+ @) w                       X( 3, 1)       0.1626000E-01        0.000000
    9 k/ |+ g, {" i* U, k                       X( 3, 2)       0.4213000E-02        0.000000  `2 K, k4 t6 D8 G0 i
                           X( 3, 3)       0.8442300E-01        0.0000005 b2 O: ~" h; w) l5 m+ f: O6 N) F
                           X( 3, 4)       0.6314420            0.000000+ M, ^! \; H1 L5 h
                           X( 3, 5)       0.3751740            0.000000
    ! @" v& q+ v# q3 `% j                       X( 3, 6)       0.5355790            0.0000008 t4 \: c1 [; Y
                           X( 4, 1)       0.4587000E-02        0.000000
    . f5 i* @8 h2 }% {) b7 N, `                       X( 4, 2)        0.000000            0.0000008 s5 g1 l- v' o5 J! |& A
                           X( 4, 3)       0.8615400E-01        0.0000009 M4 J: `" B$ e0 @& ]1 a9 ^. y
                           X( 4, 4)       0.5575220            0.000000. k5 w1 N- Z2 L$ }! _0 |3 e. L8 j5 u
                           X( 4, 5)       0.2105910            0.000000
    ) J# T2 b; \) x; f9 ?! s' H                       X( 4, 6)       0.4647210            0.000000+ s8 @' U# r9 Y2 a3 V% |
                           X( 5, 1)       0.3988000E-01        0.000000! e0 t! [. {$ v) k: a6 o: h
                           X( 5, 2)       0.1533300E-01        0.000000
    5 i0 `/ U+ j. k! X) ?4 F* Q# A                       X( 5, 3)       0.7792100E-01        0.000000) D8 L5 a/ a% j' L9 X( Z+ P
                           X( 5, 4)        1.000000            0.0000002 p) o2 p/ F% Q9 `# o. ?1 h# f
                           X( 5, 5)       0.4887600            0.000000! v; b( m% r; T& g1 ~& r2 M3 Q
                           X( 5, 6)       0.8885230            0.000000* ^' b5 H. D& M% }" `$ r2 V2 u9 _
                           X( 6, 1)       0.3290000E-02        0.000000
    2 u8 h: O+ y( ]7 Q5 I( _                       X( 6, 2)       0.2399000E-02        0.000000
    + Q3 B+ R5 a/ m2 D# |                       X( 6, 3)       0.2344100E-01        0.0000001 V2 Q+ R- P: N' ^$ U5 z
                           X( 6, 4)       0.2998440            0.0000009 n+ B$ M9 K! ]+ z4 O. ^" L
                           X( 6, 5)        1.000000            0.000000
    / X+ S, M+ I5 e% }2 h. s                       X( 6, 6)       0.2178140            0.000000
    7 J8 C7 z( U* C# n" `                       X( 7, 1)       0.2175000E-01        0.000000
    9 K. ~) j" V( v1 \4 T6 w$ {& i                       X( 7, 2)       0.2329600E-01        0.000000# V2 j; R6 k7 k
                           X( 7, 3)       0.5572700E-01        0.0000002 O; d+ ^5 J) X3 R. t0 H* H0 Y
                           X( 7, 4)       0.6350860            0.000000
    9 R. k( j6 W) R" T. G) G7 @                       X( 7, 5)       0.6575610            0.0000008 n7 N* ^; F$ e8 l  b9 S: Z7 x
                           X( 7, 6)       0.5387230            0.0000009 X0 A. }% I' @9 h# `
                           X( 8, 1)        1.000000            0.0000006 U- Q9 i) p  [) ?* X
                           X( 8, 2)        1.000000            0.000000
    4 b# e+ ~* o5 z. G0 g6 p/ r- l3 M                       X( 8, 3)        1.000000            0.000000# R! E4 ?4 t/ N$ ~2 C$ v- H
                           X( 8, 4)        0.000000            0.000000
    1 |" e; i. V  e6 O                       X( 8, 5)        0.000000            0.000000) s$ x5 N2 o0 V/ k  I; D, T
                           X( 8, 6)       0.4305890            0.000000( k9 B1 x2 T6 ?
                           X( 9, 1)       0.7096300E-01        0.000000: m+ D" X  \' V8 H3 z1 T
                           X( 9, 2)       0.2775390            0.000000
    , R: P6 w! E6 D2 r- R: I0 j                       X( 9, 3)       0.3859470            0.0000008 q- |( C/ F2 I' q
                           X( 9, 4)       0.1535660            0.000000
    3 P7 k4 i0 W3 Z. b1 d1 j                       X( 9, 5)       0.1475700E-01        0.000000
    5 b+ w( P; ?5 O9 V7 R6 ~& o; u' ~                       X( 9, 6)        0.000000            0.000000. r5 W% w. {+ h" Z6 r/ P# z1 v
                          X( 10, 1)       0.3816400E-01        0.0000001 ]- X+ T' Z' n' a" v5 Q3 \
                          X( 10, 2)       0.3837300E-01        0.0000009 b) ^$ k$ z' {0 m8 t& h! b
                          X( 10, 3)       0.2499000E-01        0.000000
    & n$ u9 t5 ^9 k$ c3 \1 a                      X( 10, 4)       0.4898490            0.000000: z; K" ]* f3 l, c. b% w9 `
                          X( 10, 5)       0.3115900E-01        0.000000
    5 M! y5 @; `$ W% R. {                      X( 10, 6)        1.000000            0.000000% x: B; w9 V* B3 x$ t
                          X( 11, 1)        0.000000            0.000000
    . z( r* U' S. e* s( I                      X( 11, 2)       0.8735000E-02        0.0000000 J. \7 ^5 V& b9 R' ]: u; T/ w+ n
                          X( 11, 3)       0.1525900E-01        0.000000
    ! R2 U8 w* U* t                      X( 11, 4)       0.2889120            0.000000
    . q2 t3 g( P3 o( N                      X( 11, 5)       0.3700900E-01        0.0000009 H/ V0 y, [- A. |
                          X( 11, 6)       0.4091800E-01        0.0000009 Q, ~! B3 J) d! G: M
                          X( 12, 1)       0.5718100E-01        0.000000
    , y6 B' ?8 F5 M4 ^& J# ~! [                      X( 12, 2)       0.3850800E-01        0.000000
    4 f8 {/ c- D$ ]" B3 i* H! y                      X( 12, 3)       0.2541090            0.000000
    3 M* I$ c- A' R  |                      X( 12, 4)       0.1842790            0.000000
    ; k- u& N" |3 M                      X( 12, 5)       0.9316000E-01        0.000000
    * t( p& r; _4 w1 I2 R) t                      X( 12, 6)       0.5224600E-01        0.000000
    : J# S5 E8 ~  E! U8 Q                      X( 13, 1)       0.4821100E-01        0.000000
    5 B* `7 \3 C+ K# g/ @5 Q) H                      X( 13, 2)       0.9641000E-02        0.000000( w1 k. @4 S8 ?' M  x
                          X( 13, 3)        0.000000            0.000000
    ) ]' a, t: u4 z5 t' A                      X( 13, 4)       0.3878190            0.000000. J3 G1 Y- ~7 q( ~  [
                          X( 13, 5)       0.1358340            0.0000009 A6 F9 x( b# \5 K1 k  y; @
                          X( 13, 6)       0.1185130            0.000000
    ' B8 v+ X6 n5 X! }4 o6 ]: `2 y9 N8 c" S
                                Row    Slack or Surplus      Dual Price( _2 n7 C! ~- Z/ f0 J. L
                                  1       0.7272518            1.000000. N- {- `% T/ F2 v
                                  2        0.000000           0.3636259; G: d9 K8 A6 t7 i  b
                                  3        1.124802            0.0000002 w) H  S; i* T9 L
                                  4        1.121677            0.0000003 V4 ~7 N& \% Q3 i
                                  5        1.162172            0.000000; W  \  }) U* D6 m2 b% |# x
                                  6        1.631867            0.000000; q8 }! k. @6 I9 o
                                  7        1.386806            0.000000$ F& @) s  F/ M3 K
                                  8        1.628048            0.000000
    1 p9 K+ e) v6 ?7 z                              9       0.8751980            0.000000
    - V* s: N4 P9 X% T                             10       0.8783228            0.000000( t) C  X) S# h- V$ F
                                 11       0.8378279            0.0000008 b- e+ o# B( `1 o4 c9 z0 _3 W  X
                                 12       0.3681325            0.000000+ a8 L4 x& Z' k% \1 ^* Q" _
                                 13       0.6131944            0.000000
    * w# Y! M/ ^, ]) [% \( j6 a                             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
    但是 我设置全局求解时：两三秒就停止了: S7 J% [4 x* H  i
    出现错误提示
    3 C4 j- J/ g: A5 r; J) i还有 上面的 全局的做大值还没有局部做大的大  为什么！！！
    回复

    使用道具 举报

    返回列表 发新帖
    高级模式
    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-8 05:38 , Processed in 0.932422 second(s), 103 queries .

    回顶部