[问题求助] 请求一个整数规划的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)8 j\" W) ^4 H* F, C& s8 F& [
global upper opt c x0 A b Aeq beq ID options;
; ]. ^$ d& _5 ~7 W/ O
if nargin<10,$ t1 w6 a; e+ T( ~' W
options =optimset({});/ [( V# V# G8 S B7 Z
options.Display='off';6 Q: D8 @, e& Y$ }
options.LargeScale='off';
' ^ k, e& L' d9 j* I
end
8 _: A* b% @% d' e1 p8 g
if nargin<9,
0 i ^\" z2 t4 p2 j
id=ones(size(f));
3 ]1 y' N, \1 @: n! S
end
$ G/ b+ ~, E1 k% k. @5 m
if nargin<8,
; I# ~6 L* {0 {
x=[];
4 a) r5 C+ ], F0 U& G! ^
end
. u- ]& B8 [$ {5 ~) h
if nargin<7 | isempty(ub),
) X: D; R2 ]\" w: T: g u3 l
ub =inf*ones(size(f));
# P& d- D' I. w0 S* B& h& z. l
end
8 P1 o$ Z2 T( x) h* O
if nargin<6 | isempty(lb),
; v h8 A9 L; [) i5 w& Y% x5 S
lb=zeros(size(f)); r, _' f1 X5 t, I4 d
end
. Q$ t; o* ^* C4 t
if nargin<5,0 {) q; U; c0 F( U* L
heq=[];
' Q8 M3 |' L7 s3 M6 ~6 ^
end2 X: Q4 a, M! e( C* F, T/ v# O
if nargin<4,0 E, ~8 ^8 R6 a8 G
Geq=[];
! y# J\" {2 d8 d' d
end4 b$ x1 G7 N8 o2 @! D
upper=inf;c=f;x0=x;- Z) `7 `: ]; }/ t
A=G;$ w: b; E& N4 `/ q2 U6 J, C
b=h;Aeq=Geq;
- T! j' H7 I0 L- T
beq=heq;ID=id;+ q* D/ U1 z n
ftemp=IntL_P(lb(:),ub(:));6 r\" z# F3 T1 I
%以下为子函数/ ~. U- v% v4 u% D; @+ s) a
function ftemp=IntL_P(vlb,vub)0 ]: \, ^4 m6 s+ b. A& }
global upper opt c x0 A b Aeq beq ID options;
+ g' V8 u, S) a8 O\" p: C3 b3 @
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);
+ ] }; Q$ e5 r+ X
if how<0& X |1 R$ @, b( j) D6 E6 R O
return;
4 Z4 l5 J) o- ^, Q {! A
end8 c! X, i5 V `3 ?9 P, G- J9 O+ d
if ftemp-upper>0.000050 i4 ?. |, R; I6 @: Q6 P+ I3 t
return;
$ ?8 U9 T9 ^1 v7 ]; M
end# j- @: i4 g: I7 Z+ g# ]
if max(abs(x.*ID-round(x.*ID)))<0.00005% w: g# n1 e, u- i) h% H b
if upper-ftemp>0.00005* l8 {5 Z2 r& {4 [) o
opt=x';5 f! _ U( y- S
upper=ftemp;& n, ], C% @' w* t5 ?* @8 G# F
return;) z, z# U3 N7 \) o
else
: Y2 L5 d% S0 Q% ?. e
opt=[opt;x'];
- y$ w! p* P) H
return;4 X% V: w; @) `& N
end
* U$ L$ p/ w# g6 A9 c+ O
end
9 m d$ w& Q/ J6 E! D\" Z
notintx=find(abs(x-round(x))<=0.00005);
: A: h' }! m* G w) j+ }! v
intx=fix(x);6 H! w& b4 w! S0 M, \/ d; s- k, v, p1 o
tempvlb=vlb;tempvub=vub;
4 M; U0 k8 s6 V% Z! O
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1
4 y% x# m& Z. y5 n6 T% c
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;
* K3 J, ^\" n\" x F( C
ftemp=IntLP(tempvlb,vub);9 O, Y9 t6 a( G# }1 s' N* b. h
end1 f: Y+ g/ @0 X( S: x3 I; D# u2 a
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)0 c: G: T* Y2 V$ G\" V
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);6 Y6 J$ |2 q1 M3 k1 _3 k& g
ftemp=IntL_P(vlb,tempvub);
1 U: t\" u2 L( b8 `8 C) S Z
end