[问题求助] 请求一个整数规划的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)
5 y. U6 [( E% A. v
global upper opt c x0 A b Aeq beq ID options;7 j+ X7 C2 ]2 g$ `
if nargin<10,# Q# K Y( f& ?; ?5 B
options =optimset({});* W: X8 F7 x/ g: K! A) K
options.Display='off';
R# U0 r- h$ ]; i5 X7 p' Y
options.LargeScale='off';: r! O, w; M% V
end
2 c& q# `) N8 S4 e- o5 j
if nargin<9,& _; V5 f% Q U( i. f- | j
id=ones(size(f));
( \& m: K2 |# ?# ?$ J- I
end
7 u7 ?& J, o3 |\" f- ~6 r% ^
if nargin<8,6 j0 j2 I3 m4 ]5 R
x=[];
y [/ Q; n6 y
end
. ~& Y ~' f2 [: O
if nargin<7 | isempty(ub),
7 r1 \1 K\" u3 x: C- z0 g B
ub =inf*ones(size(f));
1 j/ u6 Z G0 H& B9 I! l2 h% {
end2 h- j5 a# ~7 y1 g
if nargin<6 | isempty(lb),
# l\" [! q! c; X
lb=zeros(size(f));/ i. ?5 w) {. O\" e9 r% P
end! r: n& G# |1 Z\" k( L
if nargin<5,
7 e. A& q; j3 _1 S3 h% L
heq=[];
5 B+ Z* r8 R% f- q+ Q
end9 h+ Y( `' Z# c- w# B. l' \& d0 W; q
if nargin<4,
& {% z \) |\" N' |+ p) K: e
Geq=[];3 q+ v2 M7 E. A) L: F5 N
end% D, [/ s\" P7 L9 ?5 P
upper=inf;c=f;x0=x;
) l3 K) n: ?9 r( Y1 J e8 u/ h
A=G;
$ I$ L( W# ] d* ^- _& G
b=h;Aeq=Geq;
. ~; d& Y, V! h6 @
beq=heq;ID=id;. l& |1 ^! @+ [
ftemp=IntL_P(lb(:),ub(:));
( ^4 ]: P4 P1 b' a
%以下为子函数; _: s; ~/ a3 g, o
function ftemp=IntL_P(vlb,vub)/ M( H& b2 p4 A. g7 R: ]9 s( v/ X\" C
global upper opt c x0 A b Aeq beq ID options;7 ^# v6 B* K, x9 F+ {4 L1 D
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);9 l3 b, K* T8 i* o5 u5 n I; i% P
if how<0
6 s+ `9 @& R S3 ?9 i- ~6 Q: g
return;
0 u0 b3 r8 [. {: r! {& ^
end6 S. z8 ^( h, n& ]1 l: Z
if ftemp-upper>0.00005
! ^' r% l: i7 ^ H, R- P8 e: V: q
return;: \. o6 I/ k3 ?- s+ { _: g
end
6 z+ b; \2 X7 V+ }
if max(abs(x.*ID-round(x.*ID)))<0.00005
\" t\" S u+ ?% ~: T, R
if upper-ftemp>0.00005
1 ]& l a( G/ k
opt=x';8 W0 b1 S9 L1 ]1 K# r0 G: r/ }
upper=ftemp;4 J J' l& b- ~
return;% J) [7 _1 b }- w, T- t+ O% k
else- [% {, l# E) J! V* s+ |
opt=[opt;x'];4 g' E T' R9 ~& z2 G! a- ]
return;
\" w2 u3 Z) e: A- k3 Y5 z
end5 B, V: R' a8 ]) f' R+ _5 j6 C
end, O1 k7 Y+ M; H6 s5 d2 ~\" W
notintx=find(abs(x-round(x))<=0.00005);! Q/ X\" Z3 ^$ J
intx=fix(x);
% i* k; P2 d' M! J5 t! A
tempvlb=vlb;tempvub=vub; Q: F; @0 L% \5 z& n$ t
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1/ m$ c( C& |5 O6 ~+ ~
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;( R+ g4 Y u) @) W
ftemp=IntLP(tempvlb,vub);
! \% O8 P6 |! b# I7 U' s7 Z
end) c/ K; v* B. o
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)4 E& I% G Z J
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);
8 `$ V( S) K' T0 A\" X
ftemp=IntL_P(vlb,tempvub);9 S, g _\" Q S, q% M
end