[问题求助] 请求一个整数规划的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)
* m/ A; @1 w4 |# V
global upper opt c x0 A b Aeq beq ID options;) b\" D$ r7 y* {# r# {2 b v J
if nargin<10,
3 A0 ?) B' ?7 k8 o, r
options =optimset({}); l0 L; p9 z; Z* o\" T6 m V
options.Display='off';
5 s# m* h: Z$ `5 ?6 H
options.LargeScale='off';
' J4 s3 Q+ p7 I, f( k) P9 G# n( t
end
/ r' P( r% D+ [1 M, I\" p\" L f
if nargin<9,$ \* h3 W# I: F, G7 g. n
id=ones(size(f));
& x\" _- ~7 A1 q U
end
9 f( o, N+ E. N9 x, P3 O. P$ e
if nargin<8,
! X: \( X* p9 H- i
x=[];
4 ^3 g5 Q0 y6 K/ j$ }% @\" N! b
end
4 }+ L+ }& y2 V- J' D\" W
if nargin<7 | isempty(ub),
9 j- c2 M$ B# a {9 V
ub =inf*ones(size(f));& b: l- Q Y, _1 u' p' G
end
, P! I+ V- r/ Z9 X) v- i/ c& e
if nargin<6 | isempty(lb),
8 I, \7 w/ ?9 W& ]4 v5 ~
lb=zeros(size(f));! e\" H C& a0 ]# b+ ]
end
& I1 c' V4 @; B* J3 P9 m7 B1 r
if nargin<5,
' z2 p: K- W6 F3 k- j9 g# V* |
heq=[];: T0 H\" c$ }4 R3 L, q4 q$ E% I
end0 W, ^/ `# e6 a! J4 d( W
if nargin<4,\" I; I/ w2 Q8 ]
Geq=[];
3 k4 ]9 a6 x( V0 e- n
end% q3 d: y& n- u, T6 N
upper=inf;c=f;x0=x;1 E8 f- R6 ^% W) W
A=G; Q! ~2 g0 l+ |, R
b=h;Aeq=Geq;
\" j; f) p+ J0 u3 `6 Z- G8 q
beq=heq;ID=id;5 I3 X7 ^$ F* d, j: W0 {3 }8 B
ftemp=IntL_P(lb(:),ub(:));, Z6 u& Z( A4 z; F5 w5 o
%以下为子函数
7 Q0 q+ s4 d2 d2 U
function ftemp=IntL_P(vlb,vub)
6 e' v8 H8 U+ w6 U
global upper opt c x0 A b Aeq beq ID options;
' X1 J: o- }; S7 C
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);
4 A1 ?1 N/ X# @) W) ~# j4 y p
if how<0: K# l: p\" j4 Z! T# H e
return;6 y6 P* h& ]$ v' B! y0 n6 G; f- n
end
2 U# A4 v1 }, t3 Y/ `: U
if ftemp-upper>0.00005
$ h/ M l9 ^% M& F, i
return;
$ N* I$ s) T+ h) i7 z
end
! s5 ? L; B3 h2 l% n; M
if max(abs(x.*ID-round(x.*ID)))<0.00005+ Z8 J# n; {: v# w% N X% ~
if upper-ftemp>0.00005+ m2 [+ i4 n% n1 z
opt=x';
, X\" G$ g7 @$ i4 Y$ i* G
upper=ftemp;
4 I: `# F& v$ p7 b3 A\" I
return;3 h, k5 }9 n# O. N' C\" z
else, o+ f: I1 r- @\" y, N
opt=[opt;x'];
$ B* S9 o- A8 _& G w3 A' D- J% ^
return;3 |9 ^- m/ r7 B. ]6 m, `/ Q
end
! v7 R- _4 g d8 k9 v
end6 _8 E. G) S( U' S
notintx=find(abs(x-round(x))<=0.00005);
2 I+ [, s' y9 T1 |) _7 o) s
intx=fix(x);
7 [2 B/ y; \# P. K2 L1 ?
tempvlb=vlb;tempvub=vub;
( f- H6 T) F2 p/ g( X5 M8 }1 D, p
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1
3 l4 [9 v\" b B. Q2 f
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;: Q3 a\" w; _' H% ?
ftemp=IntLP(tempvlb,vub);8 x4 ~\" S3 t* Z7 F
end
0 ]\" s0 t+ F% i
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)
( H: q, v3 @8 I3 L' n! g
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);
\" V, M' M\" k5 A/ A# q
ftemp=IntL_P(vlb,tempvub);2 }- H0 P' {) Q+ a4 j! I
end