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

huiwen

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

1. function [x,y]=IntLp(f,G,h,Geq,heq,lb,ub,x,id,options)
   3 C- M5 ?) \0 B7 x/ g
2. global upper opt c x0 A b Aeq beq ID options;
   7 b% g3 [1 r$ j) b0 ]0 [0 ]; p
3. if nargin<10,3 `3 X$ m0 T% L  ]8 l& U; ^
   
4.     options =optimset({});* K( F: j* ~5 e1 W9 S
   
5.     options.Display='off';% E% W' ]) C. K4 y
   
6.     options.LargeScale='off';; ^8 w% e0 {. s: y\" E
   
7. end5 H( e9 n( ^# K' o; H& `9 @# \, q
   
8. if nargin<9,
   8 P4 |: a# t7 U% i
9.     id=ones(size(f));
   3 S9 g. q! ^( ~' M! Z
10. end
    ' w8 |# `  Q( p# K+ O! \# M$ [
11. if nargin<8,
    , @9 B8 Z6 y6 v/ ?- w' R; |
12.     x=[];0 p& \8 o9 X; y
    
13. end# r+ p$ Q* ?+ z/ o8 O6 S
    
14. if nargin<7 | isempty(ub),. p; {\" y. z0 Y9 c$ h
    
15.     ub =inf*ones(size(f));
    & P! L# E5 z: m* ^. z
16. end$ j5 i6 t- s$ \4 i
    
17. if nargin<6 | isempty(lb),
    ( s* m4 ]5 ]- G
18.     lb=zeros(size(f));
    4 z( n6 k  l% W
19. end
    ' t: A0 k3 G) l, e
20. if nargin<5,
      Y: p# k6 p\" x% D, m% T
21.     heq=[];! J. I* q. \4 ]- d( ^6 l+ Q
    
22. end
    % K2 x0 U- j, ?4 `7 \  |  g
23. if nargin<4,
    3 J' B9 t4 x% M# M* s+ C8 p
24.     Geq=[];% \\" R, S8 P9 ?1 c* ^6 N: H
    
25. end7 }' T2 _' @5 `( z- O2 I
    
26. upper=inf;c=f;x0=x;# Z1 ~* o& i  s! e. `; y
    
27. A=G;8 ^7 _$ o3 N9 {8 z( g! s\" H% X
    
28. b=h;Aeq=Geq;
      c7 v; s2 m% {# B3 `+ q) _3 h* E
29. beq=heq;ID=id;+ g( V- m; p. n: S+ j: S, ?
    
30. ftemp=IntL_P(lb(:),ub(:));2 u& `  L\" _# K' K0 e
    
31. %以下为子函数
    4 j3 N6 D4 q- d5 y
32. function ftemp=IntL_P(vlb,vub)9 a+ }) k: D0 S3 f7 k2 ~  J3 U
    
33. global upper opt c x0 A b Aeq beq ID options;
    + F2 J\" o7 ^  X9 g* d
34. [x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);
    3 t( X6 G% M; d; Y
35. if how<0
    \" H0 H5 e8 T/ M9 R1 b
36.     return;
    4 U+ f' n6 X' G- A# s9 V
37. end2 ^- P* ~: `4 I/ V% x
    
38. if ftemp-upper>0.000051 Q9 j3 g+ f8 A; |\" J) \
    
39.     return;$ j4 T1 L% @% _( J! K
    
40. end\" l+ Q$ `6 N7 t5 N3 l
    
41. if max(abs(x.*ID-round(x.*ID)))<0.00005: H, |: l8 s# p4 V# C
    
42.     if upper-ftemp>0.00005* k* W: ]* y1 Z, F: I% o) n& f, p
    
43.         opt=x';
    & a# M0 [: x* f! \  g- i8 I9 L5 A
44.         upper=ftemp;# K* `+ X! `8 x( p# B\" w6 v9 A
    
45.         return;0 M4 C\" a8 h: H- Y3 U! Z
    
46.     else
    / X% h' [: c\" F
47.         opt=[opt;x'];# G/ C$ O3 f. y2 L) L\" n7 [
    
48.         return;\" h\" p8 c+ j$ H
    
49.     end+ J/ F7 p6 c$ H8 K9 C2 f' P8 w
    
50. end4 u3 P4 o6 M# S! H: a2 O/ }
    
51. notintx=find(abs(x-round(x))<=0.00005);
    $ x  b1 K5 f3 E9 a! q1 ]# Y! \2 C
52. intx=fix(x);
    ! d- u- Y7 {: f) u; c$ ?/ \
53. tempvlb=vlb;tempvub=vub;. |' Z1 _9 f9 E0 `% S
    
54. if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1- q3 U4 {5 f/ W
    
55.     tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;  n- ?4 L2 b1 G/ d) a: ?
    
56.     ftemp=IntLP(tempvlb,vub);! ?, K7 S* h: X7 z2 E* C
    
57. end
    6 R$ ~# F4 m+ ]1 G+ ]
58. if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)( @, w9 `1 w9 b4 q) j3 T+ ^
    
59.     tempvub(notintx(1,1),1)=intx(notintx(1,1),1);6 t\" {3 o) Y4 D+ Z8 U9 e
    
60.     ftemp=IntL_P(vlb,tempvub);% h8 w, Z, F) u6 Y8 N5 o
    
61. end

复制代码

喜欢就好

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

蜗牛天使

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

huiwen

蜗牛天使 发表于 2011-8-25 17:57
0 Y# |! C2 ?  U" T% t太长了 也不是很懂 没办法了

恩 也谢谢了啊！~

huiwen

喜欢就好 发表于 2011-8-25 17:50
" z) a( i9 |5 L( S" s, y8 G来看看   楼主发的程序有点小问题   无法复制  望楼主解决

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

xueyues3

路过!!!!!!
% U: ?/ ?" {2 m: c3 a1 j


2 e$ O% z: M' X' r/ d9 [, L

+ e1 v* T$ A/ A. B5 l% |  V
, i- K, p1 }8 v
) U& u  T+ ]5 S
% U9 o' u# `% W2 Z6 M$ j


2 p! Y+ I  K; j
高密招聘网 gmzhaopin.com

大鲵2003

大鲵2003

大鲵2003

冰淇淋小芯