快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

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

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

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

7

主题

5

听众

28

积分

升级  24.21%

该用户从未签到
跳转到指定楼层
1#
发表于 2009-9-5 07:49 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
model:
$ `& d# s# a3 b, t" `. p: k" csets:
0 H# ]$ L) z& N1 X0 q1 R9 A9 i" o+ xshengshi/1..13/;
$ p; g& n7 n7 S7 D1 U( uzhibiao/1..6/:a;( _/ {8 F! E: n/ X* {9 g) v2 m% M3 [
links(shengshi,zhibiao):x;! t! x/ V  P* h0 S: _
endsets1 ?4 Q& l& S2 Y: y: j% f
data:
$ I6 X" H. Q. T+ f  f% Q1 a- mx=. |7 ]6 K+ T) }; W0 D9 J9 w) H% T
0.099599  0.104095  0.021186  0.489849  0.031419  0.399601
5 ~6 h* ^( v: K4 C0.012222  0.006320  0.035075  0.558043  0.253776  0.465170 7 b) H6 c3 S8 q. o8 W* r3 v
0.016260  0.004213  0.084423  0.631442  0.375174  0.535579
8 y8 @0 g- Z, U" Z4 K/ `0.004587  0.000000  0.086154  0.557522  0.210591  0.464721 : R% L8 T# {' W2 `8 [4 r
0.039880  0.015333  0.077921  1.000000  0.488760  0.888523
) A, R) t9 R0 S6 f) p8 S0.003290  0.002399  0.023441  0.299844  1.000000  0.217814 ( }& a% S1 Q5 Y
0.021750  0.023296  0.055727  0.635086  0.657561  0.538723
; a6 F% v3 {4 V4 z/ v1.000000  1.000000  1.000000  0.000000  0.000000  0.430589
6 f* w% H$ J) W7 J  ?: u- h! u; ?0.070963  0.277539  0.385947  0.153566  0.014757  0.000000 $ B" h  y7 ]1 {
0.038164  0.038373  0.024990  0.489849  0.031159  1.000000 . L" V1 \+ S+ w5 _
0.000000  0.008735  0.015259  0.288912  0.037009  0.040918 7 v2 N! U0 D5 l& N6 u- N: u  u2 n
0.057181  0.038508  0.254109  0.184279  0.093160  0.052246 $ R* I6 @6 J9 a: t5 ]+ [9 w
0.048211  0.009641  0.000000  0.387819  0.135834  0.118513;- n" Z& y# j* z# B7 o
enddata
5 x' D9 A1 J1 O5 R: [max=7 P8 y/ `' c5 w3 x5 G# g
(@sum(shengshi(i):2 s4 r: Z5 c1 K! H: @' Z3 u
(@sum(links(i,j):a(j)*x(i,j))-@sum(shengshi(i):
/ J% b! K9 R" W4 A+ Y6 @@sum(links(i,j):a(j)*x(i,j)))/13)^2)/12)^(1/2);
/ z: G$ E  c+ e
0 T. }4 p- l. b@sum(zhibiao(j):a(j)^2)=1;
  K0 H( o2 S' O7 ~9 T@for(zhibiao(j):a(j)>-1);) q5 N/ e  B3 p0 e+ f* R' b
@for(zhibiao(j):a(j)<1);3 C8 }3 C, m% T7 W
@for(zhibiao(j)free(a(j)));
  K1 H, l) H8 c* B+ D' X7 ]6 q) R4 dend ' B3 ?3 J  v; w+ Y5 t4 P) R
求不出全局最优解，请大侠帮帮忙！！！
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.
    . y+ [2 W" A: K% U  Objective value:                             0.7272518
    ; a. @' A- [# y% X) Y  Objective bound:                             0.7272527$ E6 U9 P  Y8 N% l/ U
      Infeasibilities:                              0.000000
    ( c# D3 v5 U2 M; i0 k  Extended solver steps:                           38843
    ) P* K5 R0 d/ ^  a9 d' Z( y! o  Total solver iterations:                       7561347
    & v" z9 s/ c0 C# o+ K' k( Z) V5 ?6 n7 s0 b
    2 b$ u8 h( L7 l! B$ y
                           Variable           Value        Reduced Cost
    ! s4 K, P9 v9 o) K+ S                          A( 1)       0.1248020            0.000000) |1 g) ~# B/ A* B' j/ S7 n/ H
                              A( 2)       0.1216772            0.000000
    " U8 f1 z- b2 M2 X                          A( 3)       0.1621721            0.000000  N. G9 j  |: [6 M+ A" x1 t& P' H$ ?5 E
                              A( 4)       0.6318675            0.000000' ]) Z" r' `# V" H. t+ _& _
                              A( 5)       0.3868056            0.000000
    ( S5 a) X! ]9 ^& _& E                          A( 6)       0.6280479            0.000000
    ! {0 F1 F9 P% F  m* P) Q0 s                       X( 1, 1)       0.9959900E-01        0.000000
    - x$ ~1 f' L! n% K  K+ E                       X( 1, 2)       0.1040950            0.000000
    . a- t9 O: f9 y0 D( J6 ?+ W& o                       X( 1, 3)       0.2118600E-01        0.0000003 t  W8 |4 P( d4 m( q' ]
                           X( 1, 4)       0.4898490            0.000000' E5 U" _9 e8 h/ k' @
                           X( 1, 5)       0.3141900E-01        0.000000
    , l3 x! L7 h1 T* ^                       X( 1, 6)       0.3996010            0.000000
    1 n, X6 z& h( N8 e4 N6 E                       X( 2, 1)       0.1222200E-01        0.000000% F0 H; e: r* R
                           X( 2, 2)       0.6320000E-02        0.000000! C' r: L4 e% u( l
                           X( 2, 3)       0.3507500E-01        0.000000
    ) q+ G: g6 Q2 j% O* n& F7 \6 ]                       X( 2, 4)       0.5580430            0.0000002 y/ a+ U9 w. m. n& Q7 x
                           X( 2, 5)       0.2537760            0.000000
    # d( p4 r/ ?- W' v+ n5 ~3 z; U                       X( 2, 6)       0.4651700            0.000000
    ' `; J1 u; u4 R3 f2 w                       X( 3, 1)       0.1626000E-01        0.000000% p; W# P% O2 {5 ^  c& }
                           X( 3, 2)       0.4213000E-02        0.000000
    8 w2 f, m3 I6 A- @$ X9 }                       X( 3, 3)       0.8442300E-01        0.0000005 J7 W9 o2 U: B1 S1 t
                           X( 3, 4)       0.6314420            0.000000
    0 }6 ]! O$ |% z" j- i: `! V6 G                       X( 3, 5)       0.3751740            0.0000006 @: a1 P8 W8 ]1 ^
                           X( 3, 6)       0.5355790            0.000000$ f: w7 D# O1 z# q7 K' j
                           X( 4, 1)       0.4587000E-02        0.000000
    $ n. n& p% p: ~) j0 F8 T  I' C                       X( 4, 2)        0.000000            0.000000
    8 y3 u9 {" R( O3 a                       X( 4, 3)       0.8615400E-01        0.0000008 K; R" U# x. s5 W& Z
                           X( 4, 4)       0.5575220            0.000000
    ; w; c* P2 q# b1 s  g: k                       X( 4, 5)       0.2105910            0.000000) f( S" X2 z  ^0 w' m# B6 w# k* X
                           X( 4, 6)       0.4647210            0.000000
    . _) g" |- }, b3 l! N, `/ H  ~                       X( 5, 1)       0.3988000E-01        0.0000004 t' S0 R9 K0 B0 |# S. ^, [
                           X( 5, 2)       0.1533300E-01        0.000000
    , A" N% c& {/ j) q                       X( 5, 3)       0.7792100E-01        0.000000: N6 n7 ^4 Z& R) z
                           X( 5, 4)        1.000000            0.000000
    ! K6 X$ G- M9 A                       X( 5, 5)       0.4887600            0.000000
    , y0 F, \. J: M7 n0 @& V" h' w                       X( 5, 6)       0.8885230            0.000000+ r+ m' `- D& l! {& [$ C) v
                           X( 6, 1)       0.3290000E-02        0.0000001 J- |- a6 Y, ~6 A- W$ u
                           X( 6, 2)       0.2399000E-02        0.000000
    1 @6 R- ?' x4 z: c0 T4 }                       X( 6, 3)       0.2344100E-01        0.000000- Y# K$ H1 C2 x" f2 q8 t
                           X( 6, 4)       0.2998440            0.000000
    0 s2 F& T- M- {3 d( b                       X( 6, 5)        1.000000            0.000000
    ' P2 o: S- o' i( x                       X( 6, 6)       0.2178140            0.000000
    / O+ [' E0 R/ o% {                       X( 7, 1)       0.2175000E-01        0.000000
    ) O4 ?# c# s8 x+ u8 [& u, O                       X( 7, 2)       0.2329600E-01        0.000000
    1 j9 }2 E1 ~$ P# c# J                       X( 7, 3)       0.5572700E-01        0.000000  _; L* i, J4 I4 V% E+ K
                           X( 7, 4)       0.6350860            0.0000000 [0 y: C7 A1 U
                           X( 7, 5)       0.6575610            0.000000
    3 Q  W( L! {$ f( G$ ?1 G                       X( 7, 6)       0.5387230            0.0000002 _( G# x6 @+ @
                           X( 8, 1)        1.000000            0.000000
    ; Q  x; z: E* [9 B                       X( 8, 2)        1.000000            0.000000% j: M; l5 Z$ V3 p6 B
                           X( 8, 3)        1.000000            0.0000009 v3 d/ ?' R5 G* v
                           X( 8, 4)        0.000000            0.000000
    - C( N; o# D3 S1 T( [                       X( 8, 5)        0.000000            0.0000000 \# p0 j( S2 o2 w, A
                           X( 8, 6)       0.4305890            0.000000
    8 Z: z7 y6 A3 r7 y/ S: O) k                       X( 9, 1)       0.7096300E-01        0.000000* P- A3 b( i* d4 I0 `
                           X( 9, 2)       0.2775390            0.000000& h! E( A8 |. x, C6 S2 x* g
                           X( 9, 3)       0.3859470            0.000000% z+ O$ ]! i: d7 B' _
                           X( 9, 4)       0.1535660            0.0000000 q9 S4 X2 i, U4 @' S; K  f
                           X( 9, 5)       0.1475700E-01        0.000000. f1 B5 N" O) Y8 {3 y) g0 u, e
                           X( 9, 6)        0.000000            0.000000
    % s1 |% w; Y( J. G5 y) k5 K                      X( 10, 1)       0.3816400E-01        0.000000
    5 }; }; x8 ~8 V# f                      X( 10, 2)       0.3837300E-01        0.000000  {6 |5 o; \9 H  x% K$ [* u2 W
                          X( 10, 3)       0.2499000E-01        0.000000
    . H% y4 T& E) i5 K" \                      X( 10, 4)       0.4898490            0.000000
    % L+ g. x5 X+ _) h                      X( 10, 5)       0.3115900E-01        0.000000
    8 C2 h# G3 D9 h# X$ _1 l( z5 z                      X( 10, 6)        1.000000            0.000000
    * i- C6 l5 N) X) i- C                      X( 11, 1)        0.000000            0.000000- p+ z$ ~7 @9 s5 H- C" I  s: P
                          X( 11, 2)       0.8735000E-02        0.000000
    - W# u/ P9 l" p9 s6 a& C3 D                      X( 11, 3)       0.1525900E-01        0.000000: F  \1 M+ C& k2 U- ]) }0 S
                          X( 11, 4)       0.2889120            0.000000
    ! u$ B+ Z1 V. e3 T                      X( 11, 5)       0.3700900E-01        0.000000
      R5 P/ L+ B+ @                      X( 11, 6)       0.4091800E-01        0.0000005 }' y0 T0 B* O2 h% A9 f: y
                          X( 12, 1)       0.5718100E-01        0.000000
    $ u8 E; d- l6 Y: R4 z                      X( 12, 2)       0.3850800E-01        0.0000002 @- T) ?! z  H  e5 M& a* q+ X
                          X( 12, 3)       0.2541090            0.000000- y1 M$ |: B7 d, r. }$ u' a1 a  B
                          X( 12, 4)       0.1842790            0.000000( F$ r# j: m3 b' f4 b/ u
                          X( 12, 5)       0.9316000E-01        0.0000001 B  U  Z' `& U5 z" }
                          X( 12, 6)       0.5224600E-01        0.0000001 K5 B" m. ]2 \5 O3 D
                          X( 13, 1)       0.4821100E-01        0.000000
    2 ]( g5 a9 w# m. q1 J# |* n                      X( 13, 2)       0.9641000E-02        0.000000
    3 D+ o. }4 |8 _6 M) ~                      X( 13, 3)        0.000000            0.000000
    ' n3 s# b5 h: x+ z  o# z3 z                      X( 13, 4)       0.3878190            0.000000  _. p  V9 F& G4 s# X
                          X( 13, 5)       0.1358340            0.000000
    6 V  o/ `& F5 H' Z$ B0 O                      X( 13, 6)       0.1185130            0.000000
    9 [' i) [7 I1 X1 ]) D/ {% ?2 ~/ ~& a! Y
                                Row    Slack or Surplus      Dual Price) p2 g3 R" W& g$ x  w0 m1 h
                                  1       0.7272518            1.000000  F; h( @& T$ j% E  S
                                  2        0.000000           0.3636259( h7 `8 z* v6 e% L; X/ P& \+ p( ]
                                  3        1.124802            0.000000
    7 a# t2 t9 T5 v2 X                              4        1.121677            0.0000007 ]; F9 o, m+ N. x
                                  5        1.162172            0.000000* N! n+ `2 p- A9 T1 K
                                  6        1.631867            0.000000
      l9 o% G  @( ?* _3 V( S                              7        1.386806            0.000000
    $ t* f7 w" K: U3 C                              8        1.628048            0.000000+ K7 t! E) m1 B, e9 \8 z6 h" W" @
                                  9       0.8751980            0.000000. O6 @# i& |: F- P* K* q% G+ I) W
                                 10       0.8783228            0.000000
    & r+ l9 j0 x7 n3 s2 d  H6 e  H9 e1 j                             11       0.8378279            0.000000# S$ _# u( {9 ~: X% J
                                 12       0.3681325            0.000000
    ) K/ o' p# r( q: c; T+ O                             13       0.6131944            0.000000- Z) M: k5 [+ w! N; _" P
                                 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
    但是 我设置全局求解时：两三秒就停止了
    2 \; J% Q* U8 @, t; S出现错误提示
    7 {- D8 y9 G  b& @+ h" q( m0 C5 z; t还有 上面的 全局的做大值还没有局部做大的大  为什么！！！
    回复

    使用道具 举报

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

    回顶部