【正文】 mn=zeros(17,5)。 function myfun1load data[m,n]=size(d)。4. 鄭阿奇，MATLAB實(shí)用教程，北京：電子工業(yè)出版社，2004年。八、參考文獻(xiàn)1. 謝季堅(jiān)、劉承平，模糊數(shù)學(xué)方法及其應(yīng)用，武漢：華中科技大學(xué)出版社，2000年2. 徐祖信等， 基于數(shù)學(xué)模型的蘇州河上游和支流水質(zhì)對(duì)干流的影響分析 ，水動(dòng)力學(xué)研究與發(fā)展，A輯第19卷第6期，733——743，2004年。建立健全的法規(guī)制度,加大對(duì)違章排污的處罰力度,尤其要從重處罰那些污水處理設(shè)備安裝后不啟用并偷排廢水的企業(yè)。對(duì)未治理或治理后未達(dá)標(biāo)的工業(yè)廢水一律不準(zhǔn)向湖泊水域直接排放。對(duì)每個(gè)污染源必須進(jìn)行有效的治理,各排污單位要根據(jù)其廢水的特點(diǎn),選擇相應(yīng)的污水處理工藝和設(shè)備對(duì)污水進(jìn)行治理,達(dá)到排放標(biāo)準(zhǔn)后才能排放。干流重點(diǎn)整治是三峽庫(kù)區(qū)段、湖北宜昌至湖南岳陽(yáng)江段沿岸、武漢至九江段、各城市江段。就重點(diǎn)治理而言,先進(jìn)行飲水源治理,后是長(zhǎng)江沿岸污染帶、長(zhǎng)江一級(jí)支流。既抓全長(zhǎng)江的整治規(guī)劃,又要突出重點(diǎn),解決緊迫的水質(zhì)問(wèn)題。這幾個(gè)江段屬于干流中水質(zhì)較差,目前污染又有加重趨勢(shì)的江段,需引起足夠重視。七、建議和意見(jiàn)1．加強(qiáng)長(zhǎng)江水質(zhì)變化監(jiān)測(cè)。2．建立穩(wěn)態(tài)一唯水質(zhì)模型時(shí)，我們假設(shè)了任意觀測(cè)點(diǎn)的污染物濃度不在斷面點(diǎn)處發(fā)生對(duì)流擴(kuò)散，即與此觀測(cè)點(diǎn)時(shí)間無(wú)關(guān)。不足之處及改進(jìn)方向：1．對(duì)未來(lái)10年污水處理量的粗略預(yù)測(cè)中，忽略了自然降解的情況，得出的結(jié)果可能比實(shí)際值偏大。4．根據(jù)前述GM（1，1）模型，利用簡(jiǎn)單的比例關(guān)系粗略地預(yù)測(cè)未來(lái)10年中每年的污水處理量。利用此模型對(duì)長(zhǎng)江流域的飲用水類(lèi)河長(zhǎng)比例進(jìn)行了預(yù)測(cè)，結(jié)果優(yōu)良。3．對(duì)未來(lái)進(jìn)行水質(zhì)預(yù)測(cè)，是治理長(zhǎng)江的迫切需要。通過(guò)修正穩(wěn)態(tài)一唯對(duì)流擴(kuò)散水質(zhì)模型，得到在自然降解情況下的任意點(diǎn)的污染物濃度計(jì)算公式，由此建立了適合本題的穩(wěn)態(tài)一唯水質(zhì)模型，得出了六段干流的污染情況。由于水質(zhì)污染程度是一個(gè)模糊的概念，用隸屬度來(lái)劃分分段界限比較合適，因此評(píng)判時(shí)用了單因數(shù)評(píng)價(jià)矩陣，結(jié)果令人滿意?；谶^(guò)去10年統(tǒng)計(jì)數(shù)據(jù)，分別以干流河長(zhǎng)和全流域河長(zhǎng)的均值作為未來(lái)10年的評(píng)價(jià)河長(zhǎng) （L干：干流評(píng)價(jià)河長(zhǎng)；L全：全流域評(píng)價(jià)河長(zhǎng)），其數(shù)學(xué)計(jì)算式可描述如下：其中：Li表示第年的污水處理量；L(1)i表示第年需要治理的Ⅳ、Ⅴ類(lèi)河長(zhǎng)，ri表示第i年Ⅳ、Ⅴ類(lèi)河長(zhǎng)的預(yù)測(cè)比例；L(2)i表示第年需要治理的Ⅴ類(lèi)河長(zhǎng)， L(2)I=L干titi表示第i年Ⅴ類(lèi)河長(zhǎng)的預(yù)測(cè)比例Pi表示第i年污水年均排放量的預(yù)測(cè)值；qi表示第年全流域非飲用水類(lèi)的預(yù)測(cè)比例。由于污水是非飲用水類(lèi)，即每年的污水主要分布在全流域Ⅳ、Ⅴ類(lèi)和劣Ⅴ類(lèi)水的河段上。（三）問(wèn)題四根據(jù)問(wèn)題三的預(yù)測(cè)，長(zhǎng)江水體的污染將日益嚴(yán)重，急待解決的是如何使未來(lái)10年內(nèi)每年的干流Ⅳ、Ⅴ類(lèi)水的比例控制在20%以?xún)?nèi)，且沒(méi)有劣Ⅴ類(lèi)水。到2014年，%，即近一半為非飲用水河段。其余相對(duì)誤差都小于9%，可以認(rèn)為使用GM（1，1）模型預(yù)測(cè)效果是令人滿意的。模型建立：1．令為GM（1，1）建模序列，表示灰導(dǎo)數(shù)其中，2．令為的AGO序列，3．令為的均值（MEAN）序列，表示白化背景值則得到GM（1，1）的灰微分方程模型為其中，經(jīng)變換后得到模型求解：=11,12…20時(shí)的（程序見(jiàn)附件6），變換后得到最后預(yù)測(cè)出長(zhǎng)江流域未來(lái)十年水文年可飲用水河長(zhǎng)比例。由于1998年長(zhǎng)江流域遭遇罕見(jiàn)的洪水災(zāi)害，數(shù)據(jù)變化很大，因此我們先將長(zhǎng)江1998年的數(shù)據(jù)作為災(zāi)變數(shù)據(jù)剔除。準(zhǔn)確的水質(zhì)預(yù)測(cè)能讓人們及早地發(fā)現(xiàn)潛伏著的自然災(zāi)害。綜上所述：長(zhǎng)江干流近一年多高錳酸鹽（CODMn）和氨氮（NH3N）的主要污染源是：湖北宜昌至湖南岳陽(yáng)江段沿岸、岷江流域的四川樂(lè)山地區(qū)、重慶朱沱至湖北宜昌江段沿岸。 ② 岷江流域兩觀測(cè)點(diǎn)污染比較 CODMn濃度比較：圖4：四川樂(lè)山與四川宜賓CODMn污染濃度比較結(jié)論：近一年中，四川樂(lè)山的CODMn污染濃度普遍高于四川宜賓，其中第11月屬于污染狀態(tài)。b. 鄱陽(yáng)湖對(duì)上游NH3N污染的凈化作用更加明顯?？梢?jiàn)洞庭湖湖畔工廠對(duì)洞庭湖水質(zhì)基本沒(méi)有NH3N污染，湖水湖畔入口流域的NH3N污染起到了凈化作用。 NH3N濃度比較：圖3：兩湖進(jìn)出口NH3N污染濃度比較 結(jié)論：a. 洞庭湖出口處NH3N濃度減小，從Ⅲ類(lèi)水變成了Ⅱ類(lèi)水。b. 鄱陽(yáng)湖出入口的水質(zhì)按照國(guó)家質(zhì)量標(biāo)準(zhǔn)，均屬于Ⅱ類(lèi)，為可飲用水?？梢?jiàn)洞庭湖湖畔工廠對(duì)洞庭湖水質(zhì)中的CODMn指標(biāo)有清度污染。為了確定具體的污染源地區(qū)，我們采用作圖法：① 兩湖區(qū)進(jìn)出口處污染比較CODMn濃度比較：圖2：兩湖進(jìn)出口CODMn污染濃度比較結(jié)論：a. 洞庭湖出口處CODMn濃度增大，從Ⅱ類(lèi)水變成了Ⅲ類(lèi)水。利用（4）式分別對(duì)高錳酸鹽指數(shù)和氨氮進(jìn)行計(jì)算（程序見(jiàn)附件5），計(jì)算結(jié)果：表4：長(zhǎng)江六個(gè)江段的污染總量（單位：萬(wàn)噸）江段122334455667CODMn江段122334455667NH3N由上表的數(shù)據(jù)可知：對(duì)于高錳酸鹽（CODMn）和氨氮（NH3N）兩種污染，湖北宜昌南津關(guān)至湖南岳陽(yáng)城陵磯江段的污染最嚴(yán)重，其次為重慶朱沱至湖北宜昌南津關(guān)，和四川攀枝花重慶朱沱。表示干流第i個(gè)觀察站（地區(qū)）的污染經(jīng)過(guò)自然降解到下一個(gè)觀測(cè)站（地區(qū)）時(shí)的剩余污染量。 ② 利用（3）式，建立反映六江段污染物總量的模型考查干流上相鄰兩點(diǎn)間的污染情況，即在考慮自然降解的情況下這段流域上支流和干流的污染狀況。s1與江水的自然降解和污染物濃度有關(guān)?；谏鲜龇治?，可利用污染物質(zhì)在河流中運(yùn)動(dòng)變化的基本模型——“穩(wěn)態(tài)一維對(duì)流擴(kuò)散模型” [2]：其中，c為整個(gè)斷面的平均濃度u為流速，E為離散系數(shù) x為河道長(zhǎng)度t為時(shí)間ss2分別反映內(nèi)部降解變化和外部影響 由于長(zhǎng)江水流狀態(tài)一般為穩(wěn)定、均勻流動(dòng)狀態(tài)，與江長(zhǎng)相比，水面寬度很小，可認(rèn)為污染物濃度在橫斷面上分布比較均勻，即對(duì)同一地點(diǎn)進(jìn)行水質(zhì)采樣在無(wú)外界環(huán)境的影響下，可認(rèn)為污染物濃度不隨時(shí)間變化。其次水流速度對(duì)污染物的濃度也有明顯的影響。所以DO的權(quán)重賦值取C i/ Si的倒數(shù)，即其中，Ii （i =1,2,3）表示第i種污染物的權(quán)重；Ci 表示第i種污染物的濃度實(shí)測(cè)均值；Si為第i種污染物的六個(gè)級(jí)別濃度標(biāo)準(zhǔn)限值的均值。這里僅需要考慮DO、CODMn，NH3N這三種污染物對(duì)水質(zhì)的影響，對(duì)它們賦予不同的權(quán)重Wi。權(quán)重是衡量因子集中某一因子對(duì)水質(zhì)污染程度影響相對(duì)大小的量,權(quán)重系數(shù)越大,則該因子對(duì)水質(zhì)的影響程度越大。以DO的監(jiān)測(cè)值為自變量x,對(duì)第j級(jí)別的隸屬度為（水質(zhì)標(biāo)準(zhǔn)見(jiàn)附錄8） 將各監(jiān)測(cè)斷面的監(jiān)測(cè)數(shù)據(jù)代入前面確定的隸屬函數(shù)中，就可以計(jì)算其隸屬度，進(jìn)而建立每個(gè)斷面的單因子模糊評(píng)價(jià)矩陣（程序見(jiàn)附錄2）。由于劣Ⅴ類(lèi)污染物指標(biāo)較極端，我們根據(jù)標(biāo)準(zhǔn)將水質(zhì)分成5級(jí)。下面建立模糊綜合評(píng)價(jià)模型來(lái)解決兩年多的水質(zhì)綜合評(píng)價(jià)問(wèn)題。模糊綜合評(píng)價(jià)是由給出的評(píng)價(jià)標(biāo)準(zhǔn)和實(shí)測(cè)值,經(jīng)過(guò)模糊變換,對(duì)待評(píng)價(jià)對(duì)象給出總的評(píng)價(jià)的一種方法。三、模型假設(shè)1． 長(zhǎng)江干流的自然凈化能力近似均勻，.2． 所有的數(shù)據(jù)都來(lái)自科學(xué)采集3． 入河排污口水量與水質(zhì)變化穩(wěn)定4． 長(zhǎng)江各觀測(cè)河段的水流水質(zhì)狀態(tài)穩(wěn)定5． 四川攀枝花上流的污染忽略。對(duì)于問(wèn)題四，長(zhǎng)江水體的污染日益嚴(yán)重，急待解決的是如何使未來(lái)10年內(nèi)每年的干流Ⅳ、Ⅴ類(lèi)水的比例控制在20%以?xún)?nèi)，且沒(méi)有劣Ⅴ類(lèi)水。然而，由于只有過(guò)去十年的統(tǒng)計(jì)數(shù)據(jù)，甚至要預(yù)測(cè)10年的污染狀況。對(duì)于問(wèn)題三，根據(jù)過(guò)去10年的主要統(tǒng)計(jì)數(shù)據(jù)，要對(duì)未來(lái)10年的水質(zhì)污染作出科學(xué)的預(yù)測(cè)，這也是治理長(zhǎng)江的迫切需要。污染物在河流中符合一級(jí)反應(yīng)動(dòng)力學(xué)，可以建立穩(wěn)定一維均勻河流