快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

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

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

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

7

主题

5

听众

28

积分

升级  24.21%

该用户从未签到
跳转到指定楼层
1#
发表于 2009-9-5 07:49 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
model:4 Y6 Q! T+ n5 `& P! ~
sets:/ e  i3 J& M7 b0 `5 r' n
shengshi/1..13/;
+ P7 P8 K  L5 u- k% rzhibiao/1..6/:a;
/ j& A/ g/ ?. U5 S, I! blinks(shengshi,zhibiao):x;
9 b% j4 u9 W. X, i$ ]( H5 P! J! H, Zendsets
$ n! r( `2 t& G# gdata:
5 @6 H7 d. i# l9 x4 ~0 E. P( @3 yx=
& Y3 Y6 p$ N% r$ \2 v2 ~! E7 I( e9 r0.099599  0.104095  0.021186  0.489849  0.031419  0.399601
& b" t+ v8 [# Y$ e! p0.012222  0.006320  0.035075  0.558043  0.253776  0.465170
6 G/ w- \( o2 k) ?3 a) d( R0.016260  0.004213  0.084423  0.631442  0.375174  0.535579
$ Y) L" I$ l, o0 p* Y; o$ c0.004587  0.000000  0.086154  0.557522  0.210591  0.464721 ( y+ D& B6 s5 r* S
0.039880  0.015333  0.077921  1.000000  0.488760  0.888523
$ j$ {! l+ K) F9 m1 Z( L: z/ ]7 x' o0.003290  0.002399  0.023441  0.299844  1.000000  0.217814
3 H& y+ ?0 K2 R4 Z) D7 R0.021750  0.023296  0.055727  0.635086  0.657561  0.538723
: x/ P+ L: v1 T; i9 c: y: ^2 l1.000000  1.000000  1.000000  0.000000  0.000000  0.430589
; I& Y3 O4 s% X% P0.070963  0.277539  0.385947  0.153566  0.014757  0.000000   g! X5 Z/ w1 ^5 ~! Z" |
0.038164  0.038373  0.024990  0.489849  0.031159  1.000000
) J3 E4 S: z5 K+ J0.000000  0.008735  0.015259  0.288912  0.037009  0.040918
9 {$ C5 ^# G) R" P5 K/ G0.057181  0.038508  0.254109  0.184279  0.093160  0.052246 ! c6 Z4 g# V( v  z
0.048211  0.009641  0.000000  0.387819  0.135834  0.118513;+ @9 l4 @1 H) g" \1 a; A
enddata; D' U# o: j/ h0 K( F% C  r( J3 |6 e
max=
. J0 N: C  Q: C/ L% F( t1 [9 y(@sum(shengshi(i):
  R- z* n* A. }: Y6 }( I7 M9 \0 }1 x(@sum(links(i,j):a(j)*x(i,j))-@sum(shengshi(i):
, h5 ^4 W, Z' N* s@sum(links(i,j):a(j)*x(i,j)))/13)^2)/12)^(1/2);
4 d  m: m; A* d
' V0 X- V% x$ y6 I- T$ _@sum(zhibiao(j):a(j)^2)=1;) b) S0 U$ h, J! ^3 p" T$ ?
@for(zhibiao(j):a(j)>-1);
& q. u! q. z7 `) A@for(zhibiao(j):a(j)<1);
9 k- |  R! a/ i5 k$ l6 r@for(zhibiao(j)free(a(j)));
, E! A" @+ |+ s0 |2 T; p- `- A; ^  W; x# Wend # c$ p% e% K7 p1 }) W( |
求不出全局最优解，请大侠帮帮忙！！！
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./ I0 q" I% |# U( {. P5 v
      Objective value:                             0.7272518+ P9 X- I& I+ s/ i/ T
      Objective bound:                             0.7272527, D* e9 P' n9 x- z5 m7 @
      Infeasibilities:                              0.000000: Z3 l- j- U1 l( h
      Extended solver steps:                           388432 }) R* y# b. v: C* _
      Total solver iterations:                       7561347) j' g$ i9 \9 V5 F. w7 T* L, `7 R

    6 I9 e; p7 z2 E& T& }) v
    0 D# e$ A) v4 _9 Z) Z8 {) f                       Variable           Value        Reduced Cost
    / L% i) G+ s8 E8 e2 K                          A( 1)       0.1248020            0.0000004 W/ n& w6 l; L
                              A( 2)       0.1216772            0.000000$ U8 G& B8 S% m! m8 c% t' a
                              A( 3)       0.1621721            0.000000
    * U: n' A& [* {                          A( 4)       0.6318675            0.0000006 ~8 q7 n5 Y' g: v
                              A( 5)       0.3868056            0.000000
    ' C0 u' j% [% I& H  ?% c- g                          A( 6)       0.6280479            0.0000008 M3 X5 f8 x; ^& H
                           X( 1, 1)       0.9959900E-01        0.000000
    . ^7 y1 Q$ k  T, V" L4 ?                       X( 1, 2)       0.1040950            0.000000" T  ]- F) ~6 J' D/ l
                           X( 1, 3)       0.2118600E-01        0.000000$ P* }1 i8 |; U! _" o, s; D
                           X( 1, 4)       0.4898490            0.000000
    * j) b1 @2 ~& Y                       X( 1, 5)       0.3141900E-01        0.000000
    1 ^: F6 |. {% O, L                       X( 1, 6)       0.3996010            0.000000
    % c- \3 E9 i8 j                       X( 2, 1)       0.1222200E-01        0.000000
    ) ^/ d5 N( F& z) O- U  n                       X( 2, 2)       0.6320000E-02        0.000000' h( {9 `0 \& G% n' `. {
                           X( 2, 3)       0.3507500E-01        0.000000* Z3 V- w* L1 k1 p/ h
                           X( 2, 4)       0.5580430            0.000000# T9 ~6 a( u* u2 s/ f0 V6 l
                           X( 2, 5)       0.2537760            0.000000( T5 ~! m/ ?% R. c: d
                           X( 2, 6)       0.4651700            0.000000
    0 M$ D! Z& {. u- _; T0 r                       X( 3, 1)       0.1626000E-01        0.000000: s9 R3 S- C( l4 @* W
                           X( 3, 2)       0.4213000E-02        0.000000
    ) T/ `1 E6 |2 R6 t! k                       X( 3, 3)       0.8442300E-01        0.000000+ B/ @3 ?# t1 X7 S; t5 t% p
                           X( 3, 4)       0.6314420            0.000000
    $ y  c/ F: |8 a$ ^                       X( 3, 5)       0.3751740            0.000000
    - J$ W7 r7 n/ m5 w6 I                       X( 3, 6)       0.5355790            0.000000
    3 T) Y9 x! [( s, V$ U* r1 o                       X( 4, 1)       0.4587000E-02        0.000000
    - q8 P# @# v' Y: U/ S4 f                       X( 4, 2)        0.000000            0.0000000 _3 i. A7 M" ]" B& W
                           X( 4, 3)       0.8615400E-01        0.000000& X& [, v9 w* r
                           X( 4, 4)       0.5575220            0.000000* q1 J1 L/ x9 i
                           X( 4, 5)       0.2105910            0.000000
    . o2 g, \5 t2 I8 n" G' {                       X( 4, 6)       0.4647210            0.000000  X( j7 ~% R3 E( k, \+ j  l! o& \* V
                           X( 5, 1)       0.3988000E-01        0.000000
    ; W& L# ?& A+ }! Q0 n                       X( 5, 2)       0.1533300E-01        0.000000
    0 Q& O3 n2 D1 Z; x+ V                       X( 5, 3)       0.7792100E-01        0.000000
    5 {- ?1 Q: d% {# @1 o                       X( 5, 4)        1.000000            0.000000
    * M  o8 \1 f0 E$ }' m/ @                       X( 5, 5)       0.4887600            0.0000002 B3 G) d& J5 A+ ^8 `: \) |2 i% u
                           X( 5, 6)       0.8885230            0.000000/ X! p- x0 v+ ?; C+ u
                           X( 6, 1)       0.3290000E-02        0.0000006 b# q* `/ c/ M# l- `/ X+ E
                           X( 6, 2)       0.2399000E-02        0.000000$ U5 U/ |* S9 E5 E" |9 V
                           X( 6, 3)       0.2344100E-01        0.0000003 R7 R6 z1 Y' H! q9 N
                           X( 6, 4)       0.2998440            0.000000; s3 P! ~( |3 n, {8 T
                           X( 6, 5)        1.000000            0.0000008 K4 ^# ?3 a4 C4 ]. U# K2 p, @
                           X( 6, 6)       0.2178140            0.000000) k8 T8 L0 O4 v: R/ H7 S$ H
                           X( 7, 1)       0.2175000E-01        0.000000
    5 Z" ~" A/ V5 e3 p4 O                       X( 7, 2)       0.2329600E-01        0.000000% t: }6 J( n; \1 s" b
                           X( 7, 3)       0.5572700E-01        0.000000
    7 C! a+ W" K- B                       X( 7, 4)       0.6350860            0.000000
    2 j* N$ }, @3 D- ]; k                       X( 7, 5)       0.6575610            0.000000) P  P3 }: U2 g8 _- @
                           X( 7, 6)       0.5387230            0.000000- l$ R; ?: A# D4 c+ J8 w
                           X( 8, 1)        1.000000            0.000000. Y3 B& X  l' w
                           X( 8, 2)        1.000000            0.000000
    ) @, h+ Z( V; X8 P- }                       X( 8, 3)        1.000000            0.000000
    0 w9 n. p3 o: A/ W' d( _) V                       X( 8, 4)        0.000000            0.000000
    * H0 K/ j# j8 h" R5 W" N                       X( 8, 5)        0.000000            0.000000
    9 U7 A* ?9 h. _1 v( q& y                       X( 8, 6)       0.4305890            0.000000
    " T$ b+ ^4 |6 d                       X( 9, 1)       0.7096300E-01        0.0000006 k7 p, s" Y1 e- a) o6 F
                           X( 9, 2)       0.2775390            0.0000007 W/ ^+ \% V- ^# N2 Z, {0 N
                           X( 9, 3)       0.3859470            0.000000
    ; W& c' w* n! ~                       X( 9, 4)       0.1535660            0.000000
    0 \* M0 _$ y3 o8 v$ x9 r                       X( 9, 5)       0.1475700E-01        0.000000; k2 g0 v. o& S( f4 [3 X" Z( n8 ]
                           X( 9, 6)        0.000000            0.0000006 |" a7 s; W0 L, I8 F
                          X( 10, 1)       0.3816400E-01        0.000000  Y3 z8 \/ Z/ r8 u
                          X( 10, 2)       0.3837300E-01        0.000000
    6 r* H4 r  k) e  N) g" [4 o) S                      X( 10, 3)       0.2499000E-01        0.000000
    8 v* G. R! O4 S+ |" o8 t$ `+ z                      X( 10, 4)       0.4898490            0.000000/ a* P' Q3 O) j; K& p+ S
                          X( 10, 5)       0.3115900E-01        0.000000) ~9 K# g, }; l; G
                          X( 10, 6)        1.000000            0.000000
    " i" K, s: ^. {" |. B" G& D! x                      X( 11, 1)        0.000000            0.000000* e' `) [6 I; @$ C+ n
                          X( 11, 2)       0.8735000E-02        0.000000. [8 f0 `0 u/ u& N( X* k
                          X( 11, 3)       0.1525900E-01        0.0000005 i2 F) D5 G' W) k
                          X( 11, 4)       0.2889120            0.000000
    7 Z( f( J: \5 T5 u, z! Q* v) U                      X( 11, 5)       0.3700900E-01        0.000000
    $ N) v9 @  J  i: ~                      X( 11, 6)       0.4091800E-01        0.000000
    0 K. C" l2 |4 I7 {                      X( 12, 1)       0.5718100E-01        0.000000
    3 N% `% _8 V7 A7 g                      X( 12, 2)       0.3850800E-01        0.000000# }6 N1 }$ l8 N5 B
                          X( 12, 3)       0.2541090            0.0000003 a4 ^9 Z3 u( Q9 u, }9 d
                          X( 12, 4)       0.1842790            0.000000
    7 h. \2 h3 o- R                      X( 12, 5)       0.9316000E-01        0.000000- s; U. k, l9 @2 j" \7 c
                          X( 12, 6)       0.5224600E-01        0.000000
    % q/ h0 ?! H- k! V. S7 I                      X( 13, 1)       0.4821100E-01        0.000000" u2 t0 @, g9 Z0 X4 h
                          X( 13, 2)       0.9641000E-02        0.000000& [6 k' U# ]* ~
                          X( 13, 3)        0.000000            0.0000000 b9 T# H0 Q! w% m7 y3 F
                          X( 13, 4)       0.3878190            0.000000  [, k0 d. F9 B  g- H3 a
                          X( 13, 5)       0.1358340            0.0000003 @; P3 u. n# K4 K2 V& o9 |
                          X( 13, 6)       0.1185130            0.000000  a' B% }# w0 P8 w. k. v5 H4 a5 U
    * t; y2 n+ E6 W+ C7 v' a/ d3 A
                                Row    Slack or Surplus      Dual Price: g0 i( X8 c- W1 D
                                  1       0.7272518            1.000000: y( J7 S4 d& s" o$ [
                                  2        0.000000           0.3636259
      S4 M# D  q; D! |2 O                              3        1.124802            0.000000
    2 ~, r  R, s6 v                              4        1.121677            0.000000' M  d+ S5 Y, c0 \- q
                                  5        1.162172            0.000000. }, \1 [2 n) `
                                  6        1.631867            0.000000* U# }" p. Z/ q- k# R* j
                                  7        1.386806            0.000000
    6 }* ~- W3 ]; ^  o                              8        1.628048            0.000000
    ; a& ~" L5 c* I# u                              9       0.8751980            0.000000
    0 ?% n( w4 ?+ B9 Q2 w/ K                             10       0.8783228            0.000000# n( j7 f) D! `5 D1 |2 v1 }
                                 11       0.8378279            0.000000
    ( C' w* s7 N- U* L' y1 Q, t                             12       0.3681325            0.000000( B0 h6 J  W7 ]: l+ G
                                 13       0.6131944            0.0000005 u8 e0 Z, a+ R/ S( y
                                 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
    但是 我设置全局求解时：两三秒就停止了# L  {; _1 f; r8 J& J% H# W
    出现错误提示, G2 I: E5 g0 B- b4 B
    还有 上面的 全局的做大值还没有局部做大的大  为什么！！！
    回复

    使用道具 举报

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

    回顶部