[问题求助] 请求一个整数规划的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)
/ e% p* e* g, G4 Z8 E4 I) a# u. u9 _
global upper opt c x0 A b Aeq beq ID options;* ?! o5 c8 f- p- L
if nargin<10,
3 |4 C2 H$ ]\" ]* f7 R9 A
options =optimset({});3 k4 g* ?2 s2 S+ L2 s+ ~) e& _( H
options.Display='off';
1 Y! g. h( `\" t$ H& a) C5 z% J
options.LargeScale='off';9 s) E& i$ B7 U- V5 S1 S
end/ M* M' r c. B! j! L6 [
if nargin<9,5 [5 r; A4 {: c0 }( c' T. N* V
id=ones(size(f));
\" t7 _' ]# [\" W; c8 N
end 2 C' `- G9 T5 I$ E
if nargin<8,0 i, W\" C; P2 x/ ]% {: [
x=[];: E( a9 e+ F( f: N, N! {
end
+ E+ d* j' p3 z/ M* `5 n4 f' c3 o9 Q
if nargin<7 | isempty(ub),+ j) ]( H/ F: O( e0 F$ T) |
ub =inf*ones(size(f));\" E+ L1 T! R/ r2 u* O6 @
end
0 F8 w% X& R4 Q; I; u
if nargin<6 | isempty(lb),
1 P N: e3 Y* w a$ ?! ]. v\" O
lb=zeros(size(f));
5 z9 ^, i1 g) f1 [8 L
end
! a2 \# J9 Y* |6 V: _
if nargin<5,0 |5 x6 D! Q5 n e {
heq=[];\" w1 w\" o, S( x1 h: x3 _\" X
end
if nargin<4,: {\" r8 F; q4 R7 ]. \4 ?- ]
Geq=[];, W$ v+ E2 H$ D9 I& v3 ?# W8 d8 J
end4 S7 d3 s# d- r6 e1 P1 U
upper=inf;c=f;x0=x;
) ~0 o [4 t4 V' d$ ?$ o7 c/ l9 }
A=G;2 S7 P6 I9 N9 }6 y. w
b=h;Aeq=Geq;, g$ U- ~* w2 \) e) o
beq=heq;ID=id;
# _; u\" \2 k, m! S! s: x6 b( i4 Q- T
ftemp=IntL_P(lb(:),ub(:));
# H& B( p, v' y- m2 J7 T
%以下为子函数2 N7 s& V/ p% X S2 D1 D\" l
function ftemp=IntL_P(vlb,vub)
, O4 s: b, {% A7 t0 i0 J' B
global upper opt c x0 A b Aeq beq ID options;+ V' X! P4 G7 ~* \$ z
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);
3 ~! q( r0 y* z: E
if how<0
$ [( w9 G6 H/ H+ C- E
return;\" ^& w; F* s O0 z7 o3 {
end
* w1 j2 D! H& P' ]$ A\" O. Y$ {
if ftemp-upper>0.00005
; @5 N6 s# w% q' c- S
return;8 f. s8 q+ [0 _& X5 U4 I; \
end
4 l. k& M( H) g
if max(abs(x.*ID-round(x.*ID)))<0.000057 P/ ?2 [, ]7 B6 x. H- k
if upper-ftemp>0.00005
# n! W\" D7 @$ G0 V\" `; b \
opt=x';
1 `' m) }\" V6 q1 w4 E) |
upper=ftemp;
) ^! E6 b* c V/ U4 t! i
return;0 s1 h$ N, ?* A5 m& { a& r
else
- S/ e' N% Q/ |& Z4 w
opt=[opt;x'];
4 z/ F* f; {# D- a3 s% ]
return;, J* ~\" Q6 s! R% `3 y
end
' ]0 h4 {1 T! h: R3 ~% S
end4 H3 z5 K+ D+ C# \6 N6 g
notintx=find(abs(x-round(x))<=0.00005);
0 S\" p/ `# r3 ^1 e
intx=fix(x);1 l9 u8 |! k6 B$ ^4 q' T Q# Z
tempvlb=vlb;tempvub=vub;
8 o: J+ N# E } @1 @
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1
* ?+ n3 S; p, ]1 P
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;
1 m6 L( U. _! s4 G0 J
ftemp=IntLP(tempvlb,vub);
* j7 b* ~# l x7 \+ I\" {
end: P3 y' }+ p\" f# [6 R; h
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)
5 z7 q- D3 d8 v) t
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);5 o5 c$ o: W+ H( t7 e' U
ftemp=IntL_P(vlb,tempvub);: P# y. q+ E( v
end