[问题求助] 请求一个整数规划的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)
: v9 g) O9 N( `: y; [4 I\" A0 }
global upper opt c x0 A b Aeq beq ID options;0 `8 Y, i2 h6 M$ {7 a; H
if nargin<10,
9 m8 A* a4 ~; K6 ^9 Q0 T& ?
options =optimset({});, X) w1 E0 X2 R- D0 X+ z
options.Display='off';
/ Z% s8 X/ D# q) v# M/ x4 |
options.LargeScale='off';
: o; \+ E* U* n- R
end% L* T+ s& X3 W/ G# b
if nargin<9,
: J1 v! o, m) g
id=ones(size(f));) D\" j# X9 C+ M3 m# B# m6 _
end
& y \0 l! z; T8 `* k+ M1 A
if nargin<8,* J9 e! U/ u/ i, c4 ?1 U9 {: ~
x=[];
' Z5 K* R5 \7 v1 r3 M
end5 T9 y- K6 q0 o1 z
if nargin<7 | isempty(ub),6 }5 Y/ i# J( `2 \
ub =inf*ones(size(f));
0 _+ h/ F7 A2 [4 t
end% R8 F( z3 ~. R- k8 r8 w2 Z
if nargin<6 | isempty(lb),
' h# A+ M# b6 I9 m# L. D+ F9 j8 s\" D
lb=zeros(size(f));8 ~6 ~3 }1 `4 O. e
end1 z: Q0 c* N7 r6 A
if nargin<5,
, q: M! R N! a
heq=[];) {; D2 ]) M. J3 `4 f# s
end
& r$ w4 h* y2 I% k* Z6 T. A
if nargin<4,4 _9 y2 @! {- @, `8 u
Geq=[];
& |) I' a3 }\" f) p1 I
end
+ z+ Q0 V' N1 a, V\" B4 n2 e
upper=inf;c=f;x0=x;
8 o& {; B+ f- ~1 a! F9 S9 A' n
A=G;
5 H% w/ S* H. J2 b# L
b=h;Aeq=Geq;; Y( i8 o4 Z- r6 ]% K
beq=heq;ID=id;
5 w, S) n* x4 ]( B
ftemp=IntL_P(lb(:),ub(:));/ J5 @/ w7 \* ~# Z% j$ Y$ S
%以下为子函数5 `8 K. d& ]1 [! ?
function ftemp=IntL_P(vlb,vub)
5 y+ C3 ]8 ~0 b5 _6 H\" X5 ]\" ]
global upper opt c x0 A b Aeq beq ID options;
0 L7 n0 a( B, c& V/ ^' m\" }: i) L
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);, b9 q; l# h; q+ |
if how<0& ]# X' _& Z5 Z6 R
return;# G$ s5 u7 i( T* a
end
% K2 J0 A3 O, S9 M$ e! s% S1 B1 ~1 d
if ftemp-upper>0.00005
, B0 Y3 q7 r( b4 R\" L0 ^- W8 ]
return;
# K; A- B I1 O& R B+ h
end. y. o/ U0 ]1 b( g
if max(abs(x.*ID-round(x.*ID)))<0.00005+ v2 U- y7 k) N; u2 |
if upper-ftemp>0.00005& A5 p/ c2 [* j, a& n9 X) x
opt=x';! w( h, e! ]& r' @. g
upper=ftemp;$ B( l, c2 Q( {4 T5 q) G1 `/ A
return;. |, l8 O. e# s& ^6 ^' B
else) X+ J1 O8 o- M) Z \# b1 U: a4 }3 g
opt=[opt;x'];
# Q; s* c\" ]; |9 B5 {1 `+ w
return;0 T0 V9 B( F( _9 W2 Z; z3 {( B' @
end+ y, P; W7 m+ f) _
end
1 a) f- T! X+ h& O( n5 @\" ~, r
notintx=find(abs(x-round(x))<=0.00005);; a8 z- ?4 U9 ` X4 y
intx=fix(x);9 S u- m! W1 [0 d
tempvlb=vlb;tempvub=vub;* g# ^3 @7 |$ J9 {1 a$ `\" W
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+19 d% Z: u1 B$ N$ [
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;
! V+ V i3 X) l6 r E
ftemp=IntLP(tempvlb,vub);- V6 k+ _% C6 e; L
end
, D) ]& h: |- ~# n# F; G
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)
# }9 W. |$ o4 [! `
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);
( n( Z8 v) G+ r& S. P\" z9 ?& j
ftemp=IntL_P(vlb,tempvub);' l\" u: i( g( m9 z) L
end