请求一个整数规划的matlab程序 - 数学建模社区-数学中国

function [x,y]=IntLp(f,G,h,Geq,heq,lb,ub,x,id,options)* e H. Q) U) P; p& h
global upper opt c x0 A b Aeq beq ID options;, o7 e, k4 I9 `$ @9 o) H
if nargin<10,) ^/ k( s+ Z( X: K4 \3 S
options =optimset({});' |: [# D- m2 @+ j& {( ?3 K5 C
options.Display='off';' H6 m2 N! Z$ G+ B' n5 A* t
options.LargeScale='off';
0 a" @. z9 @, {" V! {- p0 r1 v
end
6 e1 ?# r. P6 D/ T* f$ F
if nargin<9,% D Y: ^2 Q+ ?
id=ones(size(f));7 l6 l o$ {$ k3 h' z) h% ^
end 0 |/ T- Q3 W0 J: q2 U- P
if nargin<8,% S8 ^. D# r& k& t
x=[];
; {& ] g) ~3 @
end
% S2 o# n) I+ S5 M: W
if nargin<7 | isempty(ub),
, Z, F4 r+ @; {2 `9 B
ub =inf*ones(size(f));
q( m. y4 @9 G0 I" r5 z
end' P3 R i6 m2 I5 | D' {
if nargin<6 | isempty(lb),
3 h0 R# c" w( d( m* `
lb=zeros(size(f));
0 I' n/ K$ J& @; l% y5 T
end
' m" l2 B9 }. Z2 P! `
if nargin<5,
^- q2 R% g0 a
heq=[];& a1 H- v/ C/ B* N
end
1 D7 S6 d. Z# b! C
if nargin<4,) E$ M* N1 X# ^! B4 |
Geq=[];' |+ a4 O7 o0 x. [% m0 g+ F
end1 [/ \; a! h3 e' r5 n5 x
upper=inf;c=f;x0=x;2 E% N. h4 h# K9 R9 h0 }
A=G;) [. @$ q1 r1 D7 c$ O a
b=h;Aeq=Geq;; ?# G: j& t6 }& f$ K- _
beq=heq;ID=id;& @$ o4 _2 i. D# ?3 a) y9 Y
ftemp=IntL_P(lb(:),ub(:));
, W* H* N6 m6 N1 S/ v1 x t' t5 E
%以下为子函数! Y1 m0 ^& A: W+ Z
function ftemp=IntL_P(vlb,vub)7 L: d3 C6 m$ F9 ~
global upper opt c x0 A b Aeq beq ID options;
* @6 o. P. R" ^
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);
, z8 K1 a: w$ h+ f# N9 ?5 @
if how<0$ u+ N& p2 a0 k- d; p- w
return;
2 d$ m. q$ A. X9 Z
end* H1 v- J2 n8 x/ u+ S
if ftemp-upper>0.00005
9 z0 n' h+ a% Q3 Q* ?' G! S. t
return;9 g9 ?/ g+ T5 S. M6 q% ^7 x3 I
end
4 O5 k7 `+ o3 ~: Q% W+ ~* Z: d
if max(abs(x.*ID-round(x.*ID)))<0.00005
! J& i; L% q" ^) H/ @2 c& C: H" N
if upper-ftemp>0.00005
: f" K; |* X+ q
opt=x';& ^7 L! ]/ {7 F, i
upper=ftemp;
4 _ b7 l! }$ r6 b# L4 G3 J
return;3 }# N% L, ?- R5 x3 @
else8 O- h, r C3 B; O
opt=[opt;x'];) n9 E' k7 t1 v, K5 V+ P# O
return;- q# P" }( t6 \4 q$ a% [8 }
end
, s$ o. X7 ~. t$ z: r
end! h/ `; ^% v; t
notintx=find(abs(x-round(x))<=0.00005);/ Z/ F5 @8 _# x9 M8 n. w! M U
intx=fix(x);
6 M' a1 e/ e$ x" ]" T
tempvlb=vlb;tempvub=vub;6 d& p9 p; Q0 c% P4 O( j% p3 R
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+17 j) S c+ U2 H! ?; t6 A5 v% i, Q o
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;
% l/ s0 |, ^# }1 t& t- r
ftemp=IntLP(tempvlb,vub);1 |" P. m) O/ L# _ i
end: N* j; U4 k: z$ ^% n$ b( `5 Q
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)
& n4 d8 }) b: W* g4 h: T
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);$ R* t/ P* f: [
ftemp=IntL_P(vlb,tempvub);$ ^+ _( P* T/ @( Y7 V8 @5 a/ t
end

复制代码