【正文】 seif x1(k0,4)==5 n1(1,5)=n1(1,5)+1。y2=0。 end if s1s2 s1=s2。f2的直接運(yùn)行結(jié)果：f2s1 = ans = s1 = ans = s1 = ans = s1 = ans = s1 = ans = s1 = ans = s1 = ans = s1 = ans = s1 = ans = s1 = ans = s1 = ans = s1 = ans = s1 = ans = s1 = ans = ss = ans = %本函數(shù)存盤(pán)的函數(shù)名為“f3”a=[。sheet239。for k0=1:3 max0(1,k0)=max(x0(:,k0))。x0(:,3)=1./x0(:,3)。 elseif x1(k1,4)==5 x1(k1,5)=x1(k1,5)*y1(5,1)。for i=1:2 for j=1:7 if B(i,j)0 D(i,j)=C(j)*B(i,j)。for i=1:6 T(i+1)=(2*(E(i+1)E(i))/(F(i+1)+F(i)))/(24*60*60)。for k1=1:m if x1(k1,4)==1 x1(k1,5)=x1(k1,5)*y1(1,1)。y1=y0(:,1:3)*a。x1=zeros(m,5)。)。aa=zeros(3,1)。 y1=y0(:,1:3)*a。y0(:,2)=1./y0(:,2)。for k0=1:m if x1(k0,4)==1 n1(1,1)=n1(1,1)+1。max0=zeros(1,3)。)。201 2 3 3 2 2 4 2 4 2 2 4 302 2 2 212 3 1242503222422323223021302423222106234231603222422222因此該模型的移植性很強(qiáng)。 對(duì)于水質(zhì)污染的評(píng)價(jià)與預(yù)測(cè)，本文運(yùn)用了權(quán)重系數(shù)理論、衰減理論與曲線分析三種模型方法，針對(duì)不同問(wèn)題進(jìn)行處理。 模型檢驗(yàn)：在模型一的求解中，運(yùn)用十四次仿真計(jì)算Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、劣Ⅴ水中相鄰兩類(lèi)水的綜合評(píng)估標(biāo)準(zhǔn)值的平方差的比重和的最大值，實(shí)際上是對(duì)隨機(jī)產(chǎn)生的十四萬(wàn)組權(quán)重系數(shù)下的Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、劣Ⅴ水中相鄰兩類(lèi)水的綜合評(píng)估標(biāo)準(zhǔn)值的平方差的比重和的比較。③、加強(qiáng)處罰力度，由于我國(guó)對(duì)于企業(yè)污染排放超標(biāo)最多處罰10 萬(wàn)，因此對(duì)于企業(yè)的管理不會(huì)有效，應(yīng)該根據(jù)污染的程度加重處罰的力度。以下是具體思路和做法：(1)控制污染源根據(jù)前面的分析可以知道支流上的污染源對(duì)于水質(zhì)的影響很明顯，尤其是岷江、贛江、湘江支流對(duì)于長(zhǎng)江的污染最為嚴(yán)重所以應(yīng)該對(duì)這些流域的工廠企業(yè)的排污進(jìn)行監(jiān)控，對(duì)超標(biāo)行為進(jìn)行罰款，并通過(guò)媒體予以曝光。經(jīng)查詢可知污水的密度為：。將“長(zhǎng)江水污染附件4”中的水文年的三大類(lèi)綜合于下表八，表八：（近十年水文年的三大類(lèi)水百分比）時(shí)間類(lèi)型1995199619971998199920002001200220032004水文年ⅠIIⅢ類(lèi)干　流100支 流7475ⅣⅤ類(lèi)干　流0支 流181718劣Ⅴ類(lèi)干　流00000009支 流07先對(duì)水文年ⅠⅡⅢ類(lèi)的干支流進(jìn)行預(yù)測(cè)，取表中數(shù)據(jù)進(jìn)行分析模擬（程序見(jiàn)附錄中程序三）。長(zhǎng)江干流近一年多主要污染物氨氮的污染源主要在：湖北宜昌、湖南岳陽(yáng)、重慶朱沱。若該研究對(duì)象是屬于Ⅰ類(lèi)的，則將上一步處理結(jié)果乘以；若是Ⅱ類(lèi)的，將上一步處理結(jié)果乘以…… 模型一的解答： 本模型的MATLAB程序見(jiàn)附錄 程序一（第 頁(yè)），由“f2”進(jìn)行了十四次仿真，（必需經(jīng)過(guò)多次仿真，以作比較，去最好的）仿真結(jié)果如下： s1 第1次仿真結(jié)果： 第2次仿真結(jié)果： 第3次仿真結(jié)果： 第4次仿真結(jié)果： 第5次仿真結(jié)果： 第6次仿真結(jié)果： 第7次仿真結(jié)果： 第8次仿真結(jié)果： 第9次仿真結(jié)果： 第10次仿真結(jié)果： 第11次仿真結(jié)果： 第12次仿真結(jié)果： 第13次仿真結(jié)果： 第14次仿真結(jié)果： 由此，取仿真結(jié)果的第6組數(shù)據(jù)：= = 將上述數(shù)據(jù)帶入程序一“f3”（第 頁(yè)）求得國(guó)標(biāo)水質(zhì)標(biāo)準(zhǔn)Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、劣Ⅴ類(lèi)水的水質(zhì)綜合評(píng)估值分別為： 由此得到如下結(jié)論：①?lài)?guó)標(biāo)水質(zhì)標(biāo)準(zhǔn)及其綜合評(píng)估值，見(jiàn)下表三：表三 ：國(guó)標(biāo)水質(zhì)標(biāo)準(zhǔn)及其綜合評(píng)估值序號(hào) 分 類(lèi) 標(biāo)準(zhǔn)值項(xiàng) 目Ⅰ類(lèi)Ⅱ類(lèi)Ⅲ類(lèi)Ⅳ類(lèi)Ⅴ類(lèi)劣Ⅴ類(lèi)0綜合評(píng)估值 ≥1溶解氧(DO)通過(guò)比較不同權(quán)重系數(shù)下的Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、劣Ⅴ水中相鄰兩類(lèi)水的綜合評(píng)估標(biāo)準(zhǔn)值的平方差的比重和確定合理的一組權(quán)重系數(shù)：步驟一：將原始數(shù)據(jù)處理成DO值、CODMn值的倒數(shù)和NH3N值的倒數(shù)等三個(gè)指標(biāo)值以DO值、CODMn值的倒數(shù)和NH3N值的倒數(shù)等三個(gè)指標(biāo)值對(duì)表 給出的水質(zhì)標(biāo)準(zhǔn)數(shù)據(jù)進(jìn)行處理，結(jié)果如下： 步驟二：定義Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、劣Ⅴ類(lèi)水的水質(zhì)綜合評(píng)估標(biāo)準(zhǔn)值我們定義Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、劣Ⅴ類(lèi)水的水質(zhì)綜合評(píng)估標(biāo)準(zhǔn)值、和的計(jì)算：（、為所有研究對(duì)象第i種指標(biāo)的最大值，所所給研究對(duì)象的影響）之后用取代得到（i=1，2……6）（將其DO值、CODMn值的倒數(shù)和NH3N值的倒數(shù)等三個(gè)指標(biāo)值乘以對(duì)應(yīng)的權(quán)重，并經(jīng)一定處理）步驟三：規(guī)劃出最好的權(quán)重系數(shù)取值由于Ⅰ、Ⅱ類(lèi)水的比重占絕大部分，考慮Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、劣Ⅴ水的比重對(duì)相鄰兩類(lèi)水的綜合評(píng)估標(biāo)準(zhǔn)值的平方差的和的影響。 分 類(lèi) （4）　(單位：1/天)。 需要解決的問(wèn)題：預(yù)測(cè)模型。主要問(wèn)題為各指標(biāo)權(quán)重的確定。污染物在水環(huán)境中通過(guò)物理降解、化學(xué)降解和生物降解等使水中污染物的濃度降低，反映江河自然凈化能力的指標(biāo)稱(chēng)為降解系數(shù)。結(jié)論是：2006年廢水排放量為222億噸，2010年的為429億噸，2014年的為563億噸，可見(jiàn)廢水年排放量逐年上升，2006年至2015年的具體數(shù)據(jù)見(jiàn)表十（第17頁(yè)）、表十一(第17頁(yè))和表十二（第18頁(yè))。本文主要研究了以下四個(gè)問(wèn)題：長(zhǎng)江水質(zhì)的綜合評(píng)價(jià)、主要污染源的確定、預(yù)測(cè)問(wèn)題和污水處理問(wèn)題。我們鄭重承諾，嚴(yán)格遵守競(jìng)賽規(guī)則，以保證競(jìng)賽的公正、公平性。問(wèn)題二，我們建立模型二：針對(duì)高錳酸鹽和氨氮化合物，由于其在水中能自然降解，利用衰減理論計(jì)算各地區(qū)的排放量而得到長(zhǎng)江干流近一年多高錳酸鹽指數(shù)和氨氮的主要污染源。 長(zhǎng)江是我國(guó)第一、世界第三大河流，長(zhǎng)江水質(zhì)的污染程度日趨嚴(yán)重，已引起了相關(guān)政府部門(mén)和專(zhuān)家們的高度重視。（2）研究、分析長(zhǎng)江干流近一年多主要污染物高錳酸鹽指數(shù)和氨氮的污染源主要在哪些地區(qū)?（3）假如不采取更有效的治理措施，依照過(guò)去10年的主要統(tǒng)計(jì)數(shù)據(jù)，對(duì)長(zhǎng)江未來(lái)水質(zhì)污染的發(fā)展趨勢(shì)做出預(yù)測(cè)分析，比如研究未來(lái)10年的情況。水中氫離子濃度和溶氧量無(wú)需考慮。 （1）江水中的pH、DO、CODMn和NH3N等四項(xiàng)指標(biāo)相互獨(dú)立。將每個(gè)地區(qū)在每一月份的DO值、CODMn值的倒數(shù)和NH3N值的倒數(shù)等三個(gè)指標(biāo)值作為一個(gè)研究對(duì)象。 項(xiàng) 目Ⅰ類(lèi)Ⅱ類(lèi)Ⅲ類(lèi)Ⅳ類(lèi)Ⅴ類(lèi)劣Ⅴ類(lèi)1溶解氧(DO)水的類(lèi)別值越高，水質(zhì)越差，綜合評(píng)估值越低。這樣做是為了保證權(quán)值和為1，即化為最簡(jiǎn)形式。各觀測(cè)站點(diǎn)間水流速度取相鄰兩觀測(cè)站點(diǎn)的平均速度，即為：故各觀測(cè)站點(diǎn)間水流時(shí)間為： () 其中由降解系數(shù)（），運(yùn)用衰減理論，可得干流各觀測(cè)站點(diǎn)每秒排放CODMn()和NH3N()的濃度為： 其中，取 若 則有 否則 運(yùn)用“附錄表十八”中的數(shù)據(jù)，用Matlab程序（見(jiàn)附錄程序二）求解可得各觀測(cè)站點(diǎn)的每秒排污污染物CODMn和NH3N的質(zhì)量，見(jiàn)表五：表五：觀測(cè)站點(diǎn)一年多每秒排污污染物的質(zhì)量觀測(cè)站點(diǎn)四川攀枝花重慶朱沱湖北宜昌湖南岳陽(yáng)江西九江安徽安慶江蘇南京CODMnNH3N00CODMnNH3N0CODMn0NH3N00CODMn0NH3N00CODMn0NH3N000CODMn00NH3N000CODMn0NH3N00CODMnNH3N00CODMn0NH3N0CODMn0NH3N0CODMn0NH3N0CODMn0NH3N0CODMn0NH3N0注：（單位：10kg/s） 模型二解答： 用圖表方法將上表五的數(shù)據(jù)整理為圖一和圖二圖一：觀測(cè)站點(diǎn)排放物CODMn隨時(shí)間的變化圖二：觀測(cè)站點(diǎn)排放物隨時(shí)間的變化干流近一年多主要污染物高錳酸鹽指數(shù)和氨氮的污染源主要所在地即為其排放量的前幾名。 圖三：（未來(lái)十年廢水排放量預(yù)測(cè)圖）得到擬合二次多項(xiàng)式：然后對(duì)水文年干支流總的河長(zhǎng)進(jìn)行預(yù)測(cè)，把1