数学建模垃圾分类处理

1、垃圾分类处理与清运方案设计方案1 问题的重述在垃圾分类收集与处理中，不同类的垃圾有不同的处理方式，简述如下：1）橱余垃圾可以使用脱水干燥处理装置，处理后的干物质运送饲料加工厂做原料。不同处理规模的设备成本和运行成本（分大型和小型）见附录1 说明。2） 可 回收垃圾 将收集后分类再利用。3） 有 害垃圾 ，运送到固废处理中心集中处理。4） 其他不可回收垃圾将运送到填埋场或焚烧场处理。所有垃圾将从小区运送到附近的转运站，再运送到少数几个垃圾处理中心。显然，1）和2） 两项中， 经过处理，回收和利用，产生经济效益，而 3） 和4） 只有消耗处理费用，不产生经济效益。本项研究课题旨在为深圳市的垃圾分类

2、化进程作出贡献。为此请你们运用数学建模方法对深圳市南山区的分类化垃圾的实现做一些研究，具体的研究目标是：1） 假定现有垃圾转运站规模与位置不变条件下，给出大、小型设备（橱余垃圾）的分布设计，同时在目前的运输装备条件下给出清运路线的具体方案。以期达到最佳经济效益和环保效果。2） 假设转运站允许重新设计，请为问题1）的目标重新设计。2 基本假设1）假设各小区清运站每天的垃圾量是不变的；2）假设各小区清运站的垃圾都必须在当天清理完毕；3）不考虑运输车在行驶过程中出现的塞车、抛锚等耽误时间的情况；4）不允许运输车有超载现象；5)每个小区清运站均位于街道旁，保证运输车行驶顺畅；6)南山区人口分为不同部分

3、，每部分人口固定，每天产生垃圾量固定；7)一天只从小区清运站收一次垃圾(晚上或下午)；8)所有运输车均从垃圾转运站发车最后回到垃圾转运站；9)运输车将垃圾一起送往大型设备处和小型设备处再前往坟埋场和焚烧场；10)大型垃圾处理厂的寿命是30 年。小型垃圾处理机的寿命是10 年；(11) 建设在运输垃圾过程中没有新垃圾入站。3 符号(参数)说明(1) Xj (j=1, 2,，k)为第j个解释变量；(2) j(j=1, 2,，k)为第j个未知参数；(3) 为随机误差项；(4) S为多元线性回归模型的精度；(5) Pi(xi,yi)为第i个转运站的坐标；(6) Pj(Xj,Yj沟大型厨余垃圾处理设备建

4、在地图上的坐标；(7) cost1 为大型垃圾处理设备每日垃圾处理费用；(8) Cost2为小型垃圾处理设备每日垃圾处理费用；(9) |A| 表示A点到原点的距离，恒正(10) |B| 表示B点到原点的距离，恒正(11) |A-B|表示A,B两点之间的距离，恒正(12) Ta 表示 A 点所在地的垃圾量(13) Tb 表示A 点所在地的垃圾量(14) cost:耗油量；(15) T为规划使用年限；(16) Cik为第i座收集站运往第k座中转站单位运输量单位距离的费用(元t- 1 km- 1 );(17) Xik为第i座收集站运往第k座中转站的日运输垃圾量(t - d- 1 );(18) Lik

5、为第i座收集站运往第k座中转站运输距离(km);(19) Dk j为第k座中站运往第j座处理场单位运输量单位距离的费用(元 t- 1 km-1 )；(20) Yk j为第k座中转站运往第j座处理场日运输垃圾量(t d- 1 );(21) Sk j为第k座中转站运往第j座处理场运输距离(km);(22) Fk为规划期内待建中转站的固定投资(元)；(23) E为中转站的运行成本(元 t- 1 );(24) Qmin为中转站建设的最小控制规模(t d- 1 );(25) Qmax为中转站建设的最大控制规模(t d-1);5模型的构建与求解问题一的建模与求解5.1.1 深圳市城市生活垃圾产生量的预测表

6、一 深圳市城镇垃圾产生量历年统计表(万吨)年份20012002200320042005垃圾量296302年份20062007200820092010垃圾量321357413假定被解释变量Y,与多个解释变量Xi, X2, X3,Xk。之间具有线性关系，即Y 01X1 2X2 LkXk(8)其中Xj (j=1, 2,k)为k个解释变量，j(j=1, 2,k)为k+1个未知参数,为随机误差项被解释变量Y的期望值与解释变量Xi, X2, X3,Xk的线性方程为:(9)E Y 01X12X3 L kXk01123k对于n组观测值Yi,Xii ,X2i,Xki （i=1,2，n）,其方程组形式为:其矩阵形

7、式为其中Yi01X 1iY=XB + NXnk12X2ikX kii,(i 1,2,Ln)(10)11M1Y1Y2LLYnX112X12LkXk11X122 X22LkXk22X1n2X2nLkXknn10101Y1Y211X11X12X21X22LLXk1Xk20112MMMMMM2MMYn1X1nX2nLXknknYnY1Y2MYn为被解释变量的观测值向量；1112MX1n2122MX2nLLMLk1k2MXkn为被解释变量的观测值矩阵；(11)112 为总体回归参数向量；总体回归方程为:2为随机误差向量M(12)可采用最小二乘法对上式中的待估回归系数1, 2,L n进行估计，求得 值后，

8、即E (Y) =X0可利用多元线性回归模型进行预测了 。我们对多元线性回归分析进行数学检验，包括回归方程和回归系数的显著性检验。a.回归方程的显著性检验，采用统计量：(13)U /mQ/ n m 1式中；Un2Yj Yj为回归平方和，其自由度为m; j 1Yj2Yj为剩余平方和，其自由度为(n-m-1)。利用上式计算出F值后，再利用F分布表进行检验。给定显著性水平a ,在 吩布表 中查出自由度为m和(nm 1)的值Fa ,如果F Fa,则说明Y与X1，X2,L Xm的线性相 关密切；反之，则说明两者线性关系不密切。b,回归系数的显著性检验，采用统计量:2(14)b /Ci iiiF tQ/ n

9、 m 1式中，Gi为相关矩阵C A 1的对角线上的元素。对于给定的置信水平 ，查F分布表得Fi n m 1 ,若计算值Fi Fa ,则拒绝原 假设,即认为Xi是重要变量，反之，则认为 Xi ,变量可以剔除多元线性回归模型的精度，可以利用剩余标准差(15)S . Q/(n-m-1)来衡量。S越小，则用回归方程预测Y越精确；反之亦然。采用matlab软件编程进行城市生活垃圾量多元线性回归模型预测(预测代码见附录1)。表二为训练结束后预测值与统计值的对比表，精度达到要求后用训练好的模型来预测深圳市2011-2015年城市生活垃圾产生量，预测结果见表 311。在matlab软件中运行代码后得到生活垃圾

10、产生量的回归方程为：Y=38. 7965+0. 25178xXl+0. 10508xx2- 0. 0574xx3+O. 1292xx4-0. 0138xx5+20. 8016xx6-0. 0095xx7+0. 0066xxs 3. 1460xx9方差估计：S=25. 7642回归方程的显著性检验F统计量，F=72. 3187,所以拒绝假设，即回归模型成立。表二线性回归模型预测值与统计值对比表年份20012002200320042005预测值统计值年份20062007200820092010预测值统计值413图一线性回归模型预测值与统计值对比分析图系列1系科2从表二及图一可以看出，多元线性回归模

11、型对历史值的拟合程度较高，预测精度是 可以接受的，多元线性回归模型预测值比较接近深圳市城市生活垃圾实际产生量，稍微 偏局。表三 2011-2015年深圳市城市生活垃圾产生量多元线性回归模型预测值年份20112012201320142015预测值（万吨）5.1.2 大小型厨余垃圾设备规划5.1.2.1 模型的建立题目要求给出大、小型设备（橱余垃圾）的分布设计。由于大型厨余垃圾处理设备 处理能力为200吨/日，投资额约为4500万元，运行成本为150元/吨。而每个转运站 的垃圾数量有限，所以大型厨余垃圾处理设备必须在图上重新选址建设。小型餐厨垃圾 处理机，处理能力为200-300公斤/日，投资额约

12、为28万元，运行成本为200元/吨。所 以小型垃圾处理机可以设置在垃圾中转站内。根据表四用编程作图二（程序见附录三）表四中转站坐标名称中转站厨余垃Xy名称中转站厨余垃xy圾量圾量P1麻勘站10P19前海公园站16P2阳光站10P20玉泉站25P3书,8P21九街站20P4大石獭站30P22大冲站35P5牛城村站5P23沙河巾场站30P6长源村站5P24大新小学站30P7福光站10P25涌下村站20P8塘朗站10P26深圳大学站15P9动物园站20P27白石洲南站30P10平山村站25P28北头站15P11新围村站20P29科技园站20P12月亮湾大道站40P30南园站15P13西丽路站15P3

13、1南山村站25P14光前站20P32花果路30P15同乐村站5P33官龙村站15P1615P34望海路站30P17松坪山站25P35华侨城站70P18松坪山（二）站10P36南光站15P37南山市场25P38疏港小区站40图二中转站坐标图33 怕15 14 舞10 03ifi&101314伯从图表可知每个垃圾转运站的坐标Pi（xi,yi）,假设大型厨余垃圾处理设备建在地图上的 Pj（Xj,Yj）所以对于每个垃圾中转站来说有两种情况：（1）在站内设置垃圾处理机。（2）把垃圾运往大型厨余垃圾处理厂进行处理。从中选择最优方案，从而确定垃圾大型垃圾处理站的位置。假设大型垃圾处理厂的寿命是 30年。小型

14、垃圾处理机的寿命是 10年。大型垃圾处理设备的平均每吨耗损成本 =/ (30*365*200) =20元/吨小型垃圾处理设备的平均每吨耗损成本 =280000/ (10*365*) =256元/吨吨汽车，每车耗油20L 35L 70#气油/百公里。每升70#汽油价格为元司机月薪平均3500元。如果运往大型垃圾处理设备厂，则每日垃圾处理费用(costl)=平均每日设备耗损成本+运输费用+司机工资+垃圾处理费用。如果在垃圾转运站设置小型垃圾处理机，则每日垃圾处理费用(cost2)= 平均每日设备耗损成本+垃圾处理费用。要确定大型垃圾处理厂的位置，需要计算出选择第1种方案的点。根据以上条件建立模型：

15、Cost1=7.2 35 J(x(j) x(i)2 (y(j) y(i)2 2 (s(i) 2.5) 3500 30 (150 20) s(i)Cost2=256 200) s(i)Cost1x2, y1y2,不难看出A在B的后方，即A比B远。对于前方参考点O,要将 A,B对应垃圾点的垃圾全部取回再返回O,一共有三种方式：1. O-A-O, O-B-O单独运输。这种情况下，总的路程消费等于空载运行费用(20L/百公里)与装载时运行费用(20+6*Ta) L/百公里)的总和。于是有：Cost = 20*|A| + (20+6*Ta) *|A| + 20*|B| +(20+6*Ta) *|B|2.

16、 O-A-B-O先远点再近点，即先空载至最远处，装完 A点垃圾后再返回至B,再回O点，有：Cost = 20*|A| + (20+6*Ta) *|A-B| + (20+6*(Ta+Tb)*|B|=20*|A| + *|A|*Ta + *|B|*Tb3. O-B-A-O先近点在远点，即先装B点垃圾，然后载着B点的垃圾奔至A点，再回O点，有:Cost= 20*|B| + (20+6*Tb) *|A-B| +(20+6*(Ta+Tb) *|A|=20*|B| + (20+6*Ta)*|A| + (20+6*Tb) *|B| +(20+6*Tb) *|A-B|*2比较以上三种情况，远近点的遍历顺序，可

17、以看出，“先远后近”绝对比“先近后远在花费钱的数量上要少的多，省出(20+6*Tb) *|A-B|*2 这部分的钱主要是车载着 B点 的垃圾奔到A点再返回B点。而又注意到两者的时间花费是相等的。所以在其余同等的情况下选择“先远后近”。考虑单独运输比其余的两种运输花费的钱仍不比“先远后近”省，还多了 20*|B| ,所以一般情况下，不采用单独运输。A, B两点没有明显先后顺序。-并邻状态(如图三)图三还是一共有三种情况:1, O-A-O, O-B-O单独运输。这种情况下，跟 A,B两点有先后顺序中的情况完全相同，即有：Cost = 20*|A| +(20+6*Ta) *|A| + 20*|B|

18、+(20+6*Ta) *|B|2, O-A-B-OCost = 20*|A| + (20+6*Ta)*|A-B| + (20+6*(Ta+Tb)*|B|- 1B-A-OCost = 20*|B| +(20+6*Tb)*|A-B|+(20+6*(Ta+Tb)*|A|- 2相比之下，清晰可见并邻状态下的单独运输所花的费用最少，所以在不要求时间的情况下对于并邻两点，采用单独运输的方式最节约钱。用式与式相减，得到如下判断式：6*(Ta-Tb) *|A-B|+6*(Ta+Tb)* (|B|-|A| ) -上式 A-B-O;上式 0 时，选 O-B-A-O;上式=0时，任意选上述两路线。两点选择趋势的讨论

19、。(如图四)图四A由图中看到B,C两点没有明显的先后顺序，属于并邻点。因为当运输车载重行驶时 费用会成倍的增长，比其空载时所花费用要大的多，所以排除A-B-C或A-C-B这样的一次经过3点的往返路线，仅选择B,C中的某一点与A完成此次运输，将另一点留到 下次。那么A点选择B还是C呢不妨假设旧|C|,即B点离原点的距离比C点的更远，因为A在B,C之后，所以也 就是B点离A点更近。这样，此次的运输我们更趋向于选择 A-B因为就这三点而论，A 无论是选B还是C,三点的垃圾总要运完，所以花费的钱是一样的。 但选择A-B后，下 次运输车运C点垃圾时就无需跑的更远。综上所述，得出搜索的基本原则：1 .在两

20、点递减的情况下，不采用单独运输；2 .在其余同等的情况下选择“先远后近”；3 .不要求时间的情况下对于并邻两点，采用单独运输的方式最节约钱；一般情况下用式3作判断；4 .车在装的足够多的情况下应该直接返回中转站；5 .每一次布局和每条线路的搜索不妨由剩下未搜点中的最大值开始四.关于垃圾点的垃圾是否一次清除的讨论这里说的一次消除问题不是指一天，而是指当一辆运输车已经装载了足够多的垃 圾，不能完全清理下一个垃圾点的时候， 车在下一个站点“停还是不停”的问题。例如， 一辆运输车选择了某段路线后，当清运完前几个点后，未达到饱和，但下一个点的垃圾 量又装不完，那么此车是直接返回呢，还是继续装直至车装满为

21、止呢我们判断前者更好，就是车在装的足够多的情况下应该直接返回原点。这是因为对 于下一垃圾点（假设为A点）内的垃圾而言，无论是一次装完还是分两次装完，将它们运 回所花费用是恒定的。整体而言，两者花费的钱是相等的，但分两次装要多花装车时间, 所以选择前者。5.1.3.2模型的求解首先根据题所给的数据画出散点图图五求得总耗油为x,求解程序如附录二，运输车的最优路线如下图所示：表五:线路的站点序列和油耗站点序号油耗T线0-30-29-27-3-0二号线0-28-26-32-25-5-0三号线0-36-23-33-21-0四号线0-24-18-35-15-0五号线0-34-17-16-2-012六号线0-20-11-10-0七号线0-19-13-8-0八号线0-14-7-4-1-0九号线0-22-0州线0-12-9-08线0-31-6-06模型的评价和推广(1)模型的评价好的方面：1、用多元线性回归模型预测城市垃圾的产生量，其计算参数少、计算过程简单， 且结果通过检验都比较接近真实值，2、针对垃圾收运系统的特点，引入逆向