[问题求助] 请求一个整数规划的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)
: i/ Q2 b# @) U- O, o- k3 q
global upper opt c x0 A b Aeq beq ID options;1 B+ x+ M* ]\" _7 h
if nargin<10,
4 F3 n. G) U, Z\" w
options =optimset({});& T) ^. q. R4 l4 j0 n
options.Display='off';\" I* S( {. r4 d8 m6 k\" a! ^9 o7 m$ A
options.LargeScale='off';$ s m( V& V\" f
end9 b1 r) ^& m0 J: J
if nargin<9,
$ K6 e2 e$ ]# \( v2 S' B. R5 f
id=ones(size(f));' R( d8 Z' m8 n* c7 u6 {- d
end
0 f& E( F! B7 i! n# v, V\" `6 c
if nargin<8,
3 |- W, V& z% u7 i% B/ S! X1 w/ ~- ?
x=[];
* C% Y) s0 J; t4 O9 ], P! i
end- K8 m\" ~, B% _0 G9 i3 V2 X9 x
if nargin<7 | isempty(ub),
% E3 n, c( y8 y8 G; v8 p
ub =inf*ones(size(f));! G# p) d% o8 M, a
end) c7 U1 _& T* v2 H
if nargin<6 | isempty(lb),
$ c9 q+ e# P% `& T2 ~) s& {+ _
lb=zeros(size(f));
5 R; j! K, l4 ^% W2 I H- M' O
end
1 E: u# I, D( Z7 s( K! C! v5 a
if nargin<5,
6 J) `1 C4 X; U3 c+ ~9 z0 {- b
heq=[];
5 p9 D9 T7 E& |7 B% ?
end
, F8 J6 B! b$ y& R3 r0 U\" t2 ?
if nargin<4,
0 u4 |4 \: _\" ?. ^
Geq=[];: g$ U/ L! s i! L
end
* n+ y# D5 I+ K
upper=inf;c=f;x0=x;
' V+ ^. r# u\" }8 X( n. t& }
A=G;
7 y5 p7 H) C' } Z1 b& P- h
b=h;Aeq=Geq;
0 T! C5 k7 T6 r, c# J
beq=heq;ID=id;
) T8 X% F0 ?/ f! i' f9 q, n% R
ftemp=IntL_P(lb(:),ub(:));8 h8 E2 g& M- V
%以下为子函数1 U9 g' o& d' R0 e$ h: x( Q$ s( j
function ftemp=IntL_P(vlb,vub)
/ N; @/ y\" d0 y' a8 p- U
global upper opt c x0 A b Aeq beq ID options;
1 q7 h9 a3 ~8 F* z J. w
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);- D\" f; `0 F/ ?7 y
if how<0
4 k7 B- m9 f% c. t! G/ s5 }% d8 ~
return;) }& }0 G; A& _7 |- k& {* ]9 Y2 c ?
end& g. _7 b5 j! l/ F& C3 S
if ftemp-upper>0.000052 C3 a4 h7 U/ \ @9 e4 o
return;
; E) p# ?' M$ q0 O8 Z N$ a* Z& A
end
8 O/ V: @9 }2 ~7 Z* R* m\" i
if max(abs(x.*ID-round(x.*ID)))<0.00005
& H2 v/ ]! N+ _. v$ u! J
if upper-ftemp>0.00005
. V+ g\" F: g# ~) E8 W! [5 e
opt=x';1 _! U6 h9 S2 }5 Y r
upper=ftemp;\" J& U4 R! l, t* m7 s
return;2 G) W3 N I0 T/ U, d
else( Z y+ s7 t' S( F; H2 |
opt=[opt;x'];
# Y$ j) H- u/ z6 M% E- e
return;
5 a' {, }* T0 N: r7 y/ f( \
end
6 E# S# v. h1 B% ]\" S' r
end) z9 ]2 C( }\" o3 G! m+ n
notintx=find(abs(x-round(x))<=0.00005);
; r9 [4 P% y7 [; B# v9 u- D! R
intx=fix(x);
0 \; h0 T8 B4 q! g
tempvlb=vlb;tempvub=vub;
# j7 P+ H9 o( m9 }/ i! r+ b
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1# @( X, y+ C. [; A+ v
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;: U0 t! [$ o\" x& P1 R1 L8 t
ftemp=IntLP(tempvlb,vub);' d* A6 x9 |9 @8 {% I, Q
end) ]# a$ E3 Q- B* V9 ]0 [
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)# u7 N5 O U/ @9 ~. L; e. L
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);& J$ A b$ b. G0 M y( q
ftemp=IntL_P(vlb,tempvub);
, J. x) u. n( L
end