1、Linprogdis子程序：函數(shù)x，fval，exitflag，output，lambda=.linprogdis(ifint、f、A、b、Aeq、beq、lb、ub、x0、options)%Title:用%分支計時法求解混合整數(shù)線性規(guī)劃模型%初步完成百分比：2002年12月%最新版本： 2004年3月6日%最新說明：2004年11月20日%數(shù)據(jù)處理t1，t2=尺寸(b)；如果T2=1，b=b。%將b轉(zhuǎn)置到列向量目標%調(diào)用線性求解器x，fval，exitflag，output，lambda=linprog(f，A，b，Aeq，beq，lb，ub，x0，選項)；如果exitflag=0。%如果線

2、性規(guī)劃失敗，這個解決方案也失敗。返回目標%獲取帶有整數(shù)約束的決策變量的序列號v1=find(if int=1)；%整數(shù)變量的索引tmp=x(v1)；%整數(shù)約束的決策變量的當前值如果是空的(tmp)，%沒有整數(shù)約束，它是一個一般的線性規(guī)劃，可以直接返回返回目標v2=查找(check int(tmp)=0)；%正在尋找非整數(shù)的索引如果是空的(v2)，%如果整數(shù)約束決策變量都是整數(shù)，調(diào)用結(jié)束返回目標% kth決策變量不是整數(shù)解%注意先處理不滿足整數(shù)約束的第一個決策變量k=v1(v2(1)；%分支1:左分支tmp1=，長度(f)；%線性約束的系數(shù)向量tmp 1(k)=1；低=地板(x(k)；線性規(guī)劃不

3、等式約束的系數(shù)矩陣A實際上是以%thisA分支命名的線性規(guī)劃的不等式約束向量b實際上是在% thisb分支之后調(diào)用的if ifrowinmat(tmp1，低，A，b)=1%如果分支的約束已存在于舊的A和B中，則約束將不會更改thisA=A；thisb=b其他thisA=A；tmp 1；thisb=bthisb(end 1)=低；目標%disp(fenzhi1)，thisA，thisb%遞歸調(diào)用x1，fval1，exitflag1，output1，lambda1=.linprogdis(ifint，f，thisA，thisb，Aeq，beq，lb，ub，x0，options)；% branch

4、2:右分支tmp2=，長度(f)；%tmp1。tmp 2(k)=-1；高=-上限(x(k)；if ifrowinmat(tmp2，高，A，b)=1thisA=A；thisb=b其他thisA=A；tmp 2；thisb=bthisb(end 1)=高；目標%disp(fenzhi2)，thisA，thisbx2，fval2，exitflag2，output2，lambda2=.linprogdis(ifint，f，thisA，thisb，Aeq，beq，lb，ub，x0，options)；if is empty(v2)(exit flag 10 exit flag 2=0 f val=f va

5、l 1).|(exit flag 20 exit flag 1=0 f val=f val 2).|(exit flag 10 exit flag 20 f val=f val 1 f val=f val 2)，disp(錯誤調(diào)用)返回% is空(v2)表示所有整數(shù)目標%根據(jù)返回標志exitflag確定最終的最佳解決方案e 1如果exit flag 10 exit flag 2=0%左分支有最優(yōu)解，則右分支沒有最優(yōu)解x=x1fval=fval1。exitflag=exitflag1output=output1=1；e 2否則，如果exit flag 20 exit flag 1=0%右分支有最優(yōu)解，但左分