快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

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

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

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

7

主题

5

听众

28

积分

升级  24.21%

该用户从未签到
跳转到指定楼层
1#
发表于 2009-9-5 07:49 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
model:& {  p/ E3 Z; Y% U$ k- q  T
sets:. f4 a% Q* }6 R8 @" f2 l
shengshi/1..13/;
" V7 M9 E+ Q* C) r' dzhibiao/1..6/:a;
! g% O) i9 c2 O( Flinks(shengshi,zhibiao):x;
- a/ J# F$ j* H  j" [, h" Yendsets  X- m4 o% F5 G6 N9 c
data:# @. Y* p4 d% U' {
x=
! A6 }4 u/ h- [& d: R3 Z" V0.099599  0.104095  0.021186  0.489849  0.031419  0.399601
+ e" L6 m, X( W5 d& ^" }0.012222  0.006320  0.035075  0.558043  0.253776  0.465170 . h9 Q, W, {7 q% ~2 f& z4 f
0.016260  0.004213  0.084423  0.631442  0.375174  0.535579
, I1 `3 a* i7 S1 Q6 I0.004587  0.000000  0.086154  0.557522  0.210591  0.464721 3 v/ d& H5 o& K2 [! P0 R
0.039880  0.015333  0.077921  1.000000  0.488760  0.888523 & A  e2 H! V' {, S0 V- a9 H
0.003290  0.002399  0.023441  0.299844  1.000000  0.217814 / G* w2 N6 E4 c1 h' L0 R
0.021750  0.023296  0.055727  0.635086  0.657561  0.538723
) U7 A' _+ T8 B6 ~$ _; u. N9 V1.000000  1.000000  1.000000  0.000000  0.000000  0.430589 ) N( B& K0 J/ ~& ?3 ?
0.070963  0.277539  0.385947  0.153566  0.014757  0.000000
3 V8 L" u7 i' [0 i% {# A3 {0.038164  0.038373  0.024990  0.489849  0.031159  1.000000
+ |, k0 N+ s( e1 h+ |3 l+ z0.000000  0.008735  0.015259  0.288912  0.037009  0.040918 ) V+ D  z0 O% x. M7 e- D) a
0.057181  0.038508  0.254109  0.184279  0.093160  0.052246
$ J' ^- Z" N8 I* ]' @! U; b" k0.048211  0.009641  0.000000  0.387819  0.135834  0.118513;
9 t' ]8 X6 f/ z6 ^' cenddata
1 U4 T! X: r; O) P( L: c& }. }max=" |' O2 H" L' u) [% G! }
(@sum(shengshi(i):" t9 n+ n) J! F
(@sum(links(i,j):a(j)*x(i,j))-@sum(shengshi(i):5 S" N2 I) V% Y3 M/ n  q+ D* Y
@sum(links(i,j):a(j)*x(i,j)))/13)^2)/12)^(1/2);1 l: b3 ?. m) Q5 n
$ N4 L8 r9 j8 |- X" S
@sum(zhibiao(j):a(j)^2)=1;+ J9 N# ^$ K* D1 Y
@for(zhibiao(j):a(j)>-1);$ v$ ?$ ^" q2 W: ?' h
@for(zhibiao(j):a(j)<1);+ M* P5 i' H5 J( q0 C1 f, E
@for(zhibiao(j)free(a(j)));
* S6 C2 a% n! Pend
& u( _* m! B2 X- n8 [求不出全局最优解，请大侠帮帮忙！！！
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.7 |: v2 n3 H: J
      Objective value:                             0.7272518
    / o/ n$ O" P8 `  Objective bound:                             0.7272527# |; [; x3 I0 ^& d/ ?7 u
      Infeasibilities:                              0.000000( q" f5 z: u& N5 w
      Extended solver steps:                           38843: a3 c9 Q, o% ~6 q6 c0 d- K  }
      Total solver iterations:                       7561347
    8 g$ \1 K! @  Z. [: t6 q- M4 n. E3 j" b$ B4 w  S- |3 |

    1 Q: M7 _4 K7 k' V( Z( s                       Variable           Value        Reduced Cost# r% o+ [" {$ _% Y. o: v
                              A( 1)       0.1248020            0.000000
    : x* y# l6 S7 r* j6 @9 a                          A( 2)       0.1216772            0.000000
    , U  V0 g9 E* u0 @) U+ ]3 q                          A( 3)       0.1621721            0.000000+ Z/ M- C. Y3 ^9 Q
                              A( 4)       0.6318675            0.000000  e% K: o$ J, {) A
                              A( 5)       0.3868056            0.0000008 w! r0 z+ P0 h6 w
                              A( 6)       0.6280479            0.000000
    ' n  F) D  t. K+ J# r                       X( 1, 1)       0.9959900E-01        0.000000
    9 g: H8 I! y' g7 |/ G9 [5 |0 ?# t                       X( 1, 2)       0.1040950            0.000000: E0 [8 {% R. H* l- ]! |: e
                           X( 1, 3)       0.2118600E-01        0.000000
    " H/ h$ n, G& x0 Y7 C! I) S                       X( 1, 4)       0.4898490            0.0000005 z( O2 u/ O7 q0 h
                           X( 1, 5)       0.3141900E-01        0.0000004 ?8 g. R( t7 m- R7 c+ i! U8 _; h
                           X( 1, 6)       0.3996010            0.0000002 ]* f0 T& b/ X
                           X( 2, 1)       0.1222200E-01        0.0000004 F+ k+ |9 i1 _
                           X( 2, 2)       0.6320000E-02        0.000000
    - m1 U7 j5 I- b/ z1 C                       X( 2, 3)       0.3507500E-01        0.000000
    $ Y5 Z0 I/ A6 s                       X( 2, 4)       0.5580430            0.000000- e$ S+ ^% O( O9 s
                           X( 2, 5)       0.2537760            0.0000006 |- S$ r, q$ ~8 }5 y( s1 \3 F- Z
                           X( 2, 6)       0.4651700            0.000000
    ) k: U" T5 _1 o) l2 P6 v4 E0 B( v                       X( 3, 1)       0.1626000E-01        0.000000
    8 c5 m9 G4 P% T4 |                       X( 3, 2)       0.4213000E-02        0.0000004 n9 L$ I+ O" j/ L9 b- W+ e
                           X( 3, 3)       0.8442300E-01        0.000000
    & v  P! g, I) h* _* z                       X( 3, 4)       0.6314420            0.000000
    ( I$ L/ [; i6 d: q0 {' w/ L                       X( 3, 5)       0.3751740            0.000000+ ^! V' o7 l2 `" x* h
                           X( 3, 6)       0.5355790            0.0000003 k; n7 R7 i/ d% x
                           X( 4, 1)       0.4587000E-02        0.000000' W1 d5 L0 H' o" n
                           X( 4, 2)        0.000000            0.000000
    6 ^8 y/ ?+ c6 C& j, a' Q0 P' U                       X( 4, 3)       0.8615400E-01        0.0000007 ^! f! C1 D' N# x& X
                           X( 4, 4)       0.5575220            0.0000009 {  H" C( W* A) o# U$ t* D
                           X( 4, 5)       0.2105910            0.0000009 r) H6 {/ e$ C  M' Y! t
                           X( 4, 6)       0.4647210            0.000000
    3 T: r' ~0 W; K9 A5 A$ s                       X( 5, 1)       0.3988000E-01        0.000000
    " h( `* H' O7 @                       X( 5, 2)       0.1533300E-01        0.000000( y; V9 ~2 [( A
                           X( 5, 3)       0.7792100E-01        0.000000
    - o. m. E9 N; O1 H" R                       X( 5, 4)        1.000000            0.0000003 h$ g6 F- g4 Z% x; e' s' f0 Z# L
                           X( 5, 5)       0.4887600            0.000000
    , _7 ]: ]' `, S: U- l                       X( 5, 6)       0.8885230            0.000000
    . H, d/ {  L8 A, V                       X( 6, 1)       0.3290000E-02        0.000000
    7 t0 I( y4 N4 a8 B1 |: A8 \/ d                       X( 6, 2)       0.2399000E-02        0.0000009 f! c" D& x, g8 g+ h' v
                           X( 6, 3)       0.2344100E-01        0.000000) q2 v5 B6 S1 [9 ^- V" h
                           X( 6, 4)       0.2998440            0.000000" {$ q, G7 n6 _. }( i
                           X( 6, 5)        1.000000            0.000000
    . y3 d# Y) u( e4 P* ~5 c+ ^                       X( 6, 6)       0.2178140            0.000000
    / J# x9 k3 s% Q, Q/ |                       X( 7, 1)       0.2175000E-01        0.000000$ [$ J, W% s. A9 o3 J1 L
                           X( 7, 2)       0.2329600E-01        0.000000
    ; b, |* x8 |) \! C                       X( 7, 3)       0.5572700E-01        0.000000. ?+ c3 V  V. z. q+ X+ K
                           X( 7, 4)       0.6350860            0.000000
    ! u# m3 P9 ?: c5 M$ G" L                       X( 7, 5)       0.6575610            0.0000007 R2 J" E1 \6 c# i
                           X( 7, 6)       0.5387230            0.000000
    8 b0 M! T4 }4 ^% g3 d% ~0 ?2 r8 n5 n! u                       X( 8, 1)        1.000000            0.000000
    % \; y: ?) Z8 K8 j2 C                       X( 8, 2)        1.000000            0.000000
    9 T9 d/ j  p3 O% @% `/ d; U                       X( 8, 3)        1.000000            0.0000000 z( t; |2 m' |+ _% u
                           X( 8, 4)        0.000000            0.000000
    : m* e0 x6 G' s( D- T% |7 E4 C                       X( 8, 5)        0.000000            0.000000
    : \: }7 ~4 s  `3 F1 i6 D0 L3 k                       X( 8, 6)       0.4305890            0.000000
    1 [: e- s8 G* [8 K; G( S5 e7 N                       X( 9, 1)       0.7096300E-01        0.000000
    0 e; R! d$ J" e  j6 N4 z                       X( 9, 2)       0.2775390            0.000000) h9 P; ?% x! g& x8 x7 Y8 `) c) D
                           X( 9, 3)       0.3859470            0.000000" Z( {' N) r- \3 h: J
                           X( 9, 4)       0.1535660            0.000000, E6 C9 E% `$ r# {+ d
                           X( 9, 5)       0.1475700E-01        0.000000+ v" E  {- ]8 e7 O5 l' k
                           X( 9, 6)        0.000000            0.0000007 V6 B( n, C. M  C+ ?
                          X( 10, 1)       0.3816400E-01        0.000000
    # S: A" m* s* y- u2 u1 n                      X( 10, 2)       0.3837300E-01        0.000000
    1 w" y4 W' L# R6 w4 c3 {, W                      X( 10, 3)       0.2499000E-01        0.000000
    & z4 p- o9 o- |$ {3 @& i8 ?( L6 X                      X( 10, 4)       0.4898490            0.000000
    / T  w2 w6 V9 z+ Q                      X( 10, 5)       0.3115900E-01        0.000000
    3 `& q& U7 t9 s                      X( 10, 6)        1.000000            0.000000  X% {. p7 t; p7 m
                          X( 11, 1)        0.000000            0.000000
    9 Z8 w2 Q  R) l5 T                      X( 11, 2)       0.8735000E-02        0.0000008 [6 U) h( z# G5 w# E
                          X( 11, 3)       0.1525900E-01        0.000000
    8 }/ _4 }+ z. P0 y, Q/ ]                      X( 11, 4)       0.2889120            0.0000007 ~& |! M+ g# l9 ~5 h) H
                          X( 11, 5)       0.3700900E-01        0.0000004 P' N  I. b8 {1 ], E
                          X( 11, 6)       0.4091800E-01        0.000000& k9 S' J/ @+ M/ P
                          X( 12, 1)       0.5718100E-01        0.000000  K: A! ]% `4 h/ ^% ]: x! v8 ^
                          X( 12, 2)       0.3850800E-01        0.000000
    ! _3 G6 `* g8 I& z# j+ q                      X( 12, 3)       0.2541090            0.000000
    " T+ f+ q. \6 y9 z                      X( 12, 4)       0.1842790            0.0000006 l2 Q9 k  ?2 A. z
                          X( 12, 5)       0.9316000E-01        0.0000001 i! ]$ [0 H/ F: S) d
                          X( 12, 6)       0.5224600E-01        0.0000006 q# @' X+ o" U' P' e! A
                          X( 13, 1)       0.4821100E-01        0.000000
    * @* b6 s) e7 }                      X( 13, 2)       0.9641000E-02        0.0000009 ^/ }6 v  G& }) e/ Q5 H
                          X( 13, 3)        0.000000            0.000000
    - e3 W" z0 K- l) i+ h& ^5 k                      X( 13, 4)       0.3878190            0.000000
    2 H8 f; J! D) b8 i1 z2 L                      X( 13, 5)       0.1358340            0.000000
    5 V' E2 [: g. l5 u5 Y" ^7 K. j                      X( 13, 6)       0.1185130            0.000000
    + a5 m- J# G% o  E7 a: x% _1 a8 e- |# W+ ]4 _' U! z  r% y1 d  G* l. m2 V
                                Row    Slack or Surplus      Dual Price
    " O. H' V! C7 M- r) q" ~& X                              1       0.7272518            1.000000
    3 w* B1 k2 F& X: ~- i& E- L                              2        0.000000           0.3636259
    2 P; F* Z0 Z  G  Y4 N; B! X1 ~                              3        1.124802            0.000000
    - O  f5 b. k% p' w- [/ e) w                              4        1.121677            0.0000001 Y( p* i3 a3 S" L
                                  5        1.162172            0.000000
      p1 @& q$ ^' L" \9 v& B, X                              6        1.631867            0.000000
    3 Q  L3 \( \, I# \8 t2 m- K8 h                              7        1.386806            0.000000
    4 T# W: }4 C; J- r                              8        1.628048            0.0000001 H1 J4 N6 y- d$ h( N6 Z+ x, d
                                  9       0.8751980            0.000000( T# W) l+ g, @+ ^9 @3 e
                                 10       0.8783228            0.000000
    ( {! K& b8 N7 P. W8 M                             11       0.8378279            0.000000
    5 f; i( E3 |6 |- v                             12       0.3681325            0.000000" S: X0 U1 x2 n8 g7 S
                                 13       0.6131944            0.000000
    " z9 {8 S' \% o; y9 e; d& w                             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
    但是 我设置全局求解时：两三秒就停止了4 z. l; N8 @4 M2 d$ d
    出现错误提示
    / Q' ~! _0 F, J: q/ I4 Q6 w还有 上面的 全局的做大值还没有局部做大的大  为什么！！！
    回复

    使用道具 举报

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

    回顶部