快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

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

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

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

7

主题

5

听众

28

积分

升级  24.21%

该用户从未签到
跳转到指定楼层
1#
发表于 2009-9-5 07:49 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
model:. D% Q4 \5 t4 e! u6 e: [2 L9 P
sets:
3 G8 Y8 P2 S) O9 dshengshi/1..13/;% a/ O% V0 B+ f" Q2 V# ^% S
zhibiao/1..6/:a;" z6 x; l$ O% U1 z$ T; o
links(shengshi,zhibiao):x;8 `/ `9 ]: {1 s: s8 `
endsets
( M, u7 c/ x6 c8 Mdata:
  ^7 N; O) }1 I# H+ Sx=/ t6 x- k& L" Y) p5 Z% o7 v
0.099599  0.104095  0.021186  0.489849  0.031419  0.399601 ! K6 Y# M+ a& S/ p8 X+ \; j5 E& u
0.012222  0.006320  0.035075  0.558043  0.253776  0.465170
" q! R9 f9 G; X; }! k0.016260  0.004213  0.084423  0.631442  0.375174  0.535579
2 O6 I- x  C& r3 ?& t/ s0.004587  0.000000  0.086154  0.557522  0.210591  0.464721 2 X+ c" R% H2 c4 t1 e4 u3 _# d
0.039880  0.015333  0.077921  1.000000  0.488760  0.888523 . E) \( w: p4 T$ ]; y. G
0.003290  0.002399  0.023441  0.299844  1.000000  0.217814 7 Z0 O& d9 D- S5 k- Y
0.021750  0.023296  0.055727  0.635086  0.657561  0.538723 3 A) @: c. C9 s' _* a
1.000000  1.000000  1.000000  0.000000  0.000000  0.430589 ( z& x9 r, K5 w1 O8 }. J
0.070963  0.277539  0.385947  0.153566  0.014757  0.000000
6 u0 y3 O! ~* M+ @0.038164  0.038373  0.024990  0.489849  0.031159  1.000000
' R; U' S/ W2 l: ?& m, F; `( V8 i0.000000  0.008735  0.015259  0.288912  0.037009  0.040918
. [7 g5 z3 q0 Z- G8 @+ d/ ^; O0.057181  0.038508  0.254109  0.184279  0.093160  0.052246
9 g" m; q* b6 v  \4 ~) n0.048211  0.009641  0.000000  0.387819  0.135834  0.118513;
6 H4 U* p) b+ Q; P& A! a0 D, h' nenddata
8 F2 R$ S2 R9 w$ n# C6 w+ Nmax=+ A' S/ B6 h3 d# B/ F. Y
(@sum(shengshi(i):- Q  }% G$ M4 u8 D
(@sum(links(i,j):a(j)*x(i,j))-@sum(shengshi(i):
' H7 W$ S9 T% ~2 B  y9 M  {" O. W@sum(links(i,j):a(j)*x(i,j)))/13)^2)/12)^(1/2);1 S2 H* z9 w1 S, X

5 q- z" ]. l" d+ L; _# H3 Y@sum(zhibiao(j):a(j)^2)=1;# J) z/ l; Y3 A2 V. ^
@for(zhibiao(j):a(j)>-1);, C/ k& ^% K3 C: O8 \5 _$ a
@for(zhibiao(j):a(j)<1);
9 a( Y% N# j6 n3 C2 G@for(zhibiao(j)free(a(j)));4 e7 `+ x. I" J; J/ D$ o
end
0 H$ x/ l7 ]7 M" Z# E, b2 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.5 G9 a" q5 r; J
      Objective value:                             0.7272518
    9 h+ B# }1 g2 H" H, E. e! Z# T  Objective bound:                             0.7272527( b* E# M. E: U' f3 K8 T1 ?
      Infeasibilities:                              0.000000
    + y. X: P5 l" f" b8 J- E( X  Extended solver steps:                           38843
    . {0 {7 H" A$ |' l& \  Total solver iterations:                       7561347
    $ o) z4 }' s4 Y; c3 d/ a: u* L  |! R% ?$ {& e; ?# ?

    : p& R# |- K6 `. s3 x# A$ N- C2 m  d                       Variable           Value        Reduced Cost
    8 v  [& S6 ~, A  p                          A( 1)       0.1248020            0.000000# e1 v+ z, V4 E* W; d( R$ V
                              A( 2)       0.1216772            0.000000
    / X8 @( L" ?/ r; `9 j1 e/ c; `                          A( 3)       0.1621721            0.000000
    ; _# V$ j7 x) \/ y7 [$ P                          A( 4)       0.6318675            0.0000009 A& S& ~7 b, h2 N' I% T
                              A( 5)       0.3868056            0.0000003 K  o7 c! n7 w" Z
                              A( 6)       0.6280479            0.000000
    4 ?- E6 a0 M  b6 }$ M5 Y1 {1 Z                       X( 1, 1)       0.9959900E-01        0.000000
    7 m  s& X3 {  n4 D: d) x& M$ m                       X( 1, 2)       0.1040950            0.000000
    3 b$ c0 r% U/ Q- ~' j" C6 d                       X( 1, 3)       0.2118600E-01        0.0000009 {4 a, Z1 D0 S6 w1 c/ O, ~
                           X( 1, 4)       0.4898490            0.0000007 I* q1 O4 H, o. w# w. o/ y
                           X( 1, 5)       0.3141900E-01        0.000000) e# t1 [* \; L0 y8 s' s% K* e
                           X( 1, 6)       0.3996010            0.000000: F# c% }, y1 }3 `  l, J( \
                           X( 2, 1)       0.1222200E-01        0.000000
    * G/ }2 U& k6 X$ a6 {                       X( 2, 2)       0.6320000E-02        0.000000
    7 o7 C9 @) L3 f7 m. D' U! c                       X( 2, 3)       0.3507500E-01        0.000000( h, ^5 j8 q% q& N9 E
                           X( 2, 4)       0.5580430            0.000000- m2 k* f! X7 k2 F0 H3 o  ~& @
                           X( 2, 5)       0.2537760            0.0000005 \1 Q9 M7 Z% T/ a) _/ J
                           X( 2, 6)       0.4651700            0.0000005 K# G: H& R9 u& y3 h0 [
                           X( 3, 1)       0.1626000E-01        0.000000
    : I* V' Y9 W+ G% U% j: A                       X( 3, 2)       0.4213000E-02        0.000000" f7 M4 b. x" N0 c  @
                           X( 3, 3)       0.8442300E-01        0.000000
    ' z9 @4 L" ^- a4 S1 s& ^                       X( 3, 4)       0.6314420            0.000000
    ; ^. O: n4 v* |, n: ~& w                       X( 3, 5)       0.3751740            0.000000
    9 D" {( B4 Q" e$ E6 P2 i; C                       X( 3, 6)       0.5355790            0.000000& T6 D/ V9 h' N6 @
                           X( 4, 1)       0.4587000E-02        0.000000
    2 c+ h1 i( M3 R2 }# A. n' d! o0 g  S6 d                       X( 4, 2)        0.000000            0.000000
    6 j; D4 V! `% }2 t5 B                       X( 4, 3)       0.8615400E-01        0.000000) d! W# ?' J3 P# h' d* Z: s
                           X( 4, 4)       0.5575220            0.0000000 _- n! W* @# ~- [4 O+ _
                           X( 4, 5)       0.2105910            0.000000
    " b, H* @" k) m. l; f                       X( 4, 6)       0.4647210            0.0000003 G2 i( K: `8 F6 f2 g& [
                           X( 5, 1)       0.3988000E-01        0.000000* ?  j- U8 I- t5 i
                           X( 5, 2)       0.1533300E-01        0.000000
    - f; j( E# j6 i6 F( z0 C* u                       X( 5, 3)       0.7792100E-01        0.000000
    6 G( ]* D5 X; A/ l, G( a% n                       X( 5, 4)        1.000000            0.000000
    # g5 w, N" ], P+ g                       X( 5, 5)       0.4887600            0.000000
    " ]4 e( z7 i: S1 M0 \; R                       X( 5, 6)       0.8885230            0.000000( G1 s1 M* s9 I4 N
                           X( 6, 1)       0.3290000E-02        0.000000. r+ ]# G, J9 q# r3 w' v& P3 Z# K
                           X( 6, 2)       0.2399000E-02        0.000000+ t2 I" B8 ?" U
                           X( 6, 3)       0.2344100E-01        0.000000
    # A  ?1 G7 s  A. @. d' t" |                       X( 6, 4)       0.2998440            0.000000
    7 O1 I( U: D( k1 ]                       X( 6, 5)        1.000000            0.000000
    / S- h* c% B$ I. t1 \                       X( 6, 6)       0.2178140            0.000000: j0 z3 _- q; s
                           X( 7, 1)       0.2175000E-01        0.000000/ v) ?0 u8 h, U6 w
                           X( 7, 2)       0.2329600E-01        0.000000& W+ m: G7 K; j2 p, H/ A
                           X( 7, 3)       0.5572700E-01        0.0000000 s7 P  P3 A) Z4 V$ h- n
                           X( 7, 4)       0.6350860            0.000000
    5 ?* k2 H! q7 ^                       X( 7, 5)       0.6575610            0.000000. q) S/ C) @: b) e
                           X( 7, 6)       0.5387230            0.000000# x( d* S9 T' N% s2 y
                           X( 8, 1)        1.000000            0.000000
    : P, @% \( C% j                       X( 8, 2)        1.000000            0.000000
    2 a7 c2 V: L+ _1 s3 U                       X( 8, 3)        1.000000            0.000000
    / C* O6 B' D8 L, y) {3 s                       X( 8, 4)        0.000000            0.0000001 q9 a) H8 X9 w9 D, u7 u& N
                           X( 8, 5)        0.000000            0.0000009 [: f3 |6 h2 h3 G* ?2 X
                           X( 8, 6)       0.4305890            0.000000* f0 X) ?% ?4 Q5 O% W
                           X( 9, 1)       0.7096300E-01        0.000000
    * ~3 A7 o1 i. n" s                       X( 9, 2)       0.2775390            0.000000: ^1 p5 n; `" n" M6 h8 ^
                           X( 9, 3)       0.3859470            0.000000
    / W; Z( r  F* q% U$ H( H* `) S                       X( 9, 4)       0.1535660            0.000000
    ; E6 \' ?  L9 _" S0 q                       X( 9, 5)       0.1475700E-01        0.0000008 X+ T2 v, ]% }0 A, d5 b3 C
                           X( 9, 6)        0.000000            0.000000
    6 M' D' v8 K# G+ x5 v                      X( 10, 1)       0.3816400E-01        0.000000
    , f$ {+ c1 B; j3 t8 A                      X( 10, 2)       0.3837300E-01        0.000000
    / `* k0 y, ^  u& n* {& B2 Z, t                      X( 10, 3)       0.2499000E-01        0.000000
    $ D$ ^9 `/ E* Y# b' o: @6 D                      X( 10, 4)       0.4898490            0.000000
    4 l7 f6 f3 @% X! Y; |$ |                      X( 10, 5)       0.3115900E-01        0.000000
    ' g/ {5 \2 e, u0 P                      X( 10, 6)        1.000000            0.000000! ~: R( ^$ P! x  e% B; q" m
                          X( 11, 1)        0.000000            0.000000
      E9 h' P% M. r" d8 _/ S. Y8 e% V                      X( 11, 2)       0.8735000E-02        0.000000
    : p$ }* e6 N6 t2 F8 Y7 h  j                      X( 11, 3)       0.1525900E-01        0.000000  P2 I/ g" e: ?6 `& V8 O, L
                          X( 11, 4)       0.2889120            0.0000000 s$ D8 R3 B; `5 g7 w
                          X( 11, 5)       0.3700900E-01        0.000000
    2 X" S+ }) Z6 k$ {' G, p                      X( 11, 6)       0.4091800E-01        0.000000% d; F& x+ ^( A" P1 Y/ t1 a( b# V
                          X( 12, 1)       0.5718100E-01        0.000000
    6 A  s1 B1 w$ X& G4 b4 @                      X( 12, 2)       0.3850800E-01        0.000000( F8 w6 B: `, C/ n3 v; T
                          X( 12, 3)       0.2541090            0.000000
    5 ^. y4 n1 [% n% I2 p                      X( 12, 4)       0.1842790            0.000000
    0 x3 W8 I( F6 p: }/ V                      X( 12, 5)       0.9316000E-01        0.000000  T  H* F1 n, L3 ?/ M
                          X( 12, 6)       0.5224600E-01        0.0000007 N7 B6 P  b5 U2 x/ [  q
                          X( 13, 1)       0.4821100E-01        0.000000$ n% |' B, a  ^  \
                          X( 13, 2)       0.9641000E-02        0.000000; h% r$ m- d2 v2 n8 v% {
                          X( 13, 3)        0.000000            0.000000) z# j9 R$ @# y
                          X( 13, 4)       0.3878190            0.0000005 j9 N% O- z. C5 Q6 m1 G+ P
                          X( 13, 5)       0.1358340            0.000000! _6 b4 S; E+ q) ^$ q( U
                          X( 13, 6)       0.1185130            0.000000. k2 T; I! y6 w4 W" k

    ' \! u- D% t, l3 L                            Row    Slack or Surplus      Dual Price+ ?8 }6 [$ a6 ?6 I0 O8 Y, X
                                  1       0.7272518            1.000000
    / _; `5 l- m" s2 i; p                              2        0.000000           0.3636259
    ( l) J. W  S$ w/ n0 F% F4 A$ Q                              3        1.124802            0.000000
    5 |- Y9 ~+ _' \9 ~/ Z                              4        1.121677            0.000000
    7 w  t( f) X- O) Y2 R: v. b+ b                              5        1.162172            0.0000008 o6 x/ `7 }% F
                                  6        1.631867            0.000000
      s" x$ ]" u4 B! N& k: A                              7        1.386806            0.000000
    5 |" G1 x' g, S9 X                              8        1.628048            0.000000# S2 K3 l- a- D. N# q6 W
                                  9       0.8751980            0.000000) m; V- G9 x, `, Z
                                 10       0.8783228            0.000000
    + s7 z4 f( A- {: f" o                             11       0.8378279            0.0000008 X5 A- s8 R( D6 Q
                                 12       0.3681325            0.000000
    ; H' V1 |/ i/ S* W( H7 I                             13       0.6131944            0.0000000 f+ i- g# T' @/ N) N3 @1 Q
                                 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
    但是 我设置全局求解时：两三秒就停止了
    * I: u, [( K9 ^/ [: V出现错误提示: Z$ N) s4 f& 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-5-5 21:59 , Processed in 0.410836 second(s), 103 queries .

    回顶部