数学建模之垃圾处理.docx

资源描述

数学建模之垃圾处理.docx

《数学建模之垃圾处理.docx》由会员分享，可在线阅读，更多相关《数学建模之垃圾处理.docx（23页珍藏版）》请在冰豆网上搜索。

数学建模之垃圾处理.docx

数学建模之垃圾处理

城市生活垃圾管理问题研究

摘要

近年来，随着垃圾产量的日益增加，人们已经逐渐意识到它对生态环境及人类生存带来的极大威胁。

本文针对垃圾处理问题，先采用一元线性回归和最小二乘曲线拟合的方法，求出垃圾产量的预测模型，再采用图论法，得到垃圾最短收运路径以及最佳车辆分配方案。

对于第一问，我们根据题意找到影响垃圾产量的六个因素，查得相关数据后，运用一元线性回归的方法拟合出各个影响因素与年份的关系式，如下表：

生产总值（亿万元）

城市居民人均收入（元）

人均消费水平（元）

人口数量（万人）

人均住房面积（平方米）

城市燃气普及率（%）

接着，我们采用最小二乘曲线拟合的方法得到垃圾产量与各个影响因素的关系式，如下：

这样，在已知年份的条件下，可以通过各个影响因素的值，预测出垃圾的产量。

由于预测量考虑了实际中的各个影响因素，故具有准确性和较高的实用性。

对于第二问，我们经过数据预处理，画出以车库为原点的垃圾收集点、中转站分布图。

接着，根据题中垃圾车的最大装载量与垃圾站的分布特点将数据分成十二区域，用图论法在每个区域中找到最小生成树，为了避免垃圾收运车走重复路线，我们通过观察，将最小生成树的树叶融入树中，形成一条链，即为垃圾收运车的最短收运路线。

在得到12个区域的最短路径图后，我们将行驶时间、装车时间、卸车时间相加，得到在每个区域收集垃圾的总时间（h），如下表：

区域1

区域2

区域3

区域4

区域5

区域6

2.1054

2.035

1.9497

1.9664

1.7679

1.6508

区域7

区域8

区域9

区域10

区域11

区域12

1.6977

1.6683

1.7402

1.8511

1.8233

1.7050

由上表知：

车库每天至少需要3辆垃圾收运车收运垃圾。

故最佳的车辆分配方案为3辆垃圾收运车每辆每天前往4个区域收运垃圾。

运用以上方法得到的收运路线，不但满足题设条件（不超过垃圾车的最大装载量、日负载总量以及最多日收集点数），而且还能使垃圾的收运时间最短，另外该模型可以提出合理的车辆分配方案，提高了资源利用率。

因此，本模型具有较好的实用性和可靠性。

关键词垃圾预产量线性回归最小二乘曲线拟合图论法收运路线

1.问题的重述

由于人类生产和生活的不断发展而产生的垃圾对生态环境及人类生存带来极大的威胁已逐步成为重要的社会问题。

城市生活垃圾是居民生活、消费过程中产生的废弃物，其年增长速度达8-10%，因此导致城市垃圾的数量日益庞大，并且其组分复杂还处于不断变化中,使处理费用慢慢升高。

另一方面城市垃圾占用大量土地、污染水体、污染大气、破坏植被,严重影响城市的市容景观和居民的生活环境[1]。

城市垃圾已成为各国政府急需妥善解决的首要问题之一。

城市垃圾管理包括计划、组织、行政、金融、法律和工程等多方面，并涉及到城市生活垃圾收集、运输和处置。

而中国目前处置水平低，管理办法不多，更是急待解决的问题。

在这方面，世界许多国家在谋求解决城市生活垃圾过程中，产生出许多好的办法，并在此过程中总结了经验和教训。

一般认为城市生活垃圾的影响因素包括地理位置、人口、经济发展水平（生产总值）、居民收入以及消费水平、居民家庭能源结构等等。

城市生活垃圾产量是垃圾管理系统的关键参数，城市在建造垃圾填埋场之前,必须对该地区的未来垃圾产量进行合理预测。

若对城市垃圾产生量的估计过高,相应的填埋场库容设计必然增大,各种投资亦相应增加,将导致物力财力的巨大浪费;若对垃圾产量估计过低,在未达到填埋场设计使用年限时就将被迫关闭,不得不投资另行建造其他填埋场或对填埋场进行扩容,同样会导致物力财力的浪费[2]。

因此对未来某段时间内垃圾产量的准确预测是相关垃圾管理的部门做出管理规划的前提。

另外，城市垃圾自其产生到最终被送到处置场处理，需要环卫部门对其进行收集与运输，这一过程称为城市垃圾的收运。

收运过程可简述如下：

某城市有多个行政区，每个区内均有一个车库，假设某一车库拥有最大装载量为w的垃圾收集车k辆，并且该区的垃圾收集点（待收集垃圾的点）有n个，该城市共有垃圾中转站p座。

每天k辆垃圾车从车库出发，经过收集点收集垃圾，当垃圾负载达到最大装载量时，垃圾车运往中转站，在中转站卸下所有收运的垃圾，然后再出站收集垃圾，如此反复，直到所有收集点的垃圾都被收集完，垃圾车返回车库。

以上收运过程均在各点的工作区间之内完成。

（注：

必须在收集点的工作区间之内，垃圾车才能在该点收集垃圾。

）

请利用数学方法建立以下问题的数学模型并求解，对模型的结果做出合理分析和解释。

1.在查阅相关文献，搜集垃圾产量数据的基础上建立城市生活垃圾产量中短期预测模型，并且分析模型的准确性和实用性。

2.在收运过程已知下述

（1）

（2）（3）（4）等条件下，如何安排垃圾收运车的收运路线，使垃圾收运车的行车里程尽可能的少，或者垃圾收运时间尽可能短？

（1）车库和收集点、收集点与中转站、中转站与车库的距离；

（2）各收集点每天的垃圾产量；

（3）每辆垃圾收运车的最大载荷；

（4）垃圾收集点、车库、中转站的工作区间[a,b]。

请给出规划以上垃圾收运路线的数学模型，并设计出有效的算法，针对题中给出的数据，求解模型。

并且对模型的适用性、算法的稳定性和鲁棒性做出分析。

2.问题的分析

2.1.分析影响城市生活垃圾的因素

城市生活垃圾的影响因素决定了城市生活垃圾的总量，而这些影响因素包括城市规模、城市经济的发展程度、城市人口的多少及居民的收入、消费水平等诸多方面，故此处首先查阅相关资料[3]，了解以往在解决此类问题时的标准，再根据近几年我国居民消费水平，生活方式能源结构等方面的变化规律进行分析，最终得到对城市生活垃圾总量影响显著的几个因素如下：

（1）人口数量：

随着人口的不断增加以及城市人口流动的不断加剧，在人均垃圾产量稳定变化的基础上人口的变化必定造成城市生活垃圾产量的变化；

（2）生产总值：

生产总值的变化意味着居民生活水平的改变,居民生活水平的好坏直接关系到其生活方式的调整，不同的生活方式会对城市生活垃圾产量的大小造成影响；

（3）城市居民人均收入、人均消费水平：

是居民生活水平的指标，人均收入关系到消费水平，而消费增水平的改变导致购买力发生变化；

（4）人均住房面积、城市燃气普及率：

人均住房面积的波动关联到供暖面积，而供暖需要消耗能源，故影响到垃圾的产生量；城市燃气普及率的高低意味着煤炭等产生固体垃圾的能源消耗量的多少。

2.2.搜集城市生活垃圾产量数据

随机选定重庆市作为研究对象，在中国国家统计局网站上查到其2001~2006年人口数量、国内生产总值等各影响因素及垃圾产量的统计数据如表格2.21：

表格2.21

年份

生产总值/亿元

城市居民人均收入/元

人均消费水平/元

人口数量/万人

人均住房面积/平方米

城市燃气普及率/%

垃圾产量/万吨

2001

1765.68

6572.3

1078.06

3097.91

22.5

32.2

164.6

2002

1990.01

7238.07

1228.89

3113.83

23.9

46.6

211.7

2003

2272.82

8093.67

1415.31

3130.1

25.72

59.5

215.3

2004

2692.81

9220.96

3596

3144.23

28.25

63.52

237.2

2005

3070.49

10243.99

4782

3169.16

30.68

68.84

237.6

2006

3452.14

11569.74

5417

3198.87

31.36

75.84

243.9

2.3.整体分析得城市生活垃圾产量中短期预测模

首先画出各个影响因素与年份之间的散点图，确定为线性关系，建立一元线性回归模型，求解方程参数；再确定垃圾产量与各影响因素的近似拟合曲线，利用最小二乘原理得出精确的拟合方程；综上，将两次求得关系通过方程联立，即得最终预测模型。

2.4.构建垃圾收运路线的模型

对题目所给的数据进行预处理，画出车库、中转站、垃圾收集点分布图，以不超过每辆车的最大装载量为依据将垃圾收集点分组，接着，利用图论法画出每组的最小生成树，将树叶融入树中，形成最优路线。

根据最优路线算出从车库出发到收集完垃圾回车库所用的总时间，进一步确定合理安排车辆收运垃圾的方案。

3.模型的假设

1.假设影响垃圾产量的只有内在因素（如人口数量、居民生活水平等），不包括社会因素（如社会行为准则、社会道德规范）和个体因素（人类本身个体的行为习惯、受教育程度）[4]；

2.假设垃圾产量只包括被清运的垃圾，散落的垃圾不予统计。

查询垃圾产量统计数据时即默认为生活垃圾清运量；

4.符号说明

年份-2000

生产总值

城市居民人均收入

人均消费水平

人口数量

人均住房面积

城市燃气普及率

垃圾产量

5.模型的建立

5.1.模型一：

城市生活垃圾产量中短期预测模型

5.1.1.建立城市生活垃圾产量的各影响因素与年份的关系

（1）初步判断年份与城市生活垃圾产量的影响因素

先根据年份x与影响因素y的试验值画出散点图，根据散点图确定须配曲线的类型。

此处以生产总值为例，画出散点图，如图5.11：

图5.11

由图中散点知，数据点大致落在一条直线附近，故确定为一次函数曲线（其他几个因素均如此），建立一元线性回归模型：

其中固定的未知参数

、

称为回归系数，自变量x也称为回归变量。

将

，称为y对x的回归直线方程

（2）用n对试验值

对

、

和

作最小二乘估计

设

记

解得

或

其中

得回归方程

（3）检验回归方程

的显著性

归结为对假设

进行检验。

假设

被拒绝，则回归显著，认为y与x存在线性关系，所求的线性回归方程有意义；否则回归不显著，y与x的关系不能用一元线性回归模型来描述，所得的回归方程也无意义。

检验法：

当

成立时，

其中

（回归平方和），当

时，拒绝

，否则就接受

。

（4）判断回归系数

的置信区间

和

置信水平为1-α的置信区间分别为：

和

5.1.2.确定垃圾产量与其各影响因素的关系

（1）确定拟合函数

观测影响因素与垃圾产量的一系列数据集，考虑借助曲线拟合用一个相对简单的解析曲线去逼近所得到的数据集，但拟合的曲线往往不能完全符合给出的数据，因此需要对拟合的性能给出一个量度，这里使用最小二乘原理（极小化偏差的平方和）作为衡量曲线拟合优劣的准则，它不要求得到的曲线过所有的点（可消除误差的影响），只要求在给定点上的误差的平方和最小；并且能够尽可能表现数据的趋势，靠近原来的数据点。

这里选择形如

的式子作为拟合的曲线

式中

是线性无关的函数族

（2）根据最小二乘原理先求法方程系数

假定在

上给出一组数据

，

以及对应的一组权

，这里

为权系数，要求

使

最小，其中

实际上是关于

的多元函数，求

的最小值就是求多元函数

的极值，由极值必要条件，可得

根据内积定义引入相应带权内积记号

则可改写为

这是关于参数

的线性方程组，用矩阵表示为

上式称为法方程。

当

线性无关，且在点集

上至多只有

个不同零点，则称、

在

上满足Haar条件，此时法方程的解存在唯一，其解为

从而得到最小二乘拟合曲线

5.2.模型二：

垃圾收运路线规划模型

根据垃圾车最大装载量和数据分布将所给数据分成若干个区域，每个区域均有若干个点，通过图论生成最小生成树，找到连接所有垃圾收集点的最小路径。

通过Prim算法[5]求所分区域的最小生成树，方法如下：

（1）建立以两点间距离为权值的权邻接矩阵

，其中一两点间的Manhattan距离为权值，计算公式如下：

（2）建立初始候选边表

;

（3）从候选边中选出最短边

；

（4）调整候选边集

;

（5）重复

（2）、（3）直到

中含有

条边。

通过上述算法得到的线路即为连接各点的最小生成树，其特点是边权和最小。

但是，垃圾的收运不是简单的连接问题，而最小生成树形成的路线需要在树杈路线上回返，这样就会增加行程，故在最小生成树的基础上比较树杈与其相邻的权边，以边权增加最小为原则进行边替换，将树叶融入主干，这样便可以得到能将所分区域的所有垃圾收集点串联起来的一条“链”，且其权和最小，该“链”即为最佳的行车路线。

6.模型的求解

6.1.确定年份与城市生活垃圾产量各影响因素的关系式

6.1.1.求解影响因素与年份的一元回归模型

以生产总值因素为例，假设变量x与y满足如下关系：

，运用MATLAB中一元线性回归程序处理数据分析得：

，

，方程为

由拟合结果知：

a的置信区间为[297.7，393.4]，b的置信区间为[1144.9，1517.7]。

相关系数r2为0.9901，F的值为401.6431，与F对应的概率p为0.0000。

由于相关系数r2越接近1，说明回归方程越显著；F越大，说明回归方程越显著；与F对应的概率p小于显著性水平（缺省时为0.05）时，回归模型成立。

故回归方程及参数拟合准确。

根据上述步骤对其他指因素回归方程、相关系数r2、F的值以及F对应的概率p，如表格6.11：

表格6.11

指标

回归方程