數(shù)學(xué)建模學(xué)期論文若干工廠排污費(fèi)用的數(shù)學(xué)模型

1、選修課程結(jié)業(yè)考試答卷 課程名稱： 數(shù) 學(xué) 建 模 選課課號(hào)：（2009-2010-2)-161100200-1020060061-1 主講教師： 學(xué) 生： 學(xué) 號(hào)： 班 級(jí)： 考查時(shí)間： 考察試題： 若干工廠排污費(fèi)用的數(shù)學(xué)模型 考查評(píng)分： 若干工廠排污費(fèi)用的數(shù)學(xué)模型摘要本文首先對(duì)給定的題目進(jìn)行統(tǒng)計(jì)分析，在分析結(jié)果的基礎(chǔ)上，考慮了各種不同的污水排放的差異性，將所有污水處理信息進(jìn)行整合，為污水處理情形建立相關(guān)數(shù)據(jù)材料。基于不同情況不同特點(diǎn)，分別對(duì)問題一到問題二進(jìn)行分析、建模及求解。對(duì)問題一要求：為了使使江面所有地段得水污染達(dá)到國(guó)家標(biāo)準(zhǔn)，最少需要費(fèi)用采用方法：在一定合情理得假設(shè)下，建立一個(gè)簡(jiǎn)單普遍的

2、的模型，使之能夠清晰地表示出江水與第k個(gè)站點(diǎn)得污水混合后的污水濃度，并使?jié)舛刃∮跇?biāo)準(zhǔn)，得到普遍的約束條件。再用一個(gè)函數(shù)表示出處理費(fèi)用，得到普遍的目標(biāo)函數(shù)，再將題中數(shù)據(jù)代入公式，編寫程序，再用lindo軟件解出答案。對(duì)問題二， 要求：要求三個(gè)據(jù)名店上游的水污染達(dá)到國(guó)家標(biāo)準(zhǔn)，最少需要多少費(fèi)用 采用方法：在基于問題一得假設(shè)下，建立一個(gè)簡(jiǎn)單普遍的的模型，使之能夠清晰地表示出第k個(gè)站點(diǎn)上游的污水污水濃度，并使?jié)舛刃∮跇?biāo)準(zhǔn)，得 到普遍的約束條件。再用一個(gè)函數(shù)表示出處理費(fèi)用，得到普遍的目標(biāo)函 數(shù)，再將題中數(shù)據(jù)代入公式，編寫程序，再用lindo軟件解出答案。根據(jù)上述研究，我們認(rèn)為我們建立的數(shù)學(xué)模型易于操作，對(duì)

3、實(shí)踐有著較好的指導(dǎo)意義。關(guān)鍵詞： 若干工廠 污水處理 國(guó)家標(biāo)準(zhǔn) 最小費(fèi)用一、問題重述工廠1上游江水流量和污水濃度，國(guó)家標(biāo)準(zhǔn)規(guī)定的水的污染濃度，以及各個(gè)工廠的污水流量和污水濃度均已知道。設(shè)污水處理費(fèi)用與污水處理前后的濃度差和污水流量成正比，使每單位流量的污水下降一個(gè)濃度單位需要的處理費(fèi)用（稱處理系數(shù)）為已知。處理后的污水與江水混合，流到下一個(gè)排污口之前，自然狀態(tài)下的江水也會(huì)使污水濃度降低一個(gè)比例系數(shù)（稱自凈系數(shù)），該系數(shù)可以估計(jì)。試確定各污水處理站出口的污水濃度，使在符合國(guó)家標(biāo)準(zhǔn)規(guī)定的條件下總的處理費(fèi)用最小。推而廣之，要求處理一段河流中有若干個(gè)工廠、若干個(gè)排污處理站點(diǎn)。建立數(shù)學(xué)模型，求解。二、問

4、題的分析工廠的污水處理與排放都非常熟悉的現(xiàn)象，它以這樣或那樣的形式影響著我們的生活。例如，我們的飲水、我們周邊的生態(tài)環(huán)境。污水排放過(guò)程要經(jīng)過(guò)三個(gè)步驟：工業(yè)生產(chǎn)制造廢水、污水處理站處理廢水、國(guó)家環(huán)保部門監(jiān)測(cè)排污情況。圖1：工廠的污水處理與排放過(guò)程描述工廠的污水排放必須達(dá)到國(guó)家相應(yīng)的標(biāo)準(zhǔn)，否則將受到環(huán)境保護(hù)部門的處理，嚴(yán)重者甚至受到法律的制裁。反思問題的題設(shè)與要求，是要解決污水處理過(guò)程中的費(fèi)用問題，以求達(dá)到最好的經(jīng)濟(jì)效益。從以下三個(gè)個(gè)角度思考：環(huán)保部門的角度，工廠排放的廢水必須達(dá)到國(guó)家相關(guān)的標(biāo)準(zhǔn)，且要要在環(huán)境的自動(dòng)凈化能力范圍之內(nèi)。工廠生產(chǎn)角度，工廠為了生產(chǎn)發(fā)展，必須降低成本，盡可能地實(shí)現(xiàn)增加利潤(rùn)

5、。這就要求要盡可能的減少排污費(fèi)用。居民角度，居民要求自己的居住環(huán)境盡可能的好，這就要求周圍事物盡可能的減少環(huán)境廢物的排放。對(duì)于問題一：為了使江面所有的所有的地段達(dá)到國(guó)家的標(biāo)準(zhǔn) ，由題意自然狀態(tài)下江水可使污水濃度下降一定程度，所以如果江面各處污水濃度最大的地方達(dá)到國(guó)家標(biāo)準(zhǔn)，那么各處江面都將達(dá)到國(guó)際按標(biāo)準(zhǔn)，于是顯而易見，只要各處處理站處理后的污水與江河上游污水混合后達(dá)標(biāo)即可，順次思路，可比較簡(jiǎn)單的建立數(shù)學(xué)模型解決此問題。對(duì)于問題二：為使居民點(diǎn)上游水污染達(dá)標(biāo)，由于江面存在自凈系數(shù)，即每處處理站排出的污水與江水混合后得污水濃度乘以該站下游江面的自凈系數(shù)便可得到下游居民點(diǎn)的上游的污水濃度，由此思路建立可

6、較為簡(jiǎn)易的建立合理的數(shù)學(xué)模型。三、模型假設(shè)1、 假設(shè)題中給定的數(shù)據(jù)真實(shí)可靠，且滿足國(guó)家相關(guān)標(biāo)準(zhǔn)。2、 污水經(jīng)站點(diǎn)排除后迅速均勻分布于水中。3、 假設(shè)該江面除污水處理站外，無(wú)其他水源入江。4、 工廠與處理站的污水流出量遠(yuǎn)小于江水流量。 四、符號(hào)說(shuō)明符號(hào)符號(hào)說(shuō)明n工廠的總數(shù)、處理站點(diǎn)數(shù)（n的取值為整數(shù)）k從上游到下游計(jì)數(shù)的工廠個(gè)數(shù) 第k處理站口上游江面的流量第k處理站口下游江面的流量第k處理站口上游江面的污水濃度第k處理站口下游江面的污水濃度第k個(gè)工廠排出污水濃度第k+1個(gè)工廠排出污水濃度第k個(gè)處理站的處理系數(shù)第k+1個(gè)處理站的處理系數(shù)第k個(gè)工廠的污水經(jīng)第k個(gè)處理站后排到江面的污水濃度第k+1個(gè)工

7、廠的污水經(jīng)第k+1個(gè)處理站后排到江面的污水濃度第k段江面的自凈系數(shù)（第k段江面為第k個(gè)工廠與第k+1個(gè)工廠之間的江面）第k+1段江面的自凈系數(shù)（第k段江面為第k個(gè)工廠與第k+2個(gè)工廠之間的江面） 第k個(gè)工廠與第k個(gè)處理站的污水流出量第k+1個(gè)工廠與第k+1個(gè)處理站的污水流出量國(guó)家標(biāo)準(zhǔn)規(guī)定的合格水污染濃度s污水處理總費(fèi)用五、問題求解及結(jié)果分析l 問題一1. 問題一的模型分析問題一：為了使江面上所有低端的水污染達(dá)到國(guó)家標(biāo)準(zhǔn)，且費(fèi)用最少，并求出最少費(fèi)用。在簡(jiǎn)化為常數(shù)q的情況下，第k個(gè)處理站的污水與江水混合后的污水的濃度為=+，第k段江面的自凈后的污水濃度為=，污水處理站的處理費(fèi)用為（-），于是在江面

8、上所有地段的水污染達(dá)到國(guó)家標(biāo)準(zhǔn)的條件下，使費(fèi)用最小。2. 建立模型目標(biāo)函數(shù)： s=（-）約束條件：=+ = 3. 問題一模型求解對(duì)于該題的具體題設(shè)，q=1000l/min, =0.8mg/l, =1mg/l,=5l/min，=100 mg/l，=60 mg/l，=50mg/l,=1萬(wàn)元/（l/min）（mg/l）, =0.9,=0.6。按照以上模型，并用lingo軟件進(jìn)行計(jì)算求解。計(jì)算結(jié)果為：=40 mg/l，=20 mg/l，=50 mg/l。s=495.5萬(wàn)元4. 結(jié)果分析:站點(diǎn)一處理后污水濃度為40mg/l，站點(diǎn)二處理后污水濃度20mlg/l，站點(diǎn)3不用處理，且最小費(fèi)用為495.5萬(wàn)元，

9、其影子價(jià)格和靈敏度分析詳見附件l 問題二問題二：題設(shè)只要求三個(gè)居名點(diǎn)上游的水污染達(dá)到國(guó)家標(biāo)準(zhǔn)，就最少費(fèi)用。1. 問題二的模型分析在簡(jiǎn)化為常數(shù)q的情況下，第k個(gè)處理站的污水與江水混合后的污水的濃度為=+，第k段江面的自凈后的污水濃度為=，污水處理站的處理費(fèi)用為（-），于是在江面上所有地段的水污染達(dá)到國(guó)家標(biāo)準(zhǔn)的條件下，使費(fèi)用最小。2. 建立模型目標(biāo)函數(shù)： s=（-）約束條件：=+ = 3. 問題二模型求解對(duì)于該題的具體題設(shè)，q=1000l/min, =0.8mg/l, =1mg/l,=5l/min，=100 mg/l，=60 mg/l，=50mg/l,=1萬(wàn)元/（l/min）（mg/l）, =0.

10、9,=0.6。按照以上模型，并用lingo軟件進(jìn)行計(jì)算求解。計(jì)算結(jié)果為：=62.2 mg/l，=60 mg/l，=50 mg/l。s=189萬(wàn)元計(jì)結(jié)果分析l 第一個(gè)工廠處理后污水濃度為62.2mg/l，第二個(gè)工廠和第三個(gè)工廠無(wú)需處理污水，且費(fèi)用最少為189萬(wàn)元。其影子價(jià)格和靈敏度分析詳見附件六、模型評(píng)價(jià)及改進(jìn)、模型推廣 污染是當(dāng)今社會(huì)的一個(gè)重大問題，尤其是水污染對(duì)人的生存構(gòu)成了嚴(yán)重威脅，而河流擴(kuò)散污染的主要途徑，一方面要考慮下游居民，另一方面還要?jiǎng)?wù)必使成本降到最低，以達(dá)到最大效益。怎樣才能使河流既達(dá)到國(guó)家標(biāo)準(zhǔn)而有用資金最少？我們便考慮使用數(shù)學(xué)建模。數(shù)學(xué)建模是一種數(shù)學(xué)的思考方法，是運(yùn)用數(shù)學(xué)的語(yǔ)言

11、和方法，通過(guò)抽象、簡(jiǎn)化建立能近似刻畫并解決實(shí)際問題的一種強(qiáng)有力的數(shù)學(xué)手段。基于以上概念。而江河與處理站，江河與居民點(diǎn)可以近似的考慮成線與點(diǎn)得關(guān)系，這樣就可以得到最基本的模型，然后在做出相應(yīng)的假設(shè)之后，考慮點(diǎn)與點(diǎn)之間的相互影響，點(diǎn)與線之間的相互關(guān)系，已達(dá)到比較真實(shí)的模擬出實(shí)際情況來(lái)，于是，得出以上的模型設(shè)計(jì)，該模型再基于假設(shè)條件之內(nèi)基本能夠反映真實(shí)，在已知各段江水的自凈系數(shù)前提下，可以求出有任意k個(gè)排污點(diǎn)及處理站時(shí)得最小費(fèi)用，但是該模型存在下列問題1， 自凈系數(shù)沒有明確的得知方法，應(yīng)該考慮自凈系數(shù)與江段長(zhǎng)度的關(guān)系。2， 污水的與江水混合后求濃度沒有考慮污水自身的流量3， 江水經(jīng)過(guò)一段距離后得由于

12、蒸發(fā)和滲透沒有做考慮可以針對(duì)上述問題再設(shè)變量對(duì)模型改進(jìn)。 模型的推廣：上述模型是建立在無(wú)污水排放時(shí)濃度是符合國(guó)家標(biāo)準(zhǔn)的，那么如果上游江水本就不符合國(guó)家標(biāo)準(zhǔn)，那么就需要處理江中污水濃度，亦或，先將江水的污水濃度下降一定濃度時(shí)，再同各個(gè)處理站處理后的污水混合，成本會(huì)不會(huì)有所降低呢？而江的上游會(huì)有一段比較長(zhǎng)得無(wú)人區(qū)，可用較低的成本較簡(jiǎn)單的方法處理得到較好的效果，比如撒一些除污物質(zhì)等等可設(shè)江水處理系數(shù)。那么依然可用上述模型，那么目標(biāo)函數(shù)即變?yōu)?s=（-）+（-）約束條件=+ = 七、參考文獻(xiàn)1) 數(shù)學(xué)模型（第三版） 姜啟源 謝金星 葉?。ň幹?高等教育出版社 2003.82) 數(shù)學(xué)建模課程講義 王麗

13、 西南科技大學(xué) 3) 基礎(chǔ)運(yùn)籌學(xué)教程(配光盤） 馬良（主編） 高等教育出版社 2006.7八、附件清單、附件1) 附件清單 問題一的詳細(xì)數(shù)據(jù)、lindo運(yùn)算程序、程序運(yùn)行結(jié)果。 問題二的詳細(xì)數(shù)據(jù)、lindo運(yùn)算程序、程序運(yùn)行結(jié)果。 2) 附件 問題一的詳細(xì)數(shù)據(jù)、lindo運(yùn)算程序、程序運(yùn)行結(jié)果。a) 詳細(xì)數(shù)據(jù):q=1000l/min, =0.8mg/l, =1mg/l,=5l/min，=100 mg/l，=60 mg/l，=50mg/l,=1萬(wàn)元/（l/min）（mg/l）, =0.9,=0.6。得目標(biāo)函數(shù)s=1050-5-5-5代入公式=+得站點(diǎn)一混合后污水濃度為: =5/1000+0.8站

14、點(diǎn)二混合后污水濃度為: =4.5/1000+5/1005+0.72站點(diǎn)三混合后污水濃度為: =2.7/1000+3/1005+5/1010+0.432再代入公式并化簡(jiǎn)可得約束條件50 4.5+4.97528 2.7+2.985+4.95568代入公式得到約束條件 100 60 50b) lindo運(yùn)算程序:min -5v1-5v2-5v3st2)v1=403)4.5v1+4.975v2=2804)2.7v1+2.985v2+4.95v3=5685)v1=1006)v2=607)v3=50endc) 程序運(yùn)行結(jié)果:lp optimum found at step 3 objective func

15、tion value 1) -550.5025 variable value reduced cost v1 40.000000 0.000000 v2 20.100502 0.000000 v3 50.000000 0.000000 row slack or surplus dual prices 2) 0.000000 0.477387 3) 0.000000 1.005025 4) 152.500000 0.000000 5) 60.000000 0.000000 6) 39.899498 0.000000 7) 0.000000 5.000000 no. iterations= 3 r

16、anges in which the basis is unchanged: obj coefficient ranges variable current allowable allowable coef increase decrease v1 -5.000000 0.477387 infinity v2 -5.000000 5.000000 0.527778 v3 -5.000000 5.000000 infinity righthand side ranges row current allowable allowable rhs increase decrease 2 40.0000

17、00 22.222221 40.000000 3 280.000000 198.500000 100.000000 4 568.000000 infinity 152.500000 5 100.000000 infinity 60.000000 6 60.000000 infinity 39.899498 7 50.000000 30.808083 50.000000問題二的詳細(xì)數(shù)據(jù)、lindo運(yùn)算程序、程序運(yùn)行結(jié)果。a) 詳細(xì)數(shù)據(jù):由于居民點(diǎn)1上游污水濃度為0.8mg/l符合標(biāo)準(zhǔn)，所以不作考慮。處理站3下游沒有居民點(diǎn)，所以得： =60代入數(shù)據(jù)后得目標(biāo)函數(shù)為s=（100-）*5+（60-）*5

18、=800-5-5代入公式=+得=4.5/1000+0.72=2.7/1000+3/1005+0.432再帶入公式 并化簡(jiǎn)得到約束條件4.52802.7+2.985568代入公式得到約束條件10060 b) lindo運(yùn)算程序:min -5v1-5v2st2)4.5v1=2803)2.7v1+2.985v2=5684)v1=1005)v2=60endc) 程序運(yùn)行結(jié)果:lp optimum found at step 2 objective function value 1) -611.1111 variable value reduced cost v1 62.222221 0.000000 v2 60.000000 0.000000 row slack or surplus dual prices 2) 0.000000 1.111111 3) 220.900009 0.000000 4) 37.777779 0.000000 5) 0.000000 5.000000 no. iterations= 2 ranges in which the basis is unchanged: obj coefficient ranges variable current allowable allowable coef increase decrease v1 -5.