快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

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

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

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

7

主题

5

听众

28

积分

升级  24.21%

该用户从未签到
跳转到指定楼层
1#
发表于 2009-9-5 07:49 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
model:
) I' p( v& K5 L( m' j! {5 Isets:9 r! P7 p8 \# m5 a1 r* X
shengshi/1..13/;8 Q' \1 d5 t, }' N3 F6 f# x
zhibiao/1..6/:a;) v' G- ]5 k; s; p4 c& k
links(shengshi,zhibiao):x;- b1 L* k5 a6 ?, n1 `
endsets
  [8 h. O1 S9 H' H  r- W; Udata:4 o5 C: `7 F) d' {
x=) _1 x* E0 O2 _
0.099599  0.104095  0.021186  0.489849  0.031419  0.399601
3 G$ I% X1 Z, r5 e, c4 z0.012222  0.006320  0.035075  0.558043  0.253776  0.465170 , U2 {/ {+ M  P; M& q! C
0.016260  0.004213  0.084423  0.631442  0.375174  0.535579 ! j9 ~* T7 R: C1 T
0.004587  0.000000  0.086154  0.557522  0.210591  0.464721 * e; u* o( e; b4 j" o) ~! B
0.039880  0.015333  0.077921  1.000000  0.488760  0.888523 " Q7 U$ N! d) A1 X
0.003290  0.002399  0.023441  0.299844  1.000000  0.217814 + d! A' [% x* Z6 C5 g2 y5 w
0.021750  0.023296  0.055727  0.635086  0.657561  0.538723
2 ?5 A+ u! G+ `4 q# m- t1.000000  1.000000  1.000000  0.000000  0.000000  0.430589
7 T& P! G( A: U1 v5 e* k0.070963  0.277539  0.385947  0.153566  0.014757  0.000000
- [/ j/ P! _6 D1 u0.038164  0.038373  0.024990  0.489849  0.031159  1.000000
$ l5 Z$ @* M  f# @0.000000  0.008735  0.015259  0.288912  0.037009  0.040918 ' l+ |4 s4 U* r5 l# @9 K) k8 @* B/ B
0.057181  0.038508  0.254109  0.184279  0.093160  0.052246 ; J0 |8 e1 S7 ]2 H) u1 Y
0.048211  0.009641  0.000000  0.387819  0.135834  0.118513;
+ X) Y' |5 t* _: oenddata
8 B' n$ C/ ]" {: L  y  Smax=- `( u  D7 [8 M9 B& Q$ ^/ F/ i
(@sum(shengshi(i):9 W9 o9 I& L0 f/ |2 R
(@sum(links(i,j):a(j)*x(i,j))-@sum(shengshi(i):
: t) `; B- L1 E  L- q) Z. a@sum(links(i,j):a(j)*x(i,j)))/13)^2)/12)^(1/2);
, \: Q9 N& _, D" M# G( b4 A6 |7 c  F
@sum(zhibiao(j):a(j)^2)=1;
  h/ Z6 \7 @! `$ ^6 p@for(zhibiao(j):a(j)>-1);
) R/ s$ S! Q- r8 D@for(zhibiao(j):a(j)<1);% Y: d% v5 [. m1 ?6 f6 T  ^7 M: y
@for(zhibiao(j)free(a(j)));
! V# c* U0 ]/ {- O0 wend
" D0 p& w% k( n4 d( [4 i! m  O1 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.
    ) e' g  ^6 q4 B) c) S% B  Objective value:                             0.72725181 e& b2 t& h5 W, x
      Objective bound:                             0.7272527
    7 J7 ]- ?$ l& {4 Y$ }$ L, @3 g$ Z  Infeasibilities:                              0.0000005 d4 r" G) {- E. U" @5 q- \5 A" Y
      Extended solver steps:                           38843
    . _0 g! |& Y& w( D  Total solver iterations:                       75613472 U; b( x' S  X* R

    % E' P( K$ v( l& b0 }# A  A* y$ ]( q# V
                           Variable           Value        Reduced Cost" r: i+ m( G1 y0 x4 B6 f; E
                              A( 1)       0.1248020            0.000000" I+ L- Q' J  C! q
                              A( 2)       0.1216772            0.000000& P/ W! O  ^: [0 Q" O, ~
                              A( 3)       0.1621721            0.000000
    8 q1 E- D+ H5 o9 N6 M                          A( 4)       0.6318675            0.0000009 q/ F8 O. c! o  B
                              A( 5)       0.3868056            0.000000
    ) ^- x. _) B- _$ J! l5 D* y) M5 J                          A( 6)       0.6280479            0.000000$ e6 T- B* H- j" n# k4 g3 S; A
                           X( 1, 1)       0.9959900E-01        0.000000
    ; x  b, q6 D! W  L; N                       X( 1, 2)       0.1040950            0.000000
    7 O  t) J* u6 Y( x/ {: o# t                       X( 1, 3)       0.2118600E-01        0.000000
    " w" p" }7 m% n( }. C                       X( 1, 4)       0.4898490            0.000000
    8 ^8 I1 H+ f8 b' \+ Z  a$ q7 _                       X( 1, 5)       0.3141900E-01        0.000000
    7 {7 r& t/ L, t+ d9 a& v                       X( 1, 6)       0.3996010            0.000000
    - i- w9 q1 }% }                       X( 2, 1)       0.1222200E-01        0.000000
    7 Y+ R' T' H( @5 u                       X( 2, 2)       0.6320000E-02        0.000000; F+ k. W4 S) X5 F. Z9 s
                           X( 2, 3)       0.3507500E-01        0.000000
    ) u6 p) E( }/ a9 L$ E- E. N                       X( 2, 4)       0.5580430            0.000000
    6 d, T4 }8 p9 j5 @* h5 p( R5 v  _                       X( 2, 5)       0.2537760            0.000000) Y5 N3 u" m3 M
                           X( 2, 6)       0.4651700            0.000000
    2 H. x. B6 J9 J( K: s4 C2 D( \* P( @                       X( 3, 1)       0.1626000E-01        0.000000( |# T( X% |) S5 y7 _
                           X( 3, 2)       0.4213000E-02        0.000000
    0 ^9 k- ^- y; L2 u4 o6 _6 F; Z  k                       X( 3, 3)       0.8442300E-01        0.000000: a2 H; @% E* u& C5 B$ d
                           X( 3, 4)       0.6314420            0.000000
    # [( o# ~$ w& o' _$ j                       X( 3, 5)       0.3751740            0.000000
    0 E* ]5 w7 L+ P; k" Y3 [% u                       X( 3, 6)       0.5355790            0.000000
    ' Z: Z5 p: m/ `+ G                       X( 4, 1)       0.4587000E-02        0.000000# b* o- \, e: B2 W) B3 ?8 [4 D& O
                           X( 4, 2)        0.000000            0.000000
    8 s7 x4 }* L" s/ n                       X( 4, 3)       0.8615400E-01        0.000000: q- U) Z8 I) }8 g
                           X( 4, 4)       0.5575220            0.000000% x" i% O% ?$ M; s2 [) \8 v
                           X( 4, 5)       0.2105910            0.000000& h' n4 L# @2 E; Q0 P
                           X( 4, 6)       0.4647210            0.000000; J: @8 _% V( h7 I4 S! V
                           X( 5, 1)       0.3988000E-01        0.000000. n9 o0 k$ f/ {0 u
                           X( 5, 2)       0.1533300E-01        0.000000, \4 l$ |$ w7 M2 j9 W* _
                           X( 5, 3)       0.7792100E-01        0.000000! I/ W$ d+ ~7 D+ j0 F& M& {1 P* L: V
                           X( 5, 4)        1.000000            0.000000
    ( W  t0 B! F2 }, ^$ u                       X( 5, 5)       0.4887600            0.000000
    ) _1 w& s- Q% N8 y$ I! i% s8 x                       X( 5, 6)       0.8885230            0.0000001 A9 E# z5 H& @, y5 u
                           X( 6, 1)       0.3290000E-02        0.000000
    , ]% R! ]* g! z. Z                       X( 6, 2)       0.2399000E-02        0.000000
    7 ?$ ^3 O" \8 A- k                       X( 6, 3)       0.2344100E-01        0.000000; }& N6 r% |. B% N5 T; r" v
                           X( 6, 4)       0.2998440            0.000000! W" h: B+ [; U, n) `
                           X( 6, 5)        1.000000            0.0000007 Y' I  `- T. _
                           X( 6, 6)       0.2178140            0.000000- z: U5 e  B0 L  S& a' \0 |  a8 }
                           X( 7, 1)       0.2175000E-01        0.000000
    / e$ f6 q  N- p8 @  i) v                       X( 7, 2)       0.2329600E-01        0.000000
    6 l  g: _4 R! j+ d                       X( 7, 3)       0.5572700E-01        0.000000* z/ S; I5 n: v1 o/ |' a! ]
                           X( 7, 4)       0.6350860            0.000000
    ( Q% [* j3 {, [) }: n; M! \& l                       X( 7, 5)       0.6575610            0.0000008 I# f8 k: ?. k& ~. K
                           X( 7, 6)       0.5387230            0.0000009 D4 g( h! [  r- H
                           X( 8, 1)        1.000000            0.000000
    % a7 v( J7 ~" H. C                       X( 8, 2)        1.000000            0.000000
    ; g# S# P. k, `; t6 Z/ W2 O                       X( 8, 3)        1.000000            0.000000+ [9 [/ m3 Y% ~+ y  z6 ~. ^2 Y
                           X( 8, 4)        0.000000            0.000000% R; e6 V! [- o: q7 l
                           X( 8, 5)        0.000000            0.000000
    & W4 @7 N  Y1 d0 [; C" [( ?                       X( 8, 6)       0.4305890            0.000000) }# ?: |) C% f& j: }
                           X( 9, 1)       0.7096300E-01        0.000000
    0 V- m0 z% C0 E, ^" t                       X( 9, 2)       0.2775390            0.000000: h9 x- j. X$ g/ z! `) l' B  Z3 n
                           X( 9, 3)       0.3859470            0.000000; ]6 J+ f3 }% F' n. H; _; [; j$ _) D) Q
                           X( 9, 4)       0.1535660            0.000000+ C0 G; m0 n. R. s$ u
                           X( 9, 5)       0.1475700E-01        0.000000+ v& o; y' o) e4 L, v* ?' N
                           X( 9, 6)        0.000000            0.000000
    / W# W# q: v; x2 ~- H                      X( 10, 1)       0.3816400E-01        0.000000! v9 H+ B! i) L& q
                          X( 10, 2)       0.3837300E-01        0.0000003 C: L( q; m( r& z3 e! ~
                          X( 10, 3)       0.2499000E-01        0.000000
    0 Z& ?( ]$ n6 e( K                      X( 10, 4)       0.4898490            0.000000# h, |) x4 A& Y& w6 ~
                          X( 10, 5)       0.3115900E-01        0.000000
    3 K5 M8 I: k' p5 H- O% G+ M0 \                      X( 10, 6)        1.000000            0.0000003 ^* n) j. o% V+ `
                          X( 11, 1)        0.000000            0.000000
    3 z( b& z8 G. C$ T+ G& n& G' u                      X( 11, 2)       0.8735000E-02        0.000000( i# o% p* C* p% b0 X0 s
                          X( 11, 3)       0.1525900E-01        0.000000& P3 w% X4 {) O+ y
                          X( 11, 4)       0.2889120            0.000000/ Z+ d. ^+ q% P2 K
                          X( 11, 5)       0.3700900E-01        0.000000
    % \) G7 x  q$ n# l0 |6 b                      X( 11, 6)       0.4091800E-01        0.000000
    5 ]% u5 y# P- g' D1 E, ?                      X( 12, 1)       0.5718100E-01        0.000000
    : h) l% Z1 H# g0 |( ?: q8 \                      X( 12, 2)       0.3850800E-01        0.000000) M# l1 C! b$ ?( P
                          X( 12, 3)       0.2541090            0.000000
    1 E8 N3 n" t! J6 H                      X( 12, 4)       0.1842790            0.000000
    ' H+ D* u6 x! {7 [2 S9 _. F4 @0 s                      X( 12, 5)       0.9316000E-01        0.000000
    5 |" J2 p) ~+ \6 W                      X( 12, 6)       0.5224600E-01        0.000000
    4 V& A. g' I$ O$ \* y6 S5 U                      X( 13, 1)       0.4821100E-01        0.000000
    6 y7 p0 X7 C7 H# p0 k$ e                      X( 13, 2)       0.9641000E-02        0.000000
    0 ~+ p& k7 G* F+ s% P+ ?                      X( 13, 3)        0.000000            0.000000
    * N: [' v7 g& e! y# X0 e) O. @                      X( 13, 4)       0.3878190            0.000000
    , `. B; ^0 w4 V! W4 \                      X( 13, 5)       0.1358340            0.000000& n* t' b4 L. ^- J/ r
                          X( 13, 6)       0.1185130            0.000000
    8 \% C" c$ M7 R8 d7 k0 }0 }/ M9 e* X, Z
                                Row    Slack or Surplus      Dual Price" M5 T" R/ \2 W3 h5 v" V& J0 t4 V
                                  1       0.7272518            1.000000
    ! n  _. q. }* l, ~0 R# L                              2        0.000000           0.3636259
    . r; A# A4 _0 k( G6 H9 Y2 g                              3        1.124802            0.000000  T% B& _/ Q3 o. ?5 F6 n! c7 B8 @) h- W$ ~- o
                                  4        1.121677            0.0000006 |" W) E, M3 S' D; d
                                  5        1.162172            0.000000
    2 {% F; e+ P6 Z* G                              6        1.631867            0.0000003 k: @/ K0 O' i' q3 E# _
                                  7        1.386806            0.000000
    6 a* E& _' |9 ~3 R7 g. X0 C! T                              8        1.628048            0.000000
    ( o7 e5 y# p5 M* `9 A2 I                              9       0.8751980            0.000000
    " `6 K; @: h9 x                             10       0.8783228            0.000000
    + J6 W/ S: m4 w& F+ k                             11       0.8378279            0.000000
    + K; c& X" K3 u8 g9 c, C: ]9 Z                             12       0.3681325            0.000000
    0 U  [2 r& e* F! @1 s* i6 F4 m                             13       0.6131944            0.0000007 M3 H. v8 Z- q: Q
                                 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
    但是 我设置全局求解时：两三秒就停止了
    & k$ @5 u7 v; n# R, z出现错误提示7 S7 n6 v6 Z3 ]; O+ x; I$ n" J1 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, 2025-12-4 15:43 , Processed in 2.971535 second(s), 103 queries .

    回顶部