供應(yīng)與選址問題

實驗九學(xué)習(xí)掌握MATLAB6個建筑工地要開工，每個工地的位置（a，b表示，距離單（5,1（2,7）AB20AB各設(shè)在何處，AB兩料場分別向各工地運(yùn)送水泥總的噸千表 工地的位置（a，b）及各工地水泥日用縱坐標(biāo)d（噸

fgi(x) hj(x)＝0j＝1，2，…， xRnf(x為目標(biāo)函數(shù)g(xh(x為約束函數(shù) S＝｛x｜gi(x) i＝1，2，…，m；hj(x)＝0j＝1，2，…，k S為可行集或可行域S中的點(diǎn)稱為可行點(diǎn)。這樣（1）

f(x)，s.

x成立f(x≥f(x*x*為極小化問題的局部最優(yōu)解。=0，f(xgi(xx*gi(x（iI）x*hj(xx*連續(xù)可微，且有m和1，…，k ▽f(x*)＋▽i

(x*)＋▽j

(x*)

i

(x*)

jh

(x*)K－T；用無約束最優(yōu)化方法去求解。這類方法稱為序列無約束最小化方法SUMT（其一SUMTSUMT內(nèi)點(diǎn)法。xxconstr(funx0 求解非線性規(guī)劃模型（1；xconstr(funx0options 參數(shù)options1x=constr('fun',x0,options,vlb,vub) [x,options]=constr('fun',x0,...) 同上，同時返回參數(shù)options的值 [f,g]=fun(xf=f(x);g=[g1(x);g2(x);…;gm(x z＝1xTHx＋cT

xqp(HcAb xqp(HcAb （4；x=qp(H,c,A,b,vlb,vub) x=qp(H,c,A,b,vlb,vub,x0) 指定迭代的初始值x0；x=qp(H,cAb,vlb,vub,x0,n) nAxbn個約束條件等式約束；qplphelpqp來查閱有關(guān)

f

Ax≤b，c(x) lb≤x≤xxfmincon(funx0Ab （5xfmincon(funx0AbAeqbeq （5x=fmincon(fun,x0,A,b,Aeq,beq,lb,ub) （5x=fmincon(fun,x0,Ab,Aeq,beq,lb,ub,nonlcon 函數(shù)形式（nonlcon.m； [c,ceq]=c=c(x)；ceq=ceq(x可以在MATLAB幫助文件中查閱有關(guān)該命令的更多用法。記工地的位置為（ai，bidi，i＝1，…，6；料場位置為（xjyj，日儲量為ej，j＝1，2；從料場j向工地i的運(yùn)送量為cij。這個優(yōu)化問題的目標(biāo)函數(shù)（總 (xa)2((xa)2(yb) j1icij＝di，i jcij≤ej，j ixjyjAB兩料場往各工地的運(yùn)送量cij，問題歸結(jié)（7（6）fxjyj是非線性的，所以在求新建料場位置Mfunctiona=[1.25,8.75,0.5,5.75,3,b=[1.25,0.75,4.75,5,6.5,d=[3,5,4,7,6,e=[20,20]; f1Ax(1)x(6)A往各工地的運(yùn)送量(x(13),x(14))Afor f2Bx(7)x(12)B往各工地的運(yùn)送量,(x(15),x(16))B的位置forfor >>x0=[zeros(1,12)5127]';% >>>>工地ci>>text(1.25,1.25,'+3');text(8.75,0.75,'+5');>>text(5.75,5,'+7');text(3,6.5,'+6');>>>>+++6+4+3+5 （h）,x2800h，（108（1081.697(1/1

f(xaxbx2（元x是該季生產(chǎn)的臺數(shù)。若交貨后有剩余，可用于下季度的交貨，但