[问题求助] 请求一个整数规划的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)- Z' @& @/ Q\" y$ [
global upper opt c x0 A b Aeq beq ID options;
, K* E2 Q( U) e b
if nargin<10,; Y, w1 n; [0 e% \; E
options =optimset({});
$ D3 z& K0 H3 F- n. I i4 {) X
options.Display='off';
; H$ h3 A5 ~! R3 h( X( ~
options.LargeScale='off';, c. B7 R; ?+ c\" n; w
end, B4 A, M- E( k1 w! F4 G
if nargin<9,
* k- z0 \' B$ Y8 T0 x+ {; R
id=ones(size(f));; D V% \' z4 ?
end ( r\" }. I' r( ^' a) |2 w3 e
if nargin<8,8 P3 @: F. M/ S\" `
x=[];
: e' x9 b! c# t2 n. X
end
+ N5 C9 H. I; ~8 i) R v9 j
if nargin<7 | isempty(ub),
0 w: k U9 l* i7 I2 P
ub =inf*ones(size(f));, n2 v& d# M6 }+ `& j2 w: N+ h
end( n5 ^' T; G' V, K' J+ O- D
if nargin<6 | isempty(lb),
% O8 V) b* S0 K r\" @1 i) c
lb=zeros(size(f));
2 k! q2 |+ {! j
end+ y5 f8 `1 Z7 ~& w! Z% ^* ^
if nargin<5,
o3 D% M. \( ]( a
heq=[];
1 b% N+ u% l, T: {
end
) C# r8 _! w\" Z$ c8 Y) j+ v9 x3 ~/ ~
if nargin<4,3 F! @% ]2 u1 K% f
Geq=[];* n w* o; T/ X8 G# k6 \3 P
end
4 b6 V3 T( u! C9 O
upper=inf;c=f;x0=x;/ ]4 u9 B; T5 U! ` a' b+ K
A=G;6 Q: M; H! i4 w8 p, z6 b\" A
b=h;Aeq=Geq;- x% f+ s% O) M, ?8 k' }
beq=heq;ID=id;
8 ?& @9 y( B9 ~: d
ftemp=IntL_P(lb(:),ub(:));/ j3 D/ B3 V/ Z, r\" N3 Q. B7 s' c
%以下为子函数6 h- z\" }5 P\" j! ~& h L
function ftemp=IntL_P(vlb,vub)
\" ]. R5 L2 S. y
global upper opt c x0 A b Aeq beq ID options;
8 `4 r5 i. R; G
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);* r: U8 {( V/ [6 f3 i) {: I5 R\" \
if how<0& \/ c) Z& ]( T$ O0 M/ M
return;
5 x% P! h1 E5 y\" b
end
% y! r1 \2 K) J4 I, u2 w
if ftemp-upper>0.00005! [\" z2 V( C6 `4 ~# p7 r
return;
) b\" ]8 [3 t1 g9 ]7 c- P+ a1 x3 S
end2 \8 R2 U3 G\" _ b/ X( O$ L8 Z
if max(abs(x.*ID-round(x.*ID)))<0.000052 T* L+ i( C6 z0 b' E+ J! }% F: [
if upper-ftemp>0.000053 M/ \4 y+ f0 H: F# j& }
opt=x';% |7 g7 ^$ c5 K
upper=ftemp;
( M- X5 {) P\" s, h% b% H' Q2 u
return;
9 P! |. m X\" m5 { d, t) O
else- x$ E7 s, L+ m\" z! }/ h) g9 \1 m0 ?
opt=[opt;x'];; d* ?7 B/ F7 ~\" k! X ]
return;7 \, U4 V: ?) U m, k) h. ?7 W2 s3 [
end
X4 u$ g6 o) a/ f9 \
end
+ ]/ f/ U# E# J5 U) U: [
notintx=find(abs(x-round(x))<=0.00005);3 D7 b3 M, E' i; F7 O) r; _# K
intx=fix(x);1 @9 V1 c) ?- ~% L9 x! O* J* e- N
tempvlb=vlb;tempvub=vub;! `; v! r* `- U- c) a5 c
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1% h/ Z9 t2 R# c. \. N, J: b* h4 h9 E7 y8 Z
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;! e' S+ G# V- U( @2 T0 \8 U
ftemp=IntLP(tempvlb,vub);! N* l0 L\" n+ P\" ^1 `
end$ g$ H' @) {2 l1 o
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)
- o: V! ^7 V4 ~& D
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);
* X5 r3 y/ e# Z4 o& M, M. w
ftemp=IntL_P(vlb,tempvub);. u6 O$ B+ X) h7 C; c, ^2 s' n/ r
end