資源描述：
《MATLAB分支定界法程序.doc》由會員上傳分享，免費在線閱讀，更多相關內(nèi)容在教育資源-天天文庫。
1、源代碼如下：function[x,y]=ILp(f,G,h,Geq,heq,lb,ub,x,id,options)%整數(shù)線性規(guī)劃分支定界法，可求解純整數(shù)規(guī)劃和混合整數(shù)規(guī)劃。%y=minf’*xs.t.G*x<=hGeq*x=heqx為全整數(shù)或混合整數(shù)列向量%用法%[x,y]=ILp(f,G,h,Geq,heq,lb,ub,x,id,options)%參數(shù)說明%lb:解的下界列向量（Default:-int）%ub:解的上界列向量（Default:int）%x:迭代初值列向量%id：整數(shù)變量指標列向量，1-整數(shù)，0-實數(shù)（Default:1）global
2、upperoptcx0AbAeqbeqIDoptions;ifnargin<10,options=optimset({});options.Display='off';options.LargeScale='off';endifnargin<9,id=ones(size(f));endifnargin<8,x=[];endifnargin<7
3、isempty(ub),ub=inf*ones(size(f));endifnargin<6
4、isempty(lb),lb=zeros(size(f));endifnargin<5,heq=[];endifnar
5、gin<4,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;%下面是子函數(shù)functionftemp=ILP(vlb,vub)globalupperoptcx0AbAeqbeqIDoptions;[x,ftemp,how]=linprog(c,A,b,Aeq,beq,vlb,vub,x0,options);ifhow<=0return;end;ifftemp-upper>0.00005%inordertoavo
6、iderrorreturn;end;ifmax(abs(x.*ID-round(x.*ID)))<0.00005ifupper-ftemp>0.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(n
7、otintx(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,
8、b,[],[],[0;0;0],[inf;inf;inf])從你的程序中的如下這條語句可以看出變量ID是全局變量：globalupperoptcx0AbAeqbeqIDoptions;該全局變量（ID）是在你的函數(shù)之外被別的程序語句進行過賦值，而且也很有可能被賦成了向量/矩陣、并且它的維數(shù)大小跟向量x的維數(shù)大小并不相同，這就造成了進行運算x.*ID的時候維數(shù)不匹配的問題第二個：題目:min(4*x1+4*x2);約束條件:2*x1+5*x2<=15,2*x1-2*x2<=5,x1,x2>=0,且都為整數(shù).解這個還是很容易,算出來x1,x2都為0點幾,因
9、為題目要求是整數(shù),所以主要是這個分支定界的問題,急求一個分支定界的MATLAB算法,通用算法也可以,或者只能解這道題也可以,只要能進行計算就行,最后解出來x1,x2都為0.希望大家?guī)蛶兔?先謝謝大家了.解線性方程很容易，但主要是整數(shù)解，需要用分支定界法．希望算法盡量簡單易懂,最好有注釋,本人是初學者,謝謝了把以下程序存為ILP.m，%============================function[x,y]=ILp(f,G,h,Geq,heq,lb,ub,x,id,options)%整數(shù)線性規(guī)劃分支定界法，可求解純整數(shù)規(guī)劃和混合整數(shù)規(guī)劃。%y=
10、minf’*xs.t.G*x<=hGeq*x=heqx為全整數(shù)或混合整數(shù)列向量%用法%[x,y
