快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

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

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

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

7

主题

5

听众

28

积分

升级  24.21%

该用户从未签到
跳转到指定楼层
1#
发表于 2009-9-5 07:49 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
model:& r& |% K/ P8 c
sets:
) ?" v/ C. N3 D9 f3 m3 xshengshi/1..13/;
, X! U( T! ]! y3 w, z+ C5 wzhibiao/1..6/:a;
( Q; q1 k# r: d% N4 n" ~$ N+ nlinks(shengshi,zhibiao):x;
' i! @' |1 h& `endsets
' y3 ~- n" ]- l& Y, H. Ydata:% \% o, e. J  S% ]& W
x=- B* C1 W# |3 O5 I2 `3 o
0.099599  0.104095  0.021186  0.489849  0.031419  0.399601 3 z; D" a. c3 v1 U+ x3 Y
0.012222  0.006320  0.035075  0.558043  0.253776  0.465170
1 k/ f2 `5 _5 p: e- x* Z; L% m0.016260  0.004213  0.084423  0.631442  0.375174  0.535579 " N1 i3 U* d1 n& U' W7 I, B
0.004587  0.000000  0.086154  0.557522  0.210591  0.464721
6 ~9 k6 K8 ^3 H' `9 T7 O+ s9 l" G. m0.039880  0.015333  0.077921  1.000000  0.488760  0.888523 7 f/ b( `- t0 _. X2 B7 k1 ]
0.003290  0.002399  0.023441  0.299844  1.000000  0.217814
. L% r; ?' W9 Z  t4 p; i) {  T0.021750  0.023296  0.055727  0.635086  0.657561  0.538723
( g1 e7 ]) a7 z5 h& U" l, o1.000000  1.000000  1.000000  0.000000  0.000000  0.430589
6 l3 \% D2 c- H' K3 u0.070963  0.277539  0.385947  0.153566  0.014757  0.000000
  R( m% ]' [; r  C0.038164  0.038373  0.024990  0.489849  0.031159  1.000000
7 G& y% ~/ D) S  F0.000000  0.008735  0.015259  0.288912  0.037009  0.040918
6 A! y6 A: M* ]4 x( Y. ]3 v0.057181  0.038508  0.254109  0.184279  0.093160  0.052246 3 a9 m( j7 b, a  e
0.048211  0.009641  0.000000  0.387819  0.135834  0.118513;
' p: E$ n- {! xenddata3 Q5 I4 V  v2 i6 E" {& K
max=
! e  F) N; M# W- P: N(@sum(shengshi(i):
3 K4 E, V* D+ L$ V4 I: Y(@sum(links(i,j):a(j)*x(i,j))-@sum(shengshi(i):
1 f" D4 s* \; E/ A$ ]5 p@sum(links(i,j):a(j)*x(i,j)))/13)^2)/12)^(1/2);
- O1 b6 b$ y' u" m' y7 y/ M' V% N+ b4 m' O6 i% J
@sum(zhibiao(j):a(j)^2)=1;
' Z! B& Q# Z: y4 V/ ~@for(zhibiao(j):a(j)>-1);% f1 R, u0 N" S/ f* D( K
@for(zhibiao(j):a(j)<1);" q$ p8 }& C+ O- ~% Z6 D
@for(zhibiao(j)free(a(j)));
! C+ ~: [$ F- F" i' gend
. w) F( i: Y: e& n7 ^求不出全局最优解，请大侠帮帮忙！！！
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.
    ' A* q7 L' k7 e: L  Objective value:                             0.72725186 B6 O' i/ J9 _2 g4 M) `
      Objective bound:                             0.7272527, M$ b/ w/ |: t2 X) V
      Infeasibilities:                              0.000000& j( H2 |3 i) |: ]6 g5 s/ j3 ^
      Extended solver steps:                           38843$ Z: O5 X5 A- A4 S) ]/ h7 A
      Total solver iterations:                       7561347
    " v: V3 c  z) W
    ( e0 ^# _0 |3 l! P. F& D" w3 F5 f0 t' n) h4 ]7 O
                           Variable           Value        Reduced Cost
    : r3 `. A1 Y/ m! {& Y$ S                          A( 1)       0.1248020            0.0000006 ~  L" P5 h  H& m
                              A( 2)       0.1216772            0.000000
    # q7 Z! H: b$ u! k                          A( 3)       0.1621721            0.000000; C$ g5 W" J" c
                              A( 4)       0.6318675            0.000000
    $ T  J/ a2 p' h/ k7 U                          A( 5)       0.3868056            0.000000- ]7 e% l5 [: y& c  ^9 t  q
                              A( 6)       0.6280479            0.0000001 x) ~$ b4 ?6 T+ L  P1 p9 K
                           X( 1, 1)       0.9959900E-01        0.000000
      C: Z7 Y9 @) @& h' M                       X( 1, 2)       0.1040950            0.000000* Y# f5 h6 o% b* k
                           X( 1, 3)       0.2118600E-01        0.000000
    2 ]9 b- m! `5 W. c" s                       X( 1, 4)       0.4898490            0.000000" M8 _* m$ X) z
                           X( 1, 5)       0.3141900E-01        0.000000( E  O& l& d& W
                           X( 1, 6)       0.3996010            0.0000005 a6 U2 _; m" Z5 c
                           X( 2, 1)       0.1222200E-01        0.000000, e3 ]( k! T# |
                           X( 2, 2)       0.6320000E-02        0.000000
    , l, U: L0 b3 b# n) I, z. K                       X( 2, 3)       0.3507500E-01        0.000000/ _0 [! S' l* J) _- z
                           X( 2, 4)       0.5580430            0.000000
    " A5 u; s0 S4 _9 B! G  ]8 u                       X( 2, 5)       0.2537760            0.000000+ ~$ N( ^+ }/ N8 r8 U. H0 U. B
                           X( 2, 6)       0.4651700            0.000000
    ; @7 D& D+ Z( d7 v5 S1 w                       X( 3, 1)       0.1626000E-01        0.0000003 r) ]9 a& K% U
                           X( 3, 2)       0.4213000E-02        0.000000* ^# K- Q. ?2 }( j9 l
                           X( 3, 3)       0.8442300E-01        0.000000! x7 p9 ?1 x7 k
                           X( 3, 4)       0.6314420            0.000000
    ( H2 R0 N; u1 u: @. o5 k4 t  n                       X( 3, 5)       0.3751740            0.000000( r; @( T4 T2 E  {" G$ E9 S* m
                           X( 3, 6)       0.5355790            0.000000  c, u. H' u2 a1 }
                           X( 4, 1)       0.4587000E-02        0.000000
    ; g! {6 T7 ]. T                       X( 4, 2)        0.000000            0.000000
    $ g, f3 J# t/ v8 z" u3 C                       X( 4, 3)       0.8615400E-01        0.000000
    1 u7 K* ?3 P# P# v5 n                       X( 4, 4)       0.5575220            0.000000
    4 r: j+ W% A3 Z$ Y  v                       X( 4, 5)       0.2105910            0.000000
    3 K6 h7 a% E7 y; h* Q                       X( 4, 6)       0.4647210            0.000000& G9 n8 q4 d; G9 k2 B
                           X( 5, 1)       0.3988000E-01        0.0000009 V5 U9 M! ?$ j& k' `( |! n
                           X( 5, 2)       0.1533300E-01        0.000000( M( N. m7 y0 J/ x& c- l4 z
                           X( 5, 3)       0.7792100E-01        0.000000- k6 y3 |; v) F) W
                           X( 5, 4)        1.000000            0.000000
    / {0 X9 W* t- F  v0 k/ }                       X( 5, 5)       0.4887600            0.000000
    & R! `% a5 S* y* R% @  W) v5 V% y                       X( 5, 6)       0.8885230            0.0000002 b$ Z+ w% E, H2 |9 h
                           X( 6, 1)       0.3290000E-02        0.000000
    ; a- I: R. \! n& p5 h; Q5 c                       X( 6, 2)       0.2399000E-02        0.000000
    2 n/ [- e/ r+ w9 B" I4 g! y& {                       X( 6, 3)       0.2344100E-01        0.0000004 B$ `9 K* D2 n; e2 j; c
                           X( 6, 4)       0.2998440            0.000000. Q2 h+ @2 v5 {: o8 W' z) g' n
                           X( 6, 5)        1.000000            0.000000  g9 [5 N: w" \# v
                           X( 6, 6)       0.2178140            0.000000
    ' p) N) a, J: H# |, L/ {                       X( 7, 1)       0.2175000E-01        0.000000
      ~: o+ a+ R3 P' Q5 h* o' Z                       X( 7, 2)       0.2329600E-01        0.000000
    * [+ F' P3 S' ]+ Q7 `8 T                       X( 7, 3)       0.5572700E-01        0.000000
    . I$ W* A: p# j  ?+ Q3 F* |; B0 B                       X( 7, 4)       0.6350860            0.000000- C3 g( r/ d+ i
                           X( 7, 5)       0.6575610            0.000000, I* N$ J+ J1 P+ j- D2 P. E
                           X( 7, 6)       0.5387230            0.0000009 r9 L/ }' Y6 G
                           X( 8, 1)        1.000000            0.000000& Q/ f+ X$ P- `0 }' Q
                           X( 8, 2)        1.000000            0.0000000 @, K2 a7 o; g3 r4 ?1 Z
                           X( 8, 3)        1.000000            0.0000000 E% b$ ]+ b9 m, V: e# n
                           X( 8, 4)        0.000000            0.000000
    & Y2 L6 q7 l( V- b. w                       X( 8, 5)        0.000000            0.000000+ u4 p: z' q% o8 W; t0 I1 g
                           X( 8, 6)       0.4305890            0.000000
    - ]7 e5 n4 I5 Y7 h- g                       X( 9, 1)       0.7096300E-01        0.0000005 C8 S- n0 c* v& c7 O  v& e
                           X( 9, 2)       0.2775390            0.000000( K; D' m! |! G' {  ^9 ~' R* I  E
                           X( 9, 3)       0.3859470            0.0000006 _; p& f' c. |9 Z" a
                           X( 9, 4)       0.1535660            0.000000- l* _# |8 e- b
                           X( 9, 5)       0.1475700E-01        0.000000+ U: i9 q$ d/ y) G2 D% W1 g/ i
                           X( 9, 6)        0.000000            0.0000006 z8 m" }, D) _2 E4 Q' _3 w
                          X( 10, 1)       0.3816400E-01        0.000000# B0 }% h+ h9 e+ m# O' R) M
                          X( 10, 2)       0.3837300E-01        0.000000( m8 I- `/ D$ B2 e6 @
                          X( 10, 3)       0.2499000E-01        0.0000003 S5 Q4 v( ]$ ^9 u1 q% o9 n
                          X( 10, 4)       0.4898490            0.000000
    ' {4 z5 s* ?# x                      X( 10, 5)       0.3115900E-01        0.0000002 `/ z6 `: h: m
                          X( 10, 6)        1.000000            0.0000007 ]: z6 R8 d2 g9 `! b
                          X( 11, 1)        0.000000            0.000000+ }0 q0 G' p4 X7 L! B3 [' p4 Q' `
                          X( 11, 2)       0.8735000E-02        0.000000  q2 X/ A' D& k" _9 g
                          X( 11, 3)       0.1525900E-01        0.0000002 O" v: g; u% E5 X$ [
                          X( 11, 4)       0.2889120            0.000000/ U4 v' {+ N) g0 J  D2 _9 F3 B& l
                          X( 11, 5)       0.3700900E-01        0.000000
    9 c7 T2 x2 [3 A8 P. S2 H0 b                      X( 11, 6)       0.4091800E-01        0.000000
    . J( t. ~9 h3 m                      X( 12, 1)       0.5718100E-01        0.000000' e, W" r" q# j: J+ M
                          X( 12, 2)       0.3850800E-01        0.000000
    ( V& h- v. b7 c' n( ]$ G8 k' Y                      X( 12, 3)       0.2541090            0.000000( y  j; R5 r, v3 P
                          X( 12, 4)       0.1842790            0.000000- y0 w4 q2 e* t4 t9 i4 m) V
                          X( 12, 5)       0.9316000E-01        0.000000
      E9 c5 Y! g9 Y3 F7 J, V" K& O                      X( 12, 6)       0.5224600E-01        0.000000
    * }* p! O6 g- Z! Q5 N                      X( 13, 1)       0.4821100E-01        0.000000
    - `- U- |6 Y2 s& Y+ k$ C                      X( 13, 2)       0.9641000E-02        0.0000007 Y3 C; [; t8 N% Y6 k
                          X( 13, 3)        0.000000            0.000000. v" j8 A( p8 N. C/ r% ]+ O; v
                          X( 13, 4)       0.3878190            0.000000: s" W$ Q+ j6 }1 B" f7 |
                          X( 13, 5)       0.1358340            0.000000: ]/ h. T# R+ O# V, ~5 A
                          X( 13, 6)       0.1185130            0.000000
    " r2 T) P- R  H4 |: Y+ B' d2 I
    9 O% }7 _! c: O5 }/ x& L                            Row    Slack or Surplus      Dual Price
    ; h4 T# V- |$ T                              1       0.7272518            1.000000
    * D* {( `7 h5 d                              2        0.000000           0.3636259! t& I8 x6 E& P! y# [
                                  3        1.124802            0.000000% U. y5 r" g" C. [( v3 c
                                  4        1.121677            0.000000# [2 \3 u( G$ ~9 G
                                  5        1.162172            0.000000
    . H, b9 z+ ?9 h3 N; P' o                              6        1.631867            0.000000
    ; N( k& |6 `# N4 `. R                              7        1.386806            0.0000006 B- w4 q3 S3 x; r2 ~7 ~
                                  8        1.628048            0.000000
    % [- I* S( s8 M3 l; ?7 Y, b! p                              9       0.8751980            0.000000
    # [# e1 \1 U& k- v/ L2 Z$ l                             10       0.8783228            0.000000
    * z% L3 t/ i' Q& y! o# d* }+ X; u4 K                             11       0.8378279            0.000000$ J, e+ w% x9 T# ^$ ~
                                 12       0.3681325            0.000000
    3 n: U; p& ~- X" ~* m% d                             13       0.6131944            0.0000007 U& N# R+ E$ c" z- f4 J
                                 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
    但是 我设置全局求解时：两三秒就停止了( f8 K) K) L8 o, o0 t5 B
    出现错误提示
    2 \/ ]3 u& S( h$ |- _- ~/ ]& U* q还有 上面的 全局的做大值还没有局部做大的大  为什么！！！
    回复

    使用道具 举报

    返回列表 发新帖
    高级模式
    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, 2025-12-4 14:58 , Processed in 0.626738 second(s), 103 queries .

    回顶部