原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
实用的建议,实用类文本复习建议,环保的实用的,cf最新bug2015实用的,最安全实用的网址导航,实用的对讲机,送老人的实用礼物,最实用的家用车,最实用的手机,实用的生日礼物
%%本程序是用分枝定界法求解整数线性规划问题 %%问题的标准形式： %% min c*x %% s.t. A*x=b %% Aeq*x=beq %% x要求是整数 %%程序的编写者：快乐逍遥侠 %%程序的最后修改时间是：2007-11-24 function [y,fval]=BranchBound(c,A,b,Aeq,beq)
c=[];
A=[];
b=[];
Aeq=[];
beq=[]; NL=length(c); UB=inf; LB=-inf; FN=[0]; AA(1)={A}; BB(1)={b}; k=0; flag=0; while flag==0; [x,fval,exitFlag]=linprog(c,A,b,Aeq,beq); if (exitFlag == -2) | (fval = UB)
FN(1)=[]; if isempty(FN)==1 flag=1; else k=FN(1); A=AA{k}; b=BB{k}; end else for i=1:NL if abs(x(i)-round(x(i)))1e-7 kk=FN(end); FN=[FN,kk+1,kk+2]; temp_A=zeros(1,NL); temp_A(i)=1; temp_A1=[A;temp_A]; AA(kk+1)={temp_A1}; b1=[b;fix(x(i))]; BB(kk+1)={b1}; temp_A2=[A;-temp_A]; AA(kk+2)={temp_A2}; b2=[b;-(fix(x(i))+1)]; BB(kk+2)={b2}; FN(1)=[]; k=FN(1); A=AA{k}; b=BB{k}; break; end end if (i==NL) (abs(x(i)-round(x(i)))=1e-7) UB=fval; y=x; FN(1)=[]; if isempty(FN)==1 flag=1; else k=FN(1); A=AA{k}; b=BB{k}; end end end end y=round(y); fval=c*y;
%整数规划分支定界算法matlab通用源程序 %各参数的意义同matlab优化工具箱的线性规划函数linprog %调用前，输入参数要化成matlab的标准形式 [x,val]=kfz-f-3(n,f,a,b,aeq,beq,lb,ub) x=zeros(n,1); x1=zeros(n,1); m1=2; m2=1; [x1,val1]=linprog(f,a,b,aeq,beq,lb,ub); if (x1==0) x=x1; val=val1; elseif (round(x1)==x1) x=x1; val=val1; else e1={0,a,b,aeq,beq,lb,ub,x1,val1}; e(1,1)={e1}; zl=0; zu=-val1; while (zu~=zl) for c=1:1:m2 if (m1~=2) if (cell2mat(e{m1-1,c}(1))==1) e1={1,[],[],[],[],[],[],[],0}; e(m1,c*2-1)={e1}; e(m1,c*2)={e1}; continue; end; end; x1=cell2mat(e{m1-1,c}(8)); x2=zeros(n,1); s=0; s1=1; s2=1; lb1=cell2mat(e{m1-1,c}(6)); ub1=cell2mat(e{m1-1,c}(7)); lb2=cell2mat(e{m1-1,c}(6)); ub2=cell2mat(e{m1-1,c}(7)); for d=1:1:n if (abs((round(x1(d))-x1(d)))0.0001)(s==0) s=1; lb1(d)=fix(x1(d))+1; if (a*lb1=b) s1=0; end; ub2(d)=fix(x1(d)); if (a*lb2=b) s2=0; end; end; end; e1={s1,a,b,aeq,beq,lb1,ub1,[],0}; e2={s2,a,b,aeq,beq,lb2,ub2,[],0}; e(m1,c*2-1)={e1}; e(m1,c*2)={e2}; end; m1=m1+1; m2=m2*2; for c=1:1:m2 if (cell2mat(e{m1-1,c}(1))==0) [x1,val1]=linprog(f,cell2mat(e{m1-1,c}( 2)),cell2mat(e{m1-1,c}(3)),cell2mat(e{m1-1,c}(4)),c
文档评论（0）