【老生谈算法】MATLAB分支定界法程序源码.doc资源-CSDN文库

版权申诉

5星 · 超过95%的资源 136 浏览量 2023-05-11 21:03:25 上传评论收藏 37KB DOC 举报

资源推荐

资源详情

资源评论

源代码如下：

function [x,y]=ILp(f,G,h,Geq,heq,lb,ub,x,id,options)

%整数线性规划分支定界法，可求解纯整数规划和混合整数规划1。

%y=minf’*x s.t. G*x<=h Geq*x=heq x 为全整数或混合整数列向量

%用法

%[x,y]=ILp(f,G,h,Geq,heq,lb,ub,x,id,options)

%参数说明

%lb:解的下界列向量（Default:-int）

%ub:解的上界列向量（Default:int）

%x:迭代初值列向量

%id：整数变量指标列向量，1-整数，0-实数（Default:1）

global upper opt c x0 A b Aeq beq ID options;

if nargin<10,options=optimset({});options.Display='off';

options.LargeScale='off';end

if nargin<9,id=ones(size(f));end

if nargin<8,x=[];end

if nargin<7 |isempty(ub),ub=inf*ones(size(f));end

if nargin<6 |isempty(lb),lb=zeros(size(f));end

if nargin<5,heq=[];end

if nargin<4,Geq=[];end

upper=inf;c=f;x0=x;A=G;b=h;Aeq=Geq;beq=heq;ID=id;

ftemp=ILP(lb(:),ub(:));

x=opt;y=upper;

%下面是子函数

function ftemp=ILP(vlb,vub)

global upper opt c x0 A b Aeq beq ID options;

[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);

if how <=0

return;

end;

if ftemp-upper>0.00005 %in order to avoid error

return;

end;

if max(abs(x.*ID-round(x.*ID)))<0.00005

if upper-ftemp>0.00005 %in order to avoid error

opt=x';upper=ftemp;

return;

else

本内容试读结束，登录后可阅读更多

下载后可阅读完整内容，剩余4页未读，立即下载

内容反馈

版权申诉

suda_93

2024-01-22

感谢大佬，让我及时解决了当下的问题，解燃眉之急，必须支持！
qq_39969796

2024-04-29

果断支持这个资源，资源解决了当前遇到的问题，给了新的灵感，感谢分享~
qiyunye

2024-05-06

资源是宝藏资源，实用也是真的实用，感谢大佬分享~

阿里matlab建模师

粉丝: 3325
资源: 2784

最新资源

资源上传下载、课程学习等过程中有任何疑问或建议，欢迎提出宝贵意见哦~我们会及时处理！点击此处反馈

feedback-tip