[问题求助] 请求一个整数规划的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)
6 B- C% H t* [& u
global upper opt c x0 A b Aeq beq ID options;
: H6 H* @, ?8 a! L4 B) C: D( v+ ~
if nargin<10,2 s: A\" ^$ q J$ e+ H% U3 t
options =optimset({});( r+ @\" K6 g B
options.Display='off';
a- A! x3 d0 h6 [6 d: K2 X! G
options.LargeScale='off';2 z& I6 x7 }) S. e; K' A3 t: V
end% w; w) U8 E+ `: a3 [5 Q
if nargin<9,
0 s2 _) b% k- }: g1 p
id=ones(size(f));
! I$ s1 T\" Z, j2 J+ y4 R
end - D: ]+ ]0 P+ G\" ~' k O1 F
if nargin<8,
2 I8 b5 O2 s' D% L* ]
x=[];
4 u+ x: b/ m! a6 T2 v
end$ R1 v# s5 ]5 W: j
if nargin<7 | isempty(ub),
7 S5 [/ G; b6 \# Q0 B/ Y. ~8 k, G7 [
ub =inf*ones(size(f));3 s' i8 B: l Q# P0 ~) B8 s
end2 f' P7 z0 w$ l- y1 i
if nargin<6 | isempty(lb),
: Z0 I$ I8 P) L4 K$ n+ U4 W% U
lb=zeros(size(f));
( v/ \$ ?5 Q\" f1 J
end
* Q$ Z* e: O. @/ N+ E$ p
if nargin<5,
2 Y* }\" f* C: z
heq=[];
2 J3 _3 G5 A, _2 \. ?
end, A5 ?0 z1 f1 \/ R
if nargin<4,
; c5 [/ W! Y d% p3 L/ k+ F' V- ]! [
Geq=[];7 o/ I; U% H+ K7 |1 { W( X7 z7 q) h$ J
end
8 Y\" h7 Z% d0 A8 l
upper=inf;c=f;x0=x;& i W0 O& B6 W2 ~
A=G;
, r& ?2 ?6 x# D9 l$ k4 ^\" \
b=h;Aeq=Geq;
4 R( f4 }/ g- }# b0 i& _, G
beq=heq;ID=id;& z: p5 m8 r9 X* `7 O
ftemp=IntL_P(lb(:),ub(:));
3 E6 S' m) u3 ~9 j J: w
%以下为子函数
# k9 m8 `- f9 t: U! `+ S; s! u
function ftemp=IntL_P(vlb,vub)
3 v5 |4 r) C5 i& T+ j, P% q
global upper opt c x0 A b Aeq beq ID options;( R6 o! X8 @* C
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);( i5 y. n* \6 ] R, V3 ]! [: ]8 q
if how<0
\" A( Z6 S; D1 t$ Q' N
return;
; I$ N2 Z: ?/ t8 S; Y/ r( J
end
x3 Y) K* N! T, v. a
if ftemp-upper>0.00005
- z\" O8 \) U$ w7 F3 }
return;
' `& G5 f1 f: C8 A$ P7 w
end
9 \) M$ ?+ N\" E4 r
if max(abs(x.*ID-round(x.*ID)))<0.00005
, j7 F3 E& I4 K6 @
if upper-ftemp>0.00005
) y& L6 ?& h$ ?/ Q1 a& ^: k4 |4 e
opt=x';
/ b- v. X3 z W+ f4 p$ z7 D% |
upper=ftemp;
- Z+ w# X2 T2 p! \3 _: O
return;
6 S6 g+ Z/ f' O7 F2 }9 \4 G
else
# ]5 @' `% I1 W6 F4 E8 T( O
opt=[opt;x'];7 @; Q. ?0 b\" V) I: G& W. t3 ]* n
return;
% e1 p Q( a1 x) \3 m; L6 U( O
end% D\" D0 T# D' `! f\" W
end( v& U3 a4 h. A) ]/ o\" a0 E; B
notintx=find(abs(x-round(x))<=0.00005);1 G, O3 n5 m T- S# Y- {# P
intx=fix(x);
V2 O/ L( v3 o
tempvlb=vlb;tempvub=vub;
) r: y. u; B2 _$ Y! Z8 B; y
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1
2 p% l\" ~, K4 z$ {% w\" F
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;
8 ]9 c3 m. m5 e c( s# l2 ?
ftemp=IntLP(tempvlb,vub);: q- i8 O) t6 a. h7 v k$ ^( \
end
. N, N; p5 `: }8 Q2 l1 d
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)5 A3 q% B0 q2 j: o6 c$ g' ^
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);
' U7 e3 a& R4 A. @
ftemp=IntL_P(vlb,tempvub);
\" F! m* s+ y% l. }7 c/ c. r
end