長(zhǎng)江水質(zhì)的評(píng)價(jià)和預(yù)測(cè)摘要水是生命之源，保護(hù)水就是保護(hù)我們自己，保護(hù)水的重中之重就是保護(hù)大江大河。本文對(duì)近兩年的水質(zhì)分析，綜合評(píng)價(jià)，得出了部分地區(qū)的水質(zhì)污染情況，并根據(jù)十年的數(shù)據(jù)，對(duì)未來十年水質(zhì)污染發(fā)展趨勢(shì)做了預(yù)測(cè)，本文可以得出結(jié)論：保護(hù)母親河的行動(dòng)迫在眉睫！對(duì)于問題一，為了便于綜合評(píng)價(jià)，本文設(shè)出了綜合水質(zhì)標(biāo)識(shí)指數(shù)和單因子水質(zhì)標(biāo)識(shí)指數(shù)（具體公式計(jì)算見模型建立與求解），我們通過對(duì)單個(gè)城市28個(gè)月的綜合的評(píng)價(jià)標(biāo)識(shí)指數(shù)求平均值,數(shù)據(jù)如下（1.95222.1162.23012.41842.10192.25152.04483.54692.25092.75411.78032.8682.56282.3923.58882.44352.3802），綜合的評(píng)價(jià)標(biāo)識(shí)指數(shù)平均值越大，表示污染越嚴(yán)重。對(duì)于問題二，為了判斷主要污染源分布地區(qū)，本文采取判斷本地排放主要污染物的量，十三個(gè)月的求和取平均值來斷定主要污染源。計(jì)算數(shù)據(jù)用數(shù)列表示如下：當(dāng)為高錳酸鹽指數(shù)時(shí)，（8.986，37.1748，50.907，70.4526，58.196，59.9114，58.259）當(dāng)為氨氮時(shí)，（0.4816，3.0496，4.1418，6.3864，5.0473，5.0276，2.4794）取該數(shù)據(jù)較大的幾個(gè)為污染源，為主要污染源分布地區(qū)，結(jié)果如下：高錳酸鹽指數(shù)和氨氮的污染源主要所在地分別為：湖南岳陽(yáng)城陵磯

，江西九江河西水廠，

安徽安慶皖河口,江蘇南京林山四地；湖南岳陽(yáng)城陵磯

，江西九江河西水廠，

安徽安慶皖河口三地。對(duì)與問題三，對(duì)為來十年的排污量進(jìn)行預(yù)測(cè)時(shí)，建立了灰色系統(tǒng)模型。對(duì)這十年的預(yù)測(cè)值如下：（322．5221 343．2881 365.3912 388.9175 413.9585 440.6118 468.9812 303．123 499.1772 531.3174）對(duì)于問題四，本文根據(jù)第三問對(duì)將來十年廢水排放的預(yù)測(cè)值建立了廢水排放與IV、V類水的百分比之間的關(guān)系，Ⅳ，我們建立了百分比y與廢水派放量x之間的關(guān)系y=f(x),令y≤20,求出x的上限，則預(yù)測(cè)的廢水排放量與x的上限的差值即為需要處理的污水，從而將IV、V類水的百分比控制在20%，劣V類為0，求出了每年需要處理的污水量。對(duì)于問題五，本文參考以上問題得出的數(shù)據(jù)，并參考一些文獻(xiàn)資料，呼吁保護(hù)長(zhǎng)江人人有責(zé)，保護(hù)長(zhǎng)江一定要采取行之有效的行動(dòng)！一、問題重述水是人類賴以生存的資源，保護(hù)水資源就是保護(hù)我們自己，對(duì)于我國(guó)大江大河水資源的保護(hù)和治理應(yīng)是重中之重。專家們呼吁：“以人為本，建設(shè)文明和諧社會(huì)，改善人與自然的環(huán)境，減少污染?！遍L(zhǎng)江是我國(guó)第一、世界第三大河流，長(zhǎng)江水質(zhì)的污染程度日趨嚴(yán)重，已引起了相關(guān)政府部門和專家們的高度重視。2004年10月，由全國(guó)政協(xié)與中國(guó)發(fā)展研究院聯(lián)合組成“保護(hù)長(zhǎng)江萬里行”考察團(tuán)，從長(zhǎng)江上游宜賓到下游上海，對(duì)沿線21個(gè)重點(diǎn)城市做了實(shí)地考察，揭示了一幅長(zhǎng)江污染的真實(shí)畫面，其污染程度讓人觸目驚心。為此，專家們提出“若不及時(shí)拯救，長(zhǎng)江生態(tài)10年內(nèi)將瀕臨崩潰”（附件１），并發(fā)出了“拿什么拯救癌變長(zhǎng)江”的呼喚（附件2）。附件3給出了長(zhǎng)江沿線17個(gè)觀測(cè)站（地區(qū)）近兩年多主要水質(zhì)指標(biāo)的檢測(cè)數(shù)據(jù)，以及干流上７個(gè)觀測(cè)站近一年多的基本數(shù)據(jù)（站點(diǎn)距離、水流量和水流速）。通常認(rèn)為一個(gè)觀測(cè)站（地區(qū)）的水質(zhì)污染主要來自于本地區(qū)的排污和上游的污水。一般說來，江河自身對(duì)污染物都有一定的自然凈化能力，即污染物在水環(huán)境中通過物理降解、化學(xué)降解和生物降解等使水中污染物的濃度降低。反映江河自然凈化能力的指標(biāo)稱為降解系數(shù)。事實(shí)上，長(zhǎng)江干流的自然凈化能力可以認(rèn)為是近似均勻的，根據(jù)檢測(cè)可知，主要污染物高錳酸鹽指數(shù)和氨氮的降解系數(shù)通常介于0.1~0.5之間，比如可以考慮取0.2(單位：1/天)。附件4是“1995~2004年長(zhǎng)江流域水質(zhì)報(bào)告”給出的主要統(tǒng)計(jì)數(shù)據(jù)。下面的附表是國(guó)標(biāo)(GB3838-2002)給出的《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》中4個(gè)主要項(xiàng)目標(biāo)準(zhǔn)限值，其中Ⅰ、Ⅱ、Ⅲ類為可飲用水。請(qǐng)你們研究下列問題：（1）對(duì)長(zhǎng)江近兩年多的水質(zhì)情況做出定量的綜合評(píng)價(jià)，并分析各地區(qū)水質(zhì)的污染狀況。（2）研究、分析長(zhǎng)江干流近一年多主要污染物高錳酸鹽指數(shù)和氨氮的污染源主要在哪些地區(qū)?（3）假如不采取更有效的治理措施，依照過去10年的主要統(tǒng)計(jì)數(shù)據(jù)，對(duì)長(zhǎng)江未來水質(zhì)污染的發(fā)展趨勢(shì)做出預(yù)測(cè)分析，比如研究未來10年的情況。（4）根據(jù)你的預(yù)測(cè)分析，如果未來10年內(nèi)每年都要求長(zhǎng)江干流的Ⅳ類和Ⅴ類水的比例控制在20%以內(nèi)，且沒有劣Ⅴ類水,那么每年需要處理多少污水？（5）你對(duì)解決長(zhǎng)江水質(zhì)污染問題有什么切實(shí)可行的建議和意見。附表:《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》（GB3838—2002）中4個(gè)主要項(xiàng)目標(biāo)準(zhǔn)限值單位：mg/L序號(hào)

分類

標(biāo)準(zhǔn)值

項(xiàng)目Ⅰ類Ⅱ類Ⅲ類Ⅳ類Ⅴ類劣Ⅴ類1溶解氧(DO)

≥7.5（或飽和率90%）653202高錳酸鹽指數(shù)(CODMn)≤2461015∞3氨氮（NH3-N）

≤0.150.51.01.52.0∞4PH值（無量綱）6---9二、符號(hào)說明第個(gè)觀測(cè)站（地區(qū)）第月的綜合水質(zhì)標(biāo)識(shí)指數(shù)（=1、2、3、…、17；=1、2、3、…、28）第個(gè)觀測(cè)站（地區(qū)）第月的第種檢測(cè)項(xiàng)目的單因子水質(zhì)標(biāo)識(shí)指數(shù)（1、2、3；=1、2、3、…、17；=1、2、3、…、28）第個(gè)觀測(cè)站（地區(qū)）第月第種主要檢測(cè)項(xiàng)目指數(shù)（；=1、2、3、…、17；=1、2、3、…、28）第類第種檢測(cè)項(xiàng)目指數(shù)的上限（1、2、3；=Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、劣Ⅴ）第類第種檢測(cè)項(xiàng)目指數(shù)的下限（1、2、3；=Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、劣Ⅴ）第個(gè)觀測(cè)站（地區(qū)）28個(gè)月的綜合的評(píng)價(jià)標(biāo)識(shí)指數(shù)的平均值（=1、2、3、…、17） 降解系數(shù)（=0.1～0.5）第個(gè)觀測(cè)站（地區(qū)）第月種主要檢測(cè)項(xiàng)目的單位時(shí)間內(nèi)本地排放量（=1、2、3、…、17；2、3；=1、2、3、…、13）第月污染物在第個(gè)觀測(cè)站（地區(qū)）種主要檢測(cè)項(xiàng)目的殘留量（2、3；=1、2、3、…、17；=1、2、3、…、13）第個(gè)觀測(cè)站（地區(qū)）第月的水速（=1、2、3、…、17；=1、2、3、…、13）第個(gè)觀測(cè)站（地區(qū)）與上游觀測(cè)站（地區(qū)）的距離（=1、2、3、…、17）在第月水流從上游觀測(cè)站（地區(qū)）到第個(gè)觀測(cè)站（地區(qū)）所需要時(shí)間（=1、2、3、…、17；=1、2、3、…、13）第個(gè)觀測(cè)站（地區(qū)）第月的水流量（=1、2、3、…、17；=1、2、3、…、13）第月污染物在第個(gè)觀測(cè)站（地區(qū)）種主要檢測(cè)項(xiàng)目的13個(gè)月的平均值（2、3；=1、2、3、…、17；=1、2、3、…、13） 數(shù)據(jù)統(tǒng)計(jì)數(shù)列 數(shù)據(jù)累加數(shù)列對(duì)未來第年的Ⅳ，Ⅴ類水的所占河長(zhǎng)比例的預(yù)測(cè)值。（=1，2，…，10）劣Ⅴ類水的所占河長(zhǎng)比例的預(yù)測(cè)值。（=1，2，…，10）對(duì)未來第年總排污量值。（=1，2，…，10）個(gè)觀測(cè)站依次是：四川攀枝花，重慶朱沱

，湖北宜昌南津關(guān)，

湖南岳陽(yáng)城陵磯

，江西九江河西水廠，

安徽安慶皖河口

，江蘇南京林山

，四川樂山岷江大橋，

四川宜賓涼姜溝

，四川瀘州沱江二橋，

湖北丹江口胡家?guī)X

，湖南長(zhǎng)沙新港，

湖南岳陽(yáng)岳陽(yáng)樓

，湖北武漢宗關(guān)，

江西南昌滁槎

，江西九江蛤蟆石

，江蘇揚(yáng)州三江營(yíng)

。三、模型假設(shè)1、長(zhǎng)江干流的自然凈化能了近似均勻2、降解系數(shù)在一定時(shí)間段固定不變3、干流污染物的富集主要受上游影響四、問題說明4.1對(duì)問題一的分析理解對(duì)附件3長(zhǎng)江流域主要城市水質(zhì)檢測(cè)報(bào)告，我們對(duì)所給數(shù)據(jù)做了統(tǒng)計(jì)，有統(tǒng)計(jì)結(jié)果我們可以看出水質(zhì)分類的的依據(jù)是：?jiǎn)蝹€(gè)項(xiàng)目超過所對(duì)應(yīng)的標(biāo)準(zhǔn)限制，該項(xiàng)目所對(duì)應(yīng)的類型，即為水質(zhì)類型。我們現(xiàn)在是以單個(gè)城市為研究對(duì)象，分析水質(zhì)污染狀況，如果我們選取單個(gè)項(xiàng)目分析一個(gè)城市的污染狀況，那么四個(gè)項(xiàng)目，就要對(duì)同一個(gè)城市做四次分析，再經(jīng)過比較分析，需要得到單個(gè)城市受主要項(xiàng)目指標(biāo)的影響程度，對(duì)十七個(gè)城市，每個(gè)城市做四次分析，我們可以看到，工作量是多么巨大。并且如此得到的結(jié)果是：?jiǎn)蝹€(gè)城市近兩年多水環(huán)境質(zhì)量標(biāo)準(zhǔn)類型主要是受那個(gè)污染物的影響。由于原問題要求給出定量的綜合評(píng)定，我們?nèi)绻蒙鲜龇椒?，那么必然不符合定量分析的要求，所以我們?cè)O(shè)出了一個(gè)綜合評(píng)價(jià)的指數(shù)，于是設(shè)該指數(shù)為：綜合水質(zhì)標(biāo)識(shí)指數(shù)。在設(shè)出綜合水質(zhì)標(biāo)識(shí)指數(shù)后，我們按月份計(jì)算出每個(gè)月的綜合水質(zhì)標(biāo)識(shí)指數(shù)的數(shù)值，這樣可以以月份為橫軸，綜合水質(zhì)標(biāo)識(shí)指數(shù)為縱軸，畫出單個(gè)的城市隨月份不同，綜合水質(zhì)標(biāo)識(shí)指數(shù)變化的曲線。用這種方法，我們僅用十七個(gè)圖就綜合展現(xiàn)了各個(gè)城市水質(zhì)量污染狀況的變化趨勢(shì)，這樣更加直觀簡(jiǎn)潔，同時(shí)我們使用綜合水質(zhì)標(biāo)識(shí)指數(shù)，也確實(shí)做到了定量分析，綜合來說這種方法是很具有說服力的。4．2對(duì)問題二的分析理解對(duì)長(zhǎng)江干流近一年多主要污染物高錳酸鹽指數(shù)和氨氮的主要污染源的判斷，我們必須考慮到上游對(duì)下游的的影響。上游污染物通過水流可實(shí)現(xiàn)對(duì)下游的富集，故必須通過水流量，水流速，站點(diǎn)距離對(duì)污染物富集的影響。如果我們單純的只確定污染源，根據(jù)主要檢測(cè)項(xiàng)目值就可簡(jiǎn)單判斷，由于污染物從上游到下游，要經(jīng)過降解過程，如果下游污染物的檢測(cè)項(xiàng)目值大于上游的值，那么就可簡(jiǎn)單判斷此地是污染源。但根據(jù)題意我們必須判斷的是主要污染源，由此不能把污染源當(dāng)作主要污染源，我們必須比較的是當(dāng)?shù)氐呐盼勰芰?，才能斷定那一地點(diǎn)是主要污染源。綜合上述針對(duì)問題三，我們主要的目的就是建立模型計(jì)算出單個(gè)地點(diǎn)的排污能力，再綜合比較，得出主要污染源，另外由于降解系數(shù)（0.1—0.5）不同，我們可以得到不同的污染物富集結(jié)果，由此也得到不同的單個(gè)地點(diǎn)的排污能力。通過這種方法，我們就可方便解決判斷主要污染源的地點(diǎn)問題。4．3對(duì)問題三的分析理解首先我們對(duì)評(píng)價(jià)河長(zhǎng)，這一概念進(jìn)行了分析，通過概率論的有關(guān)知識(shí)，我們可以對(duì)評(píng)價(jià)河長(zhǎng)，理解為隨機(jī)取樣，有部分研究全體的問題。由此有評(píng)價(jià)河長(zhǎng)，得出的各類水質(zhì)所占河長(zhǎng)的百分比，可認(rèn)為是全流域的該類水質(zhì)所占河長(zhǎng)的比例。這樣就達(dá)到了對(duì)評(píng)價(jià)河長(zhǎng)的研究，實(shí)現(xiàn)了對(duì)全流域的研究。由于該問題是通過過去十年的數(shù)據(jù)實(shí)現(xiàn)對(duì)未來十年的預(yù)測(cè)，所以有必要建立一個(gè)各類水質(zhì)類型河長(zhǎng)所占全流域比例，隨年度變化的函數(shù)模型。但通過對(duì)各類水質(zhì)類型河長(zhǎng)所占全流域比例，隨年度變化的圖形，我們看到它們的在坐標(biāo)圖（以年度為橫軸，比例為縱軸）上的點(diǎn)是很離散的，我們很難得到殘差比較小，并且擬合程度很好的函數(shù)。所以我們又考慮到對(duì)后三類水質(zhì)類型河長(zhǎng)所占全流域比例之和，與年度變化建立函數(shù)模型（該模型如果建立，也可對(duì)污染程度預(yù)測(cè)，因?yàn)楹笕悶榉秋嬘盟侵匚廴舅颍?，但結(jié)果同上，依然得到很離散的圖形（該圖見圖（一）），我們同樣通過線形回歸很難得到擬合很好的函數(shù)。這時(shí)我們圖（一）圖（二）又考慮用設(shè)一個(gè)水質(zhì)年度級(jí)別來判斷，這樣每年對(duì)應(yīng)一個(gè)年度級(jí)別，綜合預(yù)測(cè)未來十年的水質(zhì)級(jí)別，但這個(gè)級(jí)別很難用一個(gè)有說服力的標(biāo)準(zhǔn)來確定，另外問題四，也要用到未來十年的Ⅳ類，Ⅴ類，劣Ⅴ水的比例數(shù)據(jù)，故設(shè)出水質(zhì)年度級(jí)別，并不方便解題，并無端增加數(shù)據(jù)統(tǒng)計(jì)上的麻煩。另外我們?cè)噲D建立各個(gè)類水的比例與排污量的關(guān)系圖（以單個(gè)類的水十年所占百分比為橫軸上的點(diǎn)（單位為：%），年排污量為縱軸上的點(diǎn)（單位為：億噸）），我們依然看到它們的圖形是很離散的（圖二為Ⅰ類水的圖形）。我們不優(yōu)先考慮此類模型。最后我們不得不考慮最初所設(shè)想的模型的可行性，建立每年排污量與年度的對(duì)應(yīng)關(guān)系函數(shù)時(shí)，我們選擇了建立了灰色系統(tǒng)模型，這樣比較精確的對(duì)問題求解。同樣年度排污量與年度的變化關(guān)系也有灰色系統(tǒng)模型求出。4．4對(duì)問題四的分析理解只要得出Ⅳ，Ⅴ類水的所占河長(zhǎng)比例的預(yù)測(cè)值、劣Ⅴ類水的所占河長(zhǎng)比例的預(yù)測(cè)值和未來第年總排污量值就可以簡(jiǎn)單的算出排污量。我們便可以得到所需模型。五、模型的建立與求解5．1根據(jù)所給兩年多主要檢測(cè)項(xiàng)目的數(shù)據(jù)，對(duì)長(zhǎng)江兩年多的水質(zhì)情況做出定量的綜合評(píng)定，并分析各個(gè)地區(qū)水質(zhì)的污染狀況。依照《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》（GB3838-2002）中的四個(gè)主要項(xiàng)目標(biāo)準(zhǔn)限值，我們根據(jù)兩年多的水質(zhì)情況表，我們判斷水質(zhì)類別是有單個(gè)項(xiàng)目超標(biāo)決定的，但為了綜合評(píng)定，我們有必要把各個(gè)項(xiàng)目都考慮在內(nèi)，由此得到綜合的評(píng)價(jià)標(biāo)識(shí)指數(shù)。我們首先建立單因子水質(zhì)標(biāo)識(shí)指數(shù)函數(shù)模型。我們可將模型分為兩類，這是因?yàn)槿萁庋?DO)的單因子水質(zhì)標(biāo)識(shí)指數(shù)判斷是有該類型的下限值比較所決定的，而高錳酸鹽指數(shù)（KMn）與氨氮（NH3-N）的單因子水質(zhì)標(biāo)識(shí)指數(shù)是有該類型的上限值比較所確定的。由于PH值在判斷水質(zhì)類別時(shí)，沒有明顯作用，我們?cè)诳紤]綜合的評(píng)價(jià)標(biāo)識(shí)指數(shù)時(shí)不在考慮它對(duì)該指數(shù)的影響。最后我們可以對(duì)單個(gè)城市28個(gè)月的綜合的評(píng)價(jià)標(biāo)識(shí)指數(shù)求平均值，作為最后水質(zhì)污染狀況的綜合評(píng)價(jià)依據(jù)，該數(shù)據(jù)越大表示污染程度越深，并設(shè)該平均值為：在這里我們求出第個(gè)觀測(cè)站（地區(qū)）第月容解氧(DO)的單因子水質(zhì)標(biāo)識(shí)指數(shù)，我們?yōu)榱伺c水質(zhì)類別建立聯(lián)系，故得到如下的計(jì)算公式：+，（=1）上式中當(dāng)=Ⅰ時(shí)=1，當(dāng)=Ⅱ時(shí)=2，當(dāng)=Ⅲ時(shí)=3，當(dāng)=Ⅳ時(shí)=4，當(dāng)=Ⅴ時(shí)=5，當(dāng)=劣Ⅴ時(shí)，=6。有上式我們先判斷種（=1表示容解氧(DO)，=2表示高錳酸鹽指數(shù)，=3表示氨氮）所屬類型對(duì)應(yīng)值的和，再進(jìn)行計(jì)算。其中當(dāng)=1，=Ⅰ時(shí)，根據(jù)附錄3給出的過去兩年多的數(shù)據(jù)，經(jīng)統(tǒng)計(jì)得出=14.4，=7.5。對(duì)于當(dāng)=1，=Ⅱ、Ⅲ、Ⅳ、Ⅴ、劣Ⅴ時(shí)，參照《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》上的數(shù)據(jù)便可得到和的值。用同樣方法我們可求出第個(gè)觀測(cè)站（地區(qū)）第月高錳酸鹽指數(shù)（KMn）與氨氮（NH3-N）的單因子水質(zhì)標(biāo)識(shí)指數(shù)，故得到如下的計(jì)算公式：+，（=2、3）上式中同樣當(dāng)=Ⅰ時(shí)=1，當(dāng)=Ⅱ時(shí)=2，當(dāng)=Ⅲ時(shí)=3，當(dāng)=Ⅳ時(shí)=4，當(dāng)=Ⅴ時(shí)=5，當(dāng)=劣Ⅴ時(shí)，=6。有上式我們先判斷種（=1表示容解氧(DO)，=2表示高錳酸鹽指數(shù)，=3表示氨氮）所屬類型對(duì)應(yīng)值的和，當(dāng)=2，=Ⅳ時(shí)，=10，=6；當(dāng)=2，=Ⅴ、劣Ⅴ時(shí)，由于根據(jù)附錄3給出的過去兩年多的數(shù)據(jù)，經(jīng)統(tǒng)計(jì)得出高錳酸鹽指數(shù)的上限值為9.9，所以在這里和不需要再給出。當(dāng)=2，=Ⅰ、Ⅱ、Ⅲ時(shí)參照《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》上的數(shù)據(jù)便可得到和的值。當(dāng)=3，=劣Ⅴ時(shí)，根據(jù)附錄3給出的過去兩年多的數(shù)據(jù)，經(jīng)統(tǒng)計(jì)得出氨氮的=24.2，=15。當(dāng)=2，=Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ時(shí)，參照《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》上的數(shù)據(jù)便可得到和的值。有以上單因子水質(zhì)標(biāo)識(shí)指數(shù)確定后，我們可得到綜合的評(píng)價(jià)標(biāo)識(shí)指數(shù)，并有計(jì)算公式可求出，公式如下：有以上公式我們用matlab編程（見附錄一）分別得出如下各個(gè)觀測(cè)站（地區(qū)）第28個(gè)月的綜合的評(píng)價(jià)標(biāo)識(shí)指數(shù)，我們列出部分城市數(shù)據(jù)及根據(jù)月份，單因子水質(zhì)標(biāo)識(shí)指數(shù)和綜合的評(píng)價(jià)標(biāo)識(shí)指數(shù)變化的曲線圖（其余數(shù)據(jù)見附錄一），如下：四川攀枝花：（1.7444 2.4599 1.7210 1.6983 2.1492 1.6394 1.8143 2.1620 2.4505 1.9345 1.9756 2.1287 1.7766 1.7758 3.0108 2.4489 1.5372 1.8444 1.6425 1.5585 1.7258 1.8533 1.6650 1.8630 1.9316 2.1237 1.8976 2.1295）重慶朱沱：（2.2703 1.9783 1.9636 1.9502 2.0670 2.0734 1.9670 2.0600 2.1903 2.2765 2.3516 2.1346 2.1778 2.2346 2.0070 2.3651 1.9288 1.9988 2.0832 2.0694 2.1787 2.2032 2.3673 2.2367 1.8434 2.1134 2.1309 2.0273） 根據(jù)matlab編程，繪圖所得結(jié)果，我們通過對(duì)單個(gè)城市28個(gè)月的綜合的評(píng)價(jià)標(biāo)識(shí)指數(shù)求平均值，公式如下：我們分析各個(gè)地區(qū)的水質(zhì)污染狀況如下：個(gè)地區(qū)的28個(gè)月的綜合的評(píng)價(jià)標(biāo)識(shí)指數(shù)求平均值依次如下：（1.95222.1162.23012.41842.10192.25152.04483.54692.25092.75411.78032.8682.56282.3923.58882.44352.3802）。越大表示該地區(qū)污染越嚴(yán)重，有該數(shù)據(jù)我們判斷四川樂山岷江大橋，江西南昌滁槎兩地污染最嚴(yán)重，重慶朱沱

、湖北宜昌南津關(guān)、湖南岳陽(yáng)城陵磯

、江西九江河西水廠、

安徽安慶皖河口

、江蘇南京林山、四川宜賓涼姜溝、四川瀘州沱江二橋、西九江蛤蟆石

、江蘇揚(yáng)州三江營(yíng)

、湖南長(zhǎng)沙新港、湖南岳陽(yáng)岳陽(yáng)樓、湖北武漢宗關(guān)屬于中度污染，其余兩地污染不是很嚴(yán)重。5.2根據(jù)近一年多的干流統(tǒng)計(jì)數(shù)據(jù)，研究、分析主要污染物高錳酸鹽指數(shù)（KMn）與氨氮（NH3-N）的主要污染源所在地。由于水從上游到下游流動(dòng)的過程中，必然會(huì)產(chǎn)生污染物的富集。又由于水系自身有自凈化能力，我們可以通過上下游主要監(jiān)測(cè)項(xiàng)目的指標(biāo)值的對(duì)比，簡(jiǎn)單判斷某一個(gè)地區(qū)是否存在污染源。但我們的目的是為了得到主要污染源，所以我們必須建立模型，通過對(duì)比各地污染物的排放量，確定污染源所在的主要地區(qū)。我們首先確定污染物第月水從上游觀測(cè)站（地區(qū)）流到第個(gè)觀測(cè)站（地區(qū)）經(jīng)過時(shí)間后，第月污染物在第個(gè)觀測(cè)站（地區(qū)）種主要檢測(cè)項(xiàng)目的殘留量（=2表示高錳酸鹽指數(shù)，=3表示氨氮）的計(jì)算公式，計(jì)算過程如下：由于（為上游觀測(cè)站（地區(qū)）第月種主要檢測(cè)項(xiàng)目的單位時(shí)間內(nèi)本地排放量，為降解系數(shù)）于是可以得到在在時(shí)間內(nèi)第個(gè)觀測(cè)站（地區(qū)）的上游污染物的降解量，再有上游觀測(cè)站（地區(qū)）第月種主要檢測(cè)項(xiàng)目的本地排放量減去降解量，便可得到殘留量的計(jì)算公式：。有以上公式，我們便可求出第個(gè)觀測(cè)站（地區(qū)）第月種主要檢測(cè)項(xiàng)目的本地排單位時(shí)間內(nèi)放量（=2表示高錳酸鹽指數(shù)，=3表示氨氮）的計(jì)算公式：（表示第個(gè)觀測(cè)站（地區(qū)）第月的水流量）為了方便分析，我們對(duì)十三個(gè)月的求和取平均值，來確定污染源的主要所在地。經(jīng)過matlab編程（見附錄二），我們可以計(jì)算出當(dāng)=0.2（=0.1～0.5）時(shí)的的數(shù)值為：當(dāng)=2時(shí)（=2表示高錳酸鹽指數(shù)）求（單位為：噸）用數(shù)列表示為：（8.986，37.1748，，50.907，70.4526，58.196，59.9114，58.259）得到長(zhǎng)江七個(gè)干流城市排污量對(duì)比圖如下：由此我們得出高錳酸鹽的主要污染源所在地為=4，5，6,7三地。即湖南岳陽(yáng)城陵磯

，江西九江河西水廠，

安徽安慶皖河口,江蘇南京林山三地。當(dāng)=3時(shí)（=3表示氨氮）求（單位為：噸）用數(shù)列表示為：（0.4816，3.0496，4.1418，6.3864，5.0473，5.0276，2.4794）得到長(zhǎng)江七個(gè)干流城市排污量對(duì)比圖如下：由此我們得出氨氮污染源的主要所在地為=4，5，6三地。即湖南岳陽(yáng)城陵磯

，江西九江河西水廠，

安徽安慶皖河口三地。另外根據(jù)取值不同，所得每個(gè)城市各個(gè)月的的大小也會(huì)不同，就可以斷定出污染源的主要所在地也不同，我們列表表示如下（具體計(jì)算數(shù)據(jù)見附錄）：當(dāng)=0.1=2時(shí)有以上圖形我們可以看出降解系數(shù)對(duì)污染源斷定影響并不是很大，經(jīng)過分析我們可以斷定污染源的主要所在地排量比積累量大的多，更加說明了模型的正確性。5.3在不采取有效治理的措施的情況下，依照過去是十年的主要統(tǒng)計(jì)數(shù)據(jù)，對(duì)長(zhǎng)江未來水質(zhì)污染發(fā)展趨勢(shì)作出預(yù)測(cè)分析，比如研究未來十年的情況。通過對(duì)《附件4：1995年～2004年長(zhǎng)江流域水質(zhì)報(bào)告》的分析，我們決定對(duì)長(zhǎng)江未來水質(zhì)污染做兩個(gè)方面的預(yù)測(cè)，這兩個(gè)方面是：一、對(duì)未來十年每年排污量進(jìn)行預(yù)測(cè)。二、對(duì)未來十年Ⅳ、Ⅴ、劣Ⅴ類水所占全流域比例作出預(yù)測(cè)。在對(duì)為來十年的排污量進(jìn)行預(yù)測(cè)時(shí)，我們建立了灰色系統(tǒng)模型。當(dāng)系統(tǒng)內(nèi)部信息和特性是部分已知的，另一部分是未知時(shí)，人們往往難以建立客觀的物理原形，內(nèi)部因素難以辨識(shí)，以及相互之間的關(guān)系較為隱蔽，難已準(zhǔn)確了解這類系統(tǒng)的行為特征，因此，對(duì)于這類問題進(jìn)行定量描述，即建立模型難度較大，我們便選擇建立灰色系統(tǒng)(GreySystem)，對(duì)問題求解?；疑到y(tǒng)分析建模方法主要是根據(jù)具體灰色系統(tǒng)的行為特征數(shù)據(jù)，充分利用數(shù)量不多的數(shù)據(jù)和信息尋求相關(guān)各因素自身和因素之間的數(shù)學(xué)關(guān)系，我們對(duì)第一方面排污量進(jìn)行預(yù)測(cè)，建立相應(yīng)的數(shù)學(xué)模型模型建立過程如下：由已知數(shù)據(jù)，對(duì)于1995年～2004年每年廢水排污量統(tǒng)計(jì)為以下數(shù)據(jù)：對(duì)作一次累加，則：，；于是得到類加數(shù)列為：取的加權(quán)均值，則，為確定參數(shù)，設(shè)，記（5.3--1）于是GM(1,1)的白化微分方程模型為（5.3--2）其中是發(fā)展灰度，是內(nèi)生控制灰度.由于，取為灰導(dǎo)數(shù)，為背景值，則將方程對(duì)（5.3--2）相應(yīng)的灰微分方程為或即矩陣形式為其中,.用最小二乘法求得參數(shù)的估計(jì)值為于是方程（5.3--2）有響應(yīng)（特解）則（5.3—3）對(duì)未來第年的Ⅳ，Ⅴ類水的所占河長(zhǎng)比例的，我們建立了百分比y與廢水派放量x之間的關(guān)系y=f(x)我們可有公式（5.3—3）求出1995～2004十年的計(jì)算數(shù)值，經(jīng)過殘差檢驗(yàn)分析（數(shù)據(jù)見附錄），我們斷定所建模型是比較合適用于解決問題的。利用matlab編程（程序見附錄三）我們得出對(duì)未來十年預(yù)測(cè)的求解數(shù)據(jù)如下（322．5221，343．2881，365.3912，388.9175，413.9585，440.6118，468.9812，303．123，499.1772，531.3174）我們對(duì)殘差分析如下圖：對(duì)未來第年的Ⅳ，Ⅴ類水的所占河長(zhǎng)比例的，我們建立了百分比y與廢水派放量x之間的關(guān)系y=f(x)由于劣V類的趨勢(shì)不明顯，因而直接利用1減去第一大類與第二大類的比例即得劣V類的比例。IV、V類水比例之和枯水期干流Y=10.64.log(x-170)-21.02.支流y=-0.003257log(x-170)+1767.豐水期干流Y=6.076log(x-170)-3.07支流y=4.1178log(x-170)+1.31水文年干流y=8.4197log(x-170)-13.03.支流Y=0.01177log(x-170)+17.22關(guān)于對(duì)未來第年的Ⅳ，Ⅴ類水的所占河長(zhǎng)比例的預(yù)測(cè)值與對(duì)未來第年總排污量值的關(guān)系函數(shù)模型建立也使用該灰色系統(tǒng)模型。5.4根據(jù)解題過程5.3中，對(duì)未來水質(zhì)污染發(fā)展趨勢(shì)做出的預(yù)測(cè)，做出未來十年內(nèi)每年長(zhǎng)江干流的Ⅳ類和Ⅴ類水的比例控制在20%，并且沒有劣Ⅴ類水的污水處理模型。有5.3中所求出數(shù)據(jù)，即Ⅳ，Ⅴ類水的所占河長(zhǎng)比例的預(yù)測(cè)值對(duì)未來第年的Ⅳ，Ⅴ類水的所占河長(zhǎng)比例的，我們建立了百分比y與廢水派放量x之間的關(guān)系y=f(x),令y≤20,求出x的上限，則預(yù)測(cè)的廢水排放量與x的上限的差值即為需要處理的污水量，便得到了年處理量。未來十年每年廢水排放量和需處理量表注：?jiǎn)挝唬簝|噸年份2005200620072008200920102011201220132014排放量308.9331.8356.3382.5410.4440.0471.3504.2538.7575.0處理量98.9121.8146.3172.5200.4230.0261.3294.2328.7365.05.5根據(jù)求解過程5.1和5.3所得結(jié)果和對(duì)長(zhǎng)江污染問題的了解，提出切實(shí)可行的建議和意見。近些年來水質(zhì)污染狀況越來約嚴(yán)重，我們賴以生存的水的現(xiàn)狀是———淮河濁臭、黃河斷流、長(zhǎng)江告急。水是萬物之靈，生命之源，如若對(duì)水質(zhì)污染置之不理，任其發(fā)展，那么子孫后代將會(huì)喝不到放心的飲用水，我們這代人將有愧與我們的民族！通過我們自己對(duì)長(zhǎng)江水系污染狀況的研究，我們深深感受到拯救長(zhǎng)江的時(shí)代責(zé)任感和緊迫感！我們深信專家所言：長(zhǎng)江已陷入深度危機(jī)，若不及時(shí)拯救，10年之內(nèi)，長(zhǎng)江水系生態(tài)將瀕臨崩潰。所以拯救長(zhǎng)江必須行動(dòng)了，我們?cè)谙挛膶?duì)解決長(zhǎng)江污染問題，提出了自己的看法和意見。第一，長(zhǎng)江沿岸地區(qū)要改變GDP增長(zhǎng)模式，提升對(duì)招商引資的技術(shù)要求。一些地區(qū)只重視GDP的增長(zhǎng)，對(duì)于招商引資項(xiàng)目缺乏管理。如江蘇某市招商引資來的都是鋼鐵廠、化工廠、造紙廠、造船廠、拆船廠這樣的重污染企業(yè)，沿幾十公里長(zhǎng)的江岸一字排開。在我國(guó)，消耗了世界1/4鋼鐵，1/3煤炭，1/2水泥才造就出世界1/30GDP。如此資源成本造成長(zhǎng)江上游嚴(yán)重水土流失，中游生態(tài)脆弱，下游污染驚人。而管理者更是缺乏對(duì)環(huán)境保護(hù)的關(guān)注，甚至有人認(rèn)為：沿江建廠，長(zhǎng)江運(yùn)輸成本低、排污方便，直排長(zhǎng)江，直沖就東海。還有寫地方重復(fù)建設(shè)嚴(yán)重，對(duì)經(jīng)濟(jì)只重視量的增長(zhǎng)，而不重視質(zhì)的提高！我國(guó)提倡大力發(fā)展經(jīng)濟(jì)，重視招商引資，但決不是不分層次，只重視量，不重視質(zhì)。所以要改變經(jīng)濟(jì)增長(zhǎng)方式，改資源密集型，勞動(dòng)密集型的企業(yè)增長(zhǎng)方式為技術(shù)密集型的企業(yè)增長(zhǎng)方式，扭轉(zhuǎn)高消耗、高污染、低產(chǎn)出的狀況，全面轉(zhuǎn)變經(jīng)濟(jì)增長(zhǎng)方式。同時(shí)對(duì)招商引資進(jìn)行管理，利用外資開發(fā)本地資源的同時(shí)，也要根據(jù)環(huán)保要求限制部分企業(yè)在沿江流域的建設(shè)。GDP增長(zhǎng)不應(yīng)以生態(tài)赤字為交換條件，否則，不是福音而是災(zāi)難。第二，完善法律，完善管理方式和制度，加強(qiáng)對(duì)長(zhǎng)江水系的管理。例如：對(duì)于POP′S，我國(guó)因無系統(tǒng)調(diào)研，也就無POP′S指標(biāo)，所以排污就無人管理。所以我國(guó)應(yīng)加快完善檢查指標(biāo)，另外一些地區(qū)的干部缺乏管理意識(shí)，忽視環(huán)保，我們應(yīng)該加強(qiáng)對(duì)對(duì)干部的環(huán)抱教育。增強(qiáng)干部環(huán)保觀念，并把環(huán)保作為考核干部指標(biāo)之一。另外就如國(guó)家防汛指揮部一樣，我們應(yīng)該成立一個(gè)同意的監(jiān)管部門，統(tǒng)一管理長(zhǎng)江水系。法不在立，而在必定，所以必須完善環(huán)保法律，絕不讓不法分子專法律的漏洞。第三，增強(qiáng)輿論向?qū)?，提高全民環(huán)保意識(shí)。很多地區(qū)除了領(lǐng)導(dǎo)干部環(huán)保意識(shí)差外，廣大群眾也是極其缺乏環(huán)保意識(shí)。嚴(yán)重污染成黑色的小河邊，有的農(nóng)民在那里釣魚。臟水里面都是劇毒物質(zhì)，魚聞起來有很重的腐臭味道，怎么能吃呢！另外根據(jù)有關(guān)資料，沿江的生活污水很多就直接排入長(zhǎng)江，長(zhǎng)江變成了隨意堆放、傾倒垃圾的“垃圾河”，進(jìn)入三峽后，垃圾帶白茫茫一片覆蓋整個(gè)江面，生活與工業(yè)垃圾、油類液體、死牛、死羊都有。如果能夠喚起群眾意識(shí)，那么必然有一個(gè)異常強(qiáng)大的環(huán)保聯(lián)盟出現(xiàn)。第四，保護(hù)長(zhǎng)江源。嘉陵江的水還有點(diǎn)藍(lán)，長(zhǎng)江的水就是蠟黃蠟黃的，像黃河一樣，含有大量泥沙。在重慶市，兩江交匯的地方形成涇渭分明的一道奇景：一藍(lán)一黃，分界線特別明顯，這說明長(zhǎng)江上游的水土流失越來越嚴(yán)重。我們應(yīng)當(dāng)采取有效措施保護(hù)長(zhǎng)江源。第五，加大對(duì)治污企業(yè)的扶持力度。有關(guān)數(shù)據(jù)顯示：沿江351座污水處理廠，絕大多數(shù)因缺乏資金或網(wǎng)管滯后而閑置。相當(dāng)數(shù)量污水處理廠心理趨同“我花很多錢治污，對(duì)岸的、隔壁的都不治污，我不污染，白不污染?！蹦呈幸幌冗M(jìn)的污水處理廠，該廠實(shí)話實(shí)說，僅治理5%，就花了4.4億元，95%便作為一般性處理了。先進(jìn)處理，1千噸需300元，一般性處理60元。所以政府部門要出臺(tái)相應(yīng)政策，引導(dǎo)該類企業(yè)發(fā)展，同時(shí)給予優(yōu)惠發(fā)展或優(yōu)先發(fā)展的措施。另外銀行部門向治理污水的企業(yè)提供優(yōu)惠的貸款政策，提高企業(yè)治理污水的積極性。第六，嚴(yán)調(diào)狠打，加重對(duì)排污企業(yè)的處罰力度?！缎旅裰芸罚簩?duì)污染排放超標(biāo)企業(yè)最多處罰10萬的力度是否過低？有這樣一個(gè)案例：有個(gè)工廠污染排放超標(biāo)，環(huán)保局為了處罰它必須取樣、化驗(yàn)……一系列取證工作就花了50萬。投入如此巨大，罰款卻寥寥無幾，環(huán)保局哪里來積極性？排污企業(yè)竊喜：有本事你就罰我！違法成本低、守法執(zhí)法成本高。致使環(huán)保局有法不依、執(zhí)法不嚴(yán)、違法不究。《水環(huán)境保護(hù)法》嚴(yán)重滯后，最多罰10萬元額度，不能有效遏制違法企業(yè)偷排、漏排、超標(biāo)排放。利益驅(qū)動(dòng)，其瞞和騙到了無顧老百姓死活程度。所以應(yīng)當(dāng)增加罰款額，對(duì)違法企業(yè)重罰！第七，因地制宜，重點(diǎn)城市重點(diǎn)治理。我們分析長(zhǎng)江近些年的水質(zhì)情況，可以知道，不同地區(qū)，水質(zhì)的污染程度不同，主要受的污染物也不同。長(zhǎng)江岸邊形成了許多污染帶，在干流21個(gè)城市中，重慶、岳陽(yáng)、武漢、南京、鎮(zhèn)江、上海六大城市累計(jì)污染帶長(zhǎng)度占長(zhǎng)江干流污染帶總長(zhǎng)的73%。我們要根據(jù)不同城市的狀況出臺(tái)有針對(duì)性的治理方案。俗話：一方水土養(yǎng)一方人。眼下長(zhǎng)江陷入一方人害一方水土，一方水土害一方人的惡性循環(huán)?！八茌d舟，亦能覆舟”，千古訓(xùn)誡?？茖W(xué)發(fā)展觀揚(yáng)載抑覆。否則，還要吞噬多少人健康與生命？我們用以上提議，七劍合壁，實(shí)現(xiàn)教育治污、法律治污、科技治污，投入治污。這樣就能有效的保護(hù)長(zhǎng)江，不讓環(huán)境污染危害人類的悲劇上演！六、誤差分析綜合水質(zhì)標(biāo)識(shí)指數(shù)較為完整的標(biāo)識(shí)了總體水質(zhì)類別，水質(zhì)情況，它急能定量評(píng)價(jià)，也能定性評(píng)價(jià)；即不會(huì)因?yàn)閭€(gè)別水質(zhì)指標(biāo)較差就否定整體水質(zhì)，又能對(duì)河流整體水質(zhì)作出合理評(píng)價(jià)；既可以在同一類別中比較水質(zhì)優(yōu)劣，也可以比較不同劣Ⅴ類水的污染嚴(yán)重程度。但是單個(gè)項(xiàng)目標(biāo)準(zhǔn)限值決定水質(zhì)分類，可以有效的解決水分類的飲用問題。綜合評(píng)價(jià)指數(shù)，對(duì)于Ⅱ、Ⅲ、Ⅳ類的計(jì)算數(shù)據(jù)上很難斷定是否為飲用水。通過對(duì)主要檢測(cè)項(xiàng)目的本地單位時(shí)間內(nèi)排放量十三個(gè)月的數(shù)值求和取平均做對(duì)比，我們分析污染源的主要所在地，此方法比較綜合的評(píng)價(jià)了一年多的某地區(qū)的排污情況，經(jīng)過對(duì)比很有說服力，但是這種方法缺少了對(duì)單個(gè)月份的變動(dòng)情況考慮，對(duì)污染源的主要所在地有可能在單個(gè)月改變?nèi)鄙兕A(yù)測(cè)?；疑到y(tǒng)模型，對(duì)于不太離散的問題情況分析較為精確，但對(duì)于離散的問題，使用該模型時(shí)，就會(huì)出現(xiàn)失真的情況。 七、模型改進(jìn)與推廣對(duì)于綜合的評(píng)價(jià)標(biāo)識(shí)指數(shù)和單因子水質(zhì)標(biāo)識(shí)指數(shù)，我們可以看到它的綜合評(píng)價(jià)作用是很優(yōu)越的，我們可以把這兩類指數(shù)推廣到其他污染評(píng)價(jià)上，甚至可以對(duì)用一定方法對(duì)作出經(jīng)濟(jì)發(fā)展指數(shù)，用于評(píng)價(jià)經(jīng)濟(jì)狀況。該指數(shù)計(jì)算方法經(jīng)過簡(jiǎn)單的變換，我們便很容易的得到社會(huì)很多方面的評(píng)價(jià)指標(biāo)?；疑到y(tǒng)分析建模方法主要是根據(jù)具體灰色系統(tǒng)的行為特征數(shù)據(jù)，充分利用數(shù)量不多的數(shù)據(jù)和信息尋求相關(guān)各因素自身和因素之間的數(shù)學(xué)關(guān)系，即建立相應(yīng)的數(shù)學(xué)模型.灰色系統(tǒng)可以推廣到工程技術(shù)、經(jīng)濟(jì)管理、氣象預(yù)報(bào)，以及政治、社會(huì)、工業(yè)、農(nóng)業(yè)等領(lǐng)域統(tǒng)計(jì)預(yù)測(cè)的方面，對(duì)于不是很離散的情況，該模型求解很有說服力。對(duì)于根據(jù)長(zhǎng)江水系的污染狀況，我們提出的治污意見，對(duì)其他水系同樣適用。八、參考文獻(xiàn)[1]何壽強(qiáng)周萍張保平，嘉陵江廣元段水質(zhì)狀況分析與評(píng)價(jià)，《四川環(huán)境》，第14卷第3期：57～61頁(yè)，1995年。[2]幸治國(guó)鐘成華，長(zhǎng)江、嘉陵江重慶干流江段－維水質(zhì)數(shù)學(xué)模型，《重慶環(huán)境科學(xué)》，第16卷第3期：19～23頁(yè)，1994年。[3]黃樹聲高冬蘭，南流江水資源保護(hù)問題不容忽視，《廣西水利水電》，第四期：24～28頁(yè)，1995年。附錄一functionzohn2X=[174,179,183,189,207,234,220.5,256,270,285];x=[174,179,183,189,207,234,220.5,256,270,285];Y=[174,179,183,189,207,234,220.5,256,270,285];y=[];Z=[];z=[];r=[];fori=2:10X(i)=X(i)+X(i-1);endforj=2:10Z(j)=0.5*X(j)+0.5*X(j-1);endZ(1)=[];Y(1)=[];B=[-Z;ones(1,9)]';C=inv(B'*B)*B'*Y';a=C(1);b=C(2);fork=1:9y(k+1)=(174-b/a)*exp(-a*k)+b/a;endy(1)=174;fort=1:9z(t+1)=y(t+1)-y(t);endz(1)=174;fors=1:10r(s)=(x(s)-z(s))/x(s);endm=[1995199619971998199920002001200220032004];errorbar(m,z,r);title('模型誤差棒棒圖');z以上是對(duì)單個(gè)城市一個(gè)月的求解程序，如果城市改變，僅需要改變輸入的X，x，Y值即可。附錄二

functionzohn2X=[174,179,183,189,207,234,220.5,256,270,285];x=[174,179,183,189,207,234,220.5,256,270,285];Y=[174,179,183,189,207,234,220.5,256,270,285];y=[];Z=[];z=[];r=[];fori=2:10X(i)=X(i)+X(i-1);endforj=2:10Z(j)=0.5*X(j)+0.5*X(j-1);endZ(1)=[];Y(1)=[];B=[-Z;ones(1,9)]';C=inv(B'*B)*B'*Y';a=C(1);b=C(2);fork=1:19y(k+1)=(174-b/a)*exp(-a*k)+b/a;endy(1)=174;fort=1:19z(t+1)=y(t+1)-y(t);endz(1)=174;s=[19951996199719981999200020012002200320042005200620072008200920102011201220132014];plot(s,z,'r');z以上是對(duì)單個(gè)城市一個(gè)月的求解程序，如果城市改變，僅需要改變輸入的X，x，Y值即可。附錄三functionzohn33A=[5.184.855.354.046