快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

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

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

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

7

主题

5

听众

28

积分

升级  24.21%

该用户从未签到
跳转到指定楼层
1#
发表于 2009-9-5 07:49 |只看该作者 |正序浏览
|招呼Ta 关注Ta
model:; T# s* T5 E2 I2 L1 _
sets:3 N. x  M" I: X
shengshi/1..13/;
$ l9 C) V7 U" S8 U% J, s2 yzhibiao/1..6/:a;
0 X5 i  ~' x, u- [% Flinks(shengshi,zhibiao):x;
4 F# O- o/ \2 h! Vendsets
& c; a" H, g7 w& Q3 C8 gdata:3 F" o- y$ s% `) Z! J3 D
x=- p+ |- W$ t5 W6 @
0.099599  0.104095  0.021186  0.489849  0.031419  0.399601 , |7 G4 D4 d  g2 E( j
0.012222  0.006320  0.035075  0.558043  0.253776  0.465170
  x0 o; a* z  M. q2 J0 _- Z# B. W! C! s0.016260  0.004213  0.084423  0.631442  0.375174  0.535579 % O1 n, Z* o  m- Y0 F
0.004587  0.000000  0.086154  0.557522  0.210591  0.464721
  Z0 ], d! ?+ D0.039880  0.015333  0.077921  1.000000  0.488760  0.888523 $ `) f8 D3 P! L% r4 W8 B& \- }  b
0.003290  0.002399  0.023441  0.299844  1.000000  0.217814
$ x/ E# `9 Q) R/ k8 S" t0.021750  0.023296  0.055727  0.635086  0.657561  0.538723 2 H2 @( ^: P( H+ L& [
1.000000  1.000000  1.000000  0.000000  0.000000  0.430589
2 @. |9 N, @; ~3 S0.070963  0.277539  0.385947  0.153566  0.014757  0.000000 2 M% |( S1 e/ @% N$ w% D
0.038164  0.038373  0.024990  0.489849  0.031159  1.000000
, _7 S; T, @0 h& d+ S0.000000  0.008735  0.015259  0.288912  0.037009  0.040918
) _; k8 T' |0 q% {3 E7 ^0.057181  0.038508  0.254109  0.184279  0.093160  0.052246
" c: o- ?/ L5 e0.048211  0.009641  0.000000  0.387819  0.135834  0.118513;. O- \& e7 O8 n
enddata7 [2 m% Z% t0 U
max=
- v# U! c! B) q(@sum(shengshi(i):& i- U) m6 x3 J* r3 q) }* H
(@sum(links(i,j):a(j)*x(i,j))-@sum(shengshi(i):5 P, D% p7 V* y0 x8 M( _% F. X
@sum(links(i,j):a(j)*x(i,j)))/13)^2)/12)^(1/2);0 L! w9 J2 C9 ^

. K6 C+ F/ R) c8 n, x& F: z$ N) u" T@sum(zhibiao(j):a(j)^2)=1;
& W7 o, m3 s2 \  z6 ~3 P8 }7 Q) v@for(zhibiao(j):a(j)>-1);9 I  f; \/ ~% F* v0 \4 O( W: H& C
@for(zhibiao(j):a(j)<1);
0 }4 z$ m0 ~8 M4 j( N/ S8 i* n# z0 e! w@for(zhibiao(j)free(a(j)));
& E* ~2 i0 n& f) r# e) p& Hend , U, z, G9 `" ~7 Y  M
求不出全局最优解，请大侠帮帮忙！！！
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
回复

使用道具 举报

jiiing        

7

主题

5

听众

28

积分

升级  24.21%

该用户从未签到
10#
发表于 2009-9-5 12:16 |只看该作者 |招呼Ta 关注Ta
但是 我设置全局求解时：两三秒就停止了
, |# Q* i& a/ U# M% u" U出现错误提示0 `2 @) p: L) {  s( K/ Z
还有 上面的 全局的做大值还没有局部做大的大  为什么！！！
回复

使用道具 举报

liulongfei 实名认证       

2

主题

2

听众

100

积分

升级  0%

该用户从未签到

群组: 数学趣味、游戏、IQ等

群组: 数学建摸协会

群组: 华中科技大学
9#
发表于 2009-9-5 09:35 |只看该作者 |招呼Ta 关注Ta
选择全局最优解，然后就是时间的问题了
回复

使用道具 举报

xuyin        

0

主题

4

听众

198

积分

升级  49%

  • TA的每日心情
    开心
    2016-3-23 08:32

    • 签到天数: 26 天

    [LV.4]偶尔看看III
    8#
    发表于 2009-9-5 09:06 |只看该作者 |招呼Ta 关注Ta
    哇，都是高手啊
    回复

    使用道具 举报

    yuan07150215        

    0

    主题

    4

    听众

    168

    积分

    升级  34%

    该用户从未签到
    7#
    发表于 2009-9-5 08:58 |只看该作者 |招呼Ta 关注Ta
    先在限制条件下得出最低值，之后运行lingo程序，当运行到20-30分钟时，人为的终止求解（按下Interrupt Solver），此时运行的解与目标函数的最优解很相近！可认为是最优解！……
    回复

    使用道具 举报

    formulaone9275        

    1

    主题

    4

    听众

    188

    积分

    升级  44%

    该用户从未签到
    6#
    发表于 2009-9-5 08:54 |只看该作者 |招呼Ta 关注Ta
    有很多时候都是求不出最优解的，局部最优就可以用了 ，我看了你的程序，其中max里的sum有对i求了两次和，你看看能不能再改改，求求试试
    回复

    使用道具 举报

    ychexingwuji        

    4

    主题

    4

    听众

    70

    积分

    升级  68.42%

    该用户从未签到
    5#
    发表于 2009-9-5 08:49 |只看该作者 |招呼Ta 关注Ta
    学习了 呵呵
    回复

    使用道具 举报

    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分钟才计算出来，给他一个初始值就好了！
    回复

    使用道具 举报

    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.
    & y) G- O6 t8 T/ n, i8 m. \& P  Objective value:                             0.7272518
    9 |  g) s. Y2 y+ N6 N7 n3 o* Y; E& P3 E: O, `  Objective bound:                             0.7272527: D+ a0 u7 b0 ?
      Infeasibilities:                              0.0000003 `2 c& o9 b) h' u9 x
      Extended solver steps:                           38843
    " P, R: E8 ?! S. b- h) {* t: S' O) T  Total solver iterations:                       7561347" T) b9 @  g* }4 @, d

    - [) p1 L1 s$ _% L' t( ?: K  y1 R  ~$ o7 U) h6 W  {
                           Variable           Value        Reduced Cost# N3 ]0 i1 y' V, @( f* X# @
                              A( 1)       0.1248020            0.000000
    ; F; [" R" I8 V% V- b8 g                          A( 2)       0.1216772            0.0000000 y8 l: ]( s4 X. K+ B  t
                              A( 3)       0.1621721            0.000000
    $ ~2 J+ C1 n$ I+ e' T- M. T                          A( 4)       0.6318675            0.000000
    $ O8 y$ _# c" ]( {                          A( 5)       0.3868056            0.000000) R; |6 V1 V. L& \  R: ^  W7 t
                              A( 6)       0.6280479            0.000000
    ( u! ~. F1 g$ O, N" X                       X( 1, 1)       0.9959900E-01        0.000000" q& a) V7 r5 B; I' E  t) n! S9 m2 I
                           X( 1, 2)       0.1040950            0.000000
    5 J: v* P; R* v) ?. |                       X( 1, 3)       0.2118600E-01        0.000000
    ) V. g9 Q7 n# P* l                       X( 1, 4)       0.4898490            0.000000' ~1 b* g' e2 w7 ?
                           X( 1, 5)       0.3141900E-01        0.000000" k! a7 h% w/ a& D  i
                           X( 1, 6)       0.3996010            0.000000
    ! W2 F4 u! O1 D- @! n) _                       X( 2, 1)       0.1222200E-01        0.000000
    4 C8 G3 F- I$ d' O8 \                       X( 2, 2)       0.6320000E-02        0.000000  Y3 U& l. l8 j  r
                           X( 2, 3)       0.3507500E-01        0.000000" O  X3 i* Z% {5 _7 U, B& F
                           X( 2, 4)       0.5580430            0.000000" z$ E" }/ Y3 d6 O
                           X( 2, 5)       0.2537760            0.000000
    & q- n" C7 m1 I% l/ t                       X( 2, 6)       0.4651700            0.000000
    3 {  @& U  M9 K4 @                       X( 3, 1)       0.1626000E-01        0.000000
    0 @- I& t1 t9 a1 M3 t, t                       X( 3, 2)       0.4213000E-02        0.000000$ [4 D1 ^9 ]( `2 C
                           X( 3, 3)       0.8442300E-01        0.0000004 g# d2 @" T8 [
                           X( 3, 4)       0.6314420            0.000000$ T/ g) a. S* H$ R# m8 u9 t8 f3 I
                           X( 3, 5)       0.3751740            0.000000
    0 i% X1 E8 V0 t                       X( 3, 6)       0.5355790            0.000000
    " G% @1 N* ^; _) p1 j) ]2 P                       X( 4, 1)       0.4587000E-02        0.0000007 z) k" C: m+ B% A+ ]% |) l
                           X( 4, 2)        0.000000            0.000000
    & e+ n% j1 P4 _7 H. b- B                       X( 4, 3)       0.8615400E-01        0.000000
    3 k& |  d: v+ u                       X( 4, 4)       0.5575220            0.000000
    ; v/ [, i4 @, n) m1 {9 p3 o5 L                       X( 4, 5)       0.2105910            0.000000+ ?7 U+ j, a5 b" l
                           X( 4, 6)       0.4647210            0.000000
    2 t% M3 z0 A+ i                       X( 5, 1)       0.3988000E-01        0.000000# T% i2 x3 N3 Y4 h7 Z+ N( x
                           X( 5, 2)       0.1533300E-01        0.0000009 k/ p8 j7 j- P* {- a6 K3 X
                           X( 5, 3)       0.7792100E-01        0.0000004 R4 }$ q4 F, J* p  W$ ]$ Z
                           X( 5, 4)        1.000000            0.000000& x. {0 K& l, w* \8 [% W* b
                           X( 5, 5)       0.4887600            0.000000
    ' c) q. x3 D3 Y+ r9 J                       X( 5, 6)       0.8885230            0.000000
    , _# z0 d/ h8 X                       X( 6, 1)       0.3290000E-02        0.0000007 T( N+ x' S$ r
                           X( 6, 2)       0.2399000E-02        0.000000( m! v, P( D( s& Z. A$ y# h: A
                           X( 6, 3)       0.2344100E-01        0.000000
    - i9 l. y- }: C4 n                       X( 6, 4)       0.2998440            0.000000
    / ]5 M- ?3 {) L/ s; ~                       X( 6, 5)        1.000000            0.000000
    6 x- d6 w0 `& }. P+ W1 V                       X( 6, 6)       0.2178140            0.000000. ]8 h6 j0 w$ L( ~* s/ x" ^! L/ T
                           X( 7, 1)       0.2175000E-01        0.0000000 x( K$ l& I2 Q
                           X( 7, 2)       0.2329600E-01        0.000000
    , V9 [' t" a, `$ ?4 l                       X( 7, 3)       0.5572700E-01        0.000000
    " U9 x( @: |2 r: W3 G: o$ p! F                       X( 7, 4)       0.6350860            0.000000
    7 a) P  Y1 j+ i4 f- H% O3 \, K                       X( 7, 5)       0.6575610            0.000000
    3 l  Q( T) Y2 B6 g$ Z1 s+ |                       X( 7, 6)       0.5387230            0.000000* Z' T# Q' @4 e- r- R: P) r/ o' A
                           X( 8, 1)        1.000000            0.000000
    ) L2 C) D! n. w( s9 k+ M& K                       X( 8, 2)        1.000000            0.0000002 A4 D% D% x. S9 v
                           X( 8, 3)        1.000000            0.000000
    ) z. Q& f/ y5 x+ j# p! u                       X( 8, 4)        0.000000            0.000000$ e/ ~0 j3 d* N) w8 x! V
                           X( 8, 5)        0.000000            0.000000% J$ v) l: v2 L# p
                           X( 8, 6)       0.4305890            0.000000
    $ S& i% }* ~. x                       X( 9, 1)       0.7096300E-01        0.000000# _' T; d% l( H) b% _$ S- E
                           X( 9, 2)       0.2775390            0.000000
    ) q$ K2 Z7 N3 a/ S6 c+ k                       X( 9, 3)       0.3859470            0.000000
    , k. b- \' L" S* z7 F                       X( 9, 4)       0.1535660            0.000000
    ; k- h! R- Y4 k6 R5 ~                       X( 9, 5)       0.1475700E-01        0.0000006 y6 d) ]& k& x: O. H
                           X( 9, 6)        0.000000            0.000000; ?1 O/ ^; o5 \: @* p4 _
                          X( 10, 1)       0.3816400E-01        0.000000
    7 }2 N& a( o  o) u                      X( 10, 2)       0.3837300E-01        0.000000
    2 c3 q; K1 a- T9 Y! \9 R                      X( 10, 3)       0.2499000E-01        0.000000
    ! v) M% A& M' g0 b4 j. I" J! Y                      X( 10, 4)       0.4898490            0.000000
    0 Q6 g1 @" z/ X9 v                      X( 10, 5)       0.3115900E-01        0.000000
    " k4 K1 w' O6 [- f& L. Q                      X( 10, 6)        1.000000            0.000000
    / G3 M9 ~. G8 g* K                      X( 11, 1)        0.000000            0.000000) ]# h: s( y2 `) Z/ l! {0 t
                          X( 11, 2)       0.8735000E-02        0.000000
    6 h& B6 r3 L: x: ^5 M1 v, `                      X( 11, 3)       0.1525900E-01        0.000000! q# ?, y4 X! A' F
                          X( 11, 4)       0.2889120            0.000000
    ( Q  u3 D- H: i" \                      X( 11, 5)       0.3700900E-01        0.000000- {, ]! B0 U5 c$ F9 t8 {
                          X( 11, 6)       0.4091800E-01        0.000000
    . C  M( G. z" A; ?1 O                      X( 12, 1)       0.5718100E-01        0.000000
      e6 `/ t4 U- S* u% N                      X( 12, 2)       0.3850800E-01        0.0000006 O  X+ s7 I. `, f- o* r. W! F0 R
                          X( 12, 3)       0.2541090            0.000000# F: e3 K, V2 h8 ?# S
                          X( 12, 4)       0.1842790            0.000000' Q; k+ y2 T7 G* e0 p6 v7 j1 i  t
                          X( 12, 5)       0.9316000E-01        0.000000
      T9 ~- a% n& d+ [2 f6 W                      X( 12, 6)       0.5224600E-01        0.000000. `1 x8 L0 `' P/ Y
                          X( 13, 1)       0.4821100E-01        0.000000
    % [6 g4 e* _2 j+ a. `                      X( 13, 2)       0.9641000E-02        0.0000001 u) |! k2 W: V/ F5 k; K8 p2 E3 y
                          X( 13, 3)        0.000000            0.000000/ j. N: L4 I! @" }) u2 X
                          X( 13, 4)       0.3878190            0.000000
    5 |% z* _: A  L1 U                      X( 13, 5)       0.1358340            0.0000000 P( T. l0 H" {4 p- h3 \
                          X( 13, 6)       0.1185130            0.000000( c; w; l2 D; S" K2 j

    ) N1 e& \! n0 E6 q' r/ h                            Row    Slack or Surplus      Dual Price+ W' F% B  K' q: M
                                  1       0.7272518            1.000000
    , n: P/ s' V. Z- ~4 A" W                              2        0.000000           0.3636259! F  e1 U3 k' F8 m
                                  3        1.124802            0.000000
    " p- }8 u! f! r6 r- J4 W                              4        1.121677            0.000000
    / G2 L- u2 N! N- D, ?- o                              5        1.162172            0.000000$ v, H: a8 a( y
                                  6        1.631867            0.0000007 j) h9 ~9 S) m& Q( N5 L' v
                                  7        1.386806            0.000000
    7 R9 ^1 W* L' X3 b                              8        1.628048            0.0000000 z# z$ A* }7 d9 z; o; z% E
                                  9       0.8751980            0.0000007 r+ \7 \- T* `& t
                                 10       0.8783228            0.000000
    " N/ Y+ c( M7 A                             11       0.8378279            0.0000009 Y' i: X4 ]) c; X
                                 12       0.3681325            0.0000006 f8 E5 |1 d* \
                                 13       0.6131944            0.000000- L7 y  X$ J" y* d% r! y3 j! R
                                 14       0.3719521            0.000000
    回复

    使用道具 举报

    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
    运算时间太长，你能否增加一个初始值就好了
    回复

    使用道具 举报

    返回列表 发新帖
    高级模式
    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-6 00:28 , Processed in 0.401247 second(s), 104 queries .

    回顶部