快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

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

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

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

7

主题

5

听众

28

积分

升级  24.21%

该用户从未签到
跳转到指定楼层
1#
发表于 2009-9-5 07:49 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
model:
# D8 e% y# s5 ^% a3 O& Q$ ~* \sets:
9 A6 Y# o% k8 \% T' ishengshi/1..13/;
9 R4 {3 V8 q3 \9 N: zzhibiao/1..6/:a;
! W# A7 X/ g/ a% ]$ g7 K" [links(shengshi,zhibiao):x;
! ~4 Z+ C- A  I; zendsets
2 A$ F6 V9 a$ Y0 X9 Wdata:
" J3 `9 `: _) O$ n3 ~! M* Bx=- I/ `" g* a8 A7 O( W: y; a+ f) `
0.099599  0.104095  0.021186  0.489849  0.031419  0.399601 2 T9 a" q! ^" [, f2 v0 y
0.012222  0.006320  0.035075  0.558043  0.253776  0.465170 - {1 S, }" q. S; Z# U7 e
0.016260  0.004213  0.084423  0.631442  0.375174  0.535579 % h8 o7 ~$ Y; f) [( j+ g1 d
0.004587  0.000000  0.086154  0.557522  0.210591  0.464721
# q& k, e6 s% H) L7 e0.039880  0.015333  0.077921  1.000000  0.488760  0.888523
  M  g$ B# w0 P/ ^4 A0.003290  0.002399  0.023441  0.299844  1.000000  0.217814
& b0 Y' w) O( m# @" n' I2 y% r0.021750  0.023296  0.055727  0.635086  0.657561  0.538723 % e8 |! l; K6 l
1.000000  1.000000  1.000000  0.000000  0.000000  0.430589
3 G2 M: b4 |. c6 Q* D* J7 M- \0.070963  0.277539  0.385947  0.153566  0.014757  0.000000 % y6 S5 @3 c8 n7 G  Y5 B" D
0.038164  0.038373  0.024990  0.489849  0.031159  1.000000
4 m6 q" u7 d9 \2 {" Z. X2 e0.000000  0.008735  0.015259  0.288912  0.037009  0.040918
) B' ^6 P0 K8 j/ M. x) y0.057181  0.038508  0.254109  0.184279  0.093160  0.052246
. d7 y4 O/ v* E% o5 {0 G0.048211  0.009641  0.000000  0.387819  0.135834  0.118513;$ f1 |' A# z! i$ J. u
enddata5 {7 b8 ?8 f% a
max=
, |+ n4 L8 @$ c$ f(@sum(shengshi(i):5 _6 o8 k: W& g0 Z3 F  _' N
(@sum(links(i,j):a(j)*x(i,j))-@sum(shengshi(i):
" ^1 M1 @8 k+ I# O+ o2 t@sum(links(i,j):a(j)*x(i,j)))/13)^2)/12)^(1/2);
7 b$ d3 q8 z( j8 D, L9 S
' Y6 k" N4 {6 G7 R! V) V2 p@sum(zhibiao(j):a(j)^2)=1;
8 O6 N- _( S. i, d@for(zhibiao(j):a(j)>-1);
0 ~. A5 }# n; k2 K, Z@for(zhibiao(j):a(j)<1);. l8 N. v- V7 n" E% R8 k
@for(zhibiao(j)free(a(j)));2 S0 ~( |/ N5 h' U8 s  U) y
end 6 {) L# F, c, J! A( T
求不出全局最优解，请大侠帮帮忙！！！
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 f9 m) ^2 `, g7 z2 e: ~, N
      Objective value:                             0.7272518. a3 I% P: ~% R0 r  I* Q, [
      Objective bound:                             0.7272527
    ! v! _0 N$ v  E( N$ w6 {1 D+ l  Infeasibilities:                              0.000000
    * V6 p' {7 p; Y. u0 S( o  Extended solver steps:                           38843
    . a0 H9 c7 M4 F  Total solver iterations:                       75613474 V! B# B+ X( U7 D8 i8 Z

    * q0 L  L4 E! ?. l4 W# J5 @5 M% C6 ]9 H; V: \5 w+ P
                           Variable           Value        Reduced Cost
    % t0 k# K3 O5 s* K" C8 p% H2 B3 q* m                          A( 1)       0.1248020            0.000000
    3 a" d! a/ l! L  N- k1 I                          A( 2)       0.1216772            0.000000  C# b1 x5 A0 X3 x& z  H
                              A( 3)       0.1621721            0.000000
    / U* {5 x! u8 l; o) b& ?/ @                          A( 4)       0.6318675            0.000000
    ) b. z& o$ W6 ^+ F4 x2 L                          A( 5)       0.3868056            0.000000# P  q+ P, [% W4 e# n, Q7 p
                              A( 6)       0.6280479            0.0000002 j$ k7 l% I' Q. W& x, ]% B. j
                           X( 1, 1)       0.9959900E-01        0.000000: F. T5 q$ e* ~# j% p* J
                           X( 1, 2)       0.1040950            0.000000
    ! D$ W7 g) t' n: C  w                       X( 1, 3)       0.2118600E-01        0.000000
    * B* l* S, C" w: x) R: A1 P                       X( 1, 4)       0.4898490            0.000000& t6 N) n( F' {1 U3 I
                           X( 1, 5)       0.3141900E-01        0.000000( E9 W7 P# [8 M1 q
                           X( 1, 6)       0.3996010            0.000000
    9 ^" @5 t* i4 \. V) B& N                       X( 2, 1)       0.1222200E-01        0.000000$ ^% W8 ~, W- g- w# E
                           X( 2, 2)       0.6320000E-02        0.0000009 v2 R1 n7 A1 k0 U, c5 D$ A
                           X( 2, 3)       0.3507500E-01        0.000000
    4 X* z, h9 A" W+ J7 ^  t4 [                       X( 2, 4)       0.5580430            0.000000
    " U6 I& U6 f, y% Y4 h7 r# f                       X( 2, 5)       0.2537760            0.000000
    5 P9 H& G6 H4 P7 v, E                       X( 2, 6)       0.4651700            0.000000: W; m; T: R3 ~6 l9 ~/ e* A( O& O( v
                           X( 3, 1)       0.1626000E-01        0.000000
    0 Z: L2 l3 [: F& T8 ^. F* |4 ~) G                       X( 3, 2)       0.4213000E-02        0.000000
    2 R4 o5 h4 O9 E2 Z                       X( 3, 3)       0.8442300E-01        0.000000
    8 n+ k0 ^/ o2 ~$ k                       X( 3, 4)       0.6314420            0.000000
    / B5 {! Z5 Z- p7 r                       X( 3, 5)       0.3751740            0.0000004 J1 T$ h- [3 [, E
                           X( 3, 6)       0.5355790            0.000000
    . Z  m$ H( g! {; H! ]3 k                       X( 4, 1)       0.4587000E-02        0.000000# \! J! P$ i6 N- m& _  i2 h
                           X( 4, 2)        0.000000            0.000000
    " g4 k& j6 y. V# ?                       X( 4, 3)       0.8615400E-01        0.000000
    ' c- h6 U. _. j& g3 f3 O- l( u                       X( 4, 4)       0.5575220            0.000000
    " z6 X0 ^' p& C1 N4 }$ O                       X( 4, 5)       0.2105910            0.000000
    * @' H+ c& u( `* _/ r3 n/ }8 ?& J                       X( 4, 6)       0.4647210            0.000000
    4 `* Z8 L% r, `' F3 m7 `                       X( 5, 1)       0.3988000E-01        0.000000
    * K9 r  `" h1 W* L+ }# U  U9 k3 ^5 i                       X( 5, 2)       0.1533300E-01        0.0000006 p8 D! f0 V1 M+ G, E& a5 v5 f
                           X( 5, 3)       0.7792100E-01        0.0000002 b, i% {& p" W
                           X( 5, 4)        1.000000            0.000000) h" q* q) A  b
                           X( 5, 5)       0.4887600            0.000000
    ; A+ j3 C8 y# `7 s8 B                       X( 5, 6)       0.8885230            0.000000
    9 u/ q$ L+ e7 {4 I( n6 h5 ~+ D                       X( 6, 1)       0.3290000E-02        0.000000
    * `+ \2 J2 l" M. b- s# }                       X( 6, 2)       0.2399000E-02        0.000000
    - A/ D2 q. `' a                       X( 6, 3)       0.2344100E-01        0.000000
    4 O5 \% K/ S. u. a                       X( 6, 4)       0.2998440            0.0000001 }7 w5 A3 |. c4 s  m; \. N5 R
                           X( 6, 5)        1.000000            0.0000003 l0 A6 t- v1 j9 E
                           X( 6, 6)       0.2178140            0.000000# L$ h1 S7 o3 e/ L8 i! X
                           X( 7, 1)       0.2175000E-01        0.000000. J  h$ Q, q. L! Q2 E/ j9 Z
                           X( 7, 2)       0.2329600E-01        0.000000* B8 F; x- A8 ?# P. ^* m' t
                           X( 7, 3)       0.5572700E-01        0.000000
    8 b: w+ u4 a8 |& |                       X( 7, 4)       0.6350860            0.0000002 M& t  I" _9 t- R, e3 v
                           X( 7, 5)       0.6575610            0.0000004 z! A- J& A+ J" Y3 W) Y
                           X( 7, 6)       0.5387230            0.0000002 d% a. @6 j7 s& |8 M
                           X( 8, 1)        1.000000            0.000000* ^0 \" `& p, c& B  z0 h
                           X( 8, 2)        1.000000            0.000000
    ; k! F* C5 g) r& ?6 g. B) |                       X( 8, 3)        1.000000            0.000000
    / H) |2 F9 n9 W/ t4 m+ `                       X( 8, 4)        0.000000            0.000000
    : E0 q$ V& S6 t1 G8 `                       X( 8, 5)        0.000000            0.000000
    0 J9 H; I( M( Z4 v  G( b5 R                       X( 8, 6)       0.4305890            0.000000
    ' l: a2 C/ Q3 B$ C7 M+ P                       X( 9, 1)       0.7096300E-01        0.000000
    ! S6 ^# H2 l4 `. l, P                       X( 9, 2)       0.2775390            0.0000003 d, _* v* n- g
                           X( 9, 3)       0.3859470            0.000000
    2 ?) ]4 {, q8 f7 r: q                       X( 9, 4)       0.1535660            0.000000/ D* k8 q* h) v8 ^% t
                           X( 9, 5)       0.1475700E-01        0.000000) J3 L: v! t' e# }. e0 k) `  [
                           X( 9, 6)        0.000000            0.000000& ?  X8 {1 J" _4 L& d+ Y% N$ |1 k
                          X( 10, 1)       0.3816400E-01        0.000000' M* j0 v. f" q
                          X( 10, 2)       0.3837300E-01        0.000000
    & c( |4 k' h# T. E% s                      X( 10, 3)       0.2499000E-01        0.000000) N' |  V' f" I  T3 L: Z: }
                          X( 10, 4)       0.4898490            0.000000
    . |$ B8 V, T- d4 e/ F. s                      X( 10, 5)       0.3115900E-01        0.000000$ B5 R8 c( A4 `4 G2 R- u$ Z, r
                          X( 10, 6)        1.000000            0.000000" Y; c. E5 d' h. y6 [; E
                          X( 11, 1)        0.000000            0.000000& V8 V3 m2 j7 C5 Y0 @
                          X( 11, 2)       0.8735000E-02        0.0000008 z' R7 y5 c2 h6 v$ @/ x3 o! @
                          X( 11, 3)       0.1525900E-01        0.000000  s. Q; q* i" d8 Q5 {
                          X( 11, 4)       0.2889120            0.000000
    / p" e/ A" A* P. V, @+ M+ k                      X( 11, 5)       0.3700900E-01        0.000000
    ) T6 n8 b' `1 ~                      X( 11, 6)       0.4091800E-01        0.000000; q+ s8 \, o3 B; Z
                          X( 12, 1)       0.5718100E-01        0.000000
    ( P0 c  R1 Y+ _1 V2 ~" e& K                      X( 12, 2)       0.3850800E-01        0.000000
    $ t3 s; r" i3 J+ ]' O- n                      X( 12, 3)       0.2541090            0.000000
    2 |1 o2 Y  H2 a1 j- E; J9 u; L                      X( 12, 4)       0.1842790            0.000000
    3 H" @4 Q. V1 C$ p! O2 P# W                      X( 12, 5)       0.9316000E-01        0.000000
    5 D. r: ?6 d! M$ @                      X( 12, 6)       0.5224600E-01        0.000000
    $ {* z# b, y2 M                      X( 13, 1)       0.4821100E-01        0.000000
    # L( d3 e! V' y7 j; }  H                      X( 13, 2)       0.9641000E-02        0.000000% a3 M; I2 T& b& V: C* y
                          X( 13, 3)        0.000000            0.000000
    2 O# F" V2 j# E! w                      X( 13, 4)       0.3878190            0.0000002 P  q0 e" c4 R0 {
                          X( 13, 5)       0.1358340            0.000000
    8 |. f5 I# w; M0 A  c                      X( 13, 6)       0.1185130            0.000000
    - [3 F& Z: |7 M; X
    $ t8 k' W/ n' p( u# o3 ?2 x* ~                            Row    Slack or Surplus      Dual Price
    , T; o+ i3 a0 q6 j! V                              1       0.7272518            1.000000
    * B1 _+ o% W  a8 c* g( W                              2        0.000000           0.3636259! k8 j4 h5 ~7 t9 p& ?
                                  3        1.124802            0.0000001 j1 [. Y7 R; q( q- f; O5 N5 q6 Y9 F
                                  4        1.121677            0.0000007 L" j/ q& _% ~& l% U6 ?
                                  5        1.162172            0.000000
    ; A6 V' j! [, x! b  _% O                              6        1.631867            0.000000" H6 m- B0 o5 A  P& b% }
                                  7        1.386806            0.000000
    5 g1 S2 t1 c( }                              8        1.628048            0.000000
    : I( M7 H9 _1 n4 K                              9       0.8751980            0.000000
    ' ^3 L8 o. W; r                             10       0.8783228            0.000000
    % z1 l7 ^9 d6 `, h4 p% v  w                             11       0.8378279            0.000000
    ( g  n% ^7 K& |9 y- N! q0 U                             12       0.3681325            0.000000* U, ^! m$ n( f5 g
                                 13       0.6131944            0.000000
    ) u+ f+ o$ e$ E0 @                             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
    但是 我设置全局求解时：两三秒就停止了
    . s4 L3 h. s2 o" u$ Q出现错误提示1 t5 u* c4 O; u2 |/ c
    还有 上面的 全局的做大值还没有局部做大的大  为什么！！！
    回复

    使用道具 举报

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

    回顶部