13的課后習(xí)題解答及電子課件07第7章目標(biāo)

1、第 7 章目標(biāo)規(guī)劃習(xí)題解答7.1解 使用 clc, clear a=1 1; b=7;求解的程序如下c1=3 1; c2=1 2;lb=zeros(2,1); ub=5 5; %定義線性規(guī)劃的下界向量和上界向量 x1,y1=linprog(-c1,a,b,lb,ub) %求解第一個線性規(guī)劃x2,y2=linprog(-c2,a,b,lb,ub) %求解第二個線性規(guī)劃go=-y1,y2; %工具箱線性規(guī)劃是求解極小值goalfun=(x)-3*x(1)-x(2);-x(1)-2*x(2);x,fval=fgoalattain(goalfun,rand(2,1),go,abs(go),a,b,lb

2、,ub); x, fval=-fval %恢復(fù)到原來的目標(biāo)函數(shù)求得的滿意解為 x1 3.7561， x2 3.2439 ，對應(yīng)的目標(biāo)函數(shù)的值 z1 14.5122 ，解 設(shè)安排商業(yè)時間 x1 小時，時間 x2 小時，音樂時間 x3 小時，7.2該問題的目標(biāo)規(guī)劃模型為min z p (d d d d d ) p d 1112232 4 d d 12,311 d d 2.4,x122x d d 1,233250 d d 36000,344, d , d , d , d , d , d , d , d 0.311223344其中36000為每天收入的上限 25060 2.4 36000求解的Ling

3、o程序如下m:sets:level/1.2/:p,z,goal; variable/1.3/:x; s_con_num/1.4/:g,dplus,dminus; s_con(s_con_num,variable):c;obj(level,s_con_num)/1 1,1 2,1 3,2 4/:wplus,wminus; endsetsdata:ctr=?; goal=?0;g=12 2.4 1 36000;c=1 1 1 1 0 0 0 1 0 250 -40 -17.5;wplus=1 1 0 0;wminus=1 1 1 1; enddata min=sum(level:p*z); p(c

4、tr)=1;。-47-for(level(i)|i#ne#ctr:p(i)=0); for(level(i):z(i)=sum(obj(i,j):wplus(i,j)*dplus(j)+wminus(i,j)*dminus(j);for(s_con_num(i):sum(variable(j):c(i,j)*x(j)+dminus(i)-dplus(i)=g(i); for(level(i)|i #lt# size(level):bnd(0,z(i),goal(i);end求得 x 2.4 ， x 1， x 8.6 ， d 35590.5 ， d 0 ，每天的純收入為1234436000 d

5、409.5。47.3解 設(shè)在允許的工作時間內(nèi)產(chǎn)品I生產(chǎn) x1 件，產(chǎn)品II生產(chǎn) x2 件；在加班時間內(nèi)產(chǎn)品I生產(chǎn) x3 件，產(chǎn)品II生產(chǎn) x4 件。建立如下的目標(biāo)規(guī)劃模型min p (d d ) pd 1122 33x 2.5x d 120,121 2.5x d 160,3342 7x d 640,10343 d 0, i 1,2,3; x 0且為整數(shù), i 1,2,3,4. ii其中第3個約束右邊的640為利潤的上界，由于無論生產(chǎn)產(chǎn)品I或II，每小時的贏利不超過4元，每周的生產(chǎn)時間不超過160小時，因而最大利潤不超過640元。計算的Lingo程序如下m:sets:level/1.2/:p,z

6、,goal; variable/1.4/:x; s_con_num/1.3/:g,dminus; s_con(s_con_num,variable):c;obj(level,s_con_num)/1 1,1 2,2 3/:wminus; endsetsdata:ctr=?; goal=?0;g=120 160 640;c=3 2.5 0 0 3 2.5 3 2.5 10 8 8.5 7;wminus=1 1 1; enddata min=sum(level:p*z); p(ctr)=1;for(level(i)|i#ne#ctr:p(i)=0); for(level(i):z(i)=sum(obj(i,j):wminus(i,j)*dminus(j);for(s_con_num(i):sum(variable(j):c(i,j)*x(j)+dminus(i)=g(i); for(level(i)|i #lt# size(level):bnd(0,z(i),goal(i); for(variable:gin(x);end-48-