供应与选址问题

实验九学习掌握MATLAB6个建筑工地要开工，每个工地的位置（a，b表示，距离单（5,1（2,7）AB20AB各设在何处，AB两料场分别向各工地运送水泥总的吨千表 工地的位置（a，b）及各工地水泥日用纵坐标d（吨

fgi(x) hj(x)＝0j＝1，2，…， xRnf(x为目标函数g(xh(x为约束函数 S＝｛x｜gi(x) i＝1，2，…，m；hj(x)＝0j＝1，2，…，k S为可行集或可行域S中的点称为可行点。这样（1）

f(x)，s.

x成立f(x≥f(x*x*为极小化问题的局部最优解。=0，f(xgi(xx*gi(x（iI）x*hj(xx*连续可微，且有m和1，…，k ▽f(x*)＋▽i

(x*)＋▽j

(x*)

i

(x*)

jh

(x*)K－T；用无约束最优化方法去求解。这类方法称为序列无约束最小化方法SUMT（其一SUMTSUMT内点法。xxconstr(funx0 求解非线性规划模型（1；xconstr(funx0options 参数options1x=constr('fun',x0,options,vlb,vub) [x,options]=constr('fun',x0,...) 同上，同时返回参数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来查阅有关

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 函数形式（nonlcon.m； [c,ceq]=c=c(x)；ceq=ceq(x可以在MATLAB帮助文件中查阅有关该命令的更多用法。记工地的位置为（ai，bidi，i＝1，…，6；料场位置为（xjyj，日储量为ej，j＝1，2；从料场j向工地i的运送量为cij。这个优化问题的目标函数（总 (xa)2((xa)2(yb) j1icij＝di，i jcij≤ej，j ixjyjAB两料场往各工地的运送量cij，问题归结（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往各工地的运送量(x(13),x(14))Afor f2Bx(7)x(12)B往各工地的运送量,(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是该季生产的台数。若交货后有剩余，可用于下季度的交货，但