[问题求助] 请求一个整数规划的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)
0 r! [, g/ ^- P# Q5 z
global upper opt c x0 A b Aeq beq ID options;
+ d+ D1 w0 F. Q( M. ~$ ]. b
if nargin<10,3 Y# k5 Y* V: N4 B- V
options =optimset({});
4 t8 j, Q3 i% _- e7 x
options.Display='off';
* b# u% k- j$ i: K% q) l
options.LargeScale='off';
* L1 j$ h B6 d' o- j- F7 G' `
end% n6 z+ [4 j$ u% J# g2 `
if nargin<9,( G- c0 R' G' ^, ~$ L9 D0 a
id=ones(size(f));9 \) \\" h* B) [1 A% l
end & f; `5 ]3 P/ i$ x! n2 E8 n5 }3 B
if nargin<8,
( r! m. C0 O7 y; b: K9 l
x=[];
1 Q$ a( N( P4 z9 u9 ?4 p
end
9 c I3 a8 y+ T9 n( c7 r
if nargin<7 | isempty(ub),\" d9 y) J/ v, P d0 C4 o
ub =inf*ones(size(f));
2 l5 J4 ^5 H\" e
end
2 C# f4 C$ I4 [$ V\" G
if nargin<6 | isempty(lb),
$ F\" b\" \8 B) b, q\" m/ Z
lb=zeros(size(f));\" w6 W! ]7 B\" \! p/ [( f1 v/ j4 t
end7 F) g$ H; \5 \
if nargin<5,
# d* A( |# E$ c4 B8 k
heq=[];2 C/ H5 w, A3 q$ ]9 C, g6 p. [
end
7 V7 Z- w1 v6 I7 {+ ~
if nargin<4,
- n# l6 B( l! j- ` W1 W
Geq=[];
1 E5 ^0 K! b( p1 n! s7 L
end
5 u7 X) N- x A; [) @, v
upper=inf;c=f;x0=x;
, D# j& `7 ?7 o+ _8 S0 }' B- N
A=G;
1 ^& }5 {6 G: `% A3 [3 V\" W
b=h;Aeq=Geq;
8 k( v: |6 V B' U
beq=heq;ID=id;+ e& u8 t. l1 j( }2 K
ftemp=IntL_P(lb(:),ub(:));\" i# M, V! v, [# k- U' h
%以下为子函数
B$ Z# T* M8 G* C0 `6 h. F( r1 Z
function ftemp=IntL_P(vlb,vub)
# F+ l) o; _4 P' B1 g\" w
global upper opt c x0 A b Aeq beq ID options;* l8 ^1 `! Z/ y
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);
4 T' i' X' {, Q& C h+ m7 ] p9 `
if how<0
3 F5 a, K( u5 h% x1 a
return;
. w& N, y8 V2 s
end* {) K# f. c\" c
if ftemp-upper>0.00005
; E7 |; c/ Q\" v! j- V5 v
return;; i p$ u* M3 h, H# l. |
end
7 k; f, K6 B6 E) u: M
if max(abs(x.*ID-round(x.*ID)))<0.00005
- G8 p, Y% C2 f6 r
if upper-ftemp>0.00005
$ I+ W) F6 c) v- Y6 s! v
opt=x';
8 w3 F- \! M0 x& V- ~% s+ T) Y8 ~3 J; C
upper=ftemp;5 d2 Y( e9 n\" P
return;
5 L& B* d1 B3 Z$ t& h
else
& G4 A% i3 G+ E, T3 h+ i z
opt=[opt;x'];
7 _2 h5 ?! B7 U% j7 ~
return;- z7 T' _. {* x6 Y: F L
end4 ?9 j$ a2 Y2 R+ B- Z- q
end: P8 |/ |, Y* z9 m0 K+ g
notintx=find(abs(x-round(x))<=0.00005);- n' f* {( ~$ a9 r; J3 A S* e5 ?0 N+ Y9 ^
intx=fix(x);5 j; ?, d* X( T6 v# a3 y
tempvlb=vlb;tempvub=vub;3 `7 \8 S; ?6 ]9 s$ d! u
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+19 r! f8 r! J# B C: b
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;8 L$ u7 T7 W. k' V8 r% {! i3 z5 h# r
ftemp=IntLP(tempvlb,vub);2 j& Q r* Y) n' W& b
end
0 X- `1 e\" j\" g2 `1 Y+ t* u' q
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1) x7 X% o# j' g: B+ m: Q
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);+ i8 N7 M X, ^- ]3 i5 g
ftemp=IntL_P(vlb,tempvub);8 G! g% p% L# V
end