[问题求助] 请求一个整数规划的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\" a( N$ z( x- D- E
global upper opt c x0 A b Aeq beq ID options;7 J/ C/ B* D; x8 O- o6 K4 a# k+ x
if nargin<10,* T% {/ S- l) Q$ o* L
options =optimset({});
3 l+ G: R' e, J' }& C
options.Display='off';8 N$ @ v( O. o
options.LargeScale='off';
$ I) O [/ \7 }; w
end) [ n6 s9 N6 t* _ e! x+ d+ F
if nargin<9,
% m2 k; Y7 z! q
id=ones(size(f));, [3 L: s }5 f+ C* b# w) {: [
end
/ g\" |# u1 l, d: ]1 c6 Q, i8 V
if nargin<8,$ E' _/ ]! T3 Q
x=[];\" @) `1 U8 ~& X9 i
end$ f3 D* F3 U7 ?$ ?, T
if nargin<7 | isempty(ub),
6 `1 B2 }, M( _
ub =inf*ones(size(f));6 D4 N0 I5 B7 g l$ l\" {4 b4 j5 Y4 B
end
( h2 a$ O& Z+ i
if nargin<6 | isempty(lb),) f% p3 ~ a. @/ ^/ I- s. T4 z
lb=zeros(size(f));4 x& W( ]& C5 d' P4 P1 {9 F. g- F
end/ B4 j4 j+ G8 a
if nargin<5,( B* D3 {0 `; B
heq=[];
. ?5 \4 V1 S8 {9 ^$ h
end6 \0 r; F8 m p\" \- w) P3 O/ S8 g H
if nargin<4,
, D8 T5 ?0 P8 e5 E- G
Geq=[];
- |$ h5 a( ^; u% `. X9 _* j
end W* o( r+ g9 {& o
upper=inf;c=f;x0=x;: ]) Q\" D' B0 m\" Y
A=G;5 z0 `- }6 T$ e( h. t% Q! e. K
b=h;Aeq=Geq;
2 }( }: X9 B$ T0 p
beq=heq;ID=id;
' ^/ ?* k8 B0 z
ftemp=IntL_P(lb(:),ub(:));
7 z& {; r1 n+ w; }) f
%以下为子函数
2 {- \# G }8 B \
function ftemp=IntL_P(vlb,vub)
$ }/ x+ \9 _5 \: ]0 T
global upper opt c x0 A b Aeq beq ID options;/ z1 g# G# ]8 t! |
[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);. ]6 f' {6 w* A3 |) r1 |\" [
if how<0
4 O5 g; L' e) r8 N% U
return;
4 h: f# o# M* V' Y7 R. [+ d: Q
end
: Z; i8 r% g+ U# @2 ^0 {. s
if ftemp-upper>0.00005
, }) h' J8 \/ ~7 [
return;
4 [6 g1 a5 ?+ D$ D2 H
end
5 \$ Y3 M& R/ ^- B r\" j! q% b. @
if max(abs(x.*ID-round(x.*ID)))<0.00005( t3 Q+ u: R7 X1 \! o
if upper-ftemp>0.00005
h4 _0 a- t; @\" W' _! O8 k4 \
opt=x';. A3 M+ s$ q( O8 ~% E
upper=ftemp;
8 v8 g! u& i; B2 L* f, K
return;1 G9 n, Y4 X6 `% k! g$ V
else
* I% p, `- J% f b
opt=[opt;x'];
' M' D( ^$ T, V: d3 D5 v; L
return;
5 b- W9 ?- Q1 [0 ~8 h% l. R% Y; v
end
7 f! J. J; m. O
end
/ m0 a0 [# `\" H
notintx=find(abs(x-round(x))<=0.00005);8 v; K. D+ E4 ~) a7 ~8 N) q* z6 [- |
intx=fix(x);9 g: j% O- i& B4 H
tempvlb=vlb;tempvub=vub;* N9 \; v+ |' ]3 q) o& ]4 \7 l
if vub(notintx(1,1),1)>=intx(notintx(1,1),1)+1
4 _7 o7 x/ R9 ]% b) H\" w4 v
tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;8 u9 j; V5 E# R2 i2 F
ftemp=IntLP(tempvlb,vub);3 l+ M& K\" S) \' I7 d3 [
end4 e6 p- z) i: Q2 ~: w
if vlb(notintx(1,1),1)<=intx(notintx(1,1),1)
7 G8 |! D; Y# ]- L9 ?; I/ Y! M
tempvub(notintx(1,1),1)=intx(notintx(1,1),1);. D* B$ E. W/ t$ Q3 p
ftemp=IntL_P(vlb,tempvub);
+ w3 n' b6 g0 |8 I: ^\" ]# ~
end