[问题求助] 请求一个整数规划的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 D) ~7 L Y$ m- L- Z8 C
global upper opt c x0 A b Aeq beq ID options;5 R! a5 V' P3 E; Z q7 |
if nargin<10,
2 Z' [) I+ F6 [
options =optimset({});
& ], d8 D3 h3 a4 c
options.Display='off';* @% |\" a- e3 X
options.LargeScale='off';
; D% ?. [4 L- g' {2 G
end& f* C! T( k% y' K- n& R: H, c H
if nargin<9,' ]' U) y7 {. [% J3 G4 |
id=ones(size(f));+ R: X; H, O# g; f4 N% C
end
V; r$ N% M p
if nargin<8,
7 I6 }+ b, q- D
x=[];
( C\" N7 Q; F* S& M& D2 X+ @- F; }
end% y7 f# n' z8 [3 Y# }0 ^3 x\" r. d
if nargin<7 | isempty(ub),
, S0 ?& y% h2 f, K6 o
ub =inf*ones(size(f));
6 f5 s3 y6 D4 e3 H
end
9 B6 K; q* S, z# x\" ?2 c9 W
if nargin<6 | isempty(lb),
' m) |. i0 ~7 e! b/ u9 ~. Y
lb=zeros(size(f));. x- x2 D& j9 j# ~
end
0 S) Y) O' Z' P$ j6 V1 v
if nargin<5,7 {# \\" R( M( q* b6 X9 o
heq=[];, i# J$ t\" X+ `( n9 A- E) W
end
4 L. g V+ w+ E$ o! F
if nargin<4,+ i' l' Q2 m- ]) K: s4 f\" d
Geq=[];
* ?( ^# P3 v0 K. J2 S9 T* j
end\" D; r( x, m+ f$ a) z1 c. Q. S
upper=inf;c=f;x0=x;! e, f, l% ?! p! p
A=G;
$ i9 E- H3 I5 G; P) o& f
b=h;Aeq=Geq;
( N4 f% k q\" E' Y ~
beq=heq;ID=id;
\" R2 | B( Z4 Z9 k! n6 M
ftemp=IntL_P(lb(:),ub(:));
( O9 R' i9 V+ }& ~, w9 _6 W
%以下为子函数
function ftemp=IntL_P(vlb,vub)
0 D\" d! s4 L2 \* Y' {( G6 g
global upper opt c x0 A b Aeq beq ID options;' e: K7 Y5 l% S/ o- A
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);. G _3 L/ w) |- l V% e U5 T
if how<0+ G, j/ q* d2 S1 q
return;
9 y5 h, }$ m- h/ e+ r
end
9 G; t, N+ F- W$ a( T
if ftemp-upper>0.000055 k0 U( o( C+ p4 |4 V7 {( V
return;
% A2 x. [, v6 i6 j
end
! M+ S/ P. Z) Z L o- z' D4 }
if max(abs(x.*ID-round(x.*ID)))<0.00005
% d( d9 t* U+ w\" ]
if upper-ftemp>0.00005( i7 M$ p4 \3 q7 W# f& _
opt=x';5 B3 E% | \5 b( V- U
upper=ftemp;6 ?8 G9 \& @* V
return;
' l s* X4 g3 K% ?! v
else
0 D9 \8 {% u; y' `! G
opt=[opt;x'];
7 t2 P8 `& {3 X+ T* v, g- U+ j. U
return;
; H. |5 g# ~) j: |% Z
end
9 c& W. R# [\" Z0 c1 k6 X' i9 K
end, [5 ]6 M0 J y( z
notintx=find(abs(x-round(x))<=0.00005);; g0 A0 T4 X0 N7 f ~\" P& L0 e
intx=fix(x);
9 C7 X8 I) z' e
tempvlb=vlb;tempvub=vub;2 b; M3 D3 N: [& K. r1 `' B
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1# {3 R: D5 r4 Z\" u5 {
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;
\" d* K, \: [# z4 b- }: a
ftemp=IntLP(tempvlb,vub);\" G+ _\" F1 q5 B
end T' a' `' b; r
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)& {0 T- Q9 W5 t1 F2 U' K
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);; ?0 S$ a( h, Y# |# T7 u% d
ftemp=IntL_P(vlb,tempvub);( L* ~. l9 U1 _! [
end