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

huiwen

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

1. function [x,y]=IntLp(f,G,h,Geq,heq,lb,ub,x,id,options)  [( c/ `4 n/ L( Z\" j  V
   
2. global upper opt c x0 A b Aeq beq ID options;
   1 o4 _) s) q: y2 m# Z0 F
3. if nargin<10,
     I$ ?3 p) Z2 W\" r( L9 b
4.     options =optimset({});
   , H7 c9 w5 o2 c( o% R2 b
5.     options.Display='off';% E1 n( x! l' k5 m* h& Z; N
   
6.     options.LargeScale='off';
   4 n, \/ P6 f1 E2 ~- ]6 R3 m
7. end
   8 u% P; ^/ R/ n! h\" I! \\" S
8. if nargin<9,* k\" `8 I# ]/ T
   
9.     id=ones(size(f));$ p- L\" I; X9 ?' ?
   
10. end
    - `0 D8 p\" R6 m! @4 Q
11. if nargin<8,9 c; B9 V$ k' R; o3 t) D; c; o
    
12.     x=[];
    ! ]8 B) Q% s4 {- D
13. end
    6 R\" a: K/ e' @\" t0 ]- l( x2 `+ j
14. if nargin<7 | isempty(ub),, B9 h) ^: T0 J- ^2 J
    
15.     ub =inf*ones(size(f));2 L& e& \6 b! Z0 ~' u) g
    
16. end
    % ^' i- z6 E3 U
17. if nargin<6 | isempty(lb),: V$ D- n+ i\" x  C3 Z% ~9 V: J
    
18.     lb=zeros(size(f));
    5 l0 b5 w/ w1 T9 q2 H2 ^5 C/ ]
19. end
    / G. J/ T5 W  e/ ]2 h
20. if nargin<5,; C' Z/ \0 }. C5 f
    
21.     heq=[];
    \" Q+ `5 I+ V7 L
22. end& G$ S/ R& o9 D0 @
    
23. if nargin<4,
    3 r0 R2 n+ l7 k) U# Q
24.     Geq=[];
    3 s7 Z9 N0 c\" H3 E
25. end
    ; F( m* n7 E0 S- H4 V
26. upper=inf;c=f;x0=x;; F% f9 Y7 k. M# ~4 ~8 ^% _# ]+ \
    
27. A=G;! w3 m9 @8 k6 S6 T# O* J
    
28. b=h;Aeq=Geq;
      |) F2 b% v: v' `
29. beq=heq;ID=id;, s/ R& I& M) ?+ z' D, B\" z
    
30. ftemp=IntL_P(lb(:),ub(:));
    8 [; B1 n% o+ P7 ?2 q\" f: V& B
31. %以下为子函数3 G5 s3 l8 o4 b  e8 F  O* B
    
32. function ftemp=IntL_P(vlb,vub)
    * p9 O& f3 s3 z  T# P
33. global upper opt c x0 A b Aeq beq ID options;# h) M! G- \: O1 J* X' Y4 o
    
34. [x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);$ u. W3 R' b$ k
    
35. if how<0* u9 Q6 Q2 |% B( o$ u7 i8 z
    
36.     return;
    \" y( Q8 h5 g( ]$ d/ G3 w7 p
37. end
    , N& w2 B\" }; ]
38. if ftemp-upper>0.00005
    5 \% |& j7 u1 w+ e) f9 A$ }
39.     return;6 w  X; O( _9 U0 c
    
40. end1 R* ~4 e) ^\" U\" k9 F. P( P/ x! S
    
41. if max(abs(x.*ID-round(x.*ID)))<0.000055 i2 o9 ?! x  u. O
    
42.     if upper-ftemp>0.00005) \- b: }* L8 w: @) Z4 U/ O/ U
    
43.         opt=x';
    8 z  h# f: m, \# x% _
44.         upper=ftemp;
    7 m4 Y2 m! y7 G7 o2 `7 r. @/ t, J) g7 C
45.         return;
    ) N% `  q# e7 M! h- K! E! Y
46.     else3 [\" Q; N( ]\" \
    
47.         opt=[opt;x'];
    # m  o. y& ?, f  S1 b4 i
48.         return;
    ! ]; j' a0 N* v, Z8 R6 z
49.     end
    1 M/ V1 b: Y7 y4 K
50. end
    + t9 b. j- O1 m: ^4 c
51. notintx=find(abs(x-round(x))<=0.00005);
    # T2 w( W9 b  `! ^; N
52. intx=fix(x);: V' F1 w- ~- i- z: ?8 f
    
53. tempvlb=vlb;tempvub=vub;5 Y6 y# q/ b  b# X. |& `
    
54. if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1
    ! N0 q: l) t3 @
55.     tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;6 X1 E1 M8 k6 G; f6 O1 A
    
56.     ftemp=IntLP(tempvlb,vub);
    : T/ N& V\" }4 b( X8 g  @
57. end7 C4 c6 t7 u0 ~% {
    
58. if vlb(notintx(1,1),1)<=intx(notintx(1,1),1); p1 p. s/ Y5 x' Q. B
    
59.     tempvub(notintx(1,1),1)=intx(notintx(1,1),1);3 y: m. ^, z0 o% ^1 ]$ K- K
    
60.     ftemp=IntL_P(vlb,tempvub);% l5 N! j% r# a/ p; d1 e, Q6 J) O. l
    
61. end

复制代码

喜欢就好

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

蜗牛天使

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

huiwen

蜗牛天使 发表于 2011-8-25 17:57
' f" x; Q  I* W" W# m$ a太长了 也不是很懂 没办法了

恩 也谢谢了啊！~

huiwen

喜欢就好 发表于 2011-8-25 17:50
来看看   楼主发的程序有点小问题   无法复制  望楼主解决

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

xueyues3

路过!!!!!!


4 @1 i9 H# a2 G4 E" {2 x: X! l

4 s) [) R# v, x9 L# l" d4 P8 h6 D

, m# D& v- S. j7 u5 C" p7 }
5 X5 @' G' H7 M+ ]$ d) X$ G- G6 X3 Z

/ b$ n" S7 N2 u


" ?5 l  T, Q7 g! @, j高密招聘网 gmzhaopin.com

大鲵2003

大鲵2003

大鲵2003

冰淇淋小芯