MATLAB分支定界法程序

源代码如下：function[x,y]=ILp(f,G,h,Geq,heq,lb,ub,x,id,options)%整数线性规划分支定界法，可求解纯整数规划和混合整数规划。%y=minf’*xs.t.G*x=hGeq*x=heqx为全整数或混合整数列向量%用法%[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）globalupperoptcx0AbAeqbeqIDoptions;ifnargin10,options=optimset({});options.Display='off';options.LargeScale='off';endifnargin9,id=ones(size(f));endifnargin8,x=[];endifnargin7|isempty(ub),ub=inf*ones(size(f));endifnargin6|isempty(lb),lb=zeros(size(f));endifnargin5,heq=[];endifnargin4,Geq=[];endupper=inf;c=f;x0=x;A=G;b=h;Aeq=Geq;beq=heq;ID=id;ftemp=ILP(lb(:),ub(:));x=opt;y=upper;%下面是子函数functionftemp=ILP(vlb,vub)globalupperoptcx0AbAeqbeqIDoptions;[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);ifhow=0return;end;ifftemp-upper0.00005%inordertoavoiderrorreturn;end;ifmax(abs(x.*ID-round(x.*ID)))0.00005ifupper-ftemp0.00005%inordertoavoiderroropt=x';upper=ftemp;return;elseopt=[opt;x'];return;end;end;notintx=find(abs(x-round(x))=0.00005);%inordertoavoiderrorintx=fix(x);tempvlb=vlb;tempvub=vub;ifvub(notintx(1,1),1)=intx(notintx(1,1),1)+1;tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;ftemp=IntLP(tempvlb,vub);end;ifvlb(notintx(1,1),1)=intx(notintx(1,1),1)tempvub(notintx(1,1),1)=intx(notintx(1,1),1);ftemp=IntLP(vlb,tempvub);end;在命令行运行c=[1,1,-4];a=[1,1,2;1,1,-1;-1,1,1];b=[9;2;4];[x,f]=ILp(c,a,b,[],[],[0;0;0],[inf;inf;inf])从你的程序中的如下这条语句可以看出变量ID是全局变量：globalupperoptcx0AbAeqbeqIDoptions;该全局变量（ID）是在你的函数之外被别的程序语句进行过赋值，而且也很有可能被赋成了向量/矩阵、并且它的维数大小跟向量x的维数大小并不相同，这就造成了进行运算x.*ID的时候维数不匹配的问题第二个：题目:min(4*x1+4*x2);约束条件:2*x1+5*x2=15,2*x1-2*x2=5,x1,x2=0,且都为整数.解这个还是很容易,算出来x1,x2都为0点几,因为题目要求是整数,所以主要是这个分支定界的问题,急求一个分支定界的MATLAB算法,通用算法也可以,或者只能解这道题也可以,只要能进行计算就行,最后解出来x1,x2都为0.希望大家帮帮忙,先谢谢大家了.解线性方程很容易，但主要是整数解，需要用分支定界法．希望算法尽量简单易懂,最好有注释,本人是初学者,谢谢了把以下程序存为ILP.m，%============================function[x,y]=ILp(f,G,h,Geq,heq,lb,ub,x,id,options)%整数线性规划分支定界法，可求解纯整数规划和混合整数规划。%y=minf’*xs.t.G*x=hGeq*x=heqx为全整数或混合整数列向量%用法%[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）globalupperoptcx0AbAeqbeqIDoptions;ifnargin10,options=optimset({});options.Display='off';options.LargeScale='off';endifnargin9,id=ones(size(f));endifnargin8,x=[];endifnargin7|isempty(ub),ub=inf*ones(size(f));endifnargin6|isempty(lb),lb=zeros(size(f));endifnargin5,heq=[];endifnargin4,Geq=[];endupper=inf;c=f;x0=x;A=G;b=h;Aeq=Geq;beq=heq;ID=id;ftemp=ILP(lb(:),ub(:));x=opt;y=upper;%下面是子函数functionftemp=ILP(vlb,vub)globalupperoptcx0AbAeqbeqIDoptions;[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);ifhow=0return;end;ifftemp-upper0.00005%inordertoavoiderrorreturn;end;ifmax(abs(x.*ID-round(x.*ID)))0.00005ifupper-ftemp0.00005%inordertoavoiderroropt=x';upper=ftemp;return;elseopt=[opt;x'];return;end;end;notintx=find(abs(x-round(x))=0.00005);%inordertoavoiderrorintx=fix(x);tempvlb=vlb;tempvub=vub;ifvub(notintx(1,1),1)=intx(notintx(1,1),1)+1;tempvlb(notintx(1,1),1)=intx(notintx(1,1),1)+1;ftemp=IntLP(tempvlb,vub);end;ifvlb(notintx(1,1),1)=intx(notintx(1,1),1)tempvub(notintx(1,1),1)=intx(notintx(1,1),1);ftemp=IntLP(vlb,tempvub);end;%====================================然后：clc;clearf=[44]A=[25;2-2]b=[15;5]Aeq=[];beq=[];LB=[00];UB=[];[xn,yn]=ILp(f,A,b,Aeq,beq,LB,UB,[11],1,[])[x,fval,exitflag]=linprog(f,A,b,Aeq,beq,LB,UB)结果：xn=00yn=0Optimizationterminated.x=1.0e-013*0.2990040786747590.503948216933779fval=3.211809182434153e-013exitflag=1matlab的整数规划功能不行，还不如EXCEL的solver。提问者评价貌似还是不能正确得出整数