【正文】 995年~2004年長江流域水質(zhì)報告中水文年干支流總的河長綜合于如下表七：表七：近十年水文年干支流總的河長時間類型1995199619971998199920002001200220032004水文年干流4456447944794529613652856012598362266341支流5469555855136429243302502728134294033228733071由表七中數(shù)據(jù)可知，干支流河長呈二次曲線上升趨勢，故可對表中數(shù)據(jù)進行二次多項式分析模擬，如圖四和圖五（程序見附錄中程序三）：圖四：未來十年水文年干流河長預測圖圖五：未來十年水文年支流河長預測圖得到中水文年干支流總的河長擬合二次多項式為：水文年干流為：水文年支流為：由于報告中對水的分類較為細致（水質(zhì)被分成了六類），而對于這樣一個細分的分類，十年的數(shù)據(jù)是很難得出一個較為明顯的規(guī)律的，從我們的初期的數(shù)據(jù)分析就很容易得到這一點結論。對污水的處理是為了滿足既定目標：未來10年內(nèi)每年都要求長江干流的Ⅳ類和Ⅴ類水的比例控制在20%以內(nèi)，且沒有劣Ⅴ類水。根據(jù)我們的數(shù)據(jù)預測如果未來十年之內(nèi)不采取保護長江的措施，十年之后長江的水域有近70%會屬于IV、V 和劣V 類，因此要采取必要的措施來防止長江遭到進一步的惡化。(4)相應措施 ①、長江流域管理機構應會同地方環(huán)保部門，根據(jù)長江水量，提出流域內(nèi)每個省（區(qū)）許的最大排污量，實行排污總量控制。應用數(shù)值方法結合蒙特卡洛法解決在矩陣處理的總體思路下穿插著非線性規(guī)劃的這么一個較為復雜的問題，是本模型求解過程中的一個創(chuàng)新點。由擬合函數(shù)預測出未來十年的長江未來水質(zhì)污染的發(fā)展趨勢。在對長江水質(zhì)的評價的過程中主要運用了權重系數(shù)理論，對長江水質(zhì)的預測方面主要采用了曲線擬合的思想，這兩種思想有很強的適應性，因此，該模型有很好的推廣性，可以用它推廣到對其它河流水質(zhì)的評價與預測，如黃河、珠江、松花江等。31241402242262323320421262325043202623263222820120292123522 6 2 2 201 3 302 3 2 201 27 2 2 2 3 2 2 )。for k0=1:3 max0(1,k0)=max(x0(:,k0))。 elseif x1(k0,4)==6 n1(1,6)=n1(1,6)+1。a0=zeros(3,1)。 a0=a。]。)。endy0=[max0(1,1) 1/max0(1,2) 1/max0(1,3) 1 2 2 6 4 3 5 6 1 4 3 10 5 2 15 2 6]。for k0=1:4 x1(:,k0)=x0(:,k0)。 elseif x1(k1,4)==6 x1(k1,5)=x1(k1,5)*y1(6,1)。 else 。x1(:,5)f3的直接運行結果：f3ans = ans = ：模型二的程序%求各觀測站點的污染物的排放量function D=fun1(E,F,A,C)T= zeros(1,7)。x1(:,5)=x1(:,5)/(sum(x1(:,5))/m)。y0(:,3)=1./y0(:,3)。max0=zeros(1,3)。39。%本函數(shù)存盤的函數(shù)名為“f2”ss=0。 a=a/sum(a)。y1=zeros(6,1)。endy0=[max0(1,1) 1/max0(1,2) 1/max0(1,3) 1 2 2 6 4 3 5 6 1 4 3 10 5 2 15 2 6]。y0=zeros(6,4)。39。2 2 201 3 3 2 2 2 3 2 3 2 2 3 2 2 1 232232322252242224220124240222302140623422232222222還可以用于對石油、礦產(chǎn)的開采進行評價與預測等。該模型沒有考慮枯水期和豐水期的數(shù)據(jù)，但水文年是指在一年內(nèi)所有檢測數(shù)據(jù)的平均值，因此所得出的結果是很強的參考性，可以為更好更合理地保護長江水質(zhì)應采取的措施提供依據(jù)。那么由此求解出的最佳權重系數(shù)可能也是不同的，進而求得的各地區(qū)在2003年7月份至2005年9月份間每個月份的綜合評估值也可能是不同的。②、調(diào)整產(chǎn)業(yè)結構和工業(yè)布局。因此，我們認為，解決水