垃圾分类模型

《垃圾分类模型》由会员分享，可在线阅读，更多相关《垃圾分类模型（24页珍藏版）》请在装配图网上搜索。

1、数学建模刘慧慧等人垃圾减量分类的量化分析摘要随着城镇化的加快和人们生活水平的提高，城市的生活垃圾对生态环境及人 类生存带来极大的威胁，这逐渐成为重要的社会问题。本文以深圳天景花园和阳 光花园为例，通过研究社会因素和个体因素对垃圾减量分析的影响，以及各类来 及之间的相关性，对垃圾减量工程进行了分析探讨，提出了垃圾减量量化模型的 和有效的预测算法。针对问题一：为了分析城市生活垃圾与社会因素以及个体因素的关系，我们 通过查阅深圳市20002010年的统计年鉴，统计了城市垃圾清运量与政府教育 支出、GDP、人均GDP、环保投入、总人口数、户籍/非户籍比例、家庭结构的数 据，运用的相关性分析法和多元线性

2、回归法，用EVIEWS软件进行分析，得到的 量化模型为：Y = 0.014799 + 0.185947X 10.47037X + 17.02541X + 0.620093X 11.52609X 2.573286X + 7.901079X1234567由于存在多重共线性，通过剔除变量，得到修正模型为：y= -0. 005174- 12.28401 +20. 1641 +0. 76265 - 13. 66949 -3.-126845 + 9. 210445 针对问题二：在相关性分析的基础上，可以直接用SPSS分析了四类垃圾之 间存在的相关性，得出了结论与原因。接着，利用第一问中构建的减量分类模型，

3、 我们将附录2, 3的各类垃圾组分的数量变化进行定性的相关性分析。最后总体 得出政府各项激励措施和减量分类效果之间存在正的相关性关系。针对问题三：基础数据项是否足够，主要是通过与一些发达地区及环保工作做得 比较好的地区之间比较以及地区之间的气候、地理等相似度,根据现有经验分析， 结果是在厨余垃圾和混合类垃圾及有害垃圾方面细分；对于颗粒度是否足够，主 要是根据附件8中的数据进行分析并参照生活垃圾采样和物理分析方法 (CJ/T313-2009)的要求判断，求解结果是足够的；通过垃圾分类的作图，在 对数据分析以及一些发达地区及地理位置相近的发达地区标准的基础上确定精 力和成本的投入方向；运用统计学上

4、的一些抽样思想来达到设置少量抽样数据来 检测一定区域内减量分类工作的效果的目的。针对问题四：通过前几问的分析我们可以看出，在自然因素，社会因素和内 在因素中，所占比重的是个体因素和社会因素，个体因素不稳定，所以政府可以 从社会因素出发。加大环保投入，注重教育来改善。我们的思路是，首先用灰色 预测法计算出政府教育支出（X ）、GDP（ X ）、人均GDP（ X ）、环保投入（X ）、1234总人口数（X ）户籍/非户籍比例（X ）家庭结构（X ）、在未来五年的预测567值，然后代入主成分分析法得到的主成分模型，就得到了城市生活垃圾清运量关键词：主成分分析灰色预测多元线形回归模型一问题的重述为了减

5、少该区垃圾占用大量土地、污染水体、污染大气、破坏植被，严重影 响地区的外容景观和居民的生活环境。该地区垃圾已成为各国政府急需妥善解决 的首要问题之一。根据该区的实际情况与环境影响则这一控制过程的研究需要进 行改良和量化分析，提出了以下问题。1、分析附件有关资料并结合你自己的经历和生活观察，考虑各项教育、督导、激励措施对居民家庭垃圾减量分类结果的影响，构建量化模型描述深圳天 景花园、阳光家园垃圾减量分类过程，模型应能以量化参数描述社会因素（如各项教育、督导、激励措施等）以及个体因素（如家庭收入水平、家 庭结构、户籍类型、生活习惯等），并在后续的进一步研究过程中通过调 整相关参数来修正模型。2、基

6、于你构建的减量分类模型，试分析试点小区四类垃圾组分本身的数量存在什么样的相关性？各项激励措施与减量分类效果存在什么相关性？原 因是什么？3、根据你构建减量分类模型的研究结果,你认为在深圳现有垃圾减量分类督导过程中，目前统计的基础数据分项及颗粒度是否足够？应该在哪些数据的获取中投放更多的成本和精力？在减量分类模式大面积推广时,如何设置少量抽样数据来检测一定区域内减量分类工作的效果？4、基于你构建的减量分类模型，指出深圳未来5年推进减量分类工作关键措施，并预测措施实施的最好与最坏结果二、问题分析2.1对于问题一分析：要分析社会因素和个体因素与城市生活垃圾的清运量之间的关系，我们就要 考虑与社会因素

7、相关的教育、激励措施，与个人因素相关的家庭收入水平、家庭 结构、户籍类型和城市生活垃圾的清运量之间的关系。我们通过查阅深圳市 20002010年的统计年鉴，统计了城市垃圾清运量与政府教育支出、GDP、人均 GDP、环保投入、总人口数、户籍/非户籍比例、家庭结构的数据，用Eviews软 件对这些数据进行相关性分析和多元线性回归分析，得出线性回归模型，也即是 这些变量对城市生活垃圾的清运量的量化关系，说明了各因素对垃圾清运量的具 体影响情况。2.2对于问题二分析：分析试点小区四类垃圾组分本身的数量存在什么样的相关性？根据附件二 由SPSS可以直接看他的相关性，由于数据的特殊性，删除了二月份的几个数

8、据, 得到跟好的结果，从结果中可以看出：总量垃圾由厨余垃圾，可回收垃圾构成， 与有害垃圾大致无关。可回收物与厨余垃圾呈弱正相关，与有害垃圾呈极弱相关, 与其他垃圾呈中等程度负相关；厨余垃圾与有害垃圾呈弱相关，与有害垃圾呈极 弱负相关；有害垃圾与其他垃圾呈极弱负相关。再基于你构建的减量分类模型可 以知道政府各项激励措施和减量分类效果之间存在正的相关性关系。2.3对于问题三分析：首先对本题中的几个名词加以分析，基础数据项：是对垃圾的具体分类。数 据颗粒度：是用于表示某数据集的组成的最小单元，根据附件8所给信息知道深 圳市对垃圾处理的颗粒度标准。基础数据分项是否足够主要通过与一些发达地区 及环保工作

9、做得比较好的地区之间比较以及根据现有经验分析；对于颗粒度是否 足够主要通过对附件8的数据分析并参照生活垃圾采样和物理分析方法(CJ/T313-2009)的要求判断；通过垃圾分类的作图，再对数据分析及根据一 些发达地区的标准的基础上确定成本和精力投入方向；运用统计学上的一些抽样 思想来达到设置少量抽样数据来检测一定区域内减量分类工作的效果的目的。2.4对于问题四分析：通过前几问的分析我们可以看出，在自然因素，社会因素和内在因素中，所占 比重的是个体因素和社会因素，个体因素不稳定，所以政府可以从社会因素出发。 加大环保投入，注重教育来改善。我们的思路是，首先用灰色预测法计算出政府 教育支出（X ）

10、、GDP（ X ）、人均GDP（ X ）、环保投入（X ）、总人口数（X ）、12345户籍/非户籍比例（X ）、家庭结构（X ）、在未来五年的预测值，然后代入主成 67分分析法得到的主成分模型，就得到了城市生活垃圾清运量。三、模型假设1. 假设题目所给数据及建模搜集数据真实可靠；2. 假设小区居民扔垃圾的量不受天气和气候的影响；3. 假设所有人都按要求对垃圾分类；4. 假设有新的变化因素引进时，或者原有的因素减少时，单一准则下, 原有因素的排序全职的比例不变；5. 各个指标因素互不相关，互不影响；6. 假设收垃圾的车不会因为故障而将收垃圾的时间延后；7. 假设建模中涉及主观分析的结论基本与事

11、实相符；8. 假设深圳队环保的投资全部用于对垃圾的减量分类过程。四、符号说明X -政府教育支出1X -GDP2X -人均 GDP3X -环保投入4X总人口数5X户籍/非户籍比例6X -家庭结构7Y-城市垃圾清运量Wi基础数据分项的抽样量权重N-总抽样量Pi各区人口P总人口五、模型建立与求解511模型的建立要分析社会因素和个体因素与城市生活垃圾的清运量之间的关系，我们就 要考虑与社会因素相关的教育、激励措施，与个人因素相关的家庭收入水平、家 庭结构、户籍类型和城市生活垃圾的清运量之间的关系。我们采用了城市垃圾清 运量与政府教育支出、GDP、人均GDP、环保投入、总人口数、户籍/非户籍比例、 家庭

12、结构的数据，对这些数据进行相关性分析和多元线性回归分析，得出最后的 量化模型。建立的线性回归模型为：Y=a +a X +a X + a X + a X + a X + a X + a X011223344556667模型的求解根据20002010年的深圳市的统计年鉴，可以得出如下的政府教育支出（X ）、GDP（ X ）、人均GDP（ X ）、环保投入（X ）、总人口数（X ）、户12345籍/非户籍比例（X ）家庭结构（X ）、城市垃圾清运量（Y）数据表如下：67年份政简教育支 岀（万元）慄圳经济发展（ciipi?!元）居民生话水 平.*灰 均EP 亿元）政简激励 探楼人万忒邕人数貝万A）总人

13、數户礬/非户礬比例家霆皱构 可平均每户 aB圾清运屋 逹万吨）20000. 02308840. 037333580. 048474460. 00040220. 077(4950. 07713490. 071949240. 089010990.053009120010. 02872190. 042369070. 051462790. 00082260. 0797012020. 07970120.073969570. 0340770. 057438720020. 04444820. 05068130. 059660530. 00168150. 0821266560. 08212670.076237

14、060. 091208790.05805010. 05529550. 061138080. 069503210. 00343470. 0856080910. 08560810.073860230. 03257140. 085197920040. 08556880. 0730840. 080154290. 00701850. 08808E3440. 08808630.086228770. 031758240. 091037420050. 06878510. 084438010. 089856570. 01434080. 0910507880. 09105080. 0508430. 0920329

15、70. 087403320060. 08556880. 0332210. 101147570. 02930370.093105550. 03310550. 100940110. 031483520. 094335120070. 10644750. 116083540.112722330. 05987520. 0347687180. 03476870. 100711650. 030M4E20. 106853120080. 1242070.133235360. 12330i：iA30. 1234330. 0364494870. 03644950. 10424708-0. 089285710. 11

16、5703T20090. 1647.3140. 139971770. 124353150. 24338450. 03SO-3345T0. 09S03350. 10635730. 089285710. 124963320100. 20492580. 16252430.139358160. 51079290. 11J3347580. 11333480. 105990540. 0S8461540. 1258277用EVIEWS对X、X、X、X、X、X、X、Y进行多元线性回归的1234567得到的结果如下表所示：Depe ndent Variable: YMethod: Least SquaresDat

17、e: 08/24/12 Time: 16:43Sample: 2000 2010In cluded observati on s: 11VariableCoefficientStd. Errort-StatisticProb.C-0.01479Q0.070929-0.2086410.8481X190.1859470.1382271.3452290.2712X2-10.470372.761677-3.7913070.0322X3/17.025414.3332873.9289830.0293X40.6200930.1837213.3751870.0432X5-11.5260Q3.008411-3.

18、8312880.0313X69-2.573280.779544-3.3010160.0457X767.9010791.8909254.1784200.0250R-squared0.995893Mean depe ndent0.09090var9Adjusted0.986310S.D. depe ndent0.02633R-squaredvar9S.E. of regressi on0.003082Akaike info-8.5713criteri on37Sum squared2.85E-05Schwarz criteri on-8.2819residLog likelihood55.1423

19、659 F-statistic103.9216 Durbi n-Watso n2.602710 Prob(F-statistic) 0.00142stat1则线性回归模型为：+ 7.901079X7Y = 0.014799 + 0.185947X -10.47037X + 17.02541X + 0.620093X - 11.52609X - 2.573286X123456513模型的结果和分析由上表可以得出：模型的可决系数为R 2=0.995893，说明模型建立的可行性较强，但变量X的P值大于0.05，不通过检验，说明X与其他各变量之间存在1 1较强的多重共线性，通过分析，将变量X剔除，进行

20、模型的修正，可以得到的 结果如下表所示：Depe ndent Variable: YMethod: Least SquaresDate: 08/24/12 Time: 16:53Sample: 2000 2010In eluded observati on s: 11VariableCoefficieStd. Error t-StatisticProb.ntC-0.00517A0.077380 -0.0668660.9499X24-12.28402.642926 -4.6478820.0097X3I20.164164.004015 5.0359870.0073X40.7626500.16456

21、7 4.6342940.0098X5-13.6694 Q2.798193 -4.8851150.0081X69-3.126840.726007 -4.3069100.0126X759.2104451.777690 5.1811320.0066R-squared0.993416Mean depe ndent0.09090var9Adjusted0.983539S.D. depe ndent0.02633R-squaredvar9S.E. of regressi on0.003379Akaike info-8.2811criteri on45Sum squared4.57E-05Schwarz c

22、riteri on-8.0279residLog likelihood52.5463039 F-statistic100.5815 Durbi n-Wats on 2.361642 Prob(F-statistic) 0.00025 stat8由上表可知：可决系数R 2=0.993416较大，模型的可行性较高，而且各变 量都通过检验，所以最终的修正模型为：y= -0. 005174- 1228401 +20. 1641 +0. 76265 - 13. 66949-3. 126845 + 9. 210445 模型的建立因素间关联性如何、关联程度如何量化等问题是系统分析的关键和起点。而 关联度和相

23、关性是把各个时刻的关联系数集中为一个平均值，利用这个我们可以 对各种问题进行因素分析。基于主成分综合模型的相关性分析，将垃圾分为四类，并将组分本身的相关 性分析相关系数与相关性关系相关性r的绝对值强相关中等程度相关弱相关极弱相关或无相关522模型的求解根据附录二给出的天景花园三个月份的垃圾回收数据，分析前三个月的垃圾收集 统计表， 绘制垃圾回收数量随时间变化如下图由上图可知：三个月中绝大部分时间各类垃圾组分回收数量在一个区间内的变化 是相对稳定，但二月初各类垃圾组分数量都有很大的波动，对数据的影响很大， 故将这几个特殊点剔除掉之后，再次利用excel画出的统计表如下：由上图可知：将几个特殊点剔

24、除掉之后，三个月中绝大部分时间各类垃圾组分回 收数量在一个相对稳定区，进而将天景小区四类垃圾组分本身的数量分别用 SPSS做相关系数计算后，所得的天景花园各类指标相关系数如下图：Correlations可回收物厨余垃圾有害垃圾苴他垃圾总星可回收物Pearson Correlation1.305.195-.524.243Sig. (2-tailed).005.076.000.026N8484040484厨余垃圾Pearson Correlation.3051.128-.660.456Sig. (2-tailed).005.245.000.000N8484343434有害垃圾Pearson Cor

25、relation.195.1281-.100.164Sig. (2-tailed).076.245.365.135N8484040484苴他垃圾Pearson Correlation-.524-.660-.1001.254Sig. (2-tailed).000.000.365.020N8484848484总量Pearson Correlation.243.456.164.2541Sig. (2-tailed).026.000.135.020N8484040484*. Correlation is 引gnificant at the 0.01 level (2-tailed).*. Correl

26、ation is significant at the 0.05 level (2-tailed).问题的结果与分析对于天景花园：1.可回收物与厨余垃圾呈弱正相关，与有害垃圾呈极弱相关，与其他垃圾呈 中等程度负相关；厨余垃圾与有害垃圾呈弱相关，与有害垃圾呈极弱负相关；有 害垃圾与其他垃圾呈极弱负相关。2从图中还可以看出，总量垃圾与厨余垃圾，可回收垃圾，其他垃圾的相关 性很大，而与有害垃圾的相关性很小，故厨余垃圾，可回收垃圾，其他垃圾才是 总量垃圾的主要成分。3其他垃圾与可回收垃圾，厨余垃圾都呈负相关性，由此可以看出：可回收 物，厨余垃圾多，则其他垃圾少，可回收物，厨余垃圾很少时，其他垃圾很多，

27、 由此可知其他垃圾中包含很多可回收垃圾和厨余垃圾，故应增强居民的垃圾分 类，明确垃圾分类的种类，使居民的垃圾分类更加彻底，合理，从源头上减少污 染，增加回收利用的效率。用同样的方法对于阳光花园进行分析可以得出以下结论：总量垃圾很大部分都由 其他垃圾构成，厨余垃圾与可回收垃圾和其他垃圾呈弱负相关，其他垃圾与可回 收垃圾呈强负相关。524激励措施的相关性分析由第一问给出的减量分类模型y= -0. 005174- 12.28401 +20. 1641 +0. 76265 - 13. 66949-3. 126845 + 9. 210445 可知，其中x4是政府鼓励措施，由于x4与y是正相关，即政府的鼓

28、励措施越多, 减量增量的效果越好，故应增加奖金，鼓励方面的措施。基础数据项分析各组分比例（）针对基础数据项，参照生活垃圾采样和物理分析方法（CJ/T313-2009），对 深圳市的垃圾组进行分析如下图：厨余类纸类橡塑类纺织类木竹类灰土类砖瓦陶瓷玻璃类金属类其他混合类图1深圳市垃圾物理成分构成成分从图中明显可看出，深圳市城市生活垃圾中厨余类有机物含量很高，平均百 分含量为35.07%,其次是塑料、混合类、纸类，含量最低的为其他类（危险固 废）平均百分含量为0.01%，其中厨余、塑料、纸类和混合所占比重较大，四类 共占90.02%，其他比重相对较小。样品中可回收物橡塑类、纸类、纺织类、木 竹类、金

29、属类、玻璃类总含量较高，为45.55%， 超过厨余的含量35.07%。样品中橡塑类多为塑料袋和一次性餐盒，纸类多为浸 湿的尿不湿、卫生巾、卫生纸、利乐包等物质，回收价值不高。可回收物中回收 价值较高的物质如塑料瓶、报纸、杂志、纸箱、金属等，在收集过程中已被拾荒 者、保洁员等拣走，含量很少。根据附件8将深圳市的垃圾物理构成与发达国家（如：美国、新加坡）以及 香港作比较如下：美国2009年生活垃圾组分百分比废纸 玻璃 金属 塑料 橡胶 纺织物 木竹 其他 厨余 庭院垃圾 渣土图2美国垃圾物理构成成分图图3新家坡垃圾物理构成成分图深圳与香港10年垃圾成分对比深圳关出图3-4深圳与香港10年垃圾成分比

30、较通过SPSS软件得到了深圳与香港的垃圾成分的误差上下限，从单纯的数据 我们并不认为深圳与香港垃圾成分有显著性差异，但是从显示的误差线可以看 到，除了纸类的误差上下限之间几乎没相差、厨余、金属与其他类的误差上下限 相差不是很大外，橡塑、纺织、木竹、玻璃等这些种类的垃圾标准偏差比较大， 即两者之间有显著性的差异。因此可以得出结论：深圳的生活垃圾组分与 香港 生活垃圾组分有较明显的差别。由对比可以得到：1、深圳市的垃圾分类中厨余类占的含量高达35%，明显高于发达，且混合 垃圾类占18%，而可回收物（纸类、橡塑、玻璃等）含量低于发达国家, 就说明深圳的基础数据分类尚不足够细，许多厨余垃圾和混合垃圾没

31、有 得到很好的分类，故仍需将这两类进行细致划分以德奥更多的分项。2、美国将灰土类细分为庭院垃圾和渣土，深圳市可以根据自身情况来细分 灰土类垃圾。3、以上2个国家、以及深圳市都没有包括有害垃圾这一项，如：废旧电池、 医药类的卫生废品、水银等，这些垃圾需要经过特殊处理，故而可以考 虑增加这一数据分项。4、深圳市与香港的垃圾组成成分由些项目具有相似之处，在这些方面可以 借鉴香港的垃圾分类及处理方法。同时也要注意到，香港与深圳的地理位 置比较接近，故某些垃圾表现出相似的性质，相似地，美国、新加坡等由 于各自特定的地理位置，气候等因素的影响，垃圾组分也会有差别，不要 盲目借鉴其他地区的经验。对数据颗粒度

32、分析针对数据颗粒度分析，从附件8可知，其统计了各垃圾成分的含水率、 灰分、可燃分以及干基热值等。在一次样品制备、二次样品制备以及混合 样品制备的过程中，将烘干后的生活垃圾样品中各种成分的粒径分级破碎至5mm以下，严格按 照生活垃圾样品物理组成的基干比例，将粒径为5mm以下的各种成分混合 均匀，缩分至500g，用研磨仪将其粒径研磨至0.5mm以下，这中程度的研 磨对计算垃圾含水率和物理组成干基比例已经足够清晰明了，又因为垃圾 含水率、物理干基和热值对分来研究占有很大的影响度，所以根据此考虑 方式颗粒度足够。533对投放成本与精力分析根据附件8中的信息及垃圾组分分类作图归纳如下：1. 深圳市的垃圾

33、分类中厨余类占得总量高达35%，混合类垃圾占18%，说明许 多厨余垃圾和混合类垃圾没有得到很好的分类，应投入精力进行有地分类。2其他的垃圾的概念定义也存在界限不明，标准模糊的问题，因此也应该在该 项数据的获取中投入更多成本和精力，可以建立一套科学规范的数据收集体系对 其进行数据统计。3.从南山垃圾焚烧厂，盐田垃圾焚烧厂，下坪垃圾填埋场，老虎坑环境园，平 湖垃圾焚烧厂，鸭湖垃圾填埋场这6个城市生活垃圾末端处理设施所统计的垃圾 含水率来看（以下坪垃圾填埋场例），各种垃圾的含水率普遍较高，部分原因是 抽样前天深圳市大部分地区下雨，而下坪垃圾填埋场未下雨，分析下坪填埋场的 垃圾含水量可得出，各项垃圾含

34、水量依然很高。故对这一部分水资源的利用是一 个值得研究的问题，因为一方面，从能量角度一般垃圾处理方式都是燃烧，在燃 烧过程中有大量水分会消耗更多的热能；另一方面水资源紧缺，若能提取出里面 大量的水资源，是对水资源的一种节约。所以从这一角度看，应投入资金和精力 探究水资源的重新利用问题。61.5557.8955.4878.61% % % o o O6 4 2 率水含分成各圾垃60.449395085559.6348.2637.947.3262.7164.58L60.6957.9551.2158.567.37混合类木竹类纺织类橡塑类纸类厨余类23样品编号平均值图3-4填埋场垃圾含水率5.3.4 对

35、设置少量抽样数据来检测垃圾分类效果分析在减量分类模式大面积推广时，为了准确的检测一定区域内减量分类工作的 效果的同时又尽量减少工作量，设置一套科学合理的抽样方法显得尤为重要。综 合分析了目前主流的几种抽样检测方法后，结合垃圾分类的工作特点，我们采用分层抽样的方式来检测一定区域内减量分类工作的效果。 1减少抽样点个数(分层抽样思想)参照生活垃圾采样和物理分析方法(CJ/T313-2009)表3-1采样点数要求人口数量/万人=200最少采样点数8162030可将深圳市按区域大t致划分为九个区，然后由各区人口 Pi占总丿1 口 P的不同各自分配出不同的抽样比例，计算出各个区的人口比重：Wi=Pi/P

36、,以此来指导采样 点的设置。2. 减少抽样的垃圾基础数据分项的量，直接以垃圾各组分的比例作为各基础数据 分项的抽样量权重Wi，假设令总抽样量为N,各基础数据分项的抽样量模型： ni=Nwi。3. 除此之外，根据历史数据统计结果，对一直比较好的地区可以减少抽样频率或 减少抽样次数。如特区外较特区内垃圾总量较少，可以相对减少抽查频率。5.4问题四工作关键措施通过前几问的分析我们可以看出，在自然因素，社会因素和内在因素中，所占比 重的是个体因素和社会因素，个体因素不稳定，所以政府可以从社会因素出发。 加大环保投入，注重教育来改善。灰色预测法灰色预测法是一种对含有不确定因素的系统进行预测的方法。灰色系

37、统是介于白色系统和 黑色系统之间的一种系统。gm(i, 1)模型设有原始数据列x(o)=(x(o)(l),x(o)(2), ，x(o)(n),n为数据个数。如果根据x(o)数据列建立GM(1,1 )来实现预测功能，则基本步骤如下：(1)原始数据累加以便弱化随机序列的波动性和随机性，得到新数据序列：x 二x(d(1),x(2),，x(n)其中,X(1) (t)中各数据表示对应前几项数据的累加。x(i) (t)二 x(o)(k)， t=1,2,，n对X(1)(t)建立下述一阶线性微分方程：即GM(1,1)模型。dx(i)+ ax(i)= u dt其中，a,u为待定系数，分别称为发展系数和灰色作用量

38、,a的有效区间是(-2,2),并记a,u构成的矩阵为灰参数a = fa。只要求出参数Iu丿a,u,就能求出x(t),进而求出x(o)的未来预测值。(3) 对累加生成数据做均值生成B与常数项向量Yn,即-丄(X(1)(1) + X(1)(2)1B= - 2(X(1)(2) + X(1)(3)1-丄(X(1)(n -1) + X (1) (n)12f x(o)(2)x(0) (3)Yn=I x(o)(n)丿(4) 用最小二乘法求解灰参数a,则Aa = (BTB)-1BTYn(5) 将灰参数a代入dx + ax(D = u，并对dx+ ax(D = u进行求解，得dtdt人uu殳(1)(t+1) =

39、 (X(0)(1) )e - at + -aa由于a是通过最小二乘法求出的近似值，所 以x (1)(t+i)是一个 近似表达式,为了与原序列X(i) (t + 1)区分开来,故记为x(t + 1)。式 中t为时间序列。对函数表达式x(t+1)及x(t)进行离散,并将二者做差以便还 原X(o)原序列，得到近似数据序列X () (t + 1)如下：X () (t + 1) = X (t + 1)- X (t)(7) 对建立的灰色模型进行检验，步骤如下：计算x() (t)与x () (t)之间的残差e()(t)和相对误差q()(t):e()(t) =x()(t)-x () (t)q()(t)二e()

40、(t)/x()(t)等等。(8) 利用模型进行预测：x () = x (0)(1), x()(2),，x()(n), x()(n+1),，x () (n+m)原数列的模拟未来数列的预测例如对政府教育支出用matlab编程我们得到，参数，同理可以得到GDP,人均GDP、环保投入、总人口数、 户籍/非户籍比例、家庭结构，的参数值，及预测值。年份政府教 育支出（万元）深圳经 济发展（GDP亿元）居民生活 水平（人均政府激励 （环保投入 万元）人 万 总数人户籍/非 户籍比例家庭结构 （平均每户 人数）GDP元）万201118549421127110.81772286100110.36133.2120

41、1222459921294611.951533016103470.36673.19201327194821487013.223043093106940.38273.17201432927901707914.636040653110530.39953.15201539869601961716.1811990920114240.41703.14我们的思路是，首先用灰色预测法计算出政府教育支出（X ）、GDP （ X ）、人1 2均GDP （ X ）、环保投入（X、总人口数（X ）、户籍/非户籍比例（X ）、3456家庭结构（X 7）、在未来五年的预测值，然后代入主成分分析法得到的主成分模 型，就得

42、到了城市生活垃圾清运结合第一问得到的减量分类模型= -0.005174- 12.28401+ 20. 1641+ 0.76265x4 - 13. 66949-3, 126845x + 9. 210445b将以上数据代入，得到城市生活垃圾清运量的预测值城市生活垃年份圾清运量（万吨）2111548.182012 594.172013 644.032014 698.072015 756.63另外我们用matlab画出了拟合值和实际值的图像量运清圾垃活生市城直接用灰色预测模型对城市生活垃圾清运量得到的预测图像年份城市生活垃圾清运量做好最坏2111550.09546.272012596.58591.76

43、2013646.35641.712014700.03696.102015762.16751.10六模型的评价与改进第一问数据都模拟的很好，可以很好的预测，但对所有的变量在取量纲的时 候不是特别合理，会影响最终的结果。第二问用SPSS分析的相关性分析的很好，很好的突出四种垃圾的关系，从 第一问也可看出政府的鼓励措施与垃圾的减量分析呈正相关。符合实际，由此可 以加强居民的垃圾分类意识，政府的鼓励措施也要加强。第四问灰色微分拟合法建立的离散拟合方程是一个近似差分方程，因而很难保证拟 合方程与待拟合系统的微分方程严格近似，也就无法保证所建立的灰色模型的固有误差为无 穷小量。七参考文献1. 谢季坚，刘承

44、平.模糊数学方法及其应用M.武汉：华中理工大学出版社,2000： 205-2112. 葛军，葛伦应.层次分析法确定水质指标权重J.当代建筑，2003, 3(1) : 22-23 3杨纶标.模糊数学原理及应用M.广州：华南理工大学出版社，2000： 67-804. 李士勇.工程模糊数学及应用M.哈尔滨：哈尔滨工业大学出版社，2004:101-1085. 宁11佟春生.系统工程的理论与方法概论M.北京：国防工业出版社，2000：185-1866姜启源.数学模型M.北京：高等教育出版社，2005： 225-2317. 晓秋.模糊数学原理与方法M.徐州：中国矿业大学出版社，2004： 203-2088

45、. 刘元高，刘耀儒.Mathematica4.0实用教程M.北京：国防工业出版社，2000：65-689李安贵,张志宏，孟艳等.模糊数学及其应用M.武汉:冶金工业出版社,2003：144-14610. 李志西，杜双奎等.实验优化设计与统计分析M.北京：科学出版社，2010:23-23911. 杨玉楠，熊运实，杨军等.固体废物的处理处置工程与管理M.北京：科学出版社，2010:45-632附录clearx0=179744223601346031430462666156535495666156828698103090012824401595355；pre_num=5;n=length(x0);di

46、sp(年份)lambda=x0(1:end-1)./x0(2:end);range=minmax(lambda)x1=cumsum(x0);z=0.5*(x1(2:end)+x1(1:end-1);Y=x0(2:end);B=-z(1:end) ones(n-1,1);u=BY; %u=inv(B*B)*B*Ya=u(1)b=u(2)x0_pre=x0(1) ones(1,n+pre_num-1);for k=1:n-1+pre_numx0_pre(k+1)=(x0(1)-b/a)*(exp(-a*k)-exp(-a*(k-1);enderr=x0 - x0_pre(1:n);epsilon=abs(err)./x0(1:n).*100disp(时间)disp(x0_pre) 入入disp(IaOIo2i)disp(epsilon)t1=2000:2010;t2=2000:2015;plot(t1,x0,d,t2,x0_pre,LineWidth,2) %O-宜箔YOeOo2船y箔Yp土 宜Ixlabel(e -Y)ylabel(3EDEuiA-.0aOEA)