MATLAB編程0-1規(guī)劃問(wèn)題

1、MATLAB語(yǔ)言應(yīng)用最優(yōu)化MATLAB編程線性規(guī)劃問(wèn)題第二章0-1規(guī)劃MATLAB的0-1規(guī)劃函數(shù)bintprog是針對(duì)下述0-1規(guī)劃：min z = / * xs.t A*xb, (2. 1) aeq *x = beqx =(X,兀， , xi = Oorl, j = 1,2, 解0-1規(guī)劃(2.1)的0-1規(guī)劃函數(shù)bintprog表述為lx,fv,exitflag,output =bintprog (f, A, b, aeq,beq)(2.2)輸入部分：f為目標(biāo)函數(shù)，實(shí)為目標(biāo)函數(shù)的系數(shù)。A為(2.1)中的不等式約朿矩陣b為(2.1)中的不等式約束向量aeq為(2.1)中的等式約朿矩陣beq

2、 (2.1)中的等式約朿向量輸出部分：x為最優(yōu)解fval為最優(yōu)值exitflag為輸出標(biāo)志exitflag二1,有最優(yōu)解exitflag二0,迭代次數(shù)超過(guò)設(shè)定次數(shù)exitflag-2,約束區(qū)域不可行exitflag二-3,問(wèn)題無(wú)解output,表明算法和迭代情況如果我們不需要了解迭代情況和存儲(chǔ)情況，可將0-1規(guī)劃函數(shù)bintprog寫(xiě)成x,fv,ex二1 inprog (f, A, b, aeq,beq)(2.3)在函數(shù)bintprog中，輸入或輸出元素的符號(hào)可以變更，如(2.3)中ex仍為輸出標(biāo)志，但元素的符 號(hào)位置不能變更。在輸岀部分，如有缺者，可用號(hào)代替。函數(shù)bintprog的使用要點(diǎn)與

3、函數(shù)linprog的使用要點(diǎn)相同。函數(shù)是為求目標(biāo)函數(shù)的最小值而設(shè)垃的，如要求函數(shù)的最大值，可先求岀(-/)的最小值，則-沖 必為/的最大值。例一 用函數(shù)bintprog求解下列0-1規(guī)劃max z = x + 2x2 + 2x3 - 6x4 一 4x5s.t 3西+2吃一尤3+兀+5x? S52x + 4x2 一 2x3 -x4 - 2x5 1x 一 x2 + 6X3 + 4x4 85x)+ 3x2 + x4 5兀=0或l(i=l,2,4)用MATLAB語(yǔ)言編程如下：f二3, 7,-1,1A二-2,1, -1, 1 ;T, 1, -6, -4;A=A;-5, -3, 0, -1b二-1, &

4、5;x, fv, ex=bintprog (f, A, b ,);求解如下 ex二 1, fv=-3, x=l, 0, 1, 1例三 用函數(shù)bintprog求解下列0-1規(guī)劃 min z = 15X| +18x2 +21x3 + 19x4 +23% +22x6 +26曲 +17xs + 16禺S/ %, + x2 + x3 = 1x4+x5+x6 = 1x7 + x8 + XK) = 1x+x4+ x7 = 1x2+x5+ 屯=1x3+x6+xk) =1兀=0 或 1J = 1,2.9用MATLAB語(yǔ)言編程如下：f二15,18, 21, 19, 23, 22, 26, 17, 16;z2=zeros(1, 2); z3=zeros(1, 3);z6二zeros (1, 6) ; o3=zeros (1, 3);q=o3, z6;z3, o3, z3;z6, o3;q=q；l, z2, 1, z2, 1, z2;q二q;0, 1, z2, 1, z2, 1, 0;q=q;z2, 1, z2, 1, z2, 1;bq 二 ones (6, 1)Lx, fv, ex=bintprog(f