请求一个整数规划的matlab程序

huiwen

最近一直在研究整数规划，可是连着看了两个相关编程，用MATLAB实践时总是出错！所以，求助于各位大侠们能否提供一个整数规划的程序。或者，能否分析一下面的程序哪里有问题？谢谢了啊!

1. function [x,y]=IntLp(f,G,h,Geq,heq,lb,ub,x,id,options)+ l* d/ m9 Z+ M5 o9 ~
   
2. global upper opt c x0 A b Aeq beq ID options;1 e! s! o( ?0 K5 ^/ ^
   
3. if nargin<10,  Z7 O2 R\" e. z
   
4.     options =optimset({});# V) x1 h+ R4 L  R0 B
   
5.     options.Display='off';/ F9 D) M5 c8 H1 r( V% W
   
6.     options.LargeScale='off';  ^+ K$ D! K* c& _1 n
   
7. end* e2 z) h) J' ~' O
   
8. if nargin<9,
   $ Y* N. B' J' R( |7 _# m, j3 f
9.     id=ones(size(f));
     J, e& o& R, @  r# \* ^
10. end
    8 v; j- [5 I  k- |' y
11. if nargin<8,
    & t( N6 ]! K$ z
12.     x=[];
      l: g* @, Z2 J2 W\" A7 O
13. end( N* w\" f8 D) p5 ]( X
    
14. if nargin<7 | isempty(ub),/ l* a& |3 n( m+ ]5 m  u: A) C
    
15.     ub =inf*ones(size(f));
    ; ]! N5 l% v  x/ g7 n3 h7 M2 L
16. end' O: u  d& p, B! L( Q
    
17. if nargin<6 | isempty(lb),* ?6 e1 Q! P: O& _: w
    
18.     lb=zeros(size(f));
    . K& i\" \7 ~3 [/ k0 H
19. end- z5 h3 _; V3 h9 P* L+ j& p' M
    
20. if nargin<5,0 f5 G- p2 C; A% z
    
21.     heq=[];
    6 j6 o' O/ g! \4 d* `, x2 R, L# S
22. end, L) p7 c0 f* L5 t( K  P0 ^% J8 j
    
23. if nargin<4,7 i  j8 [! c( U  B' }/ U& g
    
24.     Geq=[];/ W: a4 C: Z5 u5 C6 C; [7 L
    
25. end) a3 c  ^. h. E: f6 g% t' d4 P
    
26. upper=inf;c=f;x0=x;* H7 i. I8 I' R) P4 g
    
27. A=G;
    6 i* P5 |# x& u. W7 Q  c
28. b=h;Aeq=Geq;3 b1 D9 X8 U$ m( E1 z2 }
    
29. beq=heq;ID=id;; I/ [. L0 p6 ?\" l. f  u
    
30. ftemp=IntL_P(lb(:),ub(:));0 i- k. k1 E( U1 ~( a
    
31. %以下为子函数
    1 \: O) c2 a4 i8 X
32. function ftemp=IntL_P(vlb,vub)
    5 j1 M$ x& a$ M+ ^8 i$ v
33. global upper opt c x0 A b Aeq beq ID options;
    7 [  Q7 C: Z- f' ]
34. [x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);
    0 [. q- c: R/ i7 ?' B
35. if how<0& f7 Z/ l; t* k
    
36.     return;
    & [8 b5 u3 f5 F\" P/ ~7 c+ _) q
37. end7 R, {7 ]2 T7 }4 E
    
38. if ftemp-upper>0.00005
    \" S\" F) [$ e. B+ }9 n+ o! h! f
39.     return;7 O2 q. T6 Q5 _  V( \( j7 G' D! }; d
    
40. end2 u1 D, ?- D: W, m+ y( m: C
    
41. if max(abs(x.*ID-round(x.*ID)))<0.000054 H9 I* s  `\" K  Y: T
    
42.     if upper-ftemp>0.00005
    9 Q, I+ j( P$ p+ \+ C7 N9 Z
43.         opt=x';
    3 S6 T# n, l* }) t  [
44.         upper=ftemp;
    3 L% }( m) k6 M5 g2 H9 P
45.         return;5 `. S\" s+ j) x4 M! Q
    
46.     else! q: \0 M: a0 _
    
47.         opt=[opt;x'];! R\" m- B' L) b4 r
    
48.         return;+ e. i- D1 j5 \- \# q
    
49.     end5 _9 G$ N! D; ~' E1 F' d7 i) D8 m
    
50. end2 s+ G2 ]9 \& {7 B/ J
    
51. notintx=find(abs(x-round(x))<=0.00005);
    ' w* f. E- w0 B6 E$ f1 n
52. intx=fix(x);\" Z& D5 Y2 e6 N: x$ U: P6 N) l
    
53. tempvlb=vlb;tempvub=vub;
    ' E* _' z8 Y& l8 e& C/ c
54. if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1  j. s\" L, v8 x
    
55.     tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;# {8 `9 ?; \& x9 O( D
    
56.     ftemp=IntLP(tempvlb,vub);; z7 Q! c& w& g, k
    
57. end
      S0 L/ [) f  x1 x0 d# E6 n4 l  Z5 s
58. if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)
    / a# X# J) R& t; n* f! I
59.     tempvub(notintx(1,1),1)=intx(notintx(1,1),1);
    : G+ V. A\" W8 O- ~2 A: S) z3 `
60.     ftemp=IntL_P(vlb,tempvub);
    , ], X+ P6 N+ i
61. end

复制代码

喜欢就好

来看看   楼主发的程序有点小问题   无法复制  望楼主解决

蜗牛天使

太长了 也不是很懂 没办法了

huiwen

蜗牛天使 发表于 2011-8-25 17:57
/ C$ ~! I& n9 h9 p6 s) M8 ]太长了 也不是很懂 没办法了

: n; b1 T; q) p, I* j: _, h恩 也谢谢了啊！~

huiwen

喜欢就好 发表于 2011-8-25 17:50
6 d$ m9 H: C  G! p; T来看看   楼主发的程序有点小问题   无法复制  望楼主解决

我这里看的时候还是可以的啊！

xueyues3

路过!!!!!!
6 ~! b) R/ A: ~

5 |: k7 _9 n- {% K5 D( d  n# e5 {








7 O& R) }3 {* z! l, H
高密招聘网 gmzhaopin.com

大鲵2003

大鲵2003

大鲵2003

冰淇淋小芯