[问题求助] 请求一个整数规划的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)
) q8 X5 K6 F: `+ e* }2 H9 p
global upper opt c x0 A b Aeq beq ID options;$ p% T% q0 L9 o. l\" G9 p# ^3 h
if nargin<10,
' N, o9 [% F\" i1 x\" n\" E! e
options =optimset({});4 I; g\" D/ y7 e* B3 E) {
options.Display='off';
1 i7 r! m- r9 A* i0 O
options.LargeScale='off';- E3 D. t0 x4 g* f\" I
end
1 v+ k2 ^; e' o% ~
if nargin<9,
1 `! Q! L7 S4 A2 y! T2 J
id=ones(size(f));
8 c$ V8 F1 @: c\" }6 ^
end 5 o: P; _( l0 X- X5 N4 X; x C
if nargin<8,
' T$ i- U; n% ^# D2 ?
x=[];& ?: `( z, D C1 n- ^8 s
end- o9 E! O( F- a9 E( D& G
if nargin<7 | isempty(ub),
) G: P% K* U p
ub =inf*ones(size(f));
; h+ E: w4 I9 ~/ d
end
* H0 M. g( |' D8 g# u' ~, u
if nargin<6 | isempty(lb),
/ K9 n5 _, L5 j3 b7 n6 I0 y& m
lb=zeros(size(f));
0 ^, K+ k. k( @4 g
end& Q: j# N. g0 U# E+ L Q$ F% N
if nargin<5,: g2 |* @1 K6 ~, k9 a0 d% b* e
heq=[];
9 f8 |4 y+ ]+ J {8 }
end
1 H% Q* q n% L$ a\" N* O( [
if nargin<4,
1 q' j+ y6 U6 g0 ]
Geq=[];* f) S1 Z* t7 f* L/ B9 \
end5 r- D+ C5 B\" G
upper=inf;c=f;x0=x;3 E S3 _4 s' s, L& O
A=G;/ J\" e2 Z5 N8 C3 a
b=h;Aeq=Geq;3 u8 r* a c: i2 ?
beq=heq;ID=id;
e$ h2 E; U. ^0 v
ftemp=IntL_P(lb(:),ub(:));3 u3 {- H G2 o8 N3 \0 r8 g: {% T% V
%以下为子函数
1 {8 r+ C( M% R
function ftemp=IntL_P(vlb,vub)! u p2 s+ }4 V7 q9 k
global upper opt c x0 A b Aeq beq ID options;! N, P; I3 p# N/ D
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);
- W& f, Q# i' M3 p' c7 @
if how<0
, H\" c% p! j! j
return;
9 @\" A/ p$ ?& `, [. j
end1 s O8 F8 v) A, L
if ftemp-upper>0.00005
9 h: q\" z5 o: m: v
return;) ^ J- V4 _, s; L# P& C
end5 \ y, T' p/ V; V# F
if max(abs(x.*ID-round(x.*ID)))<0.00005- x$ ?2 B$ ^2 @+ v# u
if upper-ftemp>0.000055 L7 P9 R% |7 ?+ K/ N! |
opt=x'; ~9 F( F7 U: B8 C& V# p1 C4 K
upper=ftemp;. x\" K4 g$ \5 M, E
return;; W. |0 X+ v- }! E0 r' P) ^6 y$ E
else& G5 z2 M+ W, H
opt=[opt;x'];
% S& P2 Q8 d/ g2 u\" ^2 w$ }2 o
return;
% l7 W. b/ h$ C4 e4 k6 P
end
2 H; B, j5 u6 ?8 [7 j$ a1 |9 ~
end: A# Q, c& a0 a' n! ~& |
notintx=find(abs(x-round(x))<=0.00005);8 o$ _! G/ O' K+ x5 q
intx=fix(x);. ]7 J; z, U! s4 J: H3 |
tempvlb=vlb;tempvub=vub;* K# h5 G: i! T$ @\" L2 T
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1
6 a9 v+ ^) Q( H
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;4 r4 J\" V' H3 j# P\" ~
ftemp=IntLP(tempvlb,vub);
( y( K' v. g; A& ~8 a
end( q5 K7 C' G- x* t) a+ V! t
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)5 M8 Q; L% B/ w5 [; k6 b
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);9 d) L2 K( b S2 n
ftemp=IntL_P(vlb,tempvub);# S1 ^( D& A3 q) }; v3 j1 f
end