快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

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

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

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

7

主题

5

听众

28

积分

升级  24.21%

该用户从未签到
跳转到指定楼层
1#
发表于 2009-9-5 07:49 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
model:* l0 ~: J6 y, r* j  R
sets:- t, l  A* f( h% l) V6 h
shengshi/1..13/;/ O" P; I7 f2 S8 V7 v
zhibiao/1..6/:a;5 e1 l, K) h/ x3 _' v
links(shengshi,zhibiao):x;
# O' C: Q* e+ g+ ~* F6 P3 l8 G# xendsets
. _4 M" `" [! tdata:2 r# H% u3 u" c+ s$ W
x=
7 `/ H" C  ~7 ^0.099599  0.104095  0.021186  0.489849  0.031419  0.399601
" s) G) u7 i, e, ~2 A1 S7 M8 d5 F0.012222  0.006320  0.035075  0.558043  0.253776  0.465170 5 D" D( Z' f* S- h6 x
0.016260  0.004213  0.084423  0.631442  0.375174  0.535579 ) _( O$ s8 C5 j6 h  j# Y/ t4 _
0.004587  0.000000  0.086154  0.557522  0.210591  0.464721 ( p2 A+ R" Y$ \0 V# y% ]
0.039880  0.015333  0.077921  1.000000  0.488760  0.888523 " N! \# `. Y# v% B9 b4 e- q
0.003290  0.002399  0.023441  0.299844  1.000000  0.217814
; Y* m1 E+ O( v' _; f0.021750  0.023296  0.055727  0.635086  0.657561  0.538723
: ~/ D3 h% h7 C- t9 ]6 f. W1.000000  1.000000  1.000000  0.000000  0.000000  0.430589 / c3 H0 N4 z  e
0.070963  0.277539  0.385947  0.153566  0.014757  0.000000   _) V' z; a6 a1 W
0.038164  0.038373  0.024990  0.489849  0.031159  1.000000 8 }, z! g5 a+ T
0.000000  0.008735  0.015259  0.288912  0.037009  0.040918
- L1 z) q5 ?1 w; u, ~1 z0.057181  0.038508  0.254109  0.184279  0.093160  0.052246
* \( t. L' s! D0.048211  0.009641  0.000000  0.387819  0.135834  0.118513;1 O4 Z( `0 F3 S. ]- c. z8 a
enddata
, o* R- V$ }0 D; `, ?# Emax=
0 T& _: W. _7 I(@sum(shengshi(i):
* z% Q6 A; G6 a6 k$ F+ H) {0 p(@sum(links(i,j):a(j)*x(i,j))-@sum(shengshi(i):' z/ G! C9 J/ Q0 A" R( n
@sum(links(i,j):a(j)*x(i,j)))/13)^2)/12)^(1/2);1 s8 m" h' h( {9 |! e3 T
- q1 b' D! [. |+ S) ?* v1 W* \
@sum(zhibiao(j):a(j)^2)=1;0 x. i  C  i. O' o5 K. C
@for(zhibiao(j):a(j)>-1);% X: x6 C) g+ ~& L  }& {4 L
@for(zhibiao(j):a(j)<1);! ~- A" T. f: }# [4 ]8 c
@for(zhibiao(j)free(a(j)));  g4 l: o# I# V% B0 M; a0 D
end 8 i7 o  @% f8 f" w# e
求不出全局最优解，请大侠帮帮忙！！！
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.
    6 o0 V+ D1 X5 ~1 T7 [% F! S, {  Objective value:                             0.7272518
    ' J' S9 J- f4 s1 X  Objective bound:                             0.72725278 k2 N3 y0 l$ x/ s0 \: [
      Infeasibilities:                              0.000000
      G) q1 i5 i) `9 r  Extended solver steps:                           38843- Z) I$ ^9 Y: Z9 J1 {
      Total solver iterations:                       7561347
    2 C& ]! [: v9 ^" A0 k) D0 N6 G6 ~8 e* n5 q" K
    * s% Q/ c7 C, d# a
                           Variable           Value        Reduced Cost
    1 g) V! w) m' F# B. K                          A( 1)       0.1248020            0.0000000 s/ [8 D: h* ]6 w; A" X  g% H
                              A( 2)       0.1216772            0.000000  U7 T" W3 L0 ?8 M; j
                              A( 3)       0.1621721            0.000000: W6 w# a* e2 f1 w
                              A( 4)       0.6318675            0.000000% k( ~. {, t& a" A1 z5 G% S  ~
                              A( 5)       0.3868056            0.000000
    & y- I0 R" B7 _1 n% R                          A( 6)       0.6280479            0.000000
    ' X% Z5 D/ L, O7 \: j3 ?8 e                       X( 1, 1)       0.9959900E-01        0.000000
    4 r# w! I. S& b* \  T                       X( 1, 2)       0.1040950            0.000000
    5 X8 e! @2 ~9 U) a5 T- p                       X( 1, 3)       0.2118600E-01        0.000000
    ; b! i& I' F3 x; W: i3 i; O                       X( 1, 4)       0.4898490            0.000000" g; u1 H: C8 V, H# P" C$ S
                           X( 1, 5)       0.3141900E-01        0.000000
      A9 N+ C4 J5 W                       X( 1, 6)       0.3996010            0.000000  A% D: u& R% Q* J* I' J; s9 m* @9 T, _
                           X( 2, 1)       0.1222200E-01        0.000000) J( C( E: n7 d- d" F& t
                           X( 2, 2)       0.6320000E-02        0.000000! N. |6 F2 {" V. N
                           X( 2, 3)       0.3507500E-01        0.000000" ]4 d$ M  X# E; o; o& o
                           X( 2, 4)       0.5580430            0.0000005 _8 k- ~8 `: E2 V. A
                           X( 2, 5)       0.2537760            0.000000
    6 A. z# J; [5 @( T+ _4 j# i                       X( 2, 6)       0.4651700            0.000000
    5 f2 B/ I& [9 u2 r! u6 x                       X( 3, 1)       0.1626000E-01        0.000000' s' M4 p& D6 W, V
                           X( 3, 2)       0.4213000E-02        0.000000$ O. S! T7 R2 j& t% Y6 x
                           X( 3, 3)       0.8442300E-01        0.000000% Z' f, V0 s  T2 m0 |
                           X( 3, 4)       0.6314420            0.000000
    8 @; x% u. u" R2 v, m: c6 Q& e# N$ t* w                       X( 3, 5)       0.3751740            0.000000
    & x$ Q( n7 A* n0 E! I. |: O" I                       X( 3, 6)       0.5355790            0.000000$ v' }+ N* F5 x' S2 m; L, i
                           X( 4, 1)       0.4587000E-02        0.000000
    / k" S8 t: W# L4 G: q: i! h                       X( 4, 2)        0.000000            0.000000, V, S1 n* V2 ]. J9 N$ ]; h. S
                           X( 4, 3)       0.8615400E-01        0.000000& Z# X7 s. ~% r9 `
                           X( 4, 4)       0.5575220            0.000000
    ; p# v0 Y0 w$ g: E4 U: D                       X( 4, 5)       0.2105910            0.000000
    / h0 o7 U* H, c+ w; L                       X( 4, 6)       0.4647210            0.000000
    # n$ R7 O1 r; {/ S% n% k) [                       X( 5, 1)       0.3988000E-01        0.000000
    3 D' ?4 I) t/ n$ |) u                       X( 5, 2)       0.1533300E-01        0.0000000 O" y& s# Q9 P4 ]9 q
                           X( 5, 3)       0.7792100E-01        0.0000006 L1 w. ]( o7 N' p# E
                           X( 5, 4)        1.000000            0.000000
    % I- t5 @: l, o+ k8 H                       X( 5, 5)       0.4887600            0.000000* H  U& b3 o* P8 B' C
                           X( 5, 6)       0.8885230            0.0000007 A3 T  ^& V4 H) f1 p" ]: ~, X0 m
                           X( 6, 1)       0.3290000E-02        0.000000" b* w' j2 ^. ]
                           X( 6, 2)       0.2399000E-02        0.0000005 p6 t6 w& o$ q0 t3 L+ I1 E. d7 ]
                           X( 6, 3)       0.2344100E-01        0.0000004 L$ d) H$ h/ t1 {9 d. m4 a% U
                           X( 6, 4)       0.2998440            0.000000
    4 U" x- p+ Q% o, h, s7 v* h                       X( 6, 5)        1.000000            0.000000
    4 W3 I7 E( n0 s8 q% f/ I                       X( 6, 6)       0.2178140            0.000000. ~  u0 P" N7 M6 g6 R4 {
                           X( 7, 1)       0.2175000E-01        0.000000
    7 F" T" H% g# P! t                       X( 7, 2)       0.2329600E-01        0.000000
    5 d& x0 z$ v8 |& T' a8 Y3 h                       X( 7, 3)       0.5572700E-01        0.000000
    6 d( N8 D3 [& W0 e                       X( 7, 4)       0.6350860            0.000000- u3 Y4 U! k/ Y+ Y3 E* U* c( Z
                           X( 7, 5)       0.6575610            0.000000
    ! I& w+ V9 o6 t* k) ?$ w6 b0 N                       X( 7, 6)       0.5387230            0.000000
    1 u$ r0 q: h; ~3 _' z1 G                       X( 8, 1)        1.000000            0.000000
    2 V3 B% b% u/ m% [% C( k                       X( 8, 2)        1.000000            0.000000
    2 `. V+ s2 A0 P/ l2 r. w                       X( 8, 3)        1.000000            0.000000
    # t- a; ^1 s" K: d, t  E/ ~( h                       X( 8, 4)        0.000000            0.000000
    8 h) D3 e6 m7 I, h* s) e- P                       X( 8, 5)        0.000000            0.000000
    * K4 g/ }% J  N7 x                       X( 8, 6)       0.4305890            0.000000
    - C. ]4 x& N) ]/ ~$ P0 M4 A3 X                       X( 9, 1)       0.7096300E-01        0.0000006 b! H0 u6 {' ?; C2 X0 x$ [
                           X( 9, 2)       0.2775390            0.000000
    ! z4 U+ r6 V" q% s- e+ W3 ^  n+ E                       X( 9, 3)       0.3859470            0.000000
    7 U1 P! X& a% e( U. {                       X( 9, 4)       0.1535660            0.000000- S" v  C. ^0 i1 M/ K4 e
                           X( 9, 5)       0.1475700E-01        0.000000
    ) K, ]( U. `/ S0 o" \                       X( 9, 6)        0.000000            0.000000% L% R6 \  y- m
                          X( 10, 1)       0.3816400E-01        0.000000
    $ H$ j6 ~( ]) t; }* B9 z6 r                      X( 10, 2)       0.3837300E-01        0.0000005 u, s! q3 T4 X' r2 ^
                          X( 10, 3)       0.2499000E-01        0.000000
    4 }1 `: N6 H1 Z' O) ?/ D! G0 f                      X( 10, 4)       0.4898490            0.000000
    6 r( d. p& ]5 A/ v1 r                      X( 10, 5)       0.3115900E-01        0.000000
    3 d. M9 u) G! t8 R( _% u. {8 c3 N                      X( 10, 6)        1.000000            0.000000
    8 l+ G- g- v# v                      X( 11, 1)        0.000000            0.000000
    & g8 i0 b# r* y0 J' V) `                      X( 11, 2)       0.8735000E-02        0.000000
    ( w' Q! V, Q! G$ k  f7 L                      X( 11, 3)       0.1525900E-01        0.000000
    ) G) i  M+ L9 _                      X( 11, 4)       0.2889120            0.000000
    5 Z' ?/ ]$ u0 e4 E' Y: R                      X( 11, 5)       0.3700900E-01        0.000000
    7 k% B/ q3 i4 T; |, h- \( T                      X( 11, 6)       0.4091800E-01        0.000000
    & S( D3 c! @) k  P" n$ n4 E0 l                      X( 12, 1)       0.5718100E-01        0.000000
    ! L) [* s# C( S' p" ]! a8 g& ]; ]                      X( 12, 2)       0.3850800E-01        0.000000
    : X0 k1 Z; F. u, [                      X( 12, 3)       0.2541090            0.0000001 I# d: q  {' `; H; D" c1 h
                          X( 12, 4)       0.1842790            0.000000# [7 r- r" L& q- t
                          X( 12, 5)       0.9316000E-01        0.000000! D5 d* R- J4 G
                          X( 12, 6)       0.5224600E-01        0.000000) p5 V6 a  e$ r, X6 X, R" s1 S+ C
                          X( 13, 1)       0.4821100E-01        0.000000$ b4 d0 m! L0 y# k: k  J( K9 F
                          X( 13, 2)       0.9641000E-02        0.000000
    / @, ]' O7 p- n$ ?4 k6 {                      X( 13, 3)        0.000000            0.000000
    9 k9 Z  e, D1 Y                      X( 13, 4)       0.3878190            0.000000
    3 b  W; i5 Q8 d7 ?9 L                      X( 13, 5)       0.1358340            0.000000
    , W5 b% w7 {) L, g! f  v                      X( 13, 6)       0.1185130            0.000000
    . J$ f) k& }4 o8 _% q; s( I1 @( K8 m( T5 i" ]
                                Row    Slack or Surplus      Dual Price
    * Q+ G  u* \9 D( l& i# D                              1       0.7272518            1.000000
    ( _# M. M+ E; }, j+ O) n                              2        0.000000           0.3636259
    # R% C- Q% D+ G                              3        1.124802            0.0000006 u& c8 O$ ]' X# R2 r" O; s9 ]
                                  4        1.121677            0.000000: Q6 C$ l5 l7 f# E9 N
                                  5        1.162172            0.000000
    1 G: m& Q+ ^2 |$ f                              6        1.631867            0.000000" p9 A; u' z- H% z+ D$ J. F
                                  7        1.386806            0.000000
    ; y2 b8 U! ?: f                              8        1.628048            0.0000003 X+ f: c$ O7 G7 w/ l
                                  9       0.8751980            0.000000
    0 i+ }5 b7 p. |8 k0 ]                             10       0.8783228            0.000000  J. ]. c  m7 O
                                 11       0.8378279            0.000000
    0 h# Q1 F. P6 s9 c1 v                             12       0.3681325            0.000000( S! R" k- e9 n) [5 S
                                 13       0.6131944            0.000000
    . E7 q( P% ~( {# h                             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 J/ b$ ~, ?/ a" a
    出现错误提示
    5 h" {/ P* F- G还有 上面的 全局的做大值还没有局部做大的大  为什么！！！
    回复

    使用道具 举报

    返回列表 发新帖
    高级模式
    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-8-23 06:48 , Processed in 0.422784 second(s), 102 queries .

    回顶部