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

huiwen

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

1. function [x,y]=IntLp(f,G,h,Geq,heq,lb,ub,x,id,options)) x% x  v& u' G2 Z
   
2. global upper opt c x0 A b Aeq beq ID options;
   % a- j# T+ S  |0 h+ i% k( o6 L
3. if nargin<10,$ |; K# q* ~% K  ~& r& A$ ?1 y
   
4.     options =optimset({});5 z9 s5 k6 x( v
   
5.     options.Display='off';' ~- `8 ~5 M- W  Q/ u
   
6.     options.LargeScale='off';
   - l; j+ T# ]( n+ o& m
7. end
   . z\" G\" B9 P0 Q& l* t\" A( M
8. if nargin<9,( U5 l% I' k/ Y/ l
   
9.     id=ones(size(f));  A& b3 K9 c; `6 ^& s6 l
   
10. end 1 O$ K4 S2 Z+ Z4 r3 m1 _
    
11. if nargin<8,8 N\" Q3 W$ \; w! c- e8 V( D* X
    
12.     x=[];
    . Q* R2 B, ]3 s+ p\" N: q
13. end
    8 u! L8 r' Y! G# l7 q
14. if nargin<7 | isempty(ub),
    - R4 o* U3 i& R: S8 |
15.     ub =inf*ones(size(f));
    & S& n2 y' k0 Y2 B. p\" \) U7 e
16. end
    3 t; W, |5 H( |2 @. G
17. if nargin<6 | isempty(lb),
    $ C/ L6 M  [5 b; P, v% i. D
18.     lb=zeros(size(f));$ `2 i) ]0 {7 M- L) K
    
19. end/ j/ m# R9 L5 n
    
20. if nargin<5,4 Q' s& |6 J* }) `
    
21.     heq=[];
    & T1 x1 i) q6 A# h2 H% S2 ?
22. end
    + a& Z8 A  I1 ^5 l
23. if nargin<4,
    7 U5 u+ N! V2 e: J
24.     Geq=[];; w9 B; o/ h* F- t2 T
    
25. end
      j1 ]% M+ E2 T' J* R  t
26. upper=inf;c=f;x0=x;
    \" h; J/ S% j3 |' l
27. A=G;2 Y! R$ k2 @# p: e3 L\" }. M  M3 \
    
28. b=h;Aeq=Geq;
    : I# @9 l+ p\" \: V4 \. t, l
29. beq=heq;ID=id;
    ( C* }# T- A- c- ]4 g% ~8 l6 F* S, W
30. ftemp=IntL_P(lb(:),ub(:));
    3 t7 z. o+ f  P, z
31. %以下为子函数1 p9 Y- ^, M+ D+ s4 ]9 K8 r
    
32. function ftemp=IntL_P(vlb,vub)
    5 O( p1 m3 V0 @6 S
33. global upper opt c x0 A b Aeq beq ID options;
    ; L3 T% B4 p: l0 Z2 H8 G# I
34. [x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);1 Q3 J) ^2 `. H0 ]# x
    
35. if how<0
    & Y- A3 g\" T* c; \7 N) i1 E# P
36.     return;, r  b; i( F. f4 m3 b
    
37. end
    9 \+ L- ?, @5 E( |, F$ ^! N
38. if ftemp-upper>0.00005$ B+ b5 y' M' r9 f
    
39.     return;
    6 n/ i3 T4 x! I: T% N: `/ j
40. end
    3 e9 L( b' H/ U# v9 h
41. if max(abs(x.*ID-round(x.*ID)))<0.00005
    : B; x; o+ A6 Z3 E! R
42.     if upper-ftemp>0.00005
    ; U/ S) k& y- c! r' K. O8 L
43.         opt=x';7 q3 g4 `, b2 L7 @& \; m6 w4 Z
    
44.         upper=ftemp;
    5 t1 V7 J( b. G. ~3 E/ W\" \3 P
45.         return;
    , b$ W  T- @) o# A4 j
46.     else
    6 @\" A) Z) K\" Y\" m- `
47.         opt=[opt;x'];
    * `9 K- ~2 U- M8 n' ?
48.         return;1 {& j$ O7 b, v$ S* [
    
49.     end2 A, |. t: A6 {7 Z9 X  G
    
50. end
    2 I/ N) h) x5 p# e; X6 W6 W2 Z
51. notintx=find(abs(x-round(x))<=0.00005);$ B* m5 S) y2 [  Y7 `: X5 J
    
52. intx=fix(x);
    7 e0 z' x) T, d7 U
53. tempvlb=vlb;tempvub=vub;
    8 w' u$ w6 w) j1 K$ P$ w' l2 C
54. if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+13 v8 B1 x0 Y, D
    
55.     tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;
    9 G; A9 o3 J& V2 H) G
56.     ftemp=IntLP(tempvlb,vub);0 E( m2 \8 c) u& O. {
    
57. end
    ; ^( R0 D7 U1 [# c1 w3 }$ C/ S# [
58. if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)
    # F2 U; b% z\" R7 U+ o/ U
59.     tempvub(notintx(1,1),1)=intx(notintx(1,1),1);
    + V1 ^0 k+ ?\" w1 o
60.     ftemp=IntL_P(vlb,tempvub);
      P3 N1 g( J* c) J
61. end

复制代码

一束阳光

十W分感谢楼主分
" l1 T% E% D# c! \  f  s

NICE--LLL

我也是遇到同样的问题啊

冰淇淋小芯

                  

大鲵2003

大鲵2003

大鲵2003

xueyues3

路过!!!!!!

# z! U, ~. B( p6 Y; y" f5 o1 L

+ V& ?  `4 E8 s; t- m
% S. W# f; c  J, D


# `1 q. I) c$ W- O' \/ V: ?




: j6 n. y: S# u- T! N# |* O高密招聘网 gmzhaopin.com