[问题求助] 请求一个整数规划的matlab程序

[复制链接]

字体大小: 正常放大

huiwen

4 主题	3 听众	258 积分

升级 79%

TA的每日心情

	开心 2014-12-3 16:59

签到天数: 84 天

[LV.6]常住居民II

群组: Matlab讨论组

群组: 西安交大数学建模

群组: 2011年第一期数学建模

电梯直达

发表于 2011-8-25 16:51 |只看该作者 |正序浏览

|招呼Ta 关注Ta

最近一直在研究整数规划，可是连着看了两个相关编程，用MATLAB实践时总是出错！所以，求助于各位大侠们能否提供一个整数规划的程序。或者，能否分析一下面的程序哪里有问题？谢谢了啊!

function [x,y]=IntLp(f,G,h,Geq,heq,lb,ub,x,id,options); M; O s+ ^\" @1 s0 X
global upper opt c x0 A b Aeq beq ID options;
, o! {, L; \7 Z
if nargin<10,) U; B) h6 ^+ ]4 i
options =optimset({});
\" k: l/ I* i' b5 b; i7 t
options.Display='off';# S& g% ] |1 F; o J+ \/ X
options.LargeScale='off';- C0 [( s) l4 t/ C5 t1 u
end2 I4 B l4 v* l\" J/ h$ W) Y
if nargin<9,2 |9 S0 `3 P! J+ I7 _\" d+ p\" S* E5 e
id=ones(size(f));' u9 o! P9 ^# Q j
end % |: d( K- Z# o0 K/ b
if nargin<8,
) ?\" X' d$ G& f
x=[];\" e9 A$ K0 N8 C; A: ~8 T/ W
end
4 H! B3 ^0 c5 u5 G
if nargin<7 | isempty(ub),) h/ |3 J. Z9 Q. U6 c; c' n$ [
ub =inf*ones(size(f));
' u5 o* J6 \! d
end
9 ]( Q8 e; H6 Z+ e/ C
if nargin<6 | isempty(lb),
% g9 u' t; Z1 p2 T
lb=zeros(size(f)); Z: j2 p& T5 L' p$ \. P+ I
end! n* B# g( Y4 a+ |) B+ n: B4 _
if nargin<5,/ K: o% y2 U) M& |6 c1 ~# l
heq=[];
! R) u8 i: z: l3 C' @
end\" b l' I0 j\" C, K
if nargin<4,$ i: |0 Z3 W- q9 i
Geq=[];/ B$ n7 z( `2 P% V; a- W, G2 o) B5 Z
end
5 w# y4 m, ~) A/ P
upper=inf;c=f;x0=x;! C) Q$ v7 h3 U) a2 ]4 U9 j
A=G;0 g/ {) d9 i' d7 \- [6 _
b=h;Aeq=Geq;+ X0 ]7 q6 [( k
beq=heq;ID=id;- [0 O+ Y: E5 B+ C\" O# Q) a2 W% q/ N
ftemp=IntL_P(lb(:),ub(:));
; Y: u! g1 ^. d. C$ o
%以下为子函数 [5 A* ^$ h* l) t4 M* u! k
function ftemp=IntL_P(vlb,vub)# I# |1 S; d3 q3 V% W
global upper opt c x0 A b Aeq beq ID options;
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);
8 l' w\" }$ D# b. k0 J5 ?
if how<0% }0 v- O( U4 _' E# U$ [! j% }$ S1 |
return;
6 f4 w6 u* m) J- {& z1 Y
end9 s\" w3 y O4 E
if ftemp-upper>0.00005
! Z) s' V) Q0 N5 c' H
return;
9 Q( {- W* d- F1 K4 c) c( o. j: E
end- T' R' a+ h& j( I7 Z
if max(abs(x.*ID-round(x.*ID)))<0.00005 g% E# s( E2 `- G1 s v
if upper-ftemp>0.00005
: x4 \: G( p. s' d
opt=x';( S; V+ w( g) t6 r8 Z
upper=ftemp;! ^5 m% U2 v0 P7 F( p
return;$ [. h% Y% T# _1 U& q
else! `# Y. y8 A8 X+ D& d
opt=[opt;x'];
* {2 K1 s5 S3 R' p/ x
return;$ ^% a: z3 X% d1 a; I
end
4 N+ ^$ \7 f) ?2 _3 z\" q
end
0 p/ J d+ ?$ ]8 M9 N
notintx=find(abs(x-round(x))<=0.00005);
6 d$ w% i' R4 A( }4 t
intx=fix(x);\" N( z/ H. r6 [4 x ^
tempvlb=vlb;tempvub=vub;4 q, f& S( z0 S- v! f+ |; _
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1/ e4 F+ [8 M/ P
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;; v1 V2 |) a; m4 f
ftemp=IntLP(tempvlb,vub); _, Q& H9 E' `# s* Y
end3 [6 z$ M4 t9 @! d; h1 F6 l) z
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1) z: ~6 w& H& ~. Q
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);
1 D! C+ { y& _( C, c
ftemp=IntL_P(vlb,tempvub);) p$ |$ }& b# Y, ]' ^ N0 L
end