數(shù)學(xué)建模垃圾分類處理

1、數(shù)學(xué)建模垃圾分類處理 陳云中1 問題的重述在垃圾分類收集與處理中，不同類的垃圾有不同的處理方式，簡述如下：1）櫥余垃圾可以使用脫水干燥處理裝置，處理后的干物質(zhì)運送飼料加工廠做原料。不同處理規(guī)模的設(shè)備成本和運行成本（分大型和小型）見附錄1說明。2） 可回收垃圾將收集后分類再利用。3） 有害垃圾，運送到固廢處理中心集中處理。 4）其他不可回收垃圾將運送到填埋場或焚燒場處理。所有垃圾將從小區(qū)運送到附近的轉(zhuǎn)運站，再運送到少數(shù)幾個垃圾處理中心。顯然，1）和2）兩項中，經(jīng)過處理，回收和利用，產(chǎn)生經(jīng)濟效益，而3）和4）只有消耗處理費用，不產(chǎn)生經(jīng)濟效益。 1） 假定現(xiàn)有垃圾轉(zhuǎn)運站規(guī)模與位置不變條件下，給出大、

2、小型設(shè)備（櫥余垃圾）的分布設(shè)計，同時在目前的運輸裝備條件下給出清運路線的具體方案。以期達到最佳經(jīng)濟效益和環(huán)保效果。2） 假設(shè)轉(zhuǎn)運站允許重新設(shè)計，請為問題1）的目標(biāo)重新設(shè)計。2 基本假設(shè)（1）假設(shè)各小區(qū)清運站每天的垃圾量是不變的；（2）假設(shè)各小區(qū)清運站的垃圾都必須在當(dāng)天清理完畢；（3）不考慮運輸車在行駛過程中出現(xiàn)的塞車、拋錨等耽誤時間的情況；（4）不允許運輸車有超載現(xiàn)象；（5）每個小區(qū)清運站均位于街道旁，保證運輸車行駛順暢；（6）城區(qū)人口分為不同部分，每部分人口固定，每天產(chǎn)生垃圾量固定；（7）一天只從小區(qū)清運站收一次垃圾（晚上或下午）；（8）所有運輸車均從垃圾轉(zhuǎn)運站發(fā)車最后回到垃圾轉(zhuǎn)運站；（9）

3、運輸車將垃圾一起送往大型設(shè)備處和小型設(shè)備處再前往墳埋場和焚燒場；（10）大型垃圾處理廠的壽命是30年。小型垃圾處理機的壽命是10年； (11) 建設(shè)在運輸垃圾過程中沒有新垃圾入站。3 符號（參數(shù)）說明（1） (=1，2，,)為第j個解釋變量；（2） (=1，2，,) 為第j個未知參數(shù)；（3） 為隨機誤差項；（4） S為多元線性回歸模型的精度；（5） Pi(xi,yi)為第i個轉(zhuǎn)運站的坐標(biāo)；（6） Pj(Xj,Yj)為大型廚余垃圾處理設(shè)備建在地圖上的坐標(biāo)；（7） cost1為大型垃圾處理設(shè)備每日垃圾處理費用；（8） Cost2為小型垃圾處理設(shè)備每日垃圾處理費用；（9）|A| 表示A點到原點的距離

4、，恒正（10）|B| 表示B點到原點的距離，恒正（11）|A-B| 表示A,B兩點之間的距離,恒正（12）Ta 表示A點所在地的垃圾量（13）Tb 表示A點所在地的垃圾量（14）cost：耗油量; (15) T為規(guī)劃使用年限; (16) Cik為第i座收集站運往第k座中轉(zhuǎn)站單位運輸量單位距離的費用(元t- 1km- 1 ) ; (17) Xik為第i座收集站運往第k 座中轉(zhuǎn)站的日運輸垃圾量( td- 1 ) ;（18） Lik為第i座收集站運往第k座中轉(zhuǎn)站運輸距離(km) ;（19） Dk j為第k座中站運往第j座處理場單位運輸量單位距離的費用(元t- 1 km- 1 ) ; （20） Yk

5、j為第k座中轉(zhuǎn)站運往第j座處理場日運輸垃圾量( td- 1 ) ; （21） Sk j為第k座中轉(zhuǎn)站運往第j座處理場運輸距離(km)；（22） Fk 為規(guī)劃期內(nèi)待建中轉(zhuǎn)站的固定投資(元) ; （23） E為中轉(zhuǎn)站的運行成本(元t- 1 ) ; （24） Qmin為中轉(zhuǎn)站建設(shè)的最小控制規(guī)模( td- 1 ) ; （25） Qmax為中轉(zhuǎn)站建設(shè)的最大控制規(guī)模( td- 1）；. 5 模型的構(gòu)建與求解5.1問題一的建模與求解5.1.1城市生活垃圾產(chǎn)生量的預(yù)測表一 城鎮(zhèn)垃圾產(chǎn)生量歷年統(tǒng)計表（萬噸）年份20012002200320042005垃圾量281.8284.7290.4296302年份20062

6、007200820092010垃圾量321361.4357383.29413 假定被解釋變量，與多個解釋變量，,。之間具有線性關(guān)系,即 (8)其中 (=1，2，,)為k個解釋變量，(=1，2，,) 為+1個未知參數(shù)，為隨機誤差項。被解釋變量Y的期望值與解釋變量，,的線性方程為： (9)對于n組觀測值, (=1,2, ，n)，其方程組形式為： (10)即 其矩陣形式為 即 Y=X+ (11) 其中為被解釋變量的觀測值向量；=為被解釋變量的觀測值矩陣；為總體回歸參數(shù)向量；為隨機誤差向量。總體回歸方程為： E（Y）=X (12) 可采用最小二乘法對上式中的待估回歸系數(shù)進行估計，求得值后，即可利用多元

7、線性回歸模型進行預(yù)測了。我們對多元線性回歸分析進行數(shù)學(xué)檢驗，包括回歸方程和回歸系數(shù)的顯著性檢驗。a. 回歸方程的顯著性檢驗，采用統(tǒng)計量： (13) 式中；為回歸平方和，其自由度為m;為剩余平方和，其自由度為（n-m-1）。利用上式計算出F值后，再利用F分布表進行檢驗。給定顯著性水平，在F分布表中查出自由度為m和(n一m一1)的值，如果，則說明與的線性相關(guān)密切；反之，則說明兩者線性關(guān)系不密切。b回歸系數(shù)的顯著性檢驗，采用統(tǒng)計量： （14）式中，為相關(guān)矩陣的對角線上的元素。對于給定的置信水平，查分布表得，若計算值，則拒絕原假設(shè)，即認為是重要變量，反之，則認為，變量可以剔除。多元線性回歸模型的精度，

8、可以利用剩余標(biāo)準差 (15)來衡量。越小，則用回歸方程預(yù)測越精確；反之亦然。采用matlab軟件編程進行城市生活垃圾量多元線性回歸模型預(yù)測(預(yù)測代碼見附錄1)。表二為訓(xùn)練結(jié)束后預(yù)測值與統(tǒng)計值的對比表，精度達到要求后用訓(xùn)練好的模型來預(yù)測深圳市2011-2015年城市生活垃圾產(chǎn)生量，預(yù)測結(jié)果見表311。在matlab軟件中運行代碼后得到生活垃圾產(chǎn)生量的回歸方程為：Y=387965+025178xXl+010508xx200574xx3+O1292xx4-00138xx5+208016xx6-00095xx7+00066xxs一31460xx9方差估計：S=257642回歸方程的顯著性檢驗F統(tǒng)計量，

9、F=723187，所以拒絕假設(shè)，即回歸模型成立。 表二線性回歸模型預(yù)測值與統(tǒng)計值對比表年份20012002200320042005預(yù)測值280.12288.24291.97300.89308.02統(tǒng)計值281.80284.70290.40296.00302.00年份20062007200820092010預(yù)測值314.79359.58358.08390.34411.29統(tǒng)計值321.00361.40357.00383.29413 圖一 線性回歸模型預(yù)測值與統(tǒng)計值對比分析圖 從表二及圖一可以看出，多元線性回歸模型對歷史值的擬合程度較高，預(yù)測精度是可以接受的，多元線性回歸模型預(yù)測值比較接近深圳市城

10、市生活垃圾實際產(chǎn)生量，稍微偏高。 表三 2011-2015年深圳市城市生活垃圾產(chǎn)生量多元線性回歸模型預(yù)測值年份20112012201320142015預(yù)測值（萬噸） 397.3 413.0 429.0 445.3 461.8 大小型廚余垃圾設(shè)備規(guī)劃5.1.2.1模型的建立題目要求給出大、小型設(shè)備（櫥余垃圾）的分布設(shè)計。由于大型廚余垃圾處理設(shè)備處理能力為200噸/日，投資額約為4500萬元，運行成本為150元/噸。而每個轉(zhuǎn)運站的垃圾數(shù)量有限，所以大型廚余垃圾處理設(shè)備必須在圖上重新選址建設(shè)。小型餐廚垃圾處理機，處理能力為200-300公斤/日，投資額約為28萬元，運行成本為200元/噸。所以小型垃

11、圾處理機可以設(shè)置在垃圾中轉(zhuǎn)站內(nèi)。根據(jù)表四用matlab6.5編程作圖二（程序見附錄三）表四中轉(zhuǎn)站坐標(biāo)名稱中轉(zhuǎn)站廚余垃圾量Xy名稱中轉(zhuǎn)站廚余垃圾量xyP1站1109.8622.18P19站19164.8912.06P2站2108.0421.69P20站20257.5511.96P3站388.3420.92P21站21205.8611.43P4站43012.4420.39P22站22359.8811.38P5站556.9719.8P23站233011.7711.67P6站6516.0817.77P24站24305.4310.85P7站71014.3317.48P25站25205.7610.17P8

12、站81014.3417.24P26站26158.1810.65P9站92010.917.43P27站273011.0810.22P10站102510.6416.51P28站28155.529.39P11站11209.8816.37P29站29208.769.15P12站12404.8912.08P30站30155.768.96P13站13159.215.06P31站31254.58.23P14站142010.8513.9P32站32307.365.28P15站1556.1514.43P33站33159.116.85P16站161512.3513.66P34站34307.074.6P17站172

13、57.9213.12P35站357014.5311.09P18站18107.8913.12P36站36156.638.72P37站37256.259.73P38站38403.394.94圖二中轉(zhuǎn)站坐標(biāo)圖從圖表可知每個垃圾轉(zhuǎn)運站的坐標(biāo)Pi(xi,yi)，假設(shè)大型廚余垃圾處理設(shè)備建在地圖上的Pj(Xj,Yj)。所以對于每個垃圾中轉(zhuǎn)站來說有兩種情況：（1）在站內(nèi)設(shè)置垃圾處理機。(2）把垃圾運往大型廚余垃圾處理廠進行處理。從中選擇最優(yōu)方案，從而確定垃圾大型垃圾處理站的位置。假設(shè)大型垃圾處理廠的壽命是30年。小型垃圾處理機的壽命是10年。大型垃圾處理設(shè)備的平均每噸耗損成本=45000000/（30*36

14、5*200）=20元/噸小型垃圾處理設(shè)備的平均每噸耗損成本=280000/（10*365*0.3）=256元/噸2.5噸汽車，每車耗油20L35L 70#汽油/百公里。每升70#汽油價格為7.2元司機月薪平均3500元。如果運往大型垃圾處理設(shè)備廠，則每日垃圾處理費用（cost1）=平均每日設(shè)備耗損成本+運輸費用+司機工資+垃圾處理費用。如果在垃圾轉(zhuǎn)運站設(shè)置小型垃圾處理機，則每日垃圾處理費用(cost2)=平均每日設(shè)備耗損成本+垃圾處理費用。要確定大型垃圾處理廠的位置，需要計算出選擇第1種方案的點。根據(jù)以上條件建立模型：Cost1=Cost2=Cost1x2, y1y2,不難看出A在B的后方,即

15、A比B遠。對于前方參考點O,要將A,B對應(yīng)垃圾點的垃圾全部取回再返回O,一共有三種方式：1O-A-O, O-B-O單獨運輸。這種情況下，總的路程消費等于空載運行費用(20L/百公里)與裝載時運行費用（20+6*Ta）L/百公里）的總和。于是有：Cost = 20*|A| + （20+6*Ta）*|A| + 20*|B| + （20+6*Ta）*|B|2. O-A-B-O先遠點再近點，即先空載至最遠處，裝完A點垃圾后再返回至B,再回O點，有： Cost = 20*|A| + （20+6*Ta）*|A-B| +（20+6*(Ta+Tb)）*|B| = 20*|A| + 1.8*|A|*Ta + 1

16、.8*|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 比較以上三種情況，遠近點的遍歷順序，可以看出，“先遠后近”絕對比“先近后遠在花費錢的數(shù)量上要少的多，省出（20+6*Tb）*|A-B|*2 這部分的錢主要是車載著B點的垃圾奔到A點再返回B點。而又注意到兩者的時間花費是相等的。所以在其余同等的情況下選擇“先遠后近”。

17、考慮單獨運輸比其余的兩種運輸花費的錢仍不比“先遠后近”省，還多了20*|B|，所以一般情況下，不采用單獨運輸。2 A，B兩點沒有明顯先后順序。 -并鄰狀態(tài)（如圖三） 圖三還是一共有三種情況： 1O-A-O, O-B-O單獨運輸。這種情況下，跟A,B兩點有先后順序中的情況完全相同，即有：Cost = 20*|A| + （20+6*Ta）*|A| + 20*|B| + （20+6*Ta）*|B|2O-A-B-OCost = 20*|A| + (20+6*Ta)*|A-B| + (20+6*(Ta+Tb)*|B| -13.O-B-A-OCost = 20*|B| +(20+6*Tb)*|A-B|+(

18、20+6*(Ta+Tb)*|A| -2相比之下，清晰可見并鄰狀態(tài)下的單獨運輸所花的費用最少，所以在不要求時間的情況下對于并鄰兩點，采用單獨運輸?shù)姆绞阶罟?jié)約錢。用式與式相減， 得到如下判斷式：6*(Ta-Tb) *|A-B| +6*(Ta+Tb)*（|B|-|A|）-上式 A-B-O;上式 0時， 選 O-B-A-O;上式 = 0時， 任意選上述兩路線。三兩點選擇趨勢的討論。 （如圖四） 圖四由圖中看到B,C兩點沒有明顯的先后順序，屬于并鄰點。因為當(dāng)運輸車載重行駛時費用會成倍的增長，比其空載時所花費用要大的多，所以排除A-B-C或A-C-B這樣的一次經(jīng)過3點的往返路線，僅選擇B,C中的某一點與A

19、完成此次運輸，將另一點留到下次。那么A點選擇B還是C呢？不妨假設(shè)|B|C|,即B點離原點的距離比C點的更遠，因為A在B,C之后，所以也就是B點離A點更近。這樣，此次的運輸我們更趨向于選擇A-B,因為就這三點而論，A無論是選B還是C，三點的垃圾總要運完，所以花費的錢是一樣的。但選擇A-B后，下次運輸車運C點垃圾時就無需跑的更遠。綜上所述，得出搜索的基本原則：1在兩點遞減的情況下，不采用單獨運輸；2在其余同等的情況下選擇“先遠后近”；3不要求時間的情況下對于并鄰兩點，采用單獨運輸?shù)姆绞阶罟?jié)約 錢；一般情況下用式3作判斷；4車在裝的足夠多的情況下應(yīng)該直接返回中轉(zhuǎn)站；5每一次布局和每條線路的搜索不妨由

20、剩下未搜點中的最大值開始四 關(guān)于垃圾點的垃圾是否一次清除的討論這里說的一次清除問題不是指一天，而是指當(dāng)一輛運輸車已經(jīng)裝載了足夠多的垃圾，不能完全清理下一個垃圾點的時候，車在下一個站點“停還是不?！钡膯栴}。例如，一輛運輸車選擇了某段路線后，當(dāng)清運完前幾個點后，未達到飽和，但下一個點的垃圾量又裝不完，那么此車是直接返回呢，還是繼續(xù)裝直至車裝滿為止呢？我們判斷前者更好，就是車在裝的足夠多的情況下應(yīng)該直接返回原點。這是因為對于下一垃圾點(假設(shè)為A點)內(nèi)的垃圾而言，無論是一次裝完還是分兩次裝完，將它們運回所花費用是恒定的。整體而言，兩者花費的錢是相等的，但分兩次裝要多花裝車時間，所以選擇前者。5.1.3

21、.2模型的求解首先根據(jù)題所給的數(shù)據(jù)畫出散點圖 圖五求得總耗油為x，求解程序如附錄二，運輸車的最優(yōu)路線如下圖所示： 圖六表五:線路的站點序列和油耗站點序號油耗一號線0-30-29-27-3-018.4二號線0-28-26-32-25-5-017.6三號線0-36-23-33-21-016.8四號線0-24-18-35-15-013.6五號線0-34-17-16-2-012六號線0-20-11-10-011.2七號線0-19-13-8-010.8八號線0-14-7-4-1-08.8九號線0-22-08.4十號線0-12-9-08十一號線0-31-6-06.86 模型的評價和推廣（1）模型的評價好的方