快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

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

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

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

7

主题

5

听众

28

积分

升级  24.21%

该用户从未签到
跳转到指定楼层
1#
发表于 2009-9-5 07:49 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
model:
, F& @% l$ M1 l' Usets:" B9 ^! ^6 L  O) _# b
shengshi/1..13/;1 x, g4 [# x: R1 [. b  J# P9 M
zhibiao/1..6/:a;5 S: g/ x0 o1 O1 R
links(shengshi,zhibiao):x;
( S" V% R, V8 @, t5 M5 wendsets
5 r9 Z/ Q! D" T) pdata:+ r* F: j: z! ?$ Z, U8 J
x=2 P4 R# C. F7 H' }, q1 T( @% U
0.099599  0.104095  0.021186  0.489849  0.031419  0.399601 ! G9 G- L; G% _% c" |$ m
0.012222  0.006320  0.035075  0.558043  0.253776  0.465170 - b, [2 P$ k  M4 R
0.016260  0.004213  0.084423  0.631442  0.375174  0.535579 1 M7 g! i% g; j  |- J1 r
0.004587  0.000000  0.086154  0.557522  0.210591  0.464721
* c, `# B  d- Y" x$ l+ F5 l0.039880  0.015333  0.077921  1.000000  0.488760  0.888523
' |/ H  S6 e6 C$ }1 k: N( I0.003290  0.002399  0.023441  0.299844  1.000000  0.217814
) m3 r" _3 j0 `! O+ y0.021750  0.023296  0.055727  0.635086  0.657561  0.538723
  ?9 d* C& K/ C5 F1.000000  1.000000  1.000000  0.000000  0.000000  0.430589 ( h  H( U5 U" J
0.070963  0.277539  0.385947  0.153566  0.014757  0.000000 , o% f( z% {2 E
0.038164  0.038373  0.024990  0.489849  0.031159  1.000000 - Z& o" @# f* m3 b' O
0.000000  0.008735  0.015259  0.288912  0.037009  0.040918
' @: a- j0 L; K$ o3 _0 F) i0.057181  0.038508  0.254109  0.184279  0.093160  0.052246 ' {; X8 ?, F/ u& P4 q# b
0.048211  0.009641  0.000000  0.387819  0.135834  0.118513;
: L( C( ~" E+ B; Venddata
& c5 T3 ~* j5 q' m( o1 Cmax=( C9 w0 i8 S0 g: t; k
(@sum(shengshi(i):
+ [+ ]) U0 H' [- |; N, h, a(@sum(links(i,j):a(j)*x(i,j))-@sum(shengshi(i):# A) a! X. {: B4 @
@sum(links(i,j):a(j)*x(i,j)))/13)^2)/12)^(1/2);7 g' r( }& c) `1 X4 a

! ~# a& f  X0 B+ ^5 i@sum(zhibiao(j):a(j)^2)=1;
0 V5 K0 F" N- x$ E@for(zhibiao(j):a(j)>-1);
/ z/ X$ y, g+ j& \. H$ Q1 Z@for(zhibiao(j):a(j)<1);
% H% K% g" ^  X' F@for(zhibiao(j)free(a(j)));4 b5 [# v$ Y) h& b
end , ^' g! P/ g$ O" j! N1 r
求不出全局最优解，请大侠帮帮忙！！！
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.  b; ]  Y0 o# g/ F; d; Y
      Objective value:                             0.72725186 U. x8 P$ n, b2 J! x- D3 f
      Objective bound:                             0.72725271 B# ^, ]  j5 Q/ h
      Infeasibilities:                              0.000000
    " k0 t2 G$ Z$ J7 U1 \2 B  Extended solver steps:                           38843" @: k+ ]" D! N2 i4 ~/ j; p: q
      Total solver iterations:                       7561347; f; q: s6 x& ~
    $ S! q- K( t6 U

    ; B' c% M' V' `* A; B& \. p                       Variable           Value        Reduced Cost
    ! i: z" R; S) {% {                          A( 1)       0.1248020            0.000000& e9 R, l. I) Z2 [+ l# l
                              A( 2)       0.1216772            0.000000
    ( g: a* i2 ~3 k5 {$ g  v5 `  m& Z, d                          A( 3)       0.1621721            0.000000
    * `/ R3 l& r5 {$ M; d: ^! T# t/ P  z" S                          A( 4)       0.6318675            0.000000
    5 N7 T* ^3 y" g( [) D4 b                          A( 5)       0.3868056            0.000000
    6 v- y& H& A; `1 h9 U' ]5 X3 e                          A( 6)       0.6280479            0.000000
    / D0 S% d; N# p2 h' {                       X( 1, 1)       0.9959900E-01        0.000000
    * s) B/ U$ P  Q$ y+ _                       X( 1, 2)       0.1040950            0.0000000 m1 V( u; T6 ~+ l
                           X( 1, 3)       0.2118600E-01        0.000000
    $ y; n1 s  D1 @* ]                       X( 1, 4)       0.4898490            0.000000
    5 [& |: j$ @9 ^6 d                       X( 1, 5)       0.3141900E-01        0.000000! w) y4 E2 g. T! T/ `6 P
                           X( 1, 6)       0.3996010            0.0000008 L. [& Q6 [: ]0 O' ~, E
                           X( 2, 1)       0.1222200E-01        0.0000003 @4 t7 M; m! K# a( ^
                           X( 2, 2)       0.6320000E-02        0.0000007 P, \! O: R" l+ C3 V( k  B. j
                           X( 2, 3)       0.3507500E-01        0.000000$ ]* r3 E2 i( Z' Z/ {
                           X( 2, 4)       0.5580430            0.000000& r/ I& Y2 n2 r; w5 z) X
                           X( 2, 5)       0.2537760            0.000000/ A- l2 q2 c: N+ @1 k& a4 C$ P
                           X( 2, 6)       0.4651700            0.000000) E6 `# J9 o1 j3 X1 E
                           X( 3, 1)       0.1626000E-01        0.000000
    4 C$ @. e' B) q, ^                       X( 3, 2)       0.4213000E-02        0.0000002 M( W; _& w$ G( [) S7 K9 K
                           X( 3, 3)       0.8442300E-01        0.000000
    3 N, f: G; v% w' H- n/ E                       X( 3, 4)       0.6314420            0.000000. v. h% o3 s5 ^, e7 f5 x
                           X( 3, 5)       0.3751740            0.000000
    $ p, M& s2 E% u1 Q! T                       X( 3, 6)       0.5355790            0.000000( g, W+ J" a2 J- J( n% x1 _! D
                           X( 4, 1)       0.4587000E-02        0.000000
    4 k2 E: w1 t- r                       X( 4, 2)        0.000000            0.000000
    % v, s# U" K$ X2 B7 I( e" P+ c                       X( 4, 3)       0.8615400E-01        0.0000000 p0 F; ~, |& ~8 G7 e3 {( x# ^
                           X( 4, 4)       0.5575220            0.000000$ Q- h+ e: f( g* ~  l
                           X( 4, 5)       0.2105910            0.0000009 l; g2 B2 E% j4 n+ c. v
                           X( 4, 6)       0.4647210            0.000000, A& ~4 ~9 A1 q1 I8 _3 ?% T
                           X( 5, 1)       0.3988000E-01        0.000000
    ' A) {  A; @0 k$ ~# \4 X                       X( 5, 2)       0.1533300E-01        0.0000006 n; ~' X  A$ K8 R' }  r/ M
                           X( 5, 3)       0.7792100E-01        0.000000% A5 n7 l% C8 l9 w' O
                           X( 5, 4)        1.000000            0.000000
    $ U' P& n+ \6 `                       X( 5, 5)       0.4887600            0.000000$ Z' B$ \# I: E0 r5 O
                           X( 5, 6)       0.8885230            0.0000000 h8 I3 R% U) H7 B0 v1 H
                           X( 6, 1)       0.3290000E-02        0.000000; z# }2 S4 H/ P  J7 T+ L; y
                           X( 6, 2)       0.2399000E-02        0.000000
    3 h7 P1 S7 a) u# Q                       X( 6, 3)       0.2344100E-01        0.0000009 E6 U9 g" x- S" `5 T3 k7 N- S2 s; g
                           X( 6, 4)       0.2998440            0.0000000 t7 c; l2 W6 b
                           X( 6, 5)        1.000000            0.000000- Y2 c* }  a! v5 w. [) U$ T6 \7 s
                           X( 6, 6)       0.2178140            0.000000
    4 B+ g4 C$ L7 h0 i8 u' E7 X                       X( 7, 1)       0.2175000E-01        0.000000
    : F3 C& o# Z5 U                       X( 7, 2)       0.2329600E-01        0.000000
    ! s1 C- N- g, S% c" O' e9 ?$ Z/ `                       X( 7, 3)       0.5572700E-01        0.000000
    8 U1 c+ Z1 u/ Z( z1 O  z                       X( 7, 4)       0.6350860            0.000000
    4 m; k- \5 G$ X5 k                       X( 7, 5)       0.6575610            0.000000
    ; |9 c+ o: P" E' @                       X( 7, 6)       0.5387230            0.000000& f0 q  F1 ~, P, O) w
                           X( 8, 1)        1.000000            0.0000002 ~- b4 W+ V  b; n  ~
                           X( 8, 2)        1.000000            0.000000& t$ X$ M# V. m! q, R# q& K0 E, B7 X
                           X( 8, 3)        1.000000            0.0000005 Q' {0 R* i# D' @8 F% u
                           X( 8, 4)        0.000000            0.000000
    9 w# F( Y2 _' I: _8 @" D4 o                       X( 8, 5)        0.000000            0.000000
    8 b( }* K# ?. P: a1 h                       X( 8, 6)       0.4305890            0.000000
    2 V0 ]7 D# k; a3 g% _% C                       X( 9, 1)       0.7096300E-01        0.000000
    ) p$ \# ]) l( i$ v                       X( 9, 2)       0.2775390            0.000000& n; j, m1 D) A+ u. m) s
                           X( 9, 3)       0.3859470            0.000000- r$ ^& {$ e9 h/ g
                           X( 9, 4)       0.1535660            0.000000
    + }- h1 [8 d1 Z! S+ l                       X( 9, 5)       0.1475700E-01        0.000000
    : R1 b0 ]6 |) X1 n3 d9 z0 w                       X( 9, 6)        0.000000            0.000000" S$ t9 t+ ~# N" q
                          X( 10, 1)       0.3816400E-01        0.000000* a: [4 n) o( z# G7 E" i
                          X( 10, 2)       0.3837300E-01        0.000000+ n) N0 b* }3 \% p: L
                          X( 10, 3)       0.2499000E-01        0.000000
    ! X1 n/ q8 ?% ^5 F4 d8 y                      X( 10, 4)       0.4898490            0.0000002 I( L4 B) K2 _9 }* O0 t
                          X( 10, 5)       0.3115900E-01        0.0000000 U# `  T$ V4 s' a8 M
                          X( 10, 6)        1.000000            0.000000: p2 z$ ]$ L4 V: H* G% J
                          X( 11, 1)        0.000000            0.000000* K9 k* @% b" g1 ~. j2 {+ A/ N/ k
                          X( 11, 2)       0.8735000E-02        0.000000
    , i4 J$ G* C2 M  h1 T' m: z- V9 V                      X( 11, 3)       0.1525900E-01        0.0000001 E& o% v) s! r
                          X( 11, 4)       0.2889120            0.000000
    # r8 a3 S  E" m; x: y  n/ J, m" T+ P                      X( 11, 5)       0.3700900E-01        0.000000
    - j# `: \, g! u                      X( 11, 6)       0.4091800E-01        0.0000009 J$ B0 t: s% \; u. l4 A
                          X( 12, 1)       0.5718100E-01        0.000000/ Z6 B+ x9 ~+ c2 s- j
                          X( 12, 2)       0.3850800E-01        0.0000007 u, f) J" G  n  V! l7 ?% G) L
                          X( 12, 3)       0.2541090            0.0000003 M7 z; n$ {  Z$ T" _6 N+ m; H" D
                          X( 12, 4)       0.1842790            0.000000
    ; N# m% n% t' D: {1 u                      X( 12, 5)       0.9316000E-01        0.000000
    : s5 Z; U. ]: p2 F                      X( 12, 6)       0.5224600E-01        0.000000
    2 \; ^( i. a+ L# Z" R                      X( 13, 1)       0.4821100E-01        0.0000004 V9 C5 L  H3 ~  z: j  C4 o3 Z
                          X( 13, 2)       0.9641000E-02        0.000000, W3 L+ g$ o9 Y7 T3 o
                          X( 13, 3)        0.000000            0.000000/ T3 o* ^4 z9 l* M: ?7 \
                          X( 13, 4)       0.3878190            0.000000
    " N8 ]4 ^" J8 [3 B/ h                      X( 13, 5)       0.1358340            0.0000009 ?" X: l" b7 Y
                          X( 13, 6)       0.1185130            0.000000
    1 b! o9 G1 ~; C- G9 A' ?  V
    5 M- f1 F. l( r+ e. u                            Row    Slack or Surplus      Dual Price& t) ~5 u$ H' r1 S
                                  1       0.7272518            1.000000
    - M4 k/ B- Q+ n                              2        0.000000           0.36362595 l/ @+ Z8 ^) q& V3 p" p. S% y$ E# H
                                  3        1.124802            0.000000# C9 R6 z. i- m# L4 o& Q+ x! ]" ^6 F
                                  4        1.121677            0.000000
    / ~& N! c' P' P1 F9 Z5 Y9 j& m6 u                              5        1.162172            0.000000, L9 F  |% _; P; E5 o- i
                                  6        1.631867            0.000000) V" [# ^% Y- v0 f9 W2 [  j
                                  7        1.386806            0.000000  {" a/ I5 t- A* ]* b: `
                                  8        1.628048            0.000000% m) D8 N8 w) e% i& _6 [
                                  9       0.8751980            0.000000% i8 ~$ @2 T8 }" X1 J8 D
                                 10       0.8783228            0.000000
    . G+ d  o- D: l                             11       0.8378279            0.000000
    * k) F) V6 e/ ?' s4 _                             12       0.3681325            0.000000
    - A, \4 l/ e8 K/ J                             13       0.6131944            0.000000
    % e& P6 ~( n1 c$ 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
    但是 我设置全局求解时：两三秒就停止了
    # {& a$ [1 C* R- n出现错误提示
    4 c* M" C, F' }% L8 i  @* x6 h( e还有 上面的 全局的做大值还没有局部做大的大  为什么！！！
    回复

    使用道具 举报

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

    回顶部