[问题求助] 请求一个整数规划的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)# g# t6 r3 D ^2 t
global upper opt c x0 A b Aeq beq ID options;
( y9 m/ [! d9 B3 O- G' ?. s- c# h; `: Q
if nargin<10,
9 ?0 C\" C8 p! X9 t: E3 ~
options =optimset({});
. y: l4 `8 s7 x( F% i! U4 V
options.Display='off';
& m- H, x& b6 n0 V
options.LargeScale='off';/ n- W, m6 F. q. S
end1 [4 E2 M; H4 G6 ^# j
if nargin<9,
: W* _. j0 B2 ^9 _9 o
id=ones(size(f));
% o) ]\" s\" ]* N' i# L
end
! i5 _- T( L; S4 l3 T2 m
if nargin<8,
\" X; ~0 k\" l9 l; l# q$ a# V/ v6 e
x=[];
4 u# J( p' q# @2 N; g4 j1 x( ]: L, t
end; C0 R+ i9 T\" d0 ~
if nargin<7 | isempty(ub),& w' y) }6 z* H, |. w6 y2 X7 W& w
ub =inf*ones(size(f));
# U {! h7 F! `9 ~% a
end$ B: l! I9 J7 j* c, `1 @
if nargin<6 | isempty(lb),
+ |8 [8 w1 ~4 A$ x: j( w) C
lb=zeros(size(f));' W/ s( C( [/ ?: y3 t; Q
end
! x* Y# G# q8 }. i# o( y
if nargin<5,7 \! b6 j& [4 y8 a6 a. m
heq=[];4 I* C; h- \/ v d, u& b
end( @; ]; \! L7 k, ~0 W0 X
if nargin<4,
3 x\" _3 [$ d: w& F$ p( W! A* C9 D
Geq=[];3 f) ~: j& B% c m
end
; M8 i\" w) |2 H& |
upper=inf;c=f;x0=x;; i5 g& I) r6 O
A=G;! i4 ?2 r9 G, B* \- r- i
b=h;Aeq=Geq;9 z: t. T2 J8 o& _
beq=heq;ID=id;9 S: ^9 E% P/ `
ftemp=IntL_P(lb(:),ub(:));
6 D6 K4 q( S\" _ W; v
%以下为子函数
9 _1 ]$ L: y/ q5 E( |9 Y2 G, D& p
function ftemp=IntL_P(vlb,vub)
) S( s6 I L+ `2 V4 K\" H
global upper opt c x0 A b Aeq beq ID options;
{0 ~$ k/ U, d\" P6 [
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);
! l- Z8 k$ }\" H9 @
if how<0
0 E8 ] ?6 [4 [' }8 C! X$ W3 v
return;0 Z* J! U8 D$ R+ I
end1 O. W& u$ X1 p3 P6 Y, g7 L6 K
if ftemp-upper>0.00005; ]- s, E( F! T$ l: h) l( b
return;
) F: O! k. ^0 N( z! z: v\" w: ^' s
end
/ ~* D* E2 ?# ~! J\" E& [' P
if max(abs(x.*ID-round(x.*ID)))<0.00005
8 x R2 |+ K- Y* g8 B
if upper-ftemp>0.000054 f& S) E( E3 J7 i
opt=x';
! O; i- W( e5 ~% z& D
upper=ftemp;4 q4 F' T4 j1 S6 @: _/ N0 Q$ n
return;; w' H0 o4 ^; U: o
else
2 l, R# p. r/ e1 W5 }: ?
opt=[opt;x'];
6 X& ?/ \' V4 }. h' p
return;& x R5 e6 i6 r% _! G6 M
end
' v O* t9 E' g* l7 H. a
end
7 g4 Y! P8 ^0 ? N' N* s
notintx=find(abs(x-round(x))<=0.00005);0 `/ Q\" R4 S) S
intx=fix(x);
- U\" E% @+ q+ I2 e; w1 D v! b
tempvlb=vlb;tempvub=vub;
: U7 x7 g* k' E. j\" ` _& @
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1
* l! ?* m6 |% k4 o+ D8 b+ f
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;
- ~2 m) C4 S/ a! X
ftemp=IntLP(tempvlb,vub);+ m, q6 O% \! n- B% G
end1 W4 l% a% T: c/ D' ^; f K( x( U
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)8 v6 m; }6 H9 q. U
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);1 W8 b# @8 G, m* k0 W4 {
ftemp=IntL_P(vlb,tempvub);
- a, A3 z4 n# u$ N+ N
end