[问题求助] 请求一个整数规划的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)1 S2 @* h3 ?. o/ y: l+ V8 t
global upper opt c x0 A b Aeq beq ID options;
2 ]8 d! p- C$ A' o7 ?
if nargin<10,1 S: P5 n6 w4 i* V } v
options =optimset({});
# l7 X. v+ a0 A& t, [2 J% f9 }
options.Display='off';& {/ b8 ]$ T* X0 O* A# u# E& Y) g _
options.LargeScale='off';
- e( h8 N: v! F1 U
end% y4 a7 z$ c; r' Z
if nargin<9,( p( L0 j/ z; J
id=ones(size(f));
1 z3 _5 S- F/ w( s
end K2 U, h\" p# }8 h4 I/ u
if nargin<8,
* ~% m1 G7 ^3 b6 f
x=[];3 w N! w7 z/ a3 }
end) o3 Q' |7 u* S+ @
if nargin<7 | isempty(ub),
( S5 z' J9 _% }4 D7 ]; |
ub =inf*ones(size(f));
2 l7 @& D* p1 A
end# S1 @5 @3 O4 M* ]8 J6 I2 `
if nargin<6 | isempty(lb),
& s3 M2 Z; @# f! u+ _1 ?# X
lb=zeros(size(f));. g\" i+ k2 Q. c- N' h; }& t$ c
end
7 M1 H1 Y2 V# U% G
if nargin<5,
. Q- a* v' a- ?* [
heq=[];* K* T6 v1 w' j9 O: ?* X) z
end
& Y% X! V, L) c\" h5 v0 p1 M x
if nargin<4,
: U) K\" Y* X! u4 h2 f
Geq=[];
- _. u S' ~) c4 m* q) `! l' n8 {
end
5 b) G! R I) t
upper=inf;c=f;x0=x;; S9 I% n6 w/ }4 ^
A=G;
1 j j! S8 Q9 q0 v: W
b=h;Aeq=Geq;
* ]& Y$ u\" K) H1 K) H
beq=heq;ID=id;9 N9 W( q0 B\" h1 h* v' Y' ^
ftemp=IntL_P(lb(:),ub(:));! V$ _+ w/ o9 x( t: a9 V2 F
%以下为子函数
. `! u) _( F$ P8 E2 v& ~: F: |+ T
function ftemp=IntL_P(vlb,vub)$ U% b* Z) B. o8 w, T
global upper opt c x0 A b Aeq beq ID options;% G\" J) |# G: B' v5 I+ p1 x) r
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);
V: W T& V9 S2 j# Y% ] b6 u
if how<0
1 i3 Q/ X* \, D2 V8 r$ S, @; q3 _
return;9 o8 o o+ S$ i
end
$ X& u- h( J4 t4 ^
if ftemp-upper>0.00005\" C' E6 ^. f+ v# g( E
return;
4 C6 p8 Y: w( Z Y
end
6 l+ k. }2 A+ N q7 j$ S
if max(abs(x.*ID-round(x.*ID)))<0.000055 z& a' k# \3 V2 _1 n3 [8 i% I$ U
if upper-ftemp>0.00005! b( g: U; m4 G9 d6 k
opt=x';
. f5 K9 G2 z! C1 [3 j\" n( S$ K
upper=ftemp;% H5 \. e: Z7 C j1 H: l/ d
return;
/ a6 Z: o* R7 B2 ?\" a/ z( g
else8 e9 _2 b9 G/ A% Y, _
opt=[opt;x'];) u/ [5 x6 i0 H- U5 i6 w
return;
* ]; @( b9 f0 V: H
end# K0 Y; ^, e! U7 M; b
end6 S# i( S+ [8 O0 R7 _4 A l+ P( L
notintx=find(abs(x-round(x))<=0.00005);0 O* U# Y0 f' Z7 S3 L v' X; C\" {
intx=fix(x);8 Y) C* o1 u ]8 Q/ P
tempvlb=vlb;tempvub=vub;
- q: z. S$ f) x5 q/ h
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+15 |# h9 L+ W2 n\" r
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;
5 m- o0 c\" `7 g4 x. G3 R
ftemp=IntLP(tempvlb,vub);
: A; L! a: V0 v; R
end- @\" N1 a( y/ a\" R. @
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)
0 a0 R# d6 _) n2 T6 g
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);
, J4 h+ V8 l# v5 i1 a: F$ ?9 o. G1 v
ftemp=IntL_P(vlb,tempvub);
: J) y! }& f* C& u* h
end