[问题求助] 请求一个整数规划的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)
. s5 D# U+ r6 Q) j
global upper opt c x0 A b Aeq beq ID options;
2 h1 T0 J& }* L1 L
if nargin<10, l! Z1 z( t\" `) S, H
options =optimset({});2 p. l5 [) d2 U: g. y
options.Display='off';& e: W1 L z# \
options.LargeScale='off';
# W2 I* ~/ F+ k9 z9 A# x2 f0 r. }9 p( g
end
; H' h9 Q! t2 H5 A9 N
if nargin<9,
' b) H0 b* K2 n* G* ?2 O
id=ones(size(f));/ _) x9 d- @$ F4 b$ E3 i/ I, l
end
# {$ v) b$ n: b: b! H
if nargin<8,' L8 J4 j$ a0 ~) O4 c+ e/ a- t) x
x=[];# P- }8 D% J: m0 E
end, Z, m0 g6 q; z) U3 r# K6 w$ c- ~- ~
if nargin<7 | isempty(ub),* X1 ~3 a4 [5 v# w
ub =inf*ones(size(f));
' l7 T$ ]$ {0 {4 [/ L6 d- {, G! o, A
end: v& I) W* u4 l. y. A
if nargin<6 | isempty(lb),+ U0 [4 j! W6 k\" ~6 J7 m
lb=zeros(size(f));
) h/ t8 C/ n* l3 G: ]4 ?
end
) v9 R' \8 t0 Q) l% d1 O$ ]/ u
if nargin<5,! C* ~; A1 m; Z4 y h
heq=[];9 H\" e6 v5 y) z5 ]' {( Z- e. L
end
$ D: L% j1 Z% _0 |$ I3 d( y
if nargin<4,8 n I+ a l: X7 E0 c4 m+ y
Geq=[];
8 m+ g2 G& n6 E: B+ K
end
5 J2 y$ ~8 ~, y8 ~- w( ]9 S
upper=inf;c=f;x0=x;
\" D% p+ m) f( ~* l9 r
A=G;4 w. \0 ?' @% j/ {7 f& [ u5 e
b=h;Aeq=Geq;% d! Y# x, ?+ a: ]+ C* {5 q8 t
beq=heq;ID=id;3 \0 L: h. w9 t5 ?9 V
ftemp=IntL_P(lb(:),ub(:));
k8 O8 D* i! k- H, l6 g
%以下为子函数$ @1 J! A$ g; A- L
function ftemp=IntL_P(vlb,vub) H2 j% U. ]: ]3 `
global upper opt c x0 A b Aeq beq ID options;2 o: T$ ]# g( C) V
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);+ [7 L6 T) K' z0 @2 w8 m# C8 ~
if how<0
& h# p1 v3 P. L( x# I6 I0 K\" a. }5 y
return;/ D4 Q2 m% d5 _) T& @
end0 E6 [% u2 H. k6 f/ T
if ftemp-upper>0.000052 g2 @2 {8 O) @9 P9 |
return;
7 y1 t r! y3 a
end
, M& i: \5 S I V6 w
if max(abs(x.*ID-round(x.*ID)))<0.00005
6 R7 d3 |, i! w; v
if upper-ftemp>0.00005
# r& m$ Q! O+ d1 Q' p
opt=x';
* s0 h( ?2 D; ~+ L3 \9 ^7 @
upper=ftemp;
: v+ e- y- {/ X% u) ]
return;& J4 Q0 o1 w% d& c+ l% r* U8 R
else
( C: r/ S9 D: f/ j5 I& c
opt=[opt;x'];
, ?; S* V$ w8 Z% b
return;
& A# H C9 u& _* ]% g R/ Q
end- q* q# k: {0 Y) B. l
end
' H5 P3 ]; D* z7 R0 b
notintx=find(abs(x-round(x))<=0.00005); C; |1 u\" j( @' }& _
intx=fix(x);& M) w3 B/ l# L# [
tempvlb=vlb;tempvub=vub;0 U! f6 G( e8 Q }) j
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1\" _! S. b. V) Y5 D& @* @: i6 @4 X
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;
\" W. h4 j3 B- W6 ^% ]* K
ftemp=IntLP(tempvlb,vub);
; z- I# {0 Z) j
end
! a3 V0 {% [1 w
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)4 z' G7 p/ \+ I6 G( Q
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);
0 b# J; I! ]8 I; d8 z1 k9 {
ftemp=IntL_P(vlb,tempvub);3 k; m, y; b/ j+ Y* S
end