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

function [x,y]=IntLp(f,G,h,Geq,heq,lb,ub,x,id,options)
& n# J' @9 w0 x) e$ B6 Z
global upper opt c x0 A b Aeq beq ID options;
8 p* i/ o( b7 U5 y$ N
if nargin<10,# c3 m/ k2 T5 e0 y7 ~, f
options =optimset({});( E4 o! m6 k* j, {3 U
options.Display='off';0 u/ o }. g" ~$ _9 ~* P
options.LargeScale='off';
" H; l7 c! z) b2 ^: V
end# ]$ C! {! E3 |! Y9 W% W8 Q
if nargin<9,
; t4 K, z9 m8 V# a6 l# k! y
id=ones(size(f));1 c, C) x* \, I+ x: P( I. y
end % _( r( k4 e; g8 [. `
if nargin<8,* j% c1 q+ r7 I$ a# i0 ~1 p
x=[];
0 v! [7 r$ ~. r2 r2 L# {# U# Y6 p
end
; b6 ?. E3 t Y4 r; m
if nargin<7 | isempty(ub),& F5 T/ K$ U. L: b( ^/ F* o8 g
ub =inf*ones(size(f));$ m' M2 }0 H8 v3 ]& ?8 Z4 R! q" _
end
0 e; W" K8 d+ Q3 a; o
if nargin<6 | isempty(lb),
" C& ]* w0 L1 g2 R2 d9 `
lb=zeros(size(f));/ x) K& R9 f7 O4 ~: z% r$ H7 F. Y R
end
0 C- H0 }6 `" W9 {) C7 e% b
if nargin<5,
' C) K. s$ l* J2 H
heq=[];/ p% i* X- [5 r8 \
end
& x* n3 N- ^ }+ @) u4 D
if nargin<4,. x0 {" i1 d$ b% {' d7 a2 E9 I
Geq=[];
6 W. h& L$ f# A% j
end( Y' a* D5 z* E3 _8 `8 ]1 u' q
upper=inf;c=f;x0=x;0 ~* c" f, r+ M% T, [- Z+ L1 Z% f
A=G;6 U1 Q# p* W" ^$ n0 @6 f9 ^8 Y
b=h;Aeq=Geq;
# p' ~: q: u/ q* s j; J+ `
beq=heq;ID=id;, d; D4 K8 E% O
ftemp=IntL_P(lb(:),ub(:));
2 J6 n. Y/ G6 E0 Y$ i! w
%以下为子函数' Y! T' m2 q0 r8 D) h
function ftemp=IntL_P(vlb,vub); C j7 f i& s
global upper opt c x0 A b Aeq beq ID options;
$ ]8 e& S, e: M. I; B: K
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);% Q9 P6 W! D+ a7 [2 m+ T* K
if how<03 D- f4 o5 h/ t* g e% i0 `' f
return;
% K2 g/ q% Y# v, k. c
end$ n1 f* D; t( S. L
if ftemp-upper>0.00005
0 e$ P% E0 d+ n8 y( g$ k" D- ^4 C9 q
return;
7 u" }+ T" q9 d9 z
end
. d# s( [4 E8 s- F' q: T
if max(abs(x.*ID-round(x.*ID)))<0.00005- q7 w# l f; C4 r
if upper-ftemp>0.000051 j1 U1 d8 `% U
opt=x';
" Q Q' b$ K6 a* x5 W8 H: d7 V
upper=ftemp;
0 v! X+ G: U0 c, I/ j
return;
I% B" D. R' C/ K
else
. Y6 D( p: e5 n' K- W! g
opt=[opt;x'];! _5 v j% P' ?; q) s
return;4 u# T8 r( F; [' A; ]7 U: S
end
, Y# y3 ?& \# p/ x/ O. W0 |
end* w3 g e. c. {# W* m
notintx=find(abs(x-round(x))<=0.00005);+ ]+ `, A' T3 D' o8 E6 o1 B
intx=fix(x);* b9 t8 c w$ m# R
tempvlb=vlb;tempvub=vub;& ^, W; M- D& M6 X( v- y
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1% f2 t+ o: m8 ^1 @. X( `$ V
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;
# \; h) U% l# ^
ftemp=IntLP(tempvlb,vub);( t/ {+ I7 C$ m& V; A+ d
end" w* S- t) E" _# R5 y# @
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1); _) p+ W; O( u$ K$ h. d+ H% f
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);
9 S6 u" q& f9 A2 f! v$ H% V
ftemp=IntL_P(vlb,tempvub);7 }- N9 S/ e1 K3 j* a0 q5 j
end

复制代码