[问题求助] 请求一个整数规划的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)# j9 M5 @ c1 I& a& O3 @
global upper opt c x0 A b Aeq beq ID options;
. o7 i* }' W! k* V/ d' F
if nargin<10,1 V) e3 q' C4 W) h
options =optimset({});\" _7 W% [4 H8 N5 @) j+ O
options.Display='off';# T% f5 ]- m6 [8 Q. V- p
options.LargeScale='off';) l E, k1 h/ r+ w! j, R$ U5 i
end: j4 v7 r/ X' m& Z ]+ _
if nargin<9,
9 t$ M# J1 b\" J+ p
id=ones(size(f));
# ^' \; m5 L: O* r
end ) ~' r2 a6 X: C, Z% N! R
if nargin<8,
8 T) |. I6 h8 z- |9 R3 k4 l
x=[];+ G1 U% B# ?* k9 h9 E
end( I1 M/ c! E5 y8 i7 X
if nargin<7 | isempty(ub),
6 L\" J$ q$ b1 O, {
ub =inf*ones(size(f));2 I& I' U5 A: P
end, o( n; n S* Z
if nargin<6 | isempty(lb),0 Z% Q2 f: z8 q1 d
lb=zeros(size(f));
+ E9 V/ d: `* d/ k* M. _/ ~: n
end3 n4 w- v F0 F) F# D
if nargin<5,! K. d& `8 C$ ], G/ E
heq=[];
9 T$ P$ ]4 k$ L& t$ E# P
end/ ~: c' k$ W' T3 x0 L+ ?0 p' X
if nargin<4,' C1 R- R9 T0 l0 G0 @
Geq=[];% U% n- u+ l* @. M
end: [0 p' W: f+ |0 c5 ^, f4 s
upper=inf;c=f;x0=x;
2 y: j1 I: F* [
A=G;
9 O1 v' [, r9 v y. ^
b=h;Aeq=Geq;
$ y2 \4 Q) L, H0 }
beq=heq;ID=id;2 c1 I$ R0 U) H
ftemp=IntL_P(lb(:),ub(:));
; o) m# a/ s6 a$ v6 G3 y% f( ~
%以下为子函数
7 P# z' ?: @0 h+ s& y3 P
function ftemp=IntL_P(vlb,vub)) a8 [( v7 ], {- P f& U9 \/ F$ u9 S
global upper opt c x0 A b Aeq beq ID options;
; t4 D& Q- J' @9 ~$ ?0 \
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);$ Z4 @. o3 L, U( O. \7 w0 N
if how<0
0 S5 K7 |3 {+ g
return;4 n$ c9 B5 C* \) u; [; s' n* _
end, U# a$ M+ r) J' h
if ftemp-upper>0.00005
' Q z& R9 x' u\" V5 N& y, u, A
return;' h: x( y: o( q! ]9 c. ` C2 q; R
end
0 i) L$ Z F/ P\" @) Q) U$ L3 E% r
if max(abs(x.*ID-round(x.*ID)))<0.000051 [: T! T2 A\" R1 A
if upper-ftemp>0.000056 g' x$ u( `( \! o
opt=x';8 r; R0 H# p. H9 q, w- @\" X1 n
upper=ftemp;$ Y+ {+ a/ U9 ^% a0 H\" U
return;
% e& c( F& }+ O0 M* f( q+ Y
else; z {1 F4 c' S) B6 V6 i
opt=[opt;x'];
$ \4 X; u8 a: h) x0 T4 [/ o
return;
0 P% y$ Y( {8 F: P
end4 x& s' P4 ~7 O d1 Y$ V
end
9 J6 q) g Y$ U% s, B
notintx=find(abs(x-round(x))<=0.00005);. s% g4 A. w, |, }7 k0 K
intx=fix(x);- e5 D8 E5 H4 P
tempvlb=vlb;tempvub=vub;
! s( O/ {; |: q/ w
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1
3 D- i) W8 M0 S! i- l& d9 Q0 Z2 F! o
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;( s: I9 C4 e) K! m# L& x
ftemp=IntLP(tempvlb,vub);3 b' K3 v5 J3 w\" u$ a( X1 v
end: ~5 B; D, Y [
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)4 e, E( b+ V; I4 a
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);
0 D# M. [$ v a9 b. R; T4 N$ F
ftemp=IntL_P(vlb,tempvub);: x( q\" S M' z) f. a% k9 m7 o% B
end