[问题求助] 请求一个整数规划的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)
. b9 K5 B9 N3 Y3 b
global upper opt c x0 A b Aeq beq ID options;
2 ?2 ?3 N9 E. E
if nargin<10,
. y- a# N# A- T
options =optimset({});( ^* ^* Y4 T! p# V' C4 r
options.Display='off';& y3 o* d6 j( N2 U. D: ~/ e6 L
options.LargeScale='off';. }\" F I* [7 q' j/ P+ i# d9 R
end G j. L! u6 I6 w% o8 c% E
if nargin<9,0 a6 a* V* J! }' ]$ j) v k
id=ones(size(f));+ q; ^- I* Y\" c
end
. A. N4 g& G# t* q
if nargin<8,
- H* f8 A$ N1 \9 x5 B
x=[];/ @2 M) D: S. \% k1 T+ X' Z
end! U$ s+ [5 H# S+ x+ P$ @0 F
if nargin<7 | isempty(ub),
$ p( i1 y- ]+ S$ P9 Y2 _
ub =inf*ones(size(f));
q& {. c$ Q1 E3 v, `# ` ~: r
end$ Y0 ~0 G4 l) c. W
if nargin<6 | isempty(lb),- x2 h R! s6 `( Z
lb=zeros(size(f));
4 f9 x* W: s. r, Z
end
% F+ o4 q# n3 P( w% u
if nargin<5,
( t4 m) A3 M6 D1 Y$ H
heq=[];: X. t3 L- m& N9 C2 v1 |6 I% J) R
end
2 M9 z2 z# c. |! E; B
if nargin<4,
) W* d7 B9 I e) z9 ^\" s
Geq=[];
. O) r, ^\" D3 X* W
end
0 Y+ y, a5 K5 P& `2 h' A
upper=inf;c=f;x0=x;# I' j. I6 i) r4 |\" B0 A
A=G;
, W0 P4 l% P( M. s
b=h;Aeq=Geq;; U t. }7 X$ U9 o: T2 m% ^
beq=heq;ID=id;
+ m0 ?+ X+ ~; l$ t; b3 T2 N
ftemp=IntL_P(lb(:),ub(:));! K$ d- y: T6 ~: E4 i3 |0 E0 {
%以下为子函数6 E# T) N; a: O+ ^% ?! s3 P
function ftemp=IntL_P(vlb,vub)
+ a! |6 E, |\" e
global upper opt c x0 A b Aeq beq ID options;9 e% u+ [9 z) \4 |* [
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);0 F. ?* \+ W, T$ _/ z4 Z- g* M
if how<0
# ?; n/ |. C$ d! q$ k$ b) d
return;$ Y: \1 R( M, L7 g
end
% M6 X+ u4 ?* w* J2 u7 D
if ftemp-upper>0.00005
/ [4 ?. q\" [% M8 d8 F
return; o4 g& c8 ]; e( W* ~
end
6 d' S& H6 t: N3 D
if max(abs(x.*ID-round(x.*ID)))<0.000051 B\" k* w7 _) O3 z- `* z O' H
if upper-ftemp>0.00005
9 @8 q, ?! c4 w
opt=x';
. i) ]0 R6 K- d
upper=ftemp;
: y1 T$ X( Y' d+ r4 z6 G, A: x2 l
return; [+ g2 T5 A2 Y
else1 v/ E/ G0 ~2 ^1 L( ]6 B
opt=[opt;x'];; r6 Z8 ^5 U7 b1 u
return;
$ Y8 X4 p% f$ T0 {7 d& V
end
; y, j& U: [) a: F, m. L$ E
end
, N2 I$ a- e# k3 ?2 X5 h
notintx=find(abs(x-round(x))<=0.00005);
- r) F7 X; I: j( O6 P, n% ]
intx=fix(x); B\" D# O8 G4 z( ^ S
tempvlb=vlb;tempvub=vub;! {$ r& T4 C* ?) o
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1
1 s5 C8 u\" a5 {3 w2 O: X, q
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;
! K k( Z+ X\" i5 |2 y\" o
ftemp=IntLP(tempvlb,vub);
! S4 i3 c/ t8 V) j% ?. D
end
6 J% f8 L5 c0 V/ I S# I
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)
' r) H0 H/ \8 e
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);
\" m( M4 C. y) j% ~% e8 d
ftemp=IntL_P(vlb,tempvub);\" n! {- j4 [5 Y M' \3 `
end