[问题求助] 请求一个整数规划的matlab程序

[复制链接]

字体大小: 正常放大

huiwen

4 主题	3 听众	258 积分

升级 79%

TA的每日心情

	开心 2014-12-3 16:59

签到天数: 84 天

[LV.6]常住居民II

群组: Matlab讨论组

群组: 西安交大数学建模

群组: 2011年第一期数学建模

电梯直达

1^#

发表于 2011-8-25 16:51 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

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

function [x,y]=IntLp(f,G,h,Geq,heq,lb,ub,x,id,options)+ z! X: D. f5 b- H% c6 [
global upper opt c x0 A b Aeq beq ID options;
1 Q5 E: C7 D0 C6 P d9 j
if nargin<10,6 e$ U9 L\" ?$ Z\" g2 Y2 O& \5 M
options =optimset({});: B. z! J1 b7 {8 f
options.Display='off';. s: Y, G9 p9 y9 n4 S
options.LargeScale='off';- C1 s2 D3 l e4 \! w; U: @* \5 X
end
2 I' x; E5 G1 d; R6 N\" y$ t6 N' F
if nargin<9,
. W\" [) P( [+ A7 f
id=ones(size(f));
6 ?8 `9 o! R1 N) l! t) L; q
end
9 @# r- i\" N# T
if nargin<8,
) V0 J* Q( @0 B\" ?
x=[];
S, n$ `# c, m\" g
end3 v5 ?5 \* K9 i. [: C/ ^4 S4 n4 B6 {
if nargin<7 | isempty(ub),) o! r5 t& F6 q, b8 [1 ?
ub =inf*ones(size(f));
7 a\" o3 L Z) O3 V9 p' o
end
$ @' `8 f$ [7 Q4 _9 c
if nargin<6 | isempty(lb),: O1 d1 u5 U\" e# \. u' R: d3 m# J
lb=zeros(size(f));
) K) O1 _\" e- Q9 \9 G+ ~3 U8 A. ?
end
& o4 f! Z1 \; f% A! r( A
if nargin<5,+ ~7 m! M- s2 o/ i6 U
heq=[];
, C8 [6 P/ h3 ?\" V8 l% E
end
. Y/ _\" b$ O2 ~# s a
if nargin<4,8 b+ Q! j; g3 G% l
Geq=[];/ J8 m: e0 k- g8 n! ~! q
end5 Q8 A1 i8 Z3 s! o1 J5 I
upper=inf;c=f;x0=x;
' U ~) A& k% a; L* y
A=G;
: g4 T- h* ^; U$ T; s. t1 w
b=h;Aeq=Geq;\" J+ H9 N) C( E. J! K& `\" z6 _- f
beq=heq;ID=id;
8 s3 P4 r8 Z( b, V3 n- T
ftemp=IntL_P(lb(:),ub(:));( y\" X$ f; w% G* P4 q2 d
%以下为子函数
' M& q' G& J4 o9 i\" g
function ftemp=IntL_P(vlb,vub)9 i1 t( M8 I* E* Z! F$ a' ^
global upper opt c x0 A b Aeq beq ID options;# i# }4 x8 F9 {( ?4 h- _' \- z. J
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);6 T- |+ i) D' Y8 \: R0 @, x
if how<02 ~9 e- r V0 _- ~\" j
return;: P; b. E+ r6 j. I) P% m1 _# W
end
- {8 R% {' m\" |$ A( l) K
if ftemp-upper>0.00005' Y* G, W9 ~* M: G2 I( n
return;; a) f! ?) s x% J
end3 C4 u0 E8 F. ^/ F+ L/ ?9 H b\" ]
if max(abs(x.*ID-round(x.*ID)))<0.00005
! g* }# u5 W/ l# {* B& I
if upper-ftemp>0.00005
( R& Q0 J( X( e2 ?# ^7 I5 X/ k
opt=x';
/ C' _# e! N/ I/ d( C0 Y$ U
upper=ftemp;9 T, v8 ?, D% n. U1 q: X2 B2 C\" d
return; |- K, S7 l' L( s- I3 Q
else- F' ], q# X6 E
opt=[opt;x'];
. O1 h2 {+ _6 d
return; q# P( V$ K8 g8 D; m% T) G
end( x1 [9 }) t! {% W% A8 F
end
2 A7 ]: V+ o$ S\" M# T
notintx=find(abs(x-round(x))<=0.00005);
1 e6 J4 \, T& K7 m6 t8 g4 e- t
intx=fix(x);
4 z2 B1 w0 Y0 k( a: u; L
tempvlb=vlb;tempvub=vub;6 I+ c( p. A: { E, y+ k7 E9 A: X
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+12 c* T1 c5 d' k/ K' ? T
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;
( `! e. v6 G* k) D% g9 x- u
ftemp=IntLP(tempvlb,vub);1 \\" q7 C6 J( N0 e2 g
end4 r$ V7 m3 E. |4 R- J
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1) s7 T' a7 d% S8 E2 L' E+ A' s
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);
7 u9 K( J. B3 G
ftemp=IntL_P(vlb,tempvub);1 h2 d\" R+ @2 i- K( m% k\" l
end