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

function [x,y]=IntLp(f,G,h,Geq,heq,lb,ub,x,id,options)
5 @$ P& u) ~, N. i
global upper opt c x0 A b Aeq beq ID options;
2 `% R$ d9 j7 ^9 A+ b. y# w% q
if nargin<10," N% N& V+ A- Z6 _7 }% ]3 j
options =optimset({});; H( V) u9 ?1 I. ^6 C. ?
options.Display='off';$ |* B% E. a( Y& T- w
options.LargeScale='off';
+ W" `$ T) R {0 C3 i, ]
end8 F! X$ p, ? J/ e6 r. s- s
if nargin<9,& V0 K0 f3 F1 q. @' D, X# V! x
id=ones(size(f));' }+ l* {7 K0 \) K2 L
end 7 h k( S8 r( [+ r- ~9 m
if nargin<8,- }) L5 N' k/ [7 [! U/ ]# U
x=[];
4 ?( d1 b, I; t8 \7 p- g+ C$ f' g
end m. V5 ?1 T m( M' R7 m+ d W- A9 x8 q
if nargin<7 | isempty(ub),
7 m3 S3 E2 ?. l( t
ub =inf*ones(size(f));
' S) E1 e! V# j5 B
end
' b3 b% X& Y! G) L9 y
if nargin<6 | isempty(lb),
) O! H% V5 G, P3 i9 H
lb=zeros(size(f));
2 A7 u" I) ]) E, I. f
end
( l4 z* h0 ] P8 K+ w% q/ ]# y
if nargin<5,0 {" }% O- c% _$ z* \" j
heq=[];
: B3 L- J6 u- F' w. G! l
end
; ^. F! |( \( c: e+ N& ^! V( \
if nargin<4,$ e$ ~, |& B) Q* a1 W1 |
Geq=[];0 v5 x3 O6 F# k+ [" u& }% g9 S
end$ c I# L/ _' v+ V; p9 z
upper=inf;c=f;x0=x;) `2 a, y3 m4 t9 m6 W. N& T% t
A=G;
\4 e' l' J+ S0 s" w
b=h;Aeq=Geq;; J; _( V- _7 C6 u$ N
beq=heq;ID=id;; v& H0 Z2 {( q' O' g5 d, p
ftemp=IntL_P(lb(:),ub(:));+ i) x; e. E; y$ R$ s
%以下为子函数9 \0 t9 t2 Z* X' V# [ r5 K0 f8 G
function ftemp=IntL_P(vlb,vub)' k: \0 a4 W0 s' y ?! L ^. H* z
global upper opt c x0 A b Aeq beq ID options;7 q) c* J7 | z: `
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);. l. ^; u/ @* q
if how<09 Q- [& N7 `4 `! g$ z A5 [( c; h
return;
' c: u) x0 Q" L0 ^
end( I* H* ~7 K: l( i( u. n# [& {
if ftemp-upper>0.00005" k8 k( V4 u( H. U
return;
* P2 J2 X; p2 `) }' X, f
end! r5 C4 c+ m, n, R" @8 k
if max(abs(x.*ID-round(x.*ID)))<0.00005
$ p2 f9 ^; _$ W& X G4 ?1 [# U% `
if upper-ftemp>0.00005/ u8 F" ]" I' z( j4 h$ }! J
opt=x';
- l' O; k; X+ Y" E4 [" |8 I3 _3 o
upper=ftemp;
: Q$ ^+ O: B" @) |. C
return;, D, y( T5 o" _" s, r; ?% u
else* L4 M& b- h2 u% L- L
opt=[opt;x'];$ e7 l) Q, S% d% p# I# n. @- t0 M
return;) Z- Y3 f2 F/ ?7 t! Z
end
2 q! O1 k( g! O7 U/ R9 A* v
end
1 M; @- p5 r- k1 i* {
notintx=find(abs(x-round(x))<=0.00005);/ i# z, v1 z* v0 ?" l# x+ v4 D' R
intx=fix(x);
" q& V3 T* J4 r* y5 Z, j g j
tempvlb=vlb;tempvub=vub;$ S9 c7 ?/ e5 c3 K% I* A2 i
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1 [: Z2 ]/ t" `3 n4 P+ F! i
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;
/ F$ y" s# R5 O r
ftemp=IntLP(tempvlb,vub);5 Y& D' n8 o+ R+ f; C7 u
end
b+ \% \: M# X# @7 K' W
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1) A+ c. `7 ~6 N8 L: f8 o; N
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);1 Y- j h& F5 W+ ?, b
ftemp=IntL_P(vlb,tempvub);; L1 K' E/ d; n7 z# Q& t3 @
end

复制代码