[问题求助] 请求一个整数规划的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): B# x( H- j+ F0 s
global upper opt c x0 A b Aeq beq ID options;: a% N1 S+ L( h. X# O% @
if nargin<10,
5 D2 }' H; e* M- \: p5 b9 P7 [
options =optimset({});! u. K- ^+ C& _0 C2 X\" `0 o( w
options.Display='off';7 y: n( s' K% ~1 l2 h
options.LargeScale='off';
6 p/ j$ N7 s: O1 O! L
end
' S( x+ T. i, O, l* Y- q# U1 }9 t, r
if nargin<9,
/ W D4 T- }& A
id=ones(size(f));
* W0 g- D. v/ f# c- w
end ! Q, `9 U1 l- r
if nargin<8,% h4 @: m% Z6 O. E. e3 ^
x=[];
; @1 c& ~; `3 F; R
end! `# J) W- ]) u5 T. z
if nargin<7 | isempty(ub),. B6 c* K. z3 Y* r, p
ub =inf*ones(size(f));; m! Y! ~1 Z, T
end( E9 S' q) y0 o; ^9 }; v. [
if nargin<6 | isempty(lb),
! T! J5 S) K* F0 _\" F9 T
lb=zeros(size(f));, ]7 R, {- q+ f: V: N
end
S. N. B( r$ S$ U
if nargin<5,
r4 K0 s\" J! C' ~! ]+ B. W# Q( B9 p
heq=[];
& B( X2 e$ z- [
end
* v$ Q% M& \\" x, K9 _
if nargin<4,1 ~$ ^. y, W# U' y3 T( N4 `. I7 d
Geq=[];
# Q# i6 ]- x0 e1 ]
end
2 _7 N; }+ s. _! ~\" P! I
upper=inf;c=f;x0=x;! _3 x$ |6 z+ ?: K' P
A=G;, f& T: J/ r\" y, @+ X F8 j4 d
b=h;Aeq=Geq;+ ]0 w ~. [8 V1 s/ i
beq=heq;ID=id;# h7 }' u8 _3 T2 i4 s- b
ftemp=IntL_P(lb(:),ub(:));
$ B! R8 r8 Z. \; [2 s' @* s. X
%以下为子函数
) ]7 o: k' _ u0 m1 M
function ftemp=IntL_P(vlb,vub)
) C' Y) c/ w$ C& \
global upper opt c x0 A b Aeq beq ID options;
7 S\" G' t% N( X
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);; b\" e) F, p) Q% @* i2 e! |
if how<0
2 w8 L5 ~0 X4 d) x j2 ~8 u; Z
return;
5 [\" m! P1 r X
end
7 u\" E1 L9 M$ a$ i/ T5 \3 W( V- o
if ftemp-upper>0.00005
' v/ A4 g: s/ R, c
return;
% E7 W) p. o- X% S; l) y
end, ?2 O- B+ L; S7 g* s
if max(abs(x.*ID-round(x.*ID)))<0.00005+ Q' b9 V) |0 D( ?' M+ v
if upper-ftemp>0.00005
) q$ W& p9 W\" |9 z2 O% B
opt=x';
& t. m7 n1 Y: S
upper=ftemp;; y9 Z$ p( r0 j4 o
return;
- q- ~! `/ G0 V* {
else+ d% F9 D! P$ N/ ^
opt=[opt;x'];
\" J/ W+ U- w2 r8 \& [# b6 P
return;
2 ?( T& l( A$ d! b8 `1 H% ^( y& O
end5 p2 _: Y' J* E' Q$ j2 z
end
! B& Z% Y! Z6 b; s3 e
notintx=find(abs(x-round(x))<=0.00005);2 O2 j5 [1 O9 Z\" K: g
intx=fix(x);
2 v; F7 ^6 ^2 v\" }8 n& L2 p, f
tempvlb=vlb;tempvub=vub;
; O, G: t1 Y1 i1 e\" V0 \
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1
: J+ ?5 r3 G9 B: J5 @
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;! O8 c, u6 N- v$ d\" E; o' s
ftemp=IntLP(tempvlb,vub);7 X! h$ P3 k5 x; S* F6 h4 g5 F
end
& r. u( ?% O4 O9 q% o\" @
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)) g! Q/ D+ @; e$ M! S\" D2 E
tempvub(notintx(1,1),1)=intx(notintx(1,1),1); Y& W& n* |# M- V
ftemp=IntL_P(vlb,tempvub);
, p+ m. k+ E) r: n3 x6 K: e5 p
end