分支定界法的Matlab实现

1、一个调用例子：ifin t=0 1;f=10 9;a=1 0;0 1;-5 -3;b=8 10 -45;x,fval,exitflag = linprogdis(ifint,f,a,b,)fun cti onr= check int (x)%算法：如果 x(i)是整数，则返回r(i)=1;否则返回r(i)=0fun cti onr= ifrowi nmat(arow,amat)%输入参数：% arow向量，% amat 矩阵% 设计：如果arow 与矩阵amat中的某一行相等，则返回1，如果找不到现等的一行，则返回 0可以使用 ismember(arrow,amat, rows 替换 ifro

2、winmat的调用，2005-10-28 标注使用时，将下面的代码存入文件：lin progdis.mfun cti onx,fval,exitflag,output,lambda二.lin progdis(ifi nt,f,A,b,Aeq,beq,lb,ub,xO,opti ons)%Title:%分支定届法求解混合整数线性规划模型%初步完成：2002年12月%最新修订：2004-03-06%最新注释：2004-11-20%数据处理t1, t2 = size(b);if t2=1,b=b;%将b转置为列向量end%调用线性规划求解x,fval,exitflag,output,lambda=l

3、in prog(f,A,b,Aeq,beq,lb,ub,x0,opti on s);if exitflag0 & exitflag2=0& fval0 & exitflag1=0 & fval0 & exitflag20 & fval=fval1 & fval=fval2 ), disp(error call)return%isempty(v2)表示都是整数 2002-12-7非常重要确定最终的最优解end %下面分别根据返回标志exitflag%case 1& exitflag20 优解，右分支无最优解x = x1; fval = fval1; exitflag = exitflag1; o

4、utput = output1; lambda = lambda1;%case 2elseif exitflag20 & exitflag10 & exitflag20%【左、右分支均有】最优解，则比较选优if fval1fval2, %【左】分支最优(min)x = x1;fval = fvall;exitflag = exitflagl; output = outputl;lambda = lambdal;elsex = x2;,fval = fval2;exitflag = exitflag2;output = output2; lambda = Iambda2;end%fval1fva

5、l2end%【右】分支最优(min )fun cti onr= check int(x)%算法：如果x(i)是整数，则返回r(i)=1;否则返回r(i)=0%输入参数：x 向量%输出参数：r 向量fori=1:le ngth(x),ifmi n(abs(x(i)-floor(x(i),abs(x(i)-ceil(x(i)1e-03%这里用于判定是否为0的参数可以调整，如改为1e-6r(i)=1;elser(i)=0;endendfun cti onr= ifrowi nmat(arow,amat)%输入参数：% arow向量，% amat 矩阵% 设计：如果arow与矩阵amat中的某一行相等，则返回1， 如果找不到现等的一行，则返回0r = 0;rows = size(amat,1);for i=1:rows,temp= (amat(i,:)=arow);%利用 Matlab 的