快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

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

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

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

7

主题

5

听众

28

积分

升级  24.21%

该用户从未签到
跳转到指定楼层
1#
发表于 2009-9-5 07:49 |只看该作者 |正序浏览
|招呼Ta 关注Ta
model:
3 [) m2 V$ E/ l, s+ n% \sets:; Z* p( H; Q' `6 O" g2 h& a
shengshi/1..13/;
) B/ @% o; }: w5 q& bzhibiao/1..6/:a;# W. `$ J! M0 f
links(shengshi,zhibiao):x;
  N2 s' W, l  u1 J! t* C6 ]  F' lendsets
2 G, _4 i" v) A1 Rdata:' d0 ?8 y- I: {' I& L+ s
x=
6 K1 t  Q4 u6 N9 W0.099599  0.104095  0.021186  0.489849  0.031419  0.399601
- s- y" W1 p$ |( ^" F0.012222  0.006320  0.035075  0.558043  0.253776  0.465170
6 @& s. o! b* Z; k* X6 L0.016260  0.004213  0.084423  0.631442  0.375174  0.535579 1 ~* {9 K6 u, z5 R* O
0.004587  0.000000  0.086154  0.557522  0.210591  0.464721 6 S! F: E0 T. a# Y  V3 K7 S! g
0.039880  0.015333  0.077921  1.000000  0.488760  0.888523
  C, A; W$ y" E; {7 V3 p( y& w1 ?0.003290  0.002399  0.023441  0.299844  1.000000  0.217814 ( |3 B, n. D  K' {4 m) A: \; P2 n
0.021750  0.023296  0.055727  0.635086  0.657561  0.538723 7 ^5 D  Y4 N1 y& h0 y0 O
1.000000  1.000000  1.000000  0.000000  0.000000  0.430589   }( a& L+ z5 |
0.070963  0.277539  0.385947  0.153566  0.014757  0.000000 9 v6 G+ I6 d& b1 e# [+ J" x
0.038164  0.038373  0.024990  0.489849  0.031159  1.000000
4 w' [* u- b3 A. A- @- O0.000000  0.008735  0.015259  0.288912  0.037009  0.040918
% ]0 Q% u% h' O& C6 Q0.057181  0.038508  0.254109  0.184279  0.093160  0.052246 + y3 }' ]4 i3 q/ X2 b
0.048211  0.009641  0.000000  0.387819  0.135834  0.118513;
  ^3 z$ m. l6 C4 Q. d! Penddata- u7 d( A0 x3 T0 ]/ B* t9 ?) f
max=
& Y4 \, G* g$ T2 W6 Y: m(@sum(shengshi(i):' |8 h' Y; W9 h" R, V8 u& T5 L4 \5 a
(@sum(links(i,j):a(j)*x(i,j))-@sum(shengshi(i):% L- I1 k0 E* K. Z! b! l
@sum(links(i,j):a(j)*x(i,j)))/13)^2)/12)^(1/2);" N6 _( ~5 E# S! [

; x( a) c4 p+ X5 g@sum(zhibiao(j):a(j)^2)=1;" s* W. n: j7 K& q# i( R
@for(zhibiao(j):a(j)>-1);, w& A0 i/ `& a  B
@for(zhibiao(j):a(j)<1);
& _! q1 V7 F( a; B1 h5 `@for(zhibiao(j)free(a(j)));( a+ O" @& e1 X0 h. e- N
end
3 _: l2 z( t$ H! Y6 b) ^* ?求不出全局最优解，请大侠帮帮忙！！！
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
回复

使用道具 举报

jiiing        

7

主题

5

听众

28

积分

升级  24.21%

该用户从未签到
10#
发表于 2009-9-5 12:16 |只看该作者 |招呼Ta 关注Ta
但是 我设置全局求解时：两三秒就停止了
. A' w6 V& L) U8 _- r* J出现错误提示
5 o" [) o2 u" p! ]( f还有 上面的 全局的做大值还没有局部做大的大  为什么！！！
回复

使用道具 举报

liulongfei 实名认证       

2

主题

2

听众

100

积分

升级  0%

该用户从未签到

群组: 数学趣味、游戏、IQ等

群组: 数学建摸协会

群组: 华中科技大学
9#
发表于 2009-9-5 09:35 |只看该作者 |招呼Ta 关注Ta
选择全局最优解，然后就是时间的问题了
回复

使用道具 举报

xuyin        

0

主题

4

听众

198

积分

升级  49%

  • TA的每日心情
    开心
    2016-3-23 08:32

    • 签到天数: 26 天

    [LV.4]偶尔看看III
    8#
    发表于 2009-9-5 09:06 |只看该作者 |招呼Ta 关注Ta
    哇，都是高手啊
    回复

    使用道具 举报

    yuan07150215        

    0

    主题

    4

    听众

    168

    积分

    升级  34%

    该用户从未签到
    7#
    发表于 2009-9-5 08:58 |只看该作者 |招呼Ta 关注Ta
    先在限制条件下得出最低值，之后运行lingo程序，当运行到20-30分钟时，人为的终止求解（按下Interrupt Solver），此时运行的解与目标函数的最优解很相近！可认为是最优解！……
    回复

    使用道具 举报

    formulaone9275        

    1

    主题

    4

    听众

    188

    积分

    升级  44%

    该用户从未签到
    6#
    发表于 2009-9-5 08:54 |只看该作者 |招呼Ta 关注Ta
    有很多时候都是求不出最优解的，局部最优就可以用了 ，我看了你的程序，其中max里的sum有对i求了两次和，你看看能不能再改改，求求试试
    回复

    使用道具 举报

    ychexingwuji        

    4

    主题

    4

    听众

    70

    积分

    升级  68.42%

    该用户从未签到
    5#
    发表于 2009-9-5 08:49 |只看该作者 |招呼Ta 关注Ta
    学习了 呵呵
    回复

    使用道具 举报

    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分钟才计算出来，给他一个初始值就好了！
    回复

    使用道具 举报

    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 X9 M+ A6 Q; ?. J9 w( j- C# p  Objective value:                             0.7272518
    ; v3 v  W4 N' b, X/ Z- U7 c4 x; W  Objective bound:                             0.7272527
      k% B* Q: Y2 n% I  Infeasibilities:                              0.000000' u( b6 ?, S8 ~' N. j3 \7 Y
      Extended solver steps:                           38843
    % Y" D  {9 X5 [% j( Q  Total solver iterations:                       7561347
    6 R& Z6 W2 |8 n; O" T, @
    7 I* k) P1 e* Z- v' N& w5 M5 F1 u# J3 h  v
                           Variable           Value        Reduced Cost( ]& n: l: O9 Y: G9 m
                              A( 1)       0.1248020            0.0000008 I! |1 P4 ?% @" K+ D
                              A( 2)       0.1216772            0.000000( t, ~& G) H% ~9 `! F( `* X
                              A( 3)       0.1621721            0.000000, k# ]1 A2 S2 v3 d0 Q" H
                              A( 4)       0.6318675            0.000000& o/ _& V4 {) g' ~8 h% Y
                              A( 5)       0.3868056            0.000000
    , O  q2 o  F, Z. H                          A( 6)       0.6280479            0.000000
    2 ^4 L) v6 Q* W" Y9 V$ P                       X( 1, 1)       0.9959900E-01        0.000000
    / a  f3 j! b) F, h                       X( 1, 2)       0.1040950            0.000000: T# T2 A( r( t9 }% }
                           X( 1, 3)       0.2118600E-01        0.000000
    7 S8 a) M  T8 d1 W) V+ \                       X( 1, 4)       0.4898490            0.000000$ c4 G- @5 m$ U
                           X( 1, 5)       0.3141900E-01        0.000000
    + k& l. h/ t% o3 v" S9 p* z- g                       X( 1, 6)       0.3996010            0.000000
    # C# R& a: [- P: m, ?                       X( 2, 1)       0.1222200E-01        0.0000000 m( a! O/ B, |) u
                           X( 2, 2)       0.6320000E-02        0.000000- G" U! Y2 \" D9 L/ _
                           X( 2, 3)       0.3507500E-01        0.000000( L+ m' W: E, L' u
                           X( 2, 4)       0.5580430            0.000000( m, ~1 J1 N; l9 I& a
                           X( 2, 5)       0.2537760            0.000000
    - Z$ j/ B" E7 X+ l( a5 h( Q5 T2 ^                       X( 2, 6)       0.4651700            0.000000) G! x. Q& f/ c1 ]+ N- b) b" N
                           X( 3, 1)       0.1626000E-01        0.0000008 {; S& |5 R* _/ F
                           X( 3, 2)       0.4213000E-02        0.000000
    ' W2 `/ ]$ u! E4 G' X# R9 F                       X( 3, 3)       0.8442300E-01        0.000000' n- ]) f' N0 t; @' Y1 f
                           X( 3, 4)       0.6314420            0.0000004 \; |) ^4 ]) v" ]4 p7 j
                           X( 3, 5)       0.3751740            0.000000- @, D; [1 w- c, ?: I
                           X( 3, 6)       0.5355790            0.0000000 j" @' B  s, b' k: a( s/ |/ ^, d
                           X( 4, 1)       0.4587000E-02        0.000000  k$ [0 {4 ~& X
                           X( 4, 2)        0.000000            0.000000
    1 y: x0 T7 m8 h3 T5 F* S1 s0 A# a  G6 P. b7 e                       X( 4, 3)       0.8615400E-01        0.000000
    ( I. I0 S/ J3 v+ E: |# h: O  ~                       X( 4, 4)       0.5575220            0.0000001 [3 L. P# o- {& x( u
                           X( 4, 5)       0.2105910            0.0000006 m1 t; j# w) d: Y
                           X( 4, 6)       0.4647210            0.000000
    6 `; ]" X$ K2 j                       X( 5, 1)       0.3988000E-01        0.000000
    8 Q* g0 O- B4 y( G$ a1 m8 O8 @5 u. A                       X( 5, 2)       0.1533300E-01        0.000000( V1 O& Z, Z9 X. z. w. |% I' s: e
                           X( 5, 3)       0.7792100E-01        0.000000
    ) z' X3 F* _. k$ J# ]# X1 o                       X( 5, 4)        1.000000            0.000000
    ; G' u# S' q  Q* U+ ]$ I                       X( 5, 5)       0.4887600            0.000000
    % d6 K: v5 ]! H6 S! l                       X( 5, 6)       0.8885230            0.000000
    & S2 M4 F' B" u0 s/ f- I* D                       X( 6, 1)       0.3290000E-02        0.000000
    . M  ^9 O9 W2 A: `9 i& [, E                       X( 6, 2)       0.2399000E-02        0.000000# s3 O3 o, z; S+ N6 @' f/ g( d1 D
                           X( 6, 3)       0.2344100E-01        0.000000
    2 _, x% p( v) u' |- u/ |                       X( 6, 4)       0.2998440            0.000000
    3 o% ~4 q' p: K$ {                       X( 6, 5)        1.000000            0.000000. S# @4 ?) _2 l+ _; u
                           X( 6, 6)       0.2178140            0.000000, r0 B2 O2 M. m2 v* ~: k% [/ w
                           X( 7, 1)       0.2175000E-01        0.000000+ `# `* f* }* B0 h' K
                           X( 7, 2)       0.2329600E-01        0.000000# d" L0 z! N9 V: k; J
                           X( 7, 3)       0.5572700E-01        0.000000
    - K: E9 G& E- J7 W0 O- |; g                       X( 7, 4)       0.6350860            0.000000( I2 t0 `" j! f: Z
                           X( 7, 5)       0.6575610            0.000000
    ) ?" ?! q# N0 T3 W                       X( 7, 6)       0.5387230            0.000000) x9 r' j5 P9 ]4 `- @, ?! i( q
                           X( 8, 1)        1.000000            0.000000+ T$ d: G  X" R' t* u
                           X( 8, 2)        1.000000            0.000000* z2 f) b* ^: B: C; t  ^" [6 s7 a
                           X( 8, 3)        1.000000            0.000000
    ! n* C. \$ z8 H                       X( 8, 4)        0.000000            0.000000
    7 e* o" V4 |8 d- e* Y                       X( 8, 5)        0.000000            0.000000& b, d9 a2 X+ W
                           X( 8, 6)       0.4305890            0.0000004 E$ L' ?* p% ?5 o6 e: t; Y1 y. u
                           X( 9, 1)       0.7096300E-01        0.000000
    ; O% p5 ?6 z7 j6 t9 b& D' _' C% _' p                       X( 9, 2)       0.2775390            0.000000
    0 f" A# F+ Q0 _$ V7 U; b' t, S                       X( 9, 3)       0.3859470            0.000000+ y6 G+ T2 [  ^$ I- t
                           X( 9, 4)       0.1535660            0.000000
    % A4 q5 B+ ^. B+ l                       X( 9, 5)       0.1475700E-01        0.000000& }# T5 S. {. m  U1 ]2 Z
                           X( 9, 6)        0.000000            0.000000
    1 y$ s7 ~) ^; n. D- a' y  r1 v                      X( 10, 1)       0.3816400E-01        0.000000! K  i: V* o. F! `% `6 N9 u  I9 N
                          X( 10, 2)       0.3837300E-01        0.000000$ p, d7 c$ J* H, D2 \% I) c) r- W
                          X( 10, 3)       0.2499000E-01        0.000000# j& I- w. Q4 j- M: }& o4 P3 N
                          X( 10, 4)       0.4898490            0.000000& J5 Y: U! f% O
                          X( 10, 5)       0.3115900E-01        0.000000$ s/ a, M* U* [- N
                          X( 10, 6)        1.000000            0.000000
    % F; G5 R. I4 {% v* D+ E                      X( 11, 1)        0.000000            0.000000! B) ?: `2 u4 K  f8 a% D- i
                          X( 11, 2)       0.8735000E-02        0.000000
    8 ?" A; W7 R2 {                      X( 11, 3)       0.1525900E-01        0.000000
    : m  D! L% B: Y+ M) i9 K/ D/ ^9 A$ b5 n                      X( 11, 4)       0.2889120            0.000000( `. ]$ V) a! ]3 \7 j" Z
                          X( 11, 5)       0.3700900E-01        0.0000003 x# d% C+ Y6 h9 k( h
                          X( 11, 6)       0.4091800E-01        0.000000
    5 G# e. \4 x7 n; Z' ^, T* h0 _& ~                      X( 12, 1)       0.5718100E-01        0.000000& ]  x' a; D9 ?5 X& S
                          X( 12, 2)       0.3850800E-01        0.000000# y- N2 E2 W3 S( Z. y
                          X( 12, 3)       0.2541090            0.000000/ w5 D5 B2 x  I  n; ~; a
                          X( 12, 4)       0.1842790            0.000000/ A4 }. f3 Z4 l/ B! S
                          X( 12, 5)       0.9316000E-01        0.000000& F% }! D/ j  s6 z% |/ U8 ]0 `/ J; L
                          X( 12, 6)       0.5224600E-01        0.0000005 g- ]4 s" Z% n( {6 ]1 M
                          X( 13, 1)       0.4821100E-01        0.000000
    . C) w* ^. l% T9 U* P; L                      X( 13, 2)       0.9641000E-02        0.000000; D8 O3 f5 W3 {2 L
                          X( 13, 3)        0.000000            0.000000- G* D' O6 ~* X4 o4 O9 `6 M" B
                          X( 13, 4)       0.3878190            0.000000' |. Q2 Y8 i; l1 x* C
                          X( 13, 5)       0.1358340            0.000000
    2 p$ v8 m- X; ?2 b9 F" X                      X( 13, 6)       0.1185130            0.0000008 d/ D2 M3 m7 u$ \6 D

    4 p, _. M3 A; K                            Row    Slack or Surplus      Dual Price( m. O  t; V, f9 d1 i; a
                                  1       0.7272518            1.000000
    # a4 o+ w$ }& N. P                              2        0.000000           0.3636259' s2 L1 |# y# {: A
                                  3        1.124802            0.000000
    0 @5 J2 P" K6 u& T9 h1 R+ Q0 ]                              4        1.121677            0.000000
    ! i5 T, U2 o  a& c3 @+ m0 O3 q                              5        1.162172            0.000000
    # p8 f( D& O& T# r                              6        1.631867            0.000000& @5 K8 b* z5 x: O) i
                                  7        1.386806            0.000000
    * U' J% f4 ?4 R6 m( T                              8        1.628048            0.000000
    $ D/ A" O$ E# ~2 ^3 F4 H                              9       0.8751980            0.0000003 G1 }* E. B  E  L+ I& @
                                 10       0.8783228            0.0000001 ]6 ~5 D- f' B) o7 n$ d
                                 11       0.8378279            0.000000
    6 l7 x8 F7 a) o                             12       0.3681325            0.000000& T6 [+ ~+ _) p' l/ @
                                 13       0.6131944            0.000000
    / l& G4 ]9 x6 C& M# ]/ p9 A9 ?4 A: B  n6 I/ a                             14       0.3719521            0.000000
    回复

    使用道具 举报

    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
    运算时间太长，你能否增加一个初始值就好了
    回复

    使用道具 举报

    返回列表 发新帖
    高级模式
    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 04:46 , Processed in 2.357825 second(s), 104 queries .

    回顶部