[问题求助] 请求一个整数规划的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 V2 X) U! g g
global upper opt c x0 A b Aeq beq ID options;' E* p7 p! A* {. |- l. x
if nargin<10,
. S9 Z% l! F7 D6 G; @) s
options =optimset({});
, N0 M% G9 A( C0 j: Q8 G
options.Display='off';
* e' A p& z9 v2 ~2 q* O9 K! Z7 z
options.LargeScale='off';. Z$ z\" F2 E* O, g! g3 e
end# w$ p, H) ^6 K& g- q
if nargin<9,$ i: i3 q2 ~/ l/ @. i* g+ [% g
id=ones(size(f));
- @$ H+ q E4 y( s+ O( [
end
' d7 y) p! H: Z\" Y- t/ P
if nargin<8,
# \- d5 y& j( r' D( W/ s1 p8 [
x=[];4 I2 X7 X% P a! }# o1 z' |
end
. v5 o# {/ P/ j# I& R' ?# ]( o+ Z
if nargin<7 | isempty(ub),
6 g- Q U& b! V4 s: Y# N
ub =inf*ones(size(f)); i( t4 w( F\" l7 y9 B4 U3 O
end* S! u& C9 u/ O# _
if nargin<6 | isempty(lb),/ J+ ~( {7 u8 Y: |
lb=zeros(size(f));
3 g2 x) L1 L& k! E& l2 G* x
end, X3 ?( _' J) G8 o* ]3 w5 Y' C
if nargin<5,\" V4 R4 u' C\" o8 \% N# c
heq=[];4 G4 A/ q0 U, N7 |5 r$ X- I' `% _
end2 ^9 q+ a5 `: S6 B. G% W3 d
if nargin<4,
4 g4 n) E& q& O8 Y9 Z V% p9 u
Geq=[];. M! B. ~% e8 W; W) E8 I
end6 }' Q! X4 Y. x- K
upper=inf;c=f;x0=x;
j3 [1 U( Y& @
A=G;\" m& h. T2 n7 k) ?! ^. S2 L
b=h;Aeq=Geq;8 ^# M9 a6 j# [$ d# a! \
beq=heq;ID=id;5 ^! d# e6 _1 l L1 n# D4 }0 Z& B* A
ftemp=IntL_P(lb(:),ub(:));! m7 I5 [5 |: i/ ^+ Y2 V
%以下为子函数
7 Y/ \+ z* d+ e) ~- u7 h# q: n
function ftemp=IntL_P(vlb,vub)! K- L$ t$ p& D5 B: R
global upper opt c x0 A b Aeq beq ID options;* Z1 g v. o4 L6 w4 A7 x4 L
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);) V0 ?( i5 h( l\" j8 }; L
if how<0
# }! \' `9 i8 E+ R
return;; Z) Y# ]9 l2 r
end
2 O$ Q6 }) Y5 S7 J* U
if ftemp-upper>0.000054 Q7 p5 b) c' M
return;
8 r7 ~2 V2 }6 E% K: W
end% |) J' b: Q, ^1 W; U- {. B
if max(abs(x.*ID-round(x.*ID)))<0.00005% |8 j- B- X2 {2 s
if upper-ftemp>0.00005& r# l- e) f V! G# V$ W; ~3 g% {& V
opt=x';
L& C' O4 z: h7 h4 p: c
upper=ftemp;& U1 v6 S9 F. z) z6 q4 H
return;
+ ?* |3 U& h1 s7 h( B) e
else
% n9 J\" O: Z\" X! N; g* \
opt=[opt;x'];
, ?. r6 T W8 f0 c9 t
return;
/ T$ p7 V3 @8 L) e2 S
end9 s4 [# m K( M5 V5 {! p/ t
end
2 } J6 ?) ~7 m' F4 _$ y) i& ]; X
notintx=find(abs(x-round(x))<=0.00005);
) F8 _& ^1 D: k- Y* [
intx=fix(x);) X$ X, I\" r9 x& \: u- |5 |
tempvlb=vlb;tempvub=vub;
9 Q/ V4 G M' O4 ?2 a
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+10 m/ l3 K3 D0 `
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;
% Q* b; Y0 X+ r. O3 Z) ~/ P; T
ftemp=IntLP(tempvlb,vub);* Y\" |; n) H( y: W+ `6 e; V
end
4 R, a7 _* v! |; R) q
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)# M' S/ q1 w5 m8 J t
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);; z1 Q) }! k( }7 a& u/ _
ftemp=IntL_P(vlb,tempvub);' D+ |; ^+ P* [ E+ J
end