数学建模社区-数学中国
标题:
lingo编程求不出全局最优解，请大侠帮帮忙！！！
[打印本页]
作者:
jiiing
时间:
2009-9-5 07:49
标题:
lingo编程求不出全局最优解，请大侠帮帮忙！！！
model:
2 C" y) y) ]" L9 Z. u6 w
sets:
+ }4 D' `; c& `: o8 m6 {5 S. p+ E
shengshi/1..13/;
1 {# z8 l) s$ H5 t- T! Z& [3 d
zhibiao/1..6/:a;
$ i" l! H: A3 A0 \- B4 h! n" p- G u
links(shengshi,zhibiao):x;
* a) p! A: V- k& b
endsets
, d; V U9 `8 _/ W, v' B% W. \
data:
2 B. l1 C b) p& @1 ?& K
x=
+ Y2 X* G K/ B- I; B
0.099599 0.104095 0.021186 0.489849 0.031419 0.399601
. C: ]. c3 `3 E: C) [8 `1 @
0.012222 0.006320 0.035075 0.558043 0.253776 0.465170
1 n* N6 X" I" [1 h
0.016260 0.004213 0.084423 0.631442 0.375174 0.535579
1 J$ ]: w) J+ Y+ G* l
0.004587 0.000000 0.086154 0.557522 0.210591 0.464721
: s) N; P( _ w! }4 l- n) ?
0.039880 0.015333 0.077921 1.000000 0.488760 0.888523
; e& r. y+ k+ ?# i8 m- n' s/ @& {
0.003290 0.002399 0.023441 0.299844 1.000000 0.217814
5 A; X7 N/ e4 l3 K& p
0.021750 0.023296 0.055727 0.635086 0.657561 0.538723
! _. U' L p; N5 R6 a
1.000000 1.000000 1.000000 0.000000 0.000000 0.430589
+ [* g6 P2 {" x9 [
0.070963 0.277539 0.385947 0.153566 0.014757 0.000000
& v& i2 B# @/ G* ~) d( {
0.038164 0.038373 0.024990 0.489849 0.031159 1.000000
( E2 W& s% T& O+ O; g
0.000000 0.008735 0.015259 0.288912 0.037009 0.040918
1 H* m# v; |1 [/ {
0.057181 0.038508 0.254109 0.184279 0.093160 0.052246
& k. z% g* B9 B! ]8 a
0.048211 0.009641 0.000000 0.387819 0.135834 0.118513;
$ l2 P- Z$ e" ?* v) V, Z
enddata
1 C6 e( F6 J B* S* |: `0 m& z% L
max=
p! L, Y+ G( u+ H2 u
(@sum(shengshi(i):
6 b4 K1 u5 ~! J3 L8 D7 r" A2 O: E
(@sum(links(i,j):a(j)*x(i,j))-@sum(shengshi(i):
3 h+ A, \) ~9 } x% _6 i- h0 Q+ v
@sum(links(i,j):a(j)*x(i,j)))/13)^2)/12)^(1/2);
9 q# d9 q$ [0 K: G+ w+ ^
, b+ q% r- d4 t, C6 u' G) t
@sum(zhibiao(j):a(j)^2)=1;
/ k" W7 Q# \5 F$ j* i
@for(zhibiao(j):a(j)>-1);
) K! G" e7 S7 Z0 A2 w* a2 h
@for(zhibiao(j):a(j)<1);
; E6 Y4 B0 E8 T8 M, [# j6 X2 J
@for(zhibiao(j)
free(a(j)));
, f6 w7 ] B1 n) l* b9 w
end
$ [7 t- S- H5 u/ Q* H6 u
求不出全局最优解，请大侠帮帮忙！！！
作者:
kongfanqing
时间:
2009-9-5 08:32
运算时间太长，你能否增加一个初始值就好了
作者:
kongfanqing
时间:
2009-9-5 08:47
Global optimal solution found.
2 g: j. ~. G& N. }
Objective value: 0.7272518
1 P! Q5 I2 X2 z @$ p" E. t8 ~
Objective bound: 0.7272527
- F1 ^: }9 |, g
Infeasibilities: 0.000000
& c2 _3 l: K) f, F
Extended solver steps: 38843
5 Y; q+ _3 `& F* a1 _/ i% O6 N
Total solver iterations: 7561347
# j0 }/ _% K$ H9 y$ n
8 y; F* d. L/ S7 z7 T) t
! n6 N; p6 L! H" f& e, V
Variable Value Reduced Cost
! D, y3 i& a! x
A( 1) 0.1248020 0.000000
1 S+ Q1 _: A/ R% F
A( 2) 0.1216772 0.000000
+ e/ ]. {. G9 t3 P1 ~ t- s
A( 3) 0.1621721 0.000000
: H1 t8 M2 l3 L5 a/ a6 E1 A
A( 4) 0.6318675 0.000000
8 t5 h4 d+ I9 r
A( 5) 0.3868056 0.000000
0 Z7 i" j) \! x
A( 6) 0.6280479 0.000000
/ v( n) N ~9 d4 Z8 [
X( 1, 1) 0.9959900E-01 0.000000
( j, Q& b3 H. ]$ j4 f
X( 1, 2) 0.1040950 0.000000
- @/ a7 j* B0 X4 q
X( 1, 3) 0.2118600E-01 0.000000
/ \. U J% q: a8 u) H2 B. h6 b
X( 1, 4) 0.4898490 0.000000
2 x, i. U/ s& V
X( 1, 5) 0.3141900E-01 0.000000
5 y. Q" `' N1 a" L
X( 1, 6) 0.3996010 0.000000
9 S6 Y; v) A# S) Z) W
X( 2, 1) 0.1222200E-01 0.000000
: K1 k/ ~& H: m0 H" {* N
X( 2, 2) 0.6320000E-02 0.000000
' {& Z" T D! ^- \/ u8 a
X( 2, 3) 0.3507500E-01 0.000000
z4 A, I" [" m. g9 U c
X( 2, 4) 0.5580430 0.000000
* X: l' W+ Y; J7 B- O
X( 2, 5) 0.2537760 0.000000
8 r2 V# w2 U! |% B0 ]# N, n8 t
X( 2, 6) 0.4651700 0.000000
3 |, l2 \! }" e! b
X( 3, 1) 0.1626000E-01 0.000000
/ l4 o/ ^, V. R# J1 j( W# v8 R+ p
X( 3, 2) 0.4213000E-02 0.000000
; X8 p' a" T; Y2 q0 d
X( 3, 3) 0.8442300E-01 0.000000
* S( w3 O* z" _/ S$ Z7 A: m
X( 3, 4) 0.6314420 0.000000
D4 |9 ^* v5 F% j' B+ E, z
X( 3, 5) 0.3751740 0.000000
/ d' C. |* x5 j% o% A
X( 3, 6) 0.5355790 0.000000
9 z7 Y: }9 R5 \+ A4 l- j& z
X( 4, 1) 0.4587000E-02 0.000000
$ {* l4 z9 e# s4 k
X( 4, 2) 0.000000 0.000000
: J! d: d4 D$ p/ ]0 t( ?! h1 C& J- h% R
X( 4, 3) 0.8615400E-01 0.000000
& [0 G. A' u5 W" ?# n5 l
X( 4, 4) 0.5575220 0.000000
4 o/ O" q& ~4 N4 s' E
X( 4, 5) 0.2105910 0.000000
- X. N9 B9 p2 g) I- f! g
X( 4, 6) 0.4647210 0.000000
9 d% c1 _' `0 H( [- G# y1 J
X( 5, 1) 0.3988000E-01 0.000000
* D8 w" ]7 c% [$ @
X( 5, 2) 0.1533300E-01 0.000000
: j7 s4 [: E# v6 ^" m
X( 5, 3) 0.7792100E-01 0.000000
8 Z- ^7 Y% b' ]3 Q8 p7 z
X( 5, 4) 1.000000 0.000000
; V& J+ s6 N/ A$ C) Y
X( 5, 5) 0.4887600 0.000000
: `0 Q, N' K( @
X( 5, 6) 0.8885230 0.000000
) B1 k! f1 b7 p% d- L
X( 6, 1) 0.3290000E-02 0.000000
: K: b5 X) ?! ?9 J
X( 6, 2) 0.2399000E-02 0.000000
. ^9 h4 ?6 d" ^9 j4 p
X( 6, 3) 0.2344100E-01 0.000000
* u0 z5 [( l/ i# P8 [- {7 h( z, T' M4 k
X( 6, 4) 0.2998440 0.000000
1 r8 C! u) |% A3 g# O6 P
X( 6, 5) 1.000000 0.000000
' f% K6 G* t ~! k7 r
X( 6, 6) 0.2178140 0.000000
! Q- t$ n9 [' o X- q$ p* z
X( 7, 1) 0.2175000E-01 0.000000
% `- a) U* ]4 ]; X2 o$ h: ]
X( 7, 2) 0.2329600E-01 0.000000
4 r L6 H$ x& `' ]) A r0 W
X( 7, 3) 0.5572700E-01 0.000000
" J f6 G7 w! q/ [$ Q; w
X( 7, 4) 0.6350860 0.000000
. D( [, o. U0 s# D
X( 7, 5) 0.6575610 0.000000
/ X0 M5 [% V- [/ f, M, g" P
X( 7, 6) 0.5387230 0.000000
8 S7 Z( F9 [2 n! W. d7 ]
X( 8, 1) 1.000000 0.000000
: r* N% o! c! U& A
X( 8, 2) 1.000000 0.000000
& c- B2 M. ]5 a
X( 8, 3) 1.000000 0.000000
+ S+ \& _. ` v0 Y' f/ G( _
X( 8, 4) 0.000000 0.000000
& Q) \6 h' l+ x. n0 Z# Y- r7 s
X( 8, 5) 0.000000 0.000000
1 ~# L$ f c1 G; q+ ?
X( 8, 6) 0.4305890 0.000000
5 ]6 o1 m8 k( Y; ^/ n4 X1 o
X( 9, 1) 0.7096300E-01 0.000000
! U5 ]+ L$ V. `& r
X( 9, 2) 0.2775390 0.000000
# h+ F/ |2 ~: `9 T0 C- N3 X
X( 9, 3) 0.3859470 0.000000
* N6 W1 _9 t1 C+ y" E4 i( J. ~! N
X( 9, 4) 0.1535660 0.000000
1 E5 g& ?* l7 U8 V% W1 G& y+ |- T# V+ }) ~
X( 9, 5) 0.1475700E-01 0.000000
8 R% w; \* I8 ?# x! ~( D
X( 9, 6) 0.000000 0.000000
; E# p3 w7 O. m% W! m. r# [
X( 10, 1) 0.3816400E-01 0.000000
# I5 z: P6 ^; D6 d2 L5 _4 y
X( 10, 2) 0.3837300E-01 0.000000
& w. y- t* U6 d. d! d% Q0 u
X( 10, 3) 0.2499000E-01 0.000000
" E* X: C& l# c) `( E
X( 10, 4) 0.4898490 0.000000
' {0 P1 T, c: `# k# w
X( 10, 5) 0.3115900E-01 0.000000
: ^* m1 m# C# }6 r0 b2 E
X( 10, 6) 1.000000 0.000000
. {8 C5 f, g+ K& P8 i1 Y* X
X( 11, 1) 0.000000 0.000000
+ Y5 y2 \2 m0 Q
X( 11, 2) 0.8735000E-02 0.000000
# G! L5 Z' C9 r! Y4 l
X( 11, 3) 0.1525900E-01 0.000000
* j0 c6 L- I' v6 K0 c6 r, {8 U3 _
X( 11, 4) 0.2889120 0.000000
+ g, l, t, k. A* h0 c5 q" M; T5 v
X( 11, 5) 0.3700900E-01 0.000000
% F0 b$ u/ f- D/ M3 k
X( 11, 6) 0.4091800E-01 0.000000
: \4 c% y4 C+ v' ~% t2 }, w
X( 12, 1) 0.5718100E-01 0.000000
1 r- l, c: k4 L; W4 c
X( 12, 2) 0.3850800E-01 0.000000
' m6 j/ I" F$ `- h0 ~4 _$ @
X( 12, 3) 0.2541090 0.000000
$ v8 H) k- i' j6 G4 B: C& h
X( 12, 4) 0.1842790 0.000000
4 p `, a* ?0 }; V. X
X( 12, 5) 0.9316000E-01 0.000000
p( H2 s. Z" n- [7 o" z
X( 12, 6) 0.5224600E-01 0.000000
- r, ]+ a) t( P7 t- {6 ~* X
X( 13, 1) 0.4821100E-01 0.000000
9 C% p- H. M% N
X( 13, 2) 0.9641000E-02 0.000000
0 ?! _' x& |% X8 [+ R$ s0 t+ V
X( 13, 3) 0.000000 0.000000
/ X" F' e- }+ L7 Z
X( 13, 4) 0.3878190 0.000000
, k* C; b( P3 S
X( 13, 5) 0.1358340 0.000000
. P, E6 u/ o" `0 r
X( 13, 6) 0.1185130 0.000000
' c4 `- I6 k$ h; P# B# ^9 b, `
5 N% L* d+ Q f* C
Row Slack or Surplus Dual Price
9 V& h* h# [3 r/ t# V1 ?
1 0.7272518 1.000000
$ m. b2 i) \: M0 V2 [
2 0.000000 0.3636259
3 B$ }4 j/ \) Q) d0 w2 O9 C
3 1.124802 0.000000
' n* ]% k% {# P f/ {7 J" x
4 1.121677 0.000000
: x* j- I7 W0 t( u
5 1.162172 0.000000
9 x9 h J- U' Q9 x& b3 T$ S* I
6 1.631867 0.000000
. r9 o% h, B* e; k
7 1.386806 0.000000
+ ]* G0 `9 l+ N* j
8 1.628048 0.000000
+ ]- j& p3 ^( W) b8 f' Y9 {
9 0.8751980 0.000000
# t+ D! X0 c. |& E
10 0.8783228 0.000000
3 l6 n: Q6 s6 q& O
11 0.8378279 0.000000
. z6 @; k3 V9 Y. ^" J1 _
12 0.3681325 0.000000
+ ^! G' w% V+ t, q0 i/ j' h
13 0.6131944 0.000000
9 S, ~1 C# v {" C$ I9 L3 H$ q3 C
14 0.3719521 0.000000
作者:
kongfanqing
时间:
2009-9-5 08:48
运算了12分钟才计算出来，给他一个初始值就好了！
作者:
ychexingwuji
时间:
2009-9-5 08:49
学习了 呵呵
作者:
formulaone9275
时间:
2009-9-5 08:54
有很多时候都是求不出最优解的，局部最优就可以用了 ，我看了你的程序，其中max里的sum有对i求了两次和，你看看能不能再改改，求求试试
作者:
yuan07150215
时间:
2009-9-5 08:58
先在限制条件下得出最低值，之后运行lingo程序，当运行到20-30分钟时，人为的终止求解（按下Interrupt Solver），此时运行的解与目标函数的最优解很相近！可认为是最优解！……
作者:
xuyin
时间:
2009-9-5 09:06
哇，都是高手啊
作者:
liulongfei
时间:
2009-9-5 09:35
选择全局最优解，然后就是时间的问题了
作者:
jiiing
时间:
2009-9-5 12:16
但是 我设置全局求解时：两三秒就停止了
7 R% q$ Y* k( }0 [2 n. r
出现错误提示
! s9 Y' F1 Z6 B! v+ A2 B
还有 上面的 全局的做大值还没有局部做大的大 为什么！！！
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
free(a(j)));, f6 w7 ] B1 n) l* b9 w