[问题求助] 请求一个整数规划的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)/ B1 ]/ ` Q o\" X o; m
global upper opt c x0 A b Aeq beq ID options;
7 M9 W C7 s* T* Z0 B) F7 X
if nargin<10,
, D: p+ d; n! m1 F7 A
options =optimset({});
. O( v5 }$ L4 Y( P# S, ^
options.Display='off'; d+ o0 u8 q, }$ k/ ]' k' x$ Y1 g
options.LargeScale='off';3 l# g3 ?2 W# v5 n3 }2 [2 _; F
end
' V% }) n\" {, }+ e8 p0 x0 w# Y
if nargin<9,; |# m' e- H- C& I. r# x% N4 R+ L
id=ones(size(f));* L. [% U8 A/ o$ j2 \$ V
end
9 u5 q- n+ _& o! \$ T: d7 R
if nargin<8,
9 u5 Q U. B$ R2 o( x0 I8 i
x=[];
; G. S6 _$ v- s; v
end P; ]: [1 C1 J2 g# B
if nargin<7 | isempty(ub),' ]- ^3 T8 G8 R+ z# K) C
ub =inf*ones(size(f));
1 \6 f6 ^9 m: s' m% `' }
end- c( K1 j0 T\" I/ V# B+ Y a* {% q. [
if nargin<6 | isempty(lb),$ \. ~3 N4 I7 p% C
lb=zeros(size(f));
% r. [! V) ]* |% a2 _# Y# V1 Z
end; }5 T- T5 V4 O8 S& [' c' m1 V
if nargin<5, J( J3 O# x9 Z7 t9 P4 J6 d
heq=[];
6 q' z7 N5 }\" h1 o1 p5 V: N- ^
end
1 L' w; t& |, v+ t% i& O( Y3 ?! o
if nargin<4,8 C3 R% z, H, w6 i
Geq=[];
# Z) z! J! }! _; b& x
end
! Y$ p' r. [. S5 h6 N
upper=inf;c=f;x0=x;) g: B3 n4 s# G, W, @3 B1 Y4 ~
A=G;' n, U! b0 t/ o) |5 I9 Z% s/ x* D- Y
b=h;Aeq=Geq;. m3 u) ^$ \) r) S4 k\" \. J! y
beq=heq;ID=id;8 O5 A$ g8 O7 m' W Q
ftemp=IntL_P(lb(:),ub(:));
+ t# p: j; d8 m5 ]& ? k' g! \
%以下为子函数% k\" y& v1 j0 B& S: ]
function ftemp=IntL_P(vlb,vub)\" w+ o- U+ U* @5 {' ?! b+ i) J6 h5 l
global upper opt c x0 A b Aeq beq ID options;
0 A. O3 w2 ?2 h9 W2 }. z- ^' X( m: J
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);
# _5 k0 m$ P7 k9 i2 C5 c\" _
if how<02 k$ ^\" N2 l R\" ~; r
return;
* ~1 Q# k. X/ b' N
end
0 ?) q/ M* c, z' r( q# E1 W6 e
if ftemp-upper>0.00005
a& J3 k0 y' }- ^ v
return;2 o\" _7 c4 N4 i/ ^5 b
end
5 L' Z: l1 }5 y9 Q
if max(abs(x.*ID-round(x.*ID)))<0.00005- A' O8 k0 G+ [1 I* |4 r
if upper-ftemp>0.00005 L1 R' i. U! `# R' S% o( P: q
opt=x';1 y ?: z3 d3 Y! ~* U$ O$ K/ \. T
upper=ftemp;
8 c. d4 |, o+ A. [. f- a# t
return;
& r4 V! U# V9 w$ l! J
else\" A$ I% D/ `2 i, y6 o7 M9 c\" v8 [ j
opt=[opt;x'];
9 f) |, P `# {5 k3 {/ |1 h% U
return;
H; D6 U1 X! U! N- R8 l
end% V4 b, z% J4 c
end
7 e X3 g% n+ @' L
notintx=find(abs(x-round(x))<=0.00005);2 j& h+ {\" e( M1 p
intx=fix(x);
. a5 E9 q$ Y; s
tempvlb=vlb;tempvub=vub;: x# y1 W\" [\" U3 ~2 Z, b+ l6 J
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1
) O' D( P0 X# N\" ~
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;
; x\" v$ W0 O' c( o4 X, z, n) ~; ], A
ftemp=IntLP(tempvlb,vub);
) X! J c* r0 ?8 R; {7 t! M) s
end
/ X. {) X$ l, g; q
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)
4 t3 S$ N( t0 X% T
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);
3 m; U3 U; L( l, @$ b
ftemp=IntL_P(vlb,tempvub);3 ?' @( u* p! l0 P0 z; S
end