branch.m

采用分支定界算法

function [jie,fval,exflag]=branch(f,a,b,aeq,beq,lb,ub,intbegan)
[m,n]=size(a);
global optjie;
global optfval;
global time;
time=time+1;
[jie,fval,exflag]=linprog(f,a,b,aeq,beq,lb,ub);
if exflag<=0   %去掉不必要的分支（无解的分支）
    return;
else
    [t0,k,z]=panduan(jie,intbegan);
    if t0==0
        jie1=jie;
        if fval<=optfval
            optfval=fval;
            optjie=round(jie);
            optjie(1:intbegan-1)=jie(1:intbegan-1);%disp('dsfsdfsf');
            return;
        end
    elseif fval<optfval %定界，去掉不必要的分支（结果不优的分支）
        a1=a;%第一个分支
        b1=b;       
        al(m+1,:)=zeros(1,n);
        a1(m+1,k)=1;
        b1(m+1)=round(jie(k))-z;
        [jie,fval,exflag]=branch(f,a1,b1,aeq,beq,lb,ub,intbegan);
        %\\\\\\\\\\\\\\\
        a2=a;%第二个分支
        b2=b;
        a2(m+1,:)=zeros(1,n);
        a2(m+1,k)=-1;
        b2(m+1)=z-round(jie(k))-1;
        [jie,fval,exflag]=branch(f,a2,b2,aeq,beq,lb,ub,intbegan);
        return;
    end
end