[问题求助] 请求一个整数规划的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 V2 A# a8 F- J) {\" `( E
global upper opt c x0 A b Aeq beq ID options;
~4 ]5 Q0 q6 h5 o( q9 d+ V0 d
if nargin<10,3 Z8 v- G2 n3 |' v$ e' K6 }9 R
options =optimset({});. V9 c5 y7 @5 G: [: H) g
options.Display='off';
) `; Q' a3 R& y0 W
options.LargeScale='off';& J' i4 [0 k) }$ a
end1 U* Y+ g+ {. w6 y7 T& J3 M, U
if nargin<9,
: Z) e e {; A
id=ones(size(f));+ l7 O! z- K) x6 a k. ^
end $ a# `( c6 m- H\" K\" T\" S
if nargin<8,
! @2 T, ]3 \6 p! }/ s. X/ ^! G0 F6 f! F
x=[];- R1 s9 z* a6 o; U. {* I
end
4 P: G) g/ v7 I9 I& U
if nargin<7 | isempty(ub),
7 {' t\" n( J\" [
ub =inf*ones(size(f));
5 j/ X% r3 z* |4 i$ q9 @
end& C6 [, {* t2 l3 U. v. W
if nargin<6 | isempty(lb),8 y\" Q* E7 y2 Q4 r6 u
lb=zeros(size(f));% {% J$ M5 M5 W y
end8 P N( R8 @* Z8 U2 f V( j. B
if nargin<5,6 v: y6 H) t( _: S* G( E% a6 _4 y
heq=[];
: o. U& k2 Z3 Z. s0 g
end3 r/ O/ C! C$ [5 }6 W2 T% N. g7 I
if nargin<4, N, b5 r. ?! Q
Geq=[];$ ~2 g8 g3 C2 k' x- s
end) C1 n. B9 J( P5 i M8 A1 J
upper=inf;c=f;x0=x;- G4 t1 {: a3 E7 A0 L
A=G;
6 }' ? H$ `2 b+ m. `* H\" x
b=h;Aeq=Geq;
; \! H$ ~( c/ A K
beq=heq;ID=id;$ g' J+ N/ [\" }4 y: w
ftemp=IntL_P(lb(:),ub(:));
! Q! ~9 B4 J8 C: ]. B/ t
%以下为子函数7 h9 ^' l, O/ }9 P\" i6 j, c( o' z) k: M( a
function ftemp=IntL_P(vlb,vub)
4 I9 J4 ]4 s6 N7 y: _2 K* g: m
global upper opt c x0 A b Aeq beq ID options;
^! \6 q\" t\" Y$ i6 L. S! Q$ b
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);
* y; }) i. n' L; O) ?9 }
if how<01 ]( E+ p4 B) L1 {8 s+ a% U' j2 D
return;7 S% \5 I, H* v, D
end! A3 ~' P5 i+ Y
if ftemp-upper>0.00005
. E0 I l! _( d$ Y8 ^7 o& a
return;
% T6 h+ |9 W8 |4 o. O- R
end
( i, c\" n& e2 \4 W' D3 ^- g
if max(abs(x.*ID-round(x.*ID)))<0.00005
1 h. n8 j5 H- K
if upper-ftemp>0.00005
# M8 N3 n9 @; u# j1 ~8 A# T8 Y
opt=x';2 b$ A6 @9 W5 ]+ S2 m% c
upper=ftemp;
- `' z9 z9 p! S
return;8 ?1 b- N2 _; O9 q& ~( T
else8 f+ N1 v0 k9 Y1 K! }\" ?
opt=[opt;x'];3 T% `5 d4 t. m8 L! u+ z
return;
! [\" e& C) V/ m- X! _
end
2 x3 x) z# E' @# S
end
- }3 y1 g+ P\" z3 P- ~. _4 [
notintx=find(abs(x-round(x))<=0.00005);5 t c- T4 L% F; Z
intx=fix(x);3 m8 t D* B0 P+ A, g
tempvlb=vlb;tempvub=vub;. S+ A1 i: R; ]8 l
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1
( W) ?+ b6 Z4 A4 V# v+ H% I9 N1 K4 c8 \
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;- j- O- f3 n5 ?$ S
ftemp=IntLP(tempvlb,vub);8 V+ E5 g3 ?! a
end
9 W$ b2 ?* e4 x c5 g1 _7 B
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)0 `7 R& t7 O- I! m\" Y4 e
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);
4 B+ _' [1 p0 n! M) p
ftemp=IntL_P(vlb,tempvub);
) R8 ]: D0 ?- ~! Q9 g. {( Y* \
end