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

function [x,y]=IntLp(f,G,h,Geq,heq,lb,ub,x,id,options)+ a9 P4 D! W3 _+ B8 s7 x1 ~
global upper opt c x0 A b Aeq beq ID options;- ^: w5 N5 \+ p* H7 v
if nargin<10,
7 D+ r1 q! @1 r
options =optimset({});& X/ d; s( m1 O$ |1 R3 y
options.Display='off';1 o& J* ?2 _8 z; T
options.LargeScale='off';! J" D' g2 Y" _8 U
end9 r: t! Q( \6 i7 F
if nargin<9,* b7 L, b0 `; x. l
id=ones(size(f));
; Y5 W6 X+ A5 x \" h
end
4 i8 X1 f4 {0 o* b
if nargin<8,+ |$ _( L7 M( n3 P3 i
x=[];& Q5 ]! K J8 G) M, x1 d5 _
end) r) q& J+ i/ O
if nargin<7 | isempty(ub),
+ A. w* ~$ ~* ~! r
ub =inf*ones(size(f));4 h' g$ h# `" u/ ]' u6 g* j
end
$ N7 r, e' c4 P
if nargin<6 | isempty(lb),
3 V5 J3 v, r0 o- E! s
lb=zeros(size(f));
) E3 i+ f2 k3 ^
end
, n, _% r \ f1 T$ d8 p! A8 @
if nargin<5,
" w3 `$ b' V! E( h" x3 i& C
heq=[];, F l& c9 p# D! c9 m: a
end
9 Z# | q. Q; L: }0 ]% ?$ P
if nargin<4,
7 q1 q+ h9 q! |
Geq=[];) G/ L( i# v! v5 }' ?
end$ b& h: ^+ Y1 N# l9 X
upper=inf;c=f;x0=x;) I% `& L8 o" _* p6 K4 E9 o1 x
A=G;
( k3 N+ f5 ]% {+ Z/ P
b=h;Aeq=Geq;
0 ^! d3 q9 }2 d' n
beq=heq;ID=id;( M% s6 c \. u. ]. `9 u$ J
ftemp=IntL_P(lb(:),ub(:));
+ V+ c% S! ?, P
%以下为子函数2 k+ E& j. `3 u) w i, }" q5 d
function ftemp=IntL_P(vlb,vub)6 k" O" K& `# a5 w) ]5 P
global upper opt c x0 A b Aeq beq ID options;# v Z3 S2 [' r9 M2 a( i
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);8 w9 u p/ W. g3 Z
if how<0
3 O3 \6 a* p) u( b! E6 v
return;
5 N; k* r2 b& B2 f" _+ l: c) v) \
end
: c2 F% D7 T) y
if ftemp-upper>0.00005; m" ~4 {% S/ T
return;. y% n" @3 E3 q
end
- Z- y* c' U( J1 a" U
if max(abs(x.*ID-round(x.*ID)))<0.00005$ C' a5 f5 t) v5 B
if upper-ftemp>0.00005 w% f( x# T6 X
opt=x';
& h$ {( R7 y; ^, g& m
upper=ftemp;0 t& Z: h/ _7 [7 Z, Q# [5 [- }
return;
7 m- z7 K/ H! a. ]0 X9 }: k
else
4 P a5 E3 y& s/ E/ Z( H9 Q, _6 a
opt=[opt;x'];
% l7 R0 k2 [; `" L4 }2 S5 d
return;! X! p* x3 P& y% J( r
end
: U! a- |' N/ y- O$ V% z
end
5 ~" t3 ]2 K* ?7 z- z+ C o% E
notintx=find(abs(x-round(x))<=0.00005);+ n- u7 W: M1 b! n* {" h
intx=fix(x);
' _; G4 F1 M8 F3 l. W
tempvlb=vlb;tempvub=vub;
7 D4 Z) |# [; t1 D2 }
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+13 z+ L) i* r' L1 n
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;
2 {7 ^" I; `7 X a; Z( E$ ?
ftemp=IntLP(tempvlb,vub);; ]5 t% P9 o, S
end
5 U# Q- l6 o3 R1 f7 D* T, s, O
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)0 v' R+ N& v" a( o
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);2 ?/ x8 F$ Y+ o3 S* A9 Q
ftemp=IntL_P(vlb,tempvub);
7 Q2 |3 p/ I2 U6 e6 Z* U
end

复制代码