快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

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

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

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

7

主题

5

听众

28

积分

升级  24.21%

该用户从未签到
跳转到指定楼层
1#
发表于 2009-9-5 07:49 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
model:2 m/ v+ Z/ O- E2 L; q& i8 u
sets:' W& r9 P* |  y6 I! d
shengshi/1..13/;9 i1 P/ U: R4 S
zhibiao/1..6/:a;
+ f2 x( m1 W# }' e1 B% slinks(shengshi,zhibiao):x;
5 T$ |( V( I, z1 aendsets. ~7 \2 F5 B4 G& d
data:, v. F! P' p8 h  K2 [( c* z0 b3 a
x=
5 J- p: f* Z8 q0.099599  0.104095  0.021186  0.489849  0.031419  0.399601 9 {2 X1 T' r2 z8 x  Q9 ~1 z# M3 F
0.012222  0.006320  0.035075  0.558043  0.253776  0.465170
" a: r; K" d9 E: c1 g8 q8 |# {+ A0.016260  0.004213  0.084423  0.631442  0.375174  0.535579 3 F( d; D, d6 h- }; c6 E
0.004587  0.000000  0.086154  0.557522  0.210591  0.464721 0 o4 q4 }. U+ O2 O# R2 T
0.039880  0.015333  0.077921  1.000000  0.488760  0.888523 % u" W! {; K0 Z3 V: @) y# h
0.003290  0.002399  0.023441  0.299844  1.000000  0.217814
' `0 \2 e/ _% Q0 R5 K" f4 y0.021750  0.023296  0.055727  0.635086  0.657561  0.538723 ! ~! ?2 r5 h& Q: x, ?& M& H4 f
1.000000  1.000000  1.000000  0.000000  0.000000  0.430589 9 a0 `7 G5 V4 x/ b9 L& r3 z% j( V
0.070963  0.277539  0.385947  0.153566  0.014757  0.000000 ( K& g* Y8 L: z" K# D" d" ^
0.038164  0.038373  0.024990  0.489849  0.031159  1.000000
3 W! |2 g4 r5 M0.000000  0.008735  0.015259  0.288912  0.037009  0.040918 + I) _. c8 p+ m5 r
0.057181  0.038508  0.254109  0.184279  0.093160  0.052246 ' y8 p, p6 G* B# K' p1 d
0.048211  0.009641  0.000000  0.387819  0.135834  0.118513;
5 }5 M) \- s* o! B2 u% {! benddata
1 d# S6 _  C- `# G0 ?max=
$ y+ S$ c8 n& |(@sum(shengshi(i):& b9 W! x9 Z# ~- O% x0 u
(@sum(links(i,j):a(j)*x(i,j))-@sum(shengshi(i):
2 E1 i) M5 |# o/ z" Y" U@sum(links(i,j):a(j)*x(i,j)))/13)^2)/12)^(1/2);" W. P; ~/ b0 }  J

/ @8 \4 C8 ^/ G* B% x5 F@sum(zhibiao(j):a(j)^2)=1;
0 f1 W# `2 q# A6 x@for(zhibiao(j):a(j)>-1);. x- w& }% O% h/ w
@for(zhibiao(j):a(j)<1);5 L" V$ o1 Y, Q
@for(zhibiao(j)free(a(j)));, U$ V5 g/ _. ]: a+ e* _
end
" E% p  p1 i2 B. T( }7 P3 n, t! s求不出全局最优解，请大侠帮帮忙！！！
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.2 z: D1 O- o: Z4 P3 ?
      Objective value:                             0.7272518
    8 v+ d$ P* P; C" S+ _9 R4 y6 A  Objective bound:                             0.7272527
    4 B3 Z7 |: B: X  V6 V9 V  Infeasibilities:                              0.000000
    , z5 A4 a7 e' z2 @8 Z! K' C0 R  Extended solver steps:                           388433 l, F6 c; \( z4 V
      Total solver iterations:                       7561347
    + Y9 l) V# D1 m9 P5 u$ s! `" R, k8 D# U$ z8 u# j+ Q% s, H

    9 _; Y! [, R# U7 q- {/ D+ k                       Variable           Value        Reduced Cost
    ) ]; K2 a2 Y8 z                          A( 1)       0.1248020            0.000000! f; `5 J' j; T" z6 U6 V
                              A( 2)       0.1216772            0.000000# o' f9 }! i# ?* u; G. A
                              A( 3)       0.1621721            0.000000# L% ~7 {: j# x! k2 B; l
                              A( 4)       0.6318675            0.000000
    % w- z$ Y! Q9 D1 K! _& m! l                          A( 5)       0.3868056            0.000000/ ^5 h. T, x4 Q6 {/ C# \
                              A( 6)       0.6280479            0.000000  x& }- G8 H( X6 b' g6 A0 K9 A6 Q' G
                           X( 1, 1)       0.9959900E-01        0.000000
    2 n+ v6 M6 H2 x) L                       X( 1, 2)       0.1040950            0.000000; D1 {/ [2 j- s* i* b: S' a  \* x. W8 f
                           X( 1, 3)       0.2118600E-01        0.000000
    % S4 Y3 v+ m) e- Q                       X( 1, 4)       0.4898490            0.000000
    ( k$ [' E$ R! c  A                       X( 1, 5)       0.3141900E-01        0.000000
    - {1 p8 ]3 M* F: e7 L( t' z; [                       X( 1, 6)       0.3996010            0.000000
    , k% @$ Y+ K, I+ V, r- I" r                       X( 2, 1)       0.1222200E-01        0.000000' L4 P/ Y$ b: o" A8 x: z1 o7 ]/ C8 S
                           X( 2, 2)       0.6320000E-02        0.000000
    ; X. c7 W( p' w* [  A' I6 b                       X( 2, 3)       0.3507500E-01        0.0000004 i. i5 s4 X) t4 M* l2 Z9 |
                           X( 2, 4)       0.5580430            0.000000/ _+ e; @$ j- g1 c. l
                           X( 2, 5)       0.2537760            0.000000% p- `$ }0 U. Z; w" K
                           X( 2, 6)       0.4651700            0.000000" t6 E7 ^7 q; Q: Q$ Y8 B
                           X( 3, 1)       0.1626000E-01        0.000000. z$ n7 X# C: T
                           X( 3, 2)       0.4213000E-02        0.0000006 Z2 v! {: Y$ l% P, @7 G% ?# q
                           X( 3, 3)       0.8442300E-01        0.000000& y! d3 Y/ @# k" f* o0 t
                           X( 3, 4)       0.6314420            0.000000+ R3 @2 K. N5 f( s9 d) q
                           X( 3, 5)       0.3751740            0.000000
    7 T" F! O3 N  ~8 G' x) Z( I5 N                       X( 3, 6)       0.5355790            0.000000
    2 y% j0 l7 H/ l1 y                       X( 4, 1)       0.4587000E-02        0.0000003 E) g; J4 x) G9 P! \! A
                           X( 4, 2)        0.000000            0.000000
    * s  i% @9 s2 w6 L                       X( 4, 3)       0.8615400E-01        0.000000+ K( X3 \$ r2 ?
                           X( 4, 4)       0.5575220            0.0000002 C" p+ i& w( @7 I+ \* j6 w' w
                           X( 4, 5)       0.2105910            0.0000002 V  R# x2 c1 {6 }" ]; J; D3 F
                           X( 4, 6)       0.4647210            0.000000/ d" N+ G1 ^3 Q! ]/ F, t
                           X( 5, 1)       0.3988000E-01        0.000000
    ' ^& J, D6 D/ S$ O5 R                       X( 5, 2)       0.1533300E-01        0.0000001 L- e0 j* m; Y7 K) {& E
                           X( 5, 3)       0.7792100E-01        0.000000- K9 c  T/ _- K: }  k
                           X( 5, 4)        1.000000            0.000000
    9 d; @) q" A+ e                       X( 5, 5)       0.4887600            0.000000
    # S. a. V! t; S3 w                       X( 5, 6)       0.8885230            0.000000' _, Y" H' m5 y0 f5 W1 y
                           X( 6, 1)       0.3290000E-02        0.000000& \' r* ^5 k5 S: Q# t, A! m) T
                           X( 6, 2)       0.2399000E-02        0.000000
    . @0 W+ z4 ^% U7 h+ I, v+ L                       X( 6, 3)       0.2344100E-01        0.000000! `* j4 _9 p  L, C3 R0 b
                           X( 6, 4)       0.2998440            0.000000, Y! h, W$ W# l9 I6 f
                           X( 6, 5)        1.000000            0.000000) g! S3 \# z# w3 b: y1 y
                           X( 6, 6)       0.2178140            0.000000
    8 m" F, L7 w( W8 c# C* z  w                       X( 7, 1)       0.2175000E-01        0.000000
    ) F5 n, D0 w) z/ q6 l                       X( 7, 2)       0.2329600E-01        0.000000
    6 W) _8 M: c4 b7 B1 n$ `                       X( 7, 3)       0.5572700E-01        0.000000
    0 u2 O) P4 y; A. l                       X( 7, 4)       0.6350860            0.000000
    5 o4 T2 r( d$ T: R0 s: i0 A4 Q$ Z  b- I                       X( 7, 5)       0.6575610            0.000000
    " e. O% Z/ ^+ H9 c+ b                       X( 7, 6)       0.5387230            0.000000
    . S7 i$ i9 y6 Q6 E* P# q1 _2 s                       X( 8, 1)        1.000000            0.000000
    8 ~$ e* @! A6 a                       X( 8, 2)        1.000000            0.000000
    ' \' r: W& J/ B& \                       X( 8, 3)        1.000000            0.0000000 n" `. Y6 p! ]+ s
                           X( 8, 4)        0.000000            0.0000003 n) U. y& R* |  }+ [! E) @( |
                           X( 8, 5)        0.000000            0.000000
    4 F4 ?; V; ]  f: ~# s! \                       X( 8, 6)       0.4305890            0.000000% Q/ w1 E" R, a
                           X( 9, 1)       0.7096300E-01        0.0000005 r' M1 W& V: M' o% @  i4 o, [1 {
                           X( 9, 2)       0.2775390            0.000000* |. h! v- _0 Q9 k* B
                           X( 9, 3)       0.3859470            0.000000, y+ `3 [: b8 p- _; u  @7 @
                           X( 9, 4)       0.1535660            0.0000009 ]' W0 y" z4 `5 j# n
                           X( 9, 5)       0.1475700E-01        0.000000
    . h7 c7 G9 w* G                       X( 9, 6)        0.000000            0.0000001 F( v" ?7 j6 F: q
                          X( 10, 1)       0.3816400E-01        0.000000
    $ ~) [1 k6 l" \                      X( 10, 2)       0.3837300E-01        0.000000
    4 Z# J$ S3 c1 B7 [- t% N                      X( 10, 3)       0.2499000E-01        0.0000008 f  B  u. k# X' D. D& D* D
                          X( 10, 4)       0.4898490            0.000000  g2 l5 i& ~' k  O# Z" E
                          X( 10, 5)       0.3115900E-01        0.000000
    6 q; ]' y7 o) ^1 f7 |                      X( 10, 6)        1.000000            0.000000
    4 r% U$ @2 l: _4 }                      X( 11, 1)        0.000000            0.0000007 n" }/ K5 v+ s2 J- u
                          X( 11, 2)       0.8735000E-02        0.000000' d- T" Q7 N- ?+ s' D4 G- j
                          X( 11, 3)       0.1525900E-01        0.000000
    4 h8 M0 @1 u) p' W1 b" k) d3 k                      X( 11, 4)       0.2889120            0.000000
    + j- J5 k+ B: i4 j" X                      X( 11, 5)       0.3700900E-01        0.000000- x. Y$ P: e+ q' g+ n1 M+ Z
                          X( 11, 6)       0.4091800E-01        0.000000/ M% X( |+ |' ^) `# S0 Y
                          X( 12, 1)       0.5718100E-01        0.000000: p* x+ l3 v* [9 E. o
                          X( 12, 2)       0.3850800E-01        0.000000* b. t( ]) f8 N; N- o
                          X( 12, 3)       0.2541090            0.000000, X2 w" m+ B$ f/ m
                          X( 12, 4)       0.1842790            0.000000( _: P: W! W7 o7 N4 c  ]
                          X( 12, 5)       0.9316000E-01        0.000000
    * H1 H$ G7 _, ~! _# X7 B                      X( 12, 6)       0.5224600E-01        0.000000, _$ X2 F2 h4 }% u4 @' c1 Q, G( K
                          X( 13, 1)       0.4821100E-01        0.000000
    ! h9 D: U* a) R+ @                      X( 13, 2)       0.9641000E-02        0.000000( p, f/ E4 K* g2 y4 |
                          X( 13, 3)        0.000000            0.0000004 P2 g1 S1 w  L9 L
                          X( 13, 4)       0.3878190            0.0000007 e/ T+ J3 H( Z. A. @+ ^  Z; h* p
                          X( 13, 5)       0.1358340            0.000000' K  W1 Q; c9 q/ n& s& p
                          X( 13, 6)       0.1185130            0.0000009 j' Z! L4 \/ H1 H

    ' R7 U, o+ @' l: {- P  H2 x                            Row    Slack or Surplus      Dual Price/ h. F- k: q+ W  W
                                  1       0.7272518            1.000000
    ) c% d+ z* Z( A2 h$ @/ y# e                              2        0.000000           0.3636259
      R( f4 u& p9 d, ]- w                              3        1.124802            0.000000
    ) G: a1 g5 P! n  t0 n                              4        1.121677            0.0000006 s* x+ ?. p) K6 }& C
                                  5        1.162172            0.0000004 s& g9 R$ a, G* o
                                  6        1.631867            0.0000006 t% g+ L' u; b! {
                                  7        1.386806            0.0000003 X# D) M1 T: i3 r9 j* ?
                                  8        1.628048            0.000000
    $ M8 `6 `6 V$ x; M                              9       0.8751980            0.000000
    / J8 `) x! G4 C                             10       0.8783228            0.000000
    6 q& y* L: N, j% f                             11       0.8378279            0.0000004 y' m3 h% ~7 U. M' \
                                 12       0.3681325            0.0000004 [! b& z6 D: @; Y6 N6 q
                                 13       0.6131944            0.0000009 b9 l* \1 |6 w4 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
    但是 我设置全局求解时：两三秒就停止了
    # V2 g# m0 t; ]4 d出现错误提示
    9 Z7 L7 c/ C3 y" o# i7 e还有 上面的 全局的做大值还没有局部做大的大  为什么！！！
    回复

    使用道具 举报

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

    回顶部