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

huiwen

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

1. function [x,y]=IntLp(f,G,h,Geq,heq,lb,ub,x,id,options)
   + n3 V\" C0 r& r; ?+ |
2. global upper opt c x0 A b Aeq beq ID options;
   * T\" a6 j9 Z, l' P4 b* I
3. if nargin<10,
   ; y+ d0 y& J* M3 T
4.     options =optimset({});. X) T+ ?7 y; p1 B1 Y8 H
   
5.     options.Display='off';6 T2 {3 K4 |0 U/ f1 D1 z  i$ b
   
6.     options.LargeScale='off';+ k: n* ]6 w1 S6 h; K
   
7. end  x  E9 k$ y  f# q( T
   
8. if nargin<9,
   & W6 @; R  v' {
9.     id=ones(size(f));: H, j. q2 F  K& v
   
10. end
    , n' |. Q$ B. j/ [0 r
11. if nargin<8,
    5 T( a4 x9 s) u  s2 k
12.     x=[];
    7 ?# w1 E# W1 h4 n
13. end
    2 G# }6 M( }0 q8 S5 }7 O# e
14. if nargin<7 | isempty(ub),, e\" ^/ y8 x$ Q0 Z% v+ B
    
15.     ub =inf*ones(size(f));
    \" N$ }9 ^0 E$ O
16. end) P/ H0 B  s- m9 n\" q
    
17. if nargin<6 | isempty(lb),+ g8 @% D8 l3 n
    
18.     lb=zeros(size(f));
    0 l: A4 W( ~. D( ^( Z* W
19. end; q4 b' N8 V6 s
    
20. if nargin<5,
    1 A& |3 E3 z4 B5 z6 ^
21.     heq=[];0 f/ Q) B( }9 U; k. f# e! E
    
22. end
    + f) ^$ w9 k- g$ I- }; \, \( y' F
23. if nargin<4,6 A$ O6 Z' t0 c% c
    
24.     Geq=[];5 C9 o) U# Z! D/ |\" |
    
25. end* q% G) l. ]. W7 P
    
26. upper=inf;c=f;x0=x;
    9 s\" g1 _# v) D- t
27. A=G;
    ! `1 I+ Q9 J' Q
28. b=h;Aeq=Geq;
    ) N\" ?: ?# \- k+ o- Q
29. beq=heq;ID=id;
    & J+ Q' N: y, s/ r
30. ftemp=IntL_P(lb(:),ub(:));2 R4 @\" n' G7 Q: n# t9 l0 ^0 O2 Q
    
31. %以下为子函数: G  C& A# A# Z
    
32. function ftemp=IntL_P(vlb,vub)7 S( ]; x9 c3 I5 G4 w: y
    
33. global upper opt c x0 A b Aeq beq ID options;0 _3 k5 `5 d4 F% r  f1 c
    
34. [x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);
    / V9 r1 \3 v+ E
35. if how<0
      v. Y9 ^# K  Y% @+ n: l
36.     return;
    $ q2 _& \  e4 E+ _0 V9 T
37. end; Z  H* ^7 R3 t6 w2 S. _
    
38. if ftemp-upper>0.00005
    \" i, y8 j! |4 P( A
39.     return;/ b9 t- D6 K) t8 L& J
    
40. end$ {/ R1 R' n5 G9 ?
    
41. if max(abs(x.*ID-round(x.*ID)))<0.00005
    ; _, L  q; v* O% z1 E\" x
42.     if upper-ftemp>0.000056 t  s. h6 I' _( ^. y/ @; }
    
43.         opt=x';
    \" [% G3 _0 w1 z( r. Z
44.         upper=ftemp;
    6 U, W\" p9 @( E% S* m
45.         return;
    # A  J2 |) \. h: @1 e( _3 T+ k
46.     else
    - c1 M\" _: a& `  e6 M- Q; B( [
47.         opt=[opt;x'];
    / o6 s8 [' ]6 B5 W
48.         return;: S9 j* w6 a5 Q& U
    
49.     end
    8 i2 L6 d8 v- O2 t- `) E7 ]2 h5 B
50. end- t' ]. m- W  H% k4 \% i1 t\" L
    
51. notintx=find(abs(x-round(x))<=0.00005);7 \( D0 l1 H( P8 T
    
52. intx=fix(x);5 o  o0 h$ L5 f5 E1 V: U$ k
    
53. tempvlb=vlb;tempvub=vub;
    7 c4 b1 u* u6 l. k$ k; h$ X4 g
54. if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1
    ) J6 w9 W* K% h) ~% X, _
55.     tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;
    0 S0 ^, v! O+ |  ]0 _0 ^: D
56.     ftemp=IntLP(tempvlb,vub);
    8 p  ]2 p6 D4 n/ o1 [6 c: B5 d; Q8 W
57. end0 k/ W( @2 N$ C- w( e( U' ~' Z6 k
    
58. if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)! ]9 c5 p+ i' I  \- \& B/ C' ?
    
59.     tempvub(notintx(1,1),1)=intx(notintx(1,1),1);
    $ \, X' z* @3 Z  A; ^1 u( k2 v) x
60.     ftemp=IntL_P(vlb,tempvub);
    \" x+ Q! l# r\" V+ q  b) v2 u
61. end

复制代码

喜欢就好

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

蜗牛天使

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

huiwen

蜗牛天使 发表于 2011-8-25 17:57
1 l7 {7 ]7 i& U/ y% Q- m太长了 也不是很懂 没办法了

" Q) q& Q; H6 Q+ v恩 也谢谢了啊！~

huiwen

喜欢就好 发表于 2011-8-25 17:50
/ v$ [* y+ s8 ~/ N) i3 c0 `$ ^来看看   楼主发的程序有点小问题   无法复制  望楼主解决

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

xueyues3

路过!!!!!!
* ?& c$ B7 i* \9 Z0 d, C



; Z; ^: x% @1 A
. }( {5 V  g. f' {  u+ N

( H+ Y- B* A5 |: v* P

) L- W: f: T4 A0 F
" t, w0 G, a- h

高密招聘网 gmzhaopin.com

大鲵2003

大鲵2003

大鲵2003

冰淇淋小芯