快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

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

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

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

7

主题

5

听众

28

积分

升级  24.21%

该用户从未签到
跳转到指定楼层
1#
发表于 2009-9-5 07:49 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
model:
3 Y( b# Q$ U* J) L% |sets:; |3 E2 k1 s4 Q
shengshi/1..13/;
& f# Z! L' n0 f# L3 R3 M% fzhibiao/1..6/:a;4 Y3 p$ a! `5 K5 Z3 r
links(shengshi,zhibiao):x;
5 G- c# n8 E& w8 M2 {) j0 Mendsets+ u4 d+ h5 a, w4 @( W0 D5 j; F% ~
data:3 }: p5 O5 ^! z, b
x=, s  E7 S% O" u5 R
0.099599  0.104095  0.021186  0.489849  0.031419  0.399601 3 O4 p0 r$ ]" I" B) {1 }
0.012222  0.006320  0.035075  0.558043  0.253776  0.465170
8 ?5 ~1 G6 o% N3 s3 H; _' t0.016260  0.004213  0.084423  0.631442  0.375174  0.535579 & e0 h' m/ \! {" ]8 w
0.004587  0.000000  0.086154  0.557522  0.210591  0.464721 5 ]1 X& a) v0 }$ |' I6 z' f* C
0.039880  0.015333  0.077921  1.000000  0.488760  0.888523 ! _# \% M& J$ F% R- L
0.003290  0.002399  0.023441  0.299844  1.000000  0.217814 2 l8 g# B9 ~$ r  K) x# y$ j, i- f
0.021750  0.023296  0.055727  0.635086  0.657561  0.538723 . ]7 Y8 I  U* F8 w% C/ V, D
1.000000  1.000000  1.000000  0.000000  0.000000  0.430589 : Z( Z* U4 u$ U$ n8 s
0.070963  0.277539  0.385947  0.153566  0.014757  0.000000
! N' @% H, C- g# l& U0.038164  0.038373  0.024990  0.489849  0.031159  1.000000 % d0 g, C2 I. @. F4 U# ^
0.000000  0.008735  0.015259  0.288912  0.037009  0.040918
' U! X: o% v  _0.057181  0.038508  0.254109  0.184279  0.093160  0.052246
& \* b+ D# }* `7 D9 @1 M0 O! D, t0.048211  0.009641  0.000000  0.387819  0.135834  0.118513;
$ ~7 D, r9 T& h; Menddata
0 Z6 `5 H) D/ Vmax=/ J% R# M% ~6 l$ X
(@sum(shengshi(i):# k7 S" D, m; H) C* j
(@sum(links(i,j):a(j)*x(i,j))-@sum(shengshi(i):( o9 c8 ?5 E3 j9 b( U+ X
@sum(links(i,j):a(j)*x(i,j)))/13)^2)/12)^(1/2);- o* P6 k" Z- G' {
" \" |( ^  \* p* M" v% G' ~
@sum(zhibiao(j):a(j)^2)=1;: `  H3 n/ E! Q; k! d$ B3 d: I0 f
@for(zhibiao(j):a(j)>-1);
# m6 P3 A- Q+ m% B% T% o@for(zhibiao(j):a(j)<1);. z# ?9 f, m6 M* t4 o8 O
@for(zhibiao(j)free(a(j)));7 Z* O0 {/ ^: O0 e8 n, `  `
end
0 T% [* h  J0 o, 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.
    7 B) e7 ~% e, L6 s  Objective value:                             0.72725183 w. _! m. o8 J5 H& e
      Objective bound:                             0.72725273 m3 e" M1 d+ c8 q: q* Y% L; G
      Infeasibilities:                              0.000000) {9 r0 S9 l* s$ R0 C
      Extended solver steps:                           388433 R2 q% t' f; a
      Total solver iterations:                       7561347
    & s) I: i1 p- m; d+ E0 v& j' p$ _! J. [9 Y) d, ~% W$ v
    : {3 K' c- a' |
                           Variable           Value        Reduced Cost2 H# r+ p1 r! V1 U6 \/ b1 c
                              A( 1)       0.1248020            0.000000  v, o+ Q" q  @7 d
                              A( 2)       0.1216772            0.000000. E0 U) L2 |4 n2 W
                              A( 3)       0.1621721            0.000000; d4 r3 r  X- q' {% C3 @3 d
                              A( 4)       0.6318675            0.000000
      e( ]/ y7 e; d4 }$ _                          A( 5)       0.3868056            0.000000
    + h6 w) l+ e$ Y& N9 n                          A( 6)       0.6280479            0.000000
    1 B) M5 u* E7 O                       X( 1, 1)       0.9959900E-01        0.000000
    8 G* e# D( p& h) e9 N' V1 u                       X( 1, 2)       0.1040950            0.000000
    & n2 g- |' v3 K- @                       X( 1, 3)       0.2118600E-01        0.000000
    . @, w8 n  O6 M! [) u                       X( 1, 4)       0.4898490            0.0000006 [7 @, R5 f- c( w) P' J1 ~
                           X( 1, 5)       0.3141900E-01        0.000000
    : Q" K( |$ h3 S: Y8 J  c6 k                       X( 1, 6)       0.3996010            0.0000008 K+ ^9 E2 b5 [, G% M' c
                           X( 2, 1)       0.1222200E-01        0.000000' `) d6 K, B& f" b. ]- j8 G& V- ^, W" F
                           X( 2, 2)       0.6320000E-02        0.000000
    1 ^: D3 R. V- }' h                       X( 2, 3)       0.3507500E-01        0.000000
    5 X' h0 ~. I& n6 l7 M0 F1 z8 [9 |                       X( 2, 4)       0.5580430            0.0000008 Y: Y1 E) g5 }, E
                           X( 2, 5)       0.2537760            0.000000& P, N: \" \4 W8 ^% b9 D3 R
                           X( 2, 6)       0.4651700            0.000000
    + G9 g2 L$ D+ j  ~, @. ^2 K$ G                       X( 3, 1)       0.1626000E-01        0.000000; o' x+ Z2 n5 c- T( r+ M/ Q: e3 ~
                           X( 3, 2)       0.4213000E-02        0.000000
    1 i7 z# }8 `: V- P% Z5 z                       X( 3, 3)       0.8442300E-01        0.000000
    ; S4 ^( K8 P0 K. j5 y4 Y/ m3 I, T                       X( 3, 4)       0.6314420            0.000000
    6 w+ h3 h* Y7 |4 x                       X( 3, 5)       0.3751740            0.000000
    2 w( q& n1 ?# S; v' O# v0 b4 j. X                       X( 3, 6)       0.5355790            0.000000  t% L+ u1 K" X/ P' A! a; i; `
                           X( 4, 1)       0.4587000E-02        0.0000003 E% g$ U8 @  U
                           X( 4, 2)        0.000000            0.0000003 M# m: |0 C* T- o4 {4 b
                           X( 4, 3)       0.8615400E-01        0.0000002 D) |. B, K! j  p% w' z& G6 t
                           X( 4, 4)       0.5575220            0.000000
    ) b2 \& r2 k. z: H                       X( 4, 5)       0.2105910            0.000000; g" Y1 ?- K& A
                           X( 4, 6)       0.4647210            0.0000000 o3 K8 d! j! d% q% D8 U$ n* L
                           X( 5, 1)       0.3988000E-01        0.000000% \9 f- R1 f& w" A- I- K
                           X( 5, 2)       0.1533300E-01        0.000000
    ! c6 c- p1 I, e* ^  n( [+ k                       X( 5, 3)       0.7792100E-01        0.000000( Z; L0 ^( |5 b7 i4 ^5 S: t+ G
                           X( 5, 4)        1.000000            0.000000% j7 C/ a- _( e
                           X( 5, 5)       0.4887600            0.000000
    : j( q2 ]8 ?, X                       X( 5, 6)       0.8885230            0.0000008 A# E, H! s" S; t8 n
                           X( 6, 1)       0.3290000E-02        0.000000
    1 f7 Z; I+ O. G. n, p9 n                       X( 6, 2)       0.2399000E-02        0.000000, H& [, W; ~: k
                           X( 6, 3)       0.2344100E-01        0.000000
    % X- W/ W" F" F' G0 H; s; ~                       X( 6, 4)       0.2998440            0.000000
    ( a0 V( ]) O0 o" }. V6 p: L. t) r2 {                       X( 6, 5)        1.000000            0.000000
    - x$ d) `% e; M. ?9 z( d                       X( 6, 6)       0.2178140            0.000000
    + |0 F: G3 l( U2 z1 l                       X( 7, 1)       0.2175000E-01        0.000000
    : r: p! Y9 G! U3 @1 A                       X( 7, 2)       0.2329600E-01        0.000000
    : O$ K4 O% @# Z- C# O! ?; g% I                       X( 7, 3)       0.5572700E-01        0.000000
    ) }9 I. o* X# T  L                       X( 7, 4)       0.6350860            0.000000+ t3 M9 D( ?1 z5 j- }2 Y: x
                           X( 7, 5)       0.6575610            0.000000
    6 [  p# D# s8 c0 O3 Y                       X( 7, 6)       0.5387230            0.000000
    % F% Y  n6 N& u% \( n                       X( 8, 1)        1.000000            0.000000, F5 \3 g3 [& G  n4 ?# Y
                           X( 8, 2)        1.000000            0.000000
    5 D* q0 b/ j6 |; t1 R                       X( 8, 3)        1.000000            0.000000
    6 ]+ Q! z1 U. X6 J. K$ W                       X( 8, 4)        0.000000            0.000000- u7 h; b+ ]1 {. \0 |
                           X( 8, 5)        0.000000            0.000000
    - X! R2 _4 [' }1 o( [& s+ g3 ]9 }                       X( 8, 6)       0.4305890            0.000000
    6 E0 M: A9 s0 y" w. G                       X( 9, 1)       0.7096300E-01        0.000000; W- o/ T6 j4 h. d# P$ X- M; f) a
                           X( 9, 2)       0.2775390            0.000000' s  k7 K3 D" P' v6 G
                           X( 9, 3)       0.3859470            0.0000006 }  M# j: A. S- `
                           X( 9, 4)       0.1535660            0.0000003 u1 L; U6 d4 l- }+ R/ s. z
                           X( 9, 5)       0.1475700E-01        0.000000' [8 r& S: c: F: l, Q% h6 w
                           X( 9, 6)        0.000000            0.000000. c: k1 X$ Y5 V! U
                          X( 10, 1)       0.3816400E-01        0.000000+ U7 T8 L* v2 i: F3 Y
                          X( 10, 2)       0.3837300E-01        0.000000
    3 @0 n0 ?5 ~$ e* U+ c6 r                      X( 10, 3)       0.2499000E-01        0.0000009 q/ Q. ~9 ]1 c; M& N
                          X( 10, 4)       0.4898490            0.000000* i' P, C: s# e" z8 a% S9 |8 E
                          X( 10, 5)       0.3115900E-01        0.000000
    ! l3 F9 e& \2 h8 s4 @  O, d- `. _                      X( 10, 6)        1.000000            0.0000006 k3 y! ^, u& W! K( ^$ v4 M
                          X( 11, 1)        0.000000            0.000000
    - A  |  |2 ]* |2 N* O6 s                      X( 11, 2)       0.8735000E-02        0.000000
    % Z$ Y0 c, e0 w7 p3 y5 I                      X( 11, 3)       0.1525900E-01        0.0000006 ?- i1 X3 A# a$ }# r
                          X( 11, 4)       0.2889120            0.000000
    * v' Y( O* J% o3 c0 R* [                      X( 11, 5)       0.3700900E-01        0.0000006 ]7 ^! @9 Q: S* H: }- Q
                          X( 11, 6)       0.4091800E-01        0.000000
    / b. l9 l0 `* m( f0 j                      X( 12, 1)       0.5718100E-01        0.0000003 Y7 R- w; j2 E1 E
                          X( 12, 2)       0.3850800E-01        0.000000
    - {2 W6 Z. t2 C3 {& n, p5 ~                      X( 12, 3)       0.2541090            0.000000: o. G, D4 W: z, m8 V! R
                          X( 12, 4)       0.1842790            0.000000. R: B7 {: x5 l
                          X( 12, 5)       0.9316000E-01        0.000000
    % o; F- {+ K! J- g                      X( 12, 6)       0.5224600E-01        0.000000+ n8 }0 N7 m" n' B: A; u
                          X( 13, 1)       0.4821100E-01        0.000000
    % @2 o: {+ _* i0 ?! B7 {: O2 M                      X( 13, 2)       0.9641000E-02        0.000000
    ) W0 w" b+ f8 [+ s4 S6 \. W' Y9 Y3 S                      X( 13, 3)        0.000000            0.0000000 i1 e* a& G7 |' s. U
                          X( 13, 4)       0.3878190            0.000000
    % X0 Z4 f% y6 o- z+ k                      X( 13, 5)       0.1358340            0.0000009 D* _! p# b: l+ I8 F* C
                          X( 13, 6)       0.1185130            0.000000
    . q" _  p4 k9 D2 X; _( M6 a3 t$ F
                                Row    Slack or Surplus      Dual Price
    9 X* u3 o; }" h6 _' O! R* ]) c                              1       0.7272518            1.000000% b/ I1 `' j9 [5 o7 V7 v. ^+ \! B
                                  2        0.000000           0.3636259
    6 R; ?' Y% ]& h! h2 o2 z                              3        1.124802            0.000000
    . u' j* m. Y7 X. j% s                              4        1.121677            0.000000
    9 k1 ~# ^# q7 n) \                              5        1.162172            0.000000
    / D; M9 D! `  K1 P: }, x                              6        1.631867            0.000000% c  b& J% K6 L7 L; r9 u( ~
                                  7        1.386806            0.000000
    # E2 S8 P4 G4 _9 W                              8        1.628048            0.000000) v4 g( `$ W$ y" J/ B/ R
                                  9       0.8751980            0.000000) N; P7 D: Q! {: @/ c
                                 10       0.8783228            0.000000# X0 D( r, H3 |* W
                                 11       0.8378279            0.0000009 o# J6 a% P# o$ P3 o
                                 12       0.3681325            0.000000
    ; v' T7 W0 _0 ?# F! d                             13       0.6131944            0.000000
    ! B2 I  L: ?- v6 z7 R- {, e                             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
    但是 我设置全局求解时：两三秒就停止了
    & W% R) W0 N- }7 d  s6 i出现错误提示' N+ T* N7 ]7 t6 i8 E9 q9 w( H
    还有 上面的 全局的做大值还没有局部做大的大  为什么！！！
    回复

    使用道具 举报

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

    回顶部