[问题求助] 请求一个整数规划的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)- T) R& p6 Z q' I3 I6 q4 a& b5 L\" Z
global upper opt c x0 A b Aeq beq ID options;
l* s# \# g\" ]( d4 M4 Z8 U
if nargin<10,
$ l2 V) M7 Q! q. ~9 v
options =optimset({});2 k5 }% X8 W7 a/ f
options.Display='off';
( ?! H# c2 h; v) |/ `/ K4 g( S7 ^5 L
options.LargeScale='off';
- v/ e$ X8 z; x6 m. M7 e; q/ x- m- ]
end
: \( W7 R. o' }' X4 ]: w1 L
if nargin<9,9 c) v1 H6 G) M* m
id=ones(size(f));
) b/ {4 M9 t) R/ @6 ?8 m( ^
end - {7 u( G1 w4 A- Z( q; Z$ k- j2 N( W; @
if nargin<8,
- O; ]7 q6 M4 w' Y9 v9 V1 |
x=[];: i! W0 h; a. g6 ]5 N( d
end; E% s\" @* O- M; r4 z$ Z$ Q+ h7 l- b4 }
if nargin<7 | isempty(ub),
* T7 e' D9 ]6 { a8 Q
ub =inf*ones(size(f));
6 h' I. Q; j+ s
end
& w\" j3 j& S% B; G/ K/ g
if nargin<6 | isempty(lb),
' g7 l$ [( G: z; i: ~# K, o
lb=zeros(size(f));; q% E) o\" ^' r& G: V- @
end
J1 B9 |9 x- s, ^6 {
if nargin<5,
9 A8 W4 U, L4 q' C; G! L1 m/ {
heq=[];, R2 N7 x2 N# n- C: P7 B
end A. I8 |- @5 S' A% G, s
if nargin<4,
5 W9 j& j$ }& @' t
Geq=[];! C' q. p# I4 z& `
end\" b8 M5 i: Q3 n4 [% o! e
upper=inf;c=f;x0=x;9 _! u0 J8 {$ n5 X1 S$ y9 U7 \$ J
A=G; B! c7 M1 O2 Y
b=h;Aeq=Geq;
) V3 h0 E3 @, P& D3 B
beq=heq;ID=id;. z1 G, L6 a+ k6 f2 Q' B
ftemp=IntL_P(lb(:),ub(:));& z* D$ x `6 U# C6 t) v( y2 B& ?
%以下为子函数
' E5 G# r' o4 L- O
function ftemp=IntL_P(vlb,vub)
7 T1 S5 k! {5 v, X
global upper opt c x0 A b Aeq beq ID options;
\" b# p9 P1 r+ t% M6 |8 e6 k
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);
% k+ ^. P6 [& K& @# L7 y; Y
if how<0# V\" Q5 Y ]: q- G; Q& l% M3 h; U
return;+ y4 q* m! t: z+ q# w$ |: d! e- j
end0 F\" r; |6 E\" C6 y, `+ e
if ftemp-upper>0.000057 {/ T* }3 G# e! A2 a- \6 g
return;
1 ~0 C' s L8 j
end. d# ]% v3 s a! `; u4 [ s
if max(abs(x.*ID-round(x.*ID)))<0.00005
4 @5 y9 z# p( V& h/ |' S
if upper-ftemp>0.00005
6 M# Q6 H' c! o Q4 r+ K' K
opt=x';5 R8 X, a$ O/ F4 R- D; g7 A
upper=ftemp;
# R* j: C/ }) G( S8 R5 C' G# u: `0 a
return;
\" t. y: @\" Y8 I/ J+ c
else
( Y- C6 Q) m0 G
opt=[opt;x'];% J# i, P+ X- w* n; k. f0 O
return;
[1 L9 B4 [/ x/ \1 {3 l7 K
end
4 S) S L% {8 ~\" ^
end7 Y& S% X Y' q# L
notintx=find(abs(x-round(x))<=0.00005);9 ]% a! U: k& W+ g2 z3 Y4 J; X
intx=fix(x);
1 q1 Q+ n* M\" F+ S8 t
tempvlb=vlb;tempvub=vub;: g% L/ c8 D# @
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1* Z, w2 \2 b8 J0 I' ^
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;
4 ^4 w- w2 b* k8 z
ftemp=IntLP(tempvlb,vub);
8 q9 U\" m9 m9 z5 L! Q
end+ o) L3 l4 Q$ G# ~
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)
5 O$ s$ X' |: J- e- b6 M, @3 O/ d
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);
/ g; N' g3 [0 m, [( ~+ ^
ftemp=IntL_P(vlb,tempvub);7 t; \6 h9 d! [$ Z2 w
end