快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

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

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

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

7

主题

5

听众

28

积分

升级  24.21%

该用户从未签到
跳转到指定楼层
1#
发表于 2009-9-5 07:49 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
model:! R& |! {- o* P) q  f3 g- V: j* A
sets:) o+ R9 z' U6 m! R# C/ ^
shengshi/1..13/;8 @9 {  |/ y9 R" b
zhibiao/1..6/:a;
) ~- v3 h$ S. n. x  P: Dlinks(shengshi,zhibiao):x;9 ?& W: [, Z4 W4 y
endsets
. F8 k1 L* c+ {$ a* cdata:
2 D2 i7 J. p, k* K2 H1 hx=2 w  b: _5 y- v9 e& f# S
0.099599  0.104095  0.021186  0.489849  0.031419  0.399601 % u0 Y1 X6 ^8 o. |' B3 @( J! [
0.012222  0.006320  0.035075  0.558043  0.253776  0.465170 / {( g4 \, }; H
0.016260  0.004213  0.084423  0.631442  0.375174  0.535579
& |5 s1 A! Y8 Q# j6 i" T0.004587  0.000000  0.086154  0.557522  0.210591  0.464721 - A, o. m, J/ v1 e- T) p+ H
0.039880  0.015333  0.077921  1.000000  0.488760  0.888523
6 y: @6 d; T4 L5 L8 [0.003290  0.002399  0.023441  0.299844  1.000000  0.217814
# f3 A" X" w7 J& p( D0.021750  0.023296  0.055727  0.635086  0.657561  0.538723
4 X6 i9 Q, V2 n+ ~% k3 W: Q1.000000  1.000000  1.000000  0.000000  0.000000  0.430589 # ]* }$ m3 \7 r- s* q
0.070963  0.277539  0.385947  0.153566  0.014757  0.000000 0 P; @3 o0 x0 u
0.038164  0.038373  0.024990  0.489849  0.031159  1.000000
: e, T' L# f9 S, b. S1 ?1 b! m0.000000  0.008735  0.015259  0.288912  0.037009  0.040918
" E- y% [6 F+ Y4 ^$ X6 b9 @0.057181  0.038508  0.254109  0.184279  0.093160  0.052246 ! _. y6 t6 P# `5 B, R  Q3 `
0.048211  0.009641  0.000000  0.387819  0.135834  0.118513;* a" O1 D4 L- _
enddata0 `) z1 G* m" R* f- a
max=
6 M  y( w1 ]: u& U( \(@sum(shengshi(i):' Q/ |. ]7 F: k' w- x
(@sum(links(i,j):a(j)*x(i,j))-@sum(shengshi(i):
6 `; K3 c0 t" N2 E' ~@sum(links(i,j):a(j)*x(i,j)))/13)^2)/12)^(1/2);
+ H4 Y% w5 g) g' N* S& K' s; R# X! S6 E1 U" i' j$ }/ a
@sum(zhibiao(j):a(j)^2)=1;
. R+ A1 g4 Y% j3 V3 ]' n@for(zhibiao(j):a(j)>-1);: H9 F+ |/ B+ x
@for(zhibiao(j):a(j)<1);
. j$ e  [4 U1 o2 S+ i! u@for(zhibiao(j)free(a(j)));
1 ?. e" n/ I/ J% S* B  f1 Send
/ B, a  @  M) e' U! G0 R9 h求不出全局最优解，请大侠帮帮忙！！！
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.
    $ {2 U6 \  N  Z6 K5 `& @  e  Objective value:                             0.7272518& f* T) `; b) z  o  b& C$ E
      Objective bound:                             0.7272527
    - o3 q4 j5 A" _/ X1 \, H# u  Infeasibilities:                              0.000000
    ; n! l: K, |1 n" j  y  Extended solver steps:                           38843# U: r; V9 e  a% U9 f2 X* S* g
      Total solver iterations:                       7561347
    $ a. N+ S* F7 f+ h% W7 i- f" e# X3 ^6 }
    $ z% k* J9 z- G
                           Variable           Value        Reduced Cost& X: c1 D5 X( x1 j
                              A( 1)       0.1248020            0.000000; r: {% x" Q0 \; A# n% h& n  |
                              A( 2)       0.1216772            0.000000
    1 i  M2 e0 C7 ~1 q9 z/ x                          A( 3)       0.1621721            0.000000  G7 M- Q. W+ C  g  n
                              A( 4)       0.6318675            0.000000
    : Y1 ?% ^# ]- G5 ~- q. Q) H# L                          A( 5)       0.3868056            0.0000007 V- v8 `: {+ r3 k/ z6 A
                              A( 6)       0.6280479            0.000000
    ; O9 d9 I+ I' l; @2 i) I- \                       X( 1, 1)       0.9959900E-01        0.000000
    % a6 i0 g) S+ T4 R, E6 b. v; |                       X( 1, 2)       0.1040950            0.000000; M$ Q+ Z- l4 y4 R- L+ ]9 s
                           X( 1, 3)       0.2118600E-01        0.000000- [' N# G1 c9 u4 M/ X
                           X( 1, 4)       0.4898490            0.000000$ k* G2 f9 K( U6 B; `
                           X( 1, 5)       0.3141900E-01        0.000000
    6 N6 [8 D4 y+ B8 T                       X( 1, 6)       0.3996010            0.0000005 N2 K. m  f( b8 y& r2 |+ m$ z
                           X( 2, 1)       0.1222200E-01        0.0000006 o5 T: \3 M4 k  f+ J$ T
                           X( 2, 2)       0.6320000E-02        0.000000) h; Y, l' \. d& R$ Q$ N4 e8 o+ w
                           X( 2, 3)       0.3507500E-01        0.000000/ A. J# y, X5 w# o9 T1 ]: h
                           X( 2, 4)       0.5580430            0.000000# Q/ h6 ]9 r! |) {1 W( A
                           X( 2, 5)       0.2537760            0.000000
    " ~7 x. i8 d+ e* J1 m7 y/ h                       X( 2, 6)       0.4651700            0.000000
    + i3 b4 b( m1 I* M3 q3 A                       X( 3, 1)       0.1626000E-01        0.000000
    ' F' n2 t5 i- {8 ]$ |, e/ p                       X( 3, 2)       0.4213000E-02        0.000000
    2 b4 o3 H# x: ^6 x                       X( 3, 3)       0.8442300E-01        0.0000006 X8 W# o( Q# l; p
                           X( 3, 4)       0.6314420            0.000000
    / P" C  j% \9 I. d4 l6 V5 j: J, v                       X( 3, 5)       0.3751740            0.000000
    7 d2 C- ~: j9 F( b1 f                       X( 3, 6)       0.5355790            0.000000# d# [7 s8 y+ y
                           X( 4, 1)       0.4587000E-02        0.000000
    ; l0 r2 @4 S: g1 `$ Q9 S4 B) v& H( D4 {                       X( 4, 2)        0.000000            0.0000009 r) g) y9 F, F2 s  [! A+ F
                           X( 4, 3)       0.8615400E-01        0.0000003 T( n. v, c1 H/ P8 k8 _# Q6 s- I
                           X( 4, 4)       0.5575220            0.000000
    ' @8 v. T) _2 `# u                       X( 4, 5)       0.2105910            0.0000000 y$ k' ?- W1 ?- ]8 `& s7 N4 p4 b- m
                           X( 4, 6)       0.4647210            0.000000; l; o; X. ?- M2 v1 V/ m
                           X( 5, 1)       0.3988000E-01        0.0000001 I5 i. D6 c- i9 Y8 f5 E
                           X( 5, 2)       0.1533300E-01        0.000000
    6 v; m' o8 b* m                       X( 5, 3)       0.7792100E-01        0.000000' _+ a9 K: O8 w9 f2 z
                           X( 5, 4)        1.000000            0.000000+ m% S* B8 d+ u7 ]4 l
                           X( 5, 5)       0.4887600            0.000000
    4 r6 q6 N$ ~# N+ J; G5 X                       X( 5, 6)       0.8885230            0.0000008 x$ j' _" E- q7 w# I" n0 x4 ?
                           X( 6, 1)       0.3290000E-02        0.000000
      F' k. s- b, ^) j( t6 @                       X( 6, 2)       0.2399000E-02        0.000000
    # U( w$ J: @1 d6 P# D! i9 O3 ]0 M                       X( 6, 3)       0.2344100E-01        0.000000- i- c  l. \# j8 \1 e. v
                           X( 6, 4)       0.2998440            0.000000
    9 J. {- m" ~- b4 c                       X( 6, 5)        1.000000            0.000000
    + J1 N/ M( g8 ~0 D                       X( 6, 6)       0.2178140            0.0000001 a* ^+ {% Z+ _
                           X( 7, 1)       0.2175000E-01        0.000000# b. W; E) K6 ^
                           X( 7, 2)       0.2329600E-01        0.000000( S% k1 ~: X. l1 j% M) U3 |& B
                           X( 7, 3)       0.5572700E-01        0.000000
    ! k6 a( Y1 ~1 _( R7 M5 J5 m                       X( 7, 4)       0.6350860            0.000000
    $ a7 Y: e8 S" H                       X( 7, 5)       0.6575610            0.000000% ?& G8 C: E. a
                           X( 7, 6)       0.5387230            0.000000
    * S. K( o8 J1 j8 M8 k* Z* w                       X( 8, 1)        1.000000            0.000000
    : W6 N0 A7 u0 _( Y# W/ n                       X( 8, 2)        1.000000            0.000000, G) T; ?  }; F3 R0 `: G
                           X( 8, 3)        1.000000            0.000000
    " d9 U. i% V, ~% \7 m                       X( 8, 4)        0.000000            0.000000
    5 W$ c  i3 B- z" {( n% j                       X( 8, 5)        0.000000            0.0000001 `0 b" v) a: `: m% _1 |' Z' k
                           X( 8, 6)       0.4305890            0.000000/ l# B+ P9 L( k- ?+ v
                           X( 9, 1)       0.7096300E-01        0.0000009 C6 M$ U/ y8 {; g0 X+ d6 X8 m
                           X( 9, 2)       0.2775390            0.000000
    8 U  r# T( B, s! |3 e6 L0 |" q                       X( 9, 3)       0.3859470            0.000000
    7 ~' i7 Z. P9 M4 _* ~  y  b9 J                       X( 9, 4)       0.1535660            0.000000
    - e/ H* H& q* J& a6 M% C9 r                       X( 9, 5)       0.1475700E-01        0.000000/ [( F! `$ T) f$ X: \/ L) f3 a- A
                           X( 9, 6)        0.000000            0.000000
    6 _. _7 K" [4 a* h5 S                      X( 10, 1)       0.3816400E-01        0.000000
    6 Z! m/ L* M) p, m                      X( 10, 2)       0.3837300E-01        0.000000. U9 l# Q- k! M8 {5 A
                          X( 10, 3)       0.2499000E-01        0.000000
    % D! @2 C6 M: d) h: r                      X( 10, 4)       0.4898490            0.0000007 D% q$ j2 ^. n/ G" B- t2 p5 H
                          X( 10, 5)       0.3115900E-01        0.000000
    8 p# x& ^* e& w7 ]2 t3 _' {                      X( 10, 6)        1.000000            0.000000+ T9 l+ n- r: c( n. c* v
                          X( 11, 1)        0.000000            0.000000
    & j5 t3 \' a7 y4 p4 f, j                      X( 11, 2)       0.8735000E-02        0.000000
    + @. ]; n* ?0 a1 k4 K                      X( 11, 3)       0.1525900E-01        0.000000" _! p; ?& _) ^2 ?' B7 }( s5 _4 x
                          X( 11, 4)       0.2889120            0.000000
    , C7 J9 i' V3 Y" }# i                      X( 11, 5)       0.3700900E-01        0.000000* v  Z! P( D1 @* S$ w
                          X( 11, 6)       0.4091800E-01        0.000000- c4 `: \$ {( ^, I4 Q
                          X( 12, 1)       0.5718100E-01        0.000000
    + J* s  d9 i5 J0 z- }                      X( 12, 2)       0.3850800E-01        0.000000* m! L( n9 s2 p/ C" w( s+ J
                          X( 12, 3)       0.2541090            0.000000
    5 P  N0 h8 I) n4 V                      X( 12, 4)       0.1842790            0.000000& k. B8 {- [+ q8 {, D# i5 S
                          X( 12, 5)       0.9316000E-01        0.000000
    / N/ n) ~8 r( q) k: x, Y                      X( 12, 6)       0.5224600E-01        0.000000
    2 M7 U' F0 A- o1 l                      X( 13, 1)       0.4821100E-01        0.000000/ d8 b9 ^# u, {) A$ w5 m# Q
                          X( 13, 2)       0.9641000E-02        0.000000
    3 h4 m2 R* ?- j( x+ ~* \                      X( 13, 3)        0.000000            0.0000004 x+ v+ _: `. L9 e; ?
                          X( 13, 4)       0.3878190            0.000000
    % \% e% h# E9 }. ?                      X( 13, 5)       0.1358340            0.000000
    7 Q7 H; h; L- G: L2 t: V                      X( 13, 6)       0.1185130            0.000000
    " V; k: T; X3 `9 B
    : w# a+ q5 U  }8 W0 C- s; q+ [8 o                            Row    Slack or Surplus      Dual Price! ~; S+ d8 [$ a! x7 W8 _
                                  1       0.7272518            1.000000
    4 P. H  k. ]; q! q( i& l! R                              2        0.000000           0.3636259" A" M: F# ~7 v" i+ u
                                  3        1.124802            0.000000
      Q, r1 b* ~) L0 M                              4        1.121677            0.000000
    4 Y0 W7 a- }5 D9 }3 B, G                              5        1.162172            0.000000
    1 c2 i6 Y  r% V. }  o                              6        1.631867            0.000000
    : M( ^" ~7 o- @5 ?- Q! V                              7        1.386806            0.000000$ Y0 Z+ l6 L% p: a& A7 ?& b! K
                                  8        1.628048            0.000000
    ( U- f' k. R9 S- X' u6 E                              9       0.8751980            0.000000
      l' ]1 X* V2 ~! `$ n2 Q) B# a                             10       0.8783228            0.000000
    + S( V- O+ }! p                             11       0.8378279            0.0000008 Y1 J, h# `) T# R* a
                                 12       0.3681325            0.000000
    6 v, a0 Q3 a1 x# {- m8 i                             13       0.6131944            0.000000
    ( D9 ^# N0 B' z3 X                             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$ C/ W9 h/ W' L& p6 ]9 }出现错误提示- ?* G' m) N, t7 j4 B' L8 _
    还有 上面的 全局的做大值还没有局部做大的大  为什么！！！
    回复

    使用道具 举报

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

    回顶部