快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

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

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

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

7

主题

5

听众

28

积分

升级  24.21%

该用户从未签到
跳转到指定楼层
1#
发表于 2009-9-5 07:49 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
model:
) a2 O0 H2 o2 p* g4 ksets:
+ G  M- m$ ]; t3 @3 _9 p) yshengshi/1..13/;
' H  A5 H. i( h$ szhibiao/1..6/:a;
7 t% ^# G7 b# T" G/ ]& J0 zlinks(shengshi,zhibiao):x;1 r; g: ]9 K) M: [0 M' b: r  s) [
endsets
0 n; Q  j& l) A0 Ndata:( V3 `9 Y& c9 S9 {5 H1 U
x=- P& s9 E9 G% \5 L: N
0.099599  0.104095  0.021186  0.489849  0.031419  0.399601 ; p& p$ B! o+ d  p0 m) H. g
0.012222  0.006320  0.035075  0.558043  0.253776  0.465170
: ~1 x4 Q0 \4 [" g' z0.016260  0.004213  0.084423  0.631442  0.375174  0.535579
. G$ [# F% \/ ~5 J0.004587  0.000000  0.086154  0.557522  0.210591  0.464721 8 K5 W! [0 m6 M: Y5 b. L" {: ?+ {
0.039880  0.015333  0.077921  1.000000  0.488760  0.888523 - T/ c/ {$ H6 N+ ~* Y
0.003290  0.002399  0.023441  0.299844  1.000000  0.217814
6 z7 I+ ]' \  ^$ h* H0.021750  0.023296  0.055727  0.635086  0.657561  0.538723 8 y0 D. \; }& S
1.000000  1.000000  1.000000  0.000000  0.000000  0.430589
$ E% I' }4 {4 z6 \  I2 g0.070963  0.277539  0.385947  0.153566  0.014757  0.000000
4 b  t4 L' [* W1 x! {0.038164  0.038373  0.024990  0.489849  0.031159  1.000000
9 j( Z$ {; F' X! u6 ~0.000000  0.008735  0.015259  0.288912  0.037009  0.040918
7 M" p  W6 l& \% `: f0.057181  0.038508  0.254109  0.184279  0.093160  0.052246
5 r+ K, Q0 Z% h1 L1 y& o0.048211  0.009641  0.000000  0.387819  0.135834  0.118513;: L' y! e" ~( |# ]. A
enddata% H. T* q: v+ N% [
max=! i' e5 Y' |* W1 F' w
(@sum(shengshi(i):/ t6 r) a0 m: x- h: _. Z6 r
(@sum(links(i,j):a(j)*x(i,j))-@sum(shengshi(i):* c- y8 x( p3 q1 H4 I
@sum(links(i,j):a(j)*x(i,j)))/13)^2)/12)^(1/2);
+ C5 I$ W( h8 o3 N3 J# A- m* I) Z0 ^8 i3 `0 h
@sum(zhibiao(j):a(j)^2)=1;
8 q' A/ m- S( w+ `@for(zhibiao(j):a(j)>-1);% {+ U% Z2 n2 g: q- J9 `
@for(zhibiao(j):a(j)<1);" G( u) c  T2 d4 o. l! ]8 r
@for(zhibiao(j)free(a(j)));5 q6 S$ y3 C  }4 j4 y
end & E( v# P  m9 ^0 N
求不出全局最优解，请大侠帮帮忙！！！
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.
    3 M0 j. G0 p- g  Objective value:                             0.72725187 n" j5 h* x: E2 c
      Objective bound:                             0.72725275 V5 k0 R4 K% e+ @5 F: v& R3 q
      Infeasibilities:                              0.000000
    : p9 o! R% u% ~+ [  Extended solver steps:                           38843
    ; }# s; ?; z+ c, [; ?  [& A  Total solver iterations:                       7561347
    5 a9 }) V" y7 W$ `& ]5 f, u
    ! J! I% d9 ?( p* e$ O: f, v& v4 r: _; l
                           Variable           Value        Reduced Cost8 f) {# Z, N* k! g5 I
                              A( 1)       0.1248020            0.000000
    8 S; f8 c4 d& X6 p. K( `0 S                          A( 2)       0.1216772            0.000000) W' z" }9 Q( O' D; x
                              A( 3)       0.1621721            0.000000  B3 f' u' f3 D5 C/ M( N+ x$ k! B2 b
                              A( 4)       0.6318675            0.0000000 G* s4 H$ W3 S, {6 m' K# Y
                              A( 5)       0.3868056            0.000000
    0 a* I: }* k( v                          A( 6)       0.6280479            0.000000' J7 a) G- Y& T2 M
                           X( 1, 1)       0.9959900E-01        0.0000005 E7 j- K/ C- U
                           X( 1, 2)       0.1040950            0.0000001 f4 R8 ~4 X( ^, h5 g) b3 \
                           X( 1, 3)       0.2118600E-01        0.0000009 z0 m% l1 d; @; c& V( c
                           X( 1, 4)       0.4898490            0.000000
    , h  n8 @3 A! |/ p3 d# Y9 F                       X( 1, 5)       0.3141900E-01        0.000000+ l$ C3 u# h/ e$ U! m( j
                           X( 1, 6)       0.3996010            0.000000
    3 `" `+ B+ ^5 [/ q                       X( 2, 1)       0.1222200E-01        0.0000006 P* a# H, w6 ]/ X5 D; H; e
                           X( 2, 2)       0.6320000E-02        0.000000
    . k! M- q- a4 e                       X( 2, 3)       0.3507500E-01        0.000000
    " x# `: M" k7 W5 N8 h                       X( 2, 4)       0.5580430            0.000000
    % y; N  _6 C6 k                       X( 2, 5)       0.2537760            0.000000% b! `# Q7 p. a, ~9 I
                           X( 2, 6)       0.4651700            0.000000
    ( u6 y5 l9 c/ i                       X( 3, 1)       0.1626000E-01        0.000000! l. n3 ^( |( q0 W# q; h. L7 k
                           X( 3, 2)       0.4213000E-02        0.000000, m1 b3 U3 h8 }( G
                           X( 3, 3)       0.8442300E-01        0.000000
    1 x$ e/ W% K: U3 ]( ~/ v                       X( 3, 4)       0.6314420            0.000000: d0 `; v6 P& D8 W2 X$ J; G
                           X( 3, 5)       0.3751740            0.000000
    / _/ }( }/ _9 r5 @0 [0 Q) [- `0 M2 m                       X( 3, 6)       0.5355790            0.000000
    2 Z" Y' z2 a4 f# ]$ j' y9 Q  I                       X( 4, 1)       0.4587000E-02        0.000000
    ; B5 X0 h) t* R( Y                       X( 4, 2)        0.000000            0.000000
    : W# B, `; Y1 I% M+ z, o                       X( 4, 3)       0.8615400E-01        0.000000
    ; N) V6 D' T0 y6 E5 C( z5 v                       X( 4, 4)       0.5575220            0.000000
    2 V4 n" W1 M* W2 f+ }# i! r0 U                       X( 4, 5)       0.2105910            0.000000
    4 V6 @" p; i2 j! P! E6 [/ g3 r                       X( 4, 6)       0.4647210            0.0000000 D7 @/ |# |" d6 e! T! E
                           X( 5, 1)       0.3988000E-01        0.0000005 ~5 S9 [# P: C+ o) r$ ]
                           X( 5, 2)       0.1533300E-01        0.000000
    ( V# C1 g1 g& o. r! H, t, p                       X( 5, 3)       0.7792100E-01        0.000000
    ; I" N8 _4 _$ S. R1 {# L2 t                       X( 5, 4)        1.000000            0.000000
    : T* @8 M: ]0 D" c. H                       X( 5, 5)       0.4887600            0.000000
    ) O7 H! h4 _- H3 U                       X( 5, 6)       0.8885230            0.000000; c! A9 F+ L6 L$ D, {8 G. a8 k
                           X( 6, 1)       0.3290000E-02        0.0000009 x7 x1 K& G7 v9 T; c' l+ f( ]
                           X( 6, 2)       0.2399000E-02        0.0000008 S' w/ L$ w! |. V% c, V% X
                           X( 6, 3)       0.2344100E-01        0.000000
      Y' w& W# t( `9 g  @/ R- v                       X( 6, 4)       0.2998440            0.0000003 B& h7 N2 I1 @9 Q
                           X( 6, 5)        1.000000            0.000000
    : T7 [- l+ j! I                       X( 6, 6)       0.2178140            0.000000
    * V+ S2 d3 `, x                       X( 7, 1)       0.2175000E-01        0.000000
    5 a1 p0 F: y( B& U                       X( 7, 2)       0.2329600E-01        0.000000+ j  T7 r0 Q% ^& L  [
                           X( 7, 3)       0.5572700E-01        0.000000
    2 N7 I4 s) p! ^5 a  z4 Y8 {# F7 ~                       X( 7, 4)       0.6350860            0.000000
    + J( v5 l+ P! [# x, f                       X( 7, 5)       0.6575610            0.000000
    6 m* a( n3 g7 |) F* _                       X( 7, 6)       0.5387230            0.000000; z1 X) P% ~9 h- _! I/ o; u
                           X( 8, 1)        1.000000            0.000000
    $ a/ k- l4 l. f+ y7 Z                       X( 8, 2)        1.000000            0.000000
    $ x& }/ {5 r4 v  U                       X( 8, 3)        1.000000            0.000000
    . w9 x& j3 I8 K) q, C$ x, _* c0 I                       X( 8, 4)        0.000000            0.000000) `# o; k& S2 c9 x' L8 [2 t
                           X( 8, 5)        0.000000            0.000000' j' x9 W8 g* M3 s% l7 C
                           X( 8, 6)       0.4305890            0.000000
      g. m. R4 G  o                       X( 9, 1)       0.7096300E-01        0.000000; [( c, R  E9 s8 r. D  D7 H- O
                           X( 9, 2)       0.2775390            0.0000003 @+ D* E0 U+ K. y$ y) n9 p
                           X( 9, 3)       0.3859470            0.000000
    3 o, y/ @1 w& v  x  X                       X( 9, 4)       0.1535660            0.000000
      e8 i: U' a2 N( O# O, [. W* S                       X( 9, 5)       0.1475700E-01        0.000000
    . R9 e/ y2 g" A& b" l" L9 [                       X( 9, 6)        0.000000            0.000000" m' |2 |5 N; f2 s' b/ T# ^; f
                          X( 10, 1)       0.3816400E-01        0.000000
    % P4 ]) \% W9 s  M                      X( 10, 2)       0.3837300E-01        0.0000001 w6 S- b8 v0 H3 @8 x: I. x# G
                          X( 10, 3)       0.2499000E-01        0.000000+ T1 U, T) y9 L, D6 u+ K
                          X( 10, 4)       0.4898490            0.000000
    5 o2 K* b# N  t                      X( 10, 5)       0.3115900E-01        0.000000
    6 T/ K0 v  `! j" y/ ]% [8 D                      X( 10, 6)        1.000000            0.0000006 ^* F2 C3 V" M, K% M  x, A2 L! X
                          X( 11, 1)        0.000000            0.000000( {$ G- f* k4 z; k0 U2 w
                          X( 11, 2)       0.8735000E-02        0.000000/ ]) N  c2 q) u
                          X( 11, 3)       0.1525900E-01        0.000000
    1 y/ C# U8 z  P$ E( T& s% j, u                      X( 11, 4)       0.2889120            0.000000
    6 f; Z; W: k* [& q2 }" W# U                      X( 11, 5)       0.3700900E-01        0.000000
    & K, Z8 T' b" \4 Z$ z                      X( 11, 6)       0.4091800E-01        0.000000
    - e) Q8 A& {2 o, {7 {                      X( 12, 1)       0.5718100E-01        0.000000
    * l, a+ V" @! D                      X( 12, 2)       0.3850800E-01        0.000000
    / O9 l4 m' o) l) B- }: E* ^                      X( 12, 3)       0.2541090            0.000000+ K% p4 C% I3 Z- n- v4 ?/ {$ _
                          X( 12, 4)       0.1842790            0.000000
    4 H' M3 W: ]: A% {                      X( 12, 5)       0.9316000E-01        0.000000
    ) B- Y6 M( V9 R) m+ j, u                      X( 12, 6)       0.5224600E-01        0.000000  N$ M4 s5 `# [9 r/ U7 e- f
                          X( 13, 1)       0.4821100E-01        0.000000
    8 K# y: Y9 y- C* V' a                      X( 13, 2)       0.9641000E-02        0.000000: x7 D& y3 I: |" A: `8 h
                          X( 13, 3)        0.000000            0.000000
    0 }- n, H; d* a                      X( 13, 4)       0.3878190            0.000000& r+ p! B, |. \% v$ d$ Y
                          X( 13, 5)       0.1358340            0.000000
    , H( ]( ?* D. e                      X( 13, 6)       0.1185130            0.0000006 y( M; c, v4 [( a& B2 M9 K' Z: @/ u

    - v/ M* `) e1 ~& c. M% i                            Row    Slack or Surplus      Dual Price# Y- Y" P0 G3 m/ n. B* p7 c. d4 v
                                  1       0.7272518            1.000000. W! o  t6 J  A3 Y5 k6 D
                                  2        0.000000           0.3636259
    . Z# e) H8 E! {' R1 C% ~6 }                              3        1.124802            0.0000008 S" `, S! `5 M, A/ x
                                  4        1.121677            0.0000001 N  a9 a0 j  r" {5 z
                                  5        1.162172            0.000000; y! v+ p$ e8 Z: x. M$ F# V
                                  6        1.631867            0.000000
    2 _& a, X* Y2 C: y! l                              7        1.386806            0.000000% I3 d" `# v1 ?  P; x
                                  8        1.628048            0.000000" r) V2 F# C# f: ?5 Q! m
                                  9       0.8751980            0.000000
    ! [* C4 W1 }, R+ E- O                             10       0.8783228            0.000000
    7 N* T* {3 \- R5 b0 @) r- R+ N- e                             11       0.8378279            0.000000: h$ I9 s' T; i& Q
                                 12       0.3681325            0.000000; U! o4 F6 _2 C1 D: M
                                 13       0.6131944            0.000000, n3 U+ G4 o0 I. f% a7 U; W# d
                                 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
    但是 我设置全局求解时：两三秒就停止了
    0 Y% b' x/ q/ R" I  |' j出现错误提示
    7 m- x4 u; z  r6 F# k# P还有 上面的 全局的做大值还没有局部做大的大  为什么！！！
    回复

    使用道具 举报

    返回列表 发新帖
    高级模式
    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, 2025-8-22 03:29 , Processed in 0.921353 second(s), 105 queries .

    回顶部