[问题求助] 请求一个整数规划的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)
# x. `) t a- E
global upper opt c x0 A b Aeq beq ID options;
1 u& q/ C& K6 P$ O/ N h
if nargin<10,
7 i0 c& f9 [! C0 f. h5 Y
options =optimset({});
4 Y: C\" Y& d3 j% s: v- e
options.Display='off';9 E- l; H( h( W& m% U- _
options.LargeScale='off';
6 ]3 N; H C( c3 z2 H) W* _
end( R8 U# r! t+ ^- e) F
if nargin<9,
. j9 A# m% |! t& v
id=ones(size(f));6 t4 u! D$ O/ x; b+ C$ O
end 8 ]! R$ Z. |. F' p- S( a
if nargin<8,
7 J: |& s$ `* A
x=[];6 z3 m1 `! |$ p/ w+ Z
end
) i+ E& S }, A n X7 l4 n9 U
if nargin<7 | isempty(ub),% b* g! @( d+ r- W
ub =inf*ones(size(f));4 n+ j: |% p. O: c- E
end8 Y' C3 D$ _4 G, u
if nargin<6 | isempty(lb),
0 v% n$ H4 F( C! Y! T
lb=zeros(size(f));
7 j% w. N6 F: O1 c* V8 O
end
. z; v8 \4 N9 X\" \
if nargin<5,
) ?5 ?# C6 B4 n5 }, t
heq=[];8 S9 k+ t/ d! } Z* h; I9 e( u
end
8 l5 s9 r0 _3 D* n' \( W
if nargin<4,: S' S7 B* L: X1 n) m
Geq=[];
- I; R8 h. J) U# I: |
end' F( g* r\" Z5 G/ v6 X
upper=inf;c=f;x0=x;
6 z- I; V; [4 \( q* J7 y8 f
A=G;\" p4 H% j/ P/ A+ _, l8 z. G. H( Y
b=h;Aeq=Geq;, L\" |: @! F$ {5 h- |5 V% q+ Y
beq=heq;ID=id;
: q0 f4 X3 }6 y$ d9 D+ n
ftemp=IntL_P(lb(:),ub(:));( Z* Z' u# b4 w6 N) E
%以下为子函数! A0 j5 l2 q: K, `' X% B
function ftemp=IntL_P(vlb,vub)
' Q9 P4 Z0 v G+ j
global upper opt c x0 A b Aeq beq ID options;
# D2 H# v6 R9 w; y
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);. g2 I6 b) R v( r% x
if how<03 }\" m9 Q( ~4 o' B
return;
1 y# @1 Y6 Z2 ^; p8 Y0 P7 n
end. b, a- V8 H( H) `2 _
if ftemp-upper>0.00005* {3 \! q7 L) x$ p* y, j3 F
return;) \( f4 Y3 P& l8 g) D
end& B$ G1 b0 c$ J/ e3 {
if max(abs(x.*ID-round(x.*ID)))<0.00005
6 Z% j5 F5 O% e0 j+ T; I
if upper-ftemp>0.000056 t+ m5 a) h% x+ j) N0 d
opt=x';
0 o1 X, T\" O7 ?9 v2 E
upper=ftemp;
% K$ q! x1 O7 W! S/ S7 }
return;, p: q) u+ K6 W& p. A6 g }
else7 {4 i d( e |- _) E1 x& D
opt=[opt;x'];5 [ X6 ^9 g/ e! h
return;9 S7 C; b$ C) T! k( R. f
end' k |/ b\" [2 k
end
: R$ M, {0 z, y
notintx=find(abs(x-round(x))<=0.00005);! I9 e7 h3 K4 L$ G k
intx=fix(x);
. |- g* Z6 w2 {/ S, I4 g/ `& A
tempvlb=vlb;tempvub=vub;' E( n8 b: n9 M
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1
, Z, f; H9 e/ J) T
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;* N\" i) R0 C5 I5 P
ftemp=IntLP(tempvlb,vub);3 J$ u1 B r$ A! w% n/ r
end) Z+ L6 W/ w4 z u2 a, ~! |2 {
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)
- c# j. Z$ I& {. P5 K! y3 ~$ t1 O1 j% Y
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);
' u$ R( o8 Y, L+ n4 I
ftemp=IntL_P(vlb,tempvub);
9 ]. g; f* N/ v
end