[问题求助] 请求一个整数规划的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)
2 N; F7 p9 d. G3 w8 v: [, r9 P
global upper opt c x0 A b Aeq beq ID options;
- i5 ?5 r2 b3 K( F. w, c8 G( N3 p
if nargin<10,
`! b; S/ b5 e( X8 k
options =optimset({});
\" h% u9 C6 U2 L
options.Display='off';2 _4 y. O$ j, ~0 o& p
options.LargeScale='off';' H Q9 T3 a& E' ~; Y; v/ u$ B
end% ?$ w\" K2 l) f8 A
if nargin<9,; T- x& W\" t- O9 j# i r2 `
id=ones(size(f));1 r6 l! D3 q\" U: C6 m( d
end
: G: P' @* S9 E/ _1 h1 a# j( G) Z( u$ e
if nargin<8,
7 `+ |% r4 V# G
x=[];* x' J. U, D: B
end
9 t7 c4 l6 J- {8 z
if nargin<7 | isempty(ub),. B1 \$ B( U# a0 N. `4 ]( i) |& u0 t
ub =inf*ones(size(f));* ^% u1 c+ t/ m% _
end% A) n# A; a( c. `% p0 l9 {
if nargin<6 | isempty(lb),
$ Y2 s, @% A, S
lb=zeros(size(f));, j2 H: z7 B+ X
end& U! D# m\" A7 v\" i* @& y# h
if nargin<5,$ j2 U& w9 l# O
heq=[]; c% Q8 N+ _! Z! H* w$ X: R# I: D
end
5 ^$ P& R! a R( w3 ?
if nargin<4,) E! i f0 `0 A# [! c0 D
Geq=[];
+ @* |$ S9 ^* {) z R\" p/ D
end
\" b% k/ p0 Q0 ]& {& }
upper=inf;c=f;x0=x;. f) x\" J7 {0 |0 d
A=G;
+ u! L/ U9 h6 W
b=h;Aeq=Geq;9 f2 S1 f D, V7 q+ l% Y
beq=heq;ID=id;/ F2 H& z0 ?* \3 q# @
ftemp=IntL_P(lb(:),ub(:));
% B' y# [9 i6 ?6 Q/ x1 c+ j0 F
%以下为子函数' L- S0 B( Q5 V/ Z
function ftemp=IntL_P(vlb,vub)
4 A& T1 W% d0 F/ Q+ G+ P
global upper opt c x0 A b Aeq beq ID options;
( C F8 ]- X1 l6 {& I7 |5 ]
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);
. q\" c0 M& d0 A% i: W, @
if how<05 T* a\" d( i! D* D
return;* ^: S/ C* y1 b D4 v$ H. ]( U4 o- A
end/ w8 |( w) b- t6 q& q% w* r- c3 ]3 e
if ftemp-upper>0.00005
8 A: t5 m& a6 C9 C* x) S* e
return;
6 |0 Q3 O Z\" A5 N0 O' j
end( i P5 U2 T- ?/ t A
if max(abs(x.*ID-round(x.*ID)))<0.00005+ B6 V; H9 }% H
if upper-ftemp>0.00005\" @0 F, c* n0 Y3 w$ x
opt=x';
& p- T' d\" ]& e; a- d- r( Y5 ?, w
upper=ftemp;- H. }- s\" l4 P5 q: y, L4 ?
return;
\" u1 {* a( f0 M% n\" o n% k
else
5 t, M+ D+ V2 K
opt=[opt;x'];+ b2 e- c\" ~; C+ \9 ^6 P6 L7 m( I4 ^
return;
4 P. P0 Q! D4 [1 W2 p
end( I) U) E% z a- F8 n% Z
end
8 {- R; a7 s) V
notintx=find(abs(x-round(x))<=0.00005);: q% D9 ?( \8 c
intx=fix(x);
9 o8 q7 Z+ D6 i& X9 G; `
tempvlb=vlb;tempvub=vub;+ W! @. d: C f+ H8 t, V\" _/ Q
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1& L$ A3 X( r6 P+ Z: L; ~
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;
6 R. g1 w$ z2 M* u& w2 B
ftemp=IntLP(tempvlb,vub);4 l$ W) g( w- q+ M) s% {
end
- S4 _8 M. i. }5 I2 n8 n& \2 [
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)
% |! b7 s! |( Z- x7 J
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);
6 U, j& A, x% D4 E! P4 m* W\" q( f, G
ftemp=IntL_P(vlb,tempvub);& L! I# l, g2 G% q: w- ~( U
end