[数学建模资料集]选址问题.doc

其它优化问题： 例 某公司有6个建筑工地要开工，每个工地的位置（用平面坐标系表示，距离单位：千米 ）及水泥日用量(吨)由下表给出。 1 2 3 4 5 6 a 1.25 8.75 0.5 5.75 3 7.25 b 1.25 0.75 4.75 5 6.5 7.25 d 3 5 4 7 6 11 目前有两个临时料场位于A(5,1)，B(2,7)，日储量各有20吨。假设从料场到工地之间均有直线道路相连。 问题（1）：试制定每天的供应计划，即从A，B两料场分别向各工地运送多少吨水泥，使总的吨千米数最小； 问题（2）：为了进一步减少吨千米数，打算舍弃两个临时料场，改建两个新的，日储量各为20吨，问应建在何处，节省的吨千米数有多大？ 分别用matlab和Lingo写出对应的求解程序。 解 记工地的位置为(ai, bi)，水泥日用量为di，i=1,…,6; 料场位置为(xj, yj)，日储量为ej，j=1,2；从料场j向工地i的运送量为Xij。即 工地 料场 i=1 i=2 i=3 i=4 i=5 i=6 j=1 X11 X21 X31 X41 X51 X61 j=2 X12 X22 X32 X42 X52 X62 d 3 5 4 7 6 11 则模型 设X11=X1，X21= X 2，X31= X3，X41= X4，X51= X5, X61= X6；X12= X7，X22= X8，X32= X9，X42= X10，X52= X11，X62= X12， 而工地位置，i=1,…,6和料场位置，j=1,2均已知。 问题一的Matlab编程如下： % 场地运输问题的线性规划模型 % clear all clc tic %----------------------------------- % 构造目标函数系数向量与不等式和等式约束矩阵等 % 构造目标函数系数向量 a=[1.25 8.75 0.5 5.75 3 7.25]; b=[1.25 0.75 4.75 5 6.5 7.75]; d=[3 5 4 7 6 11]; x=[5 2]; y=[1 7]; e=[20 20]; for i=1:6 for j=1:2 aa(i,j)=sqrt((x(j)-a(i))^2+(y(j)-b(i))^2); end end CC=[aa(:,1); aa(:,2)]; % 构造不等式和等式约束矩阵 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=[d(1);d(2);d(3);d(4);d(5);d(6)]; % 变量的上下界 VLB=[0 0 0 0 0 0 0 0 0 0 0 0]; VUB=[]; x0=[1 2 3 0 1 0 0 1 0 1 0 1]; % 参数设置---options options=optimset(LargeScale,off); % 函数调用 [xx,fval,EXITFLAG,output]=linprog(CC,A,B,Aeq,beq,VLB,VUB,x0,options); % 结果输出 x=reshape(xx,6,2), fval, EXITFLAG, output % 注意 T=toc 问题二的模型是形式一样，所不同的是，其中料场位置，j=1,2也是变量。 设X11=X1，X21= X 2，X31= X3，X41= X4，X51= X5, X61= X6；X12= X7，X22= X8，X32= X9，X42= X10，X52= X11，X62= X12，和料场位置变量 x1=X13, y1=X14； x2=X15, y2=X16 问题二的Matlab编程如下： 1）函数调用主程序 % 场地运输问题的非线性规划模型 % clear all clc tic %-------------------------- % 初始可行解与替换 % x0=[3 5 0 7 0 1 0 0 4 0 6 10 5 1 2 7]; % x0=[ 3.0000 5.0000 0.0707 7.0000 0 0.9293