《Matlab供应与选址问题附详细编程PPT学习教案》由会员分享，可在线阅读，更多相关《Matlab供应与选址问题附详细编程PPT学习教案（10页珍藏版）》请在装配图网上搜索。
1、会计学1Matlab供应与选址问题附详细编程供应与选址问题附详细编程（一）、建立模型 记工地的位置为(ai，bi)，水泥日用量为di，i=1,6;料场位置为(xj，yj)，日储量为ej，j=1,2；从料场j向工地i的运送量为Xij。目标函数为：216122)()(minjiijijijbyaxXf 约束条件为：2 , 1 ,6 , 2 , 1 ,6121jeXidXjiijijij 当用临时料场时决策变量为：Xij，当不用临时料场时决策变量为：Xij，xj，yj。第1页/共10页（二）使用临时料场的情形 使用两个临时料场A(5,1)，B(2,7).求从料场j向工地i的运送量为Xij，在各工地用
2、量必须满足和各料场运送量不超过日储量的条件下，使总的吨千米数最小，这是线性规划问题. 线性规划模型为：2161),(minjiijXjiaaf2 , 1 , 6 , 2 , 1 , s.t.6121jeXidXjiijijij其中 22)()(),(ijijbyaxjiaa，i=1,2,6,j=1,2,为常数。 设X11=X1, X21= X 2, X31= X 3, X41= X 4, X51= X 5, X61= X 6X21= X 7, X22= X 8, X32= X 9, X42= X 10, X52= X 11, X62= X 12 编写程序第2页/共10页其详细的程序为：c1=s
3、qrt(5-1.25)2+(1-1.25)2);c2=sqrt(5-8.75)2+(1-0.75)2);c3=sqrt(5-0.5)2+(1-4.75)2);c4=sqrt(5-5.75)2+(1-5)2);c5=sqrt(5-3)2+(1-6.5)2);c6=sqrt(5-7.25)2+(1-7.25)2);c7=sqrt(2-1.25)2+(7-1.25)2);c8=sqrt(2-8.75)2+(7-0.75)2);c9=sqrt(2-0.5)2+(7-4.75)2);c10=sqrt(2-5.75)2+(7-5)2);c11=sqrt(2-3)2+(7-6.5)2);c12=sqrt(2
4、-7.25)2+(7-7.25)2);c=c1;c2;c3;c4;c5;c6;c7;c8;c9;c10;c11;c12;A=1 1 1 1 1 1 0 0 0 0 0 0;0 0 0 0 0 0 1 1 1 1 1 1;B=20;20;Aeq=1 0 0 0 0 0 1 0 0 0 0 0;0 1 0 0 0 0 0 1 0 0 0 0;0 0 1 0 0 0 0 0 1 0 0 0; 0 0 0 1 0 0 0 0 0 1 0 0;0 0 0 0 1 0 0 0 0 0 1 0;0 0 0 0 0 1 0 0 0 0 0 1;Beq=3;5;4;7;6;11;vlb=0;0;0;0;0;0;
5、0;0;0;0;0;0;x,f=linprog(c,A,B,Aeq,Beq,vlb)第3页/共10页计算结果为：x = 3.0000 5.0000 0.0000 7.0000 0.0000 1.0000 0.0000 0.0000 4.0000 0.0000 6.0000 10.0000fval = 第4页/共10页（三）改建两个新料场的情形 改建两个新料场，要同时确定料场的位置(xj,yj)和运送量Xij，在同样条件下使总吨千米数最小。这是非线性规划问题。非线性规划模型为：216122)()(minjiijijijbyaxXf2 , 1 , 6 , 2 , 1 , . .6121jeXidX
6、tsjiijijij第5页/共10页设 X11=X1, X21= X 2, X31= X 3, X41= X 4, X51= X 5, X61= X 6 X21= X 7, X22= X 8, X32= X 9, X42= X 10, X52= X 11, X62= X 12 x1=X13, y1=X14, x2=X15, y2=X16 （1）先编写M文件定义目标函数。function y=liaoch(x)y=x(1)*sqrt(x(13)-1.25)2+(x(14)-1.25)2)+x(2)*sqrt(x(13)-8.75)2+(x(14)-0.75)2)+x(3)*sqrt(x(13)-
7、0.5)2+(x(14)-4.75)2)+x(4)*sqrt(x(13)-5.75)2+(x(14)-5)2)+x(5)*sqrt(x(13)-3)2+(x(14)-6.5)2)+x(6)*sqrt(x(13)-7.25)2+(x(14)-7.25)2)+x(7)*sqrt(x(15)-1.25)2+(x(16)-1.25)2)+x(8)*sqrt(x(15)-8.75)2+(x(16)-0.75)2)+x(9)*sqrt(x(15)-0.5)2+(x(16)-4.75)2)+x(10)*sqrt(x(15)-5.75)2+(x(16)-5)2)+x(11)*sqrt(x(15)-3)2+(x
8、(16)-6.5)2)+x(12)*sqrt(x(15)-7.25)2+(x(16)-7.25)2);(2) 取初值为线性规划的计算结果及临时料场的坐标: x0=3 5 0 7 0 1 0 0 4 0 6 10 5 1 2 7;编写主程序gying2.m.第6页/共10页(3) 计算结果为：x=3.0000 5.0000 4.0000 7.0000 1.0000 0 0 0 0 0 5.0000 11.0000 5.6959 4.9285 7.2500 7.7500注意初始点的选取第7页/共10页 某厂向用户提供发动机，合同规定，第一、二、三季度末分别交货40台、60台、80台每季度的生产费用为 （元），其中x是该季生产的台数若交货后有剩余，可用于下季度交货，但需支付存储费，每台每季度c元已知工厂每季度最大生产能力为100台，第一季度开始时无存货，设a=50、c=4，问工厂应如何安排生产计划，才能既满足合同又使总费用最低讨论a、b、c变化对计划的影响，并作出合理的解释练习 2bxaxxf第8页/共10页文件文件.即时commandwidow命令文件运行结果如下：第9页/共10页
Matlab供应与选址问题附详细编程PPT学习教案
