【正文】 for i=1:2M1(i)=A1(i)*A1(i+2)。WP1=0。W1=A1(:,1)。return。A1矩陣輸入錯誤，請重新輸入A1矩陣!39。)。switch ncase 2disp(39。)。A1=input(39。請輸入判斷矩陣A1(一級判斷)39。 將此程序代碼直接輸入MATLAB的命令窗口中，或新建一個M文件，再調(diào)用該M文件的文件名，都可執(zhí)行以上的程序，可分別得到各個區(qū)縣的評價結(jié)果。)。%計算評估分值大于等于500enddisp(39。%計算評估分值大于等于400但小于500otherwisedengji=39。%計算評估分值大于等于300但小于400case {4} dengji=39。%計算評估分值大于等于200但小于300case {3} dengji=39。%計算評估分值小于200case {2} dengji=39。switch fix(KP/100)case {0,1} dengji=39。)。disp(39。P=D*Z39。)。end%計算評估分值disp(39。for p=1:jB=W1(1).*M1。B=M1(:)。[J,j]=size(M1)。)。Y=input(39。請輸入社會因子相對比例標(biāo)度39。)。X=input(39。請輸入自然因子相應(yīng)比例標(biāo)度39。M2=WT2239。W1=WT139。名山縣39。雨城區(qū) 39。滎經(jīng)縣 39。天全縣 39。蘆山縣 39。寶興縣 39。漢源縣 39。石棉縣 39。return。A22矩陣的一致性不滿足要求，請返回重新評分!39。)。disp(39。if CRCI。RI=[0,0,9/10,]。disp(S)。AW22=39。for i=1:nS=S+AW22(i)./(n*WT22(i))。endWT22=W22/WP。endW22=M22.^(1/n)。S=0。WT22=A22(:,1)。%A22矩陣輸入錯誤!endM22=A22(:,1)。)。otherwisedisp(39。A22矩陣輸入正確，可繼續(xù)進行下一步39。%判斷矩陣定義[n,n]=size(A22)。A22=39。)。enddisp(39。)。elsedisp(39。A21矩陣的一致性可以接受!39。CR。CR=CI/RI(n)。disp(WT21)%以下是一致性檢驗CI=(Sn)/(n1)。)。enddisp(39。AW21=A21*WT21。for i=1:nWP=WP+W21(i)。for i=1:6M21(i)=A21(i)*A21(i+6)*A21(i+12)*A21(i+18)*A21(i+24)*A21(i+30)。WP=0。W21=A21(:,1)。return。A21矩陣輸入錯誤，請重新輸入A21矩陣!39。)。switch ncase 6disp(39。)。A21=input(39。請輸入判斷矩陣A21(二級判斷自然因子)39。)。disp(WT1)%以下是一致性檢驗disp(39。)。enddisp(39。AW1=A1*WT1。for i=1:nWP1=WP1+W1(i)。for i=1:2M1(i)=A1(i)*A1(i+2)。WP1=0。W1=A1(:,1)。return。A1矩陣輸入錯誤，請重新輸入A1矩陣!39。)。switch ncase 2disp(39。)。A1=input(39。請輸入判斷矩陣A1(一級判斷)39。 MATLAB工具編程 依據(jù)層次分析和模糊綜合評判評價雅安水環(huán)境生態(tài)風(fēng)險的步驟，為避免復(fù)雜的人工評價方法，在此可運用MATLAB編程，以快速、準(zhǔn)確的得到各個區(qū)縣的風(fēng)險等級。(5) 評價矩陣建立在建立了各級因素集、各級權(quán)重集之后，又有了針對各種不同風(fēng)險等級的隸屬度，就可構(gòu)建相應(yīng)的評價矩陣。隸屬度的確定，需要專家通過經(jīng)驗確定。建立因素子集——社會因子的因素子集：其中，城市居民用水、農(nóng)村居民用水均指人均生活用水，與已收集資料對應(yīng)。各因子的相對權(quán)重值設(shè)為，則各個區(qū)縣的評估分值為： （其中：i取1～8）則可得出各個區(qū)縣的評過分值。表10 自然因子權(quán)重　供水總量（M1）多年平均蒸發(fā)量（M2）多年平均徑流深（M3）多年平均降水量（M4）水資源總量（M5）生態(tài)環(huán)境用水（M6）石棉縣533511漢源縣335531寶興縣533511蘆山縣113351天全縣511311滎經(jīng)縣531311雨城區(qū)511151名山縣535351表11 社會因子權(quán)重城市居民人均生活用水（M7）農(nóng)村居民人均生活用水（M8）生產(chǎn)用水（M9）廢污水年排放量（M10）石棉縣5315漢源縣1351寶興縣1111蘆山縣3513天全縣5515滎經(jīng)縣3113雨城區(qū)3339名山縣3155(3) 計算評估分值綜上可知，層次分析法評價雅安水環(huán)境生態(tài)風(fēng)險的指標(biāo)體系及各因素兩兩之間的權(quán)重值已得出，見表12所示。經(jīng)專家的綜合評價，最后取其平均值即可作為相應(yīng)因子的權(quán)重。結(jié)果如下：一級判斷：二級判斷：④ 判斷矩陣的一致性檢驗依據(jù)前述方法，檢驗一級判斷和二級判斷的特征向量，經(jīng)MATLAB編程檢驗，可知以上的結(jié)果滿足要求。根據(jù)相關(guān)專家意見，一級判斷和二級判斷中的自然因子和社會因子的結(jié)果分別見表7～9。 表5 自然因子數(shù)據(jù)統(tǒng)計NMinimumMaximumMeanStd. Deviation多年平均降水量81240多年平均蒸發(fā)量88391573多年平均徑流深87371702水資源總量8供水總量8生態(tài)環(huán)境用水8Valid N (listwise)8表6 社會因子數(shù)據(jù)統(tǒng)計NMinimumMaximumMeanStd. Deviation城市居民人均生活用水8農(nóng)村居民人均生活用水8生產(chǎn)用水8廢污水年排放量8Valid N (listwise)8 風(fēng)險分析與評價 層次分析法 層次分析法的步驟如下：(1) 生態(tài)風(fēng)險程度及相應(yīng)的比例標(biāo)度① 專家打分進行層次分析評價，要確定各個同層次的評價因子之間的相對比例標(biāo)度，評價的區(qū)域不同，各種因子的選取也就不同。其中：N代表樣本個數(shù)，Minimum、Maximum、Mean、Std. Deviation分別代表樣本的最小值、最大值、平均值和方差。多年平均降水量多年平均蒸發(fā)量多年平均徑流深廢污水年排放量生產(chǎn)用水自然因子社會因子雅安水環(huán)境生態(tài)風(fēng)險城市人均生活用水水資源總量供水總量農(nóng)村人均生活用水生態(tài)環(huán)境用水量圖4 雅安水環(huán)境生態(tài)風(fēng)險模型 獲取基本數(shù)據(jù)通過資料的收集與校核，得到自然因子和社會因子的已知數(shù)據(jù)，分別見表3和表4。經(jīng)分析可知，生產(chǎn)用水、廢污水年排放量、城市居民人均生活用水、農(nóng)村居民人均生活用水的數(shù)值愈大時，對水環(huán)境的生態(tài)風(fēng)險影響也愈大。我國是一個農(nóng)業(yè)大國，農(nóng)村居民生活質(zhì)量的好壞直接關(guān)系到社會經(jīng)濟的全面發(fā)展，因此在農(nóng)村居民人均生活用水不滿足要求時，對水環(huán)境也會產(chǎn)生一定的負面影響。因此，通常用城市居民人均生活用水來衡量城市正常運行的標(biāo)準(zhǔn)之一。(3) 城市居民人均生活用水城市用水是城市生產(chǎn)、生活、消防和市政管理等活動所需用水的統(tǒng)稱。生產(chǎn)污水需要進行凈化處理；生產(chǎn)廢水不需要凈化處理或僅需做簡單的處理，如冷卻處理。工業(yè)廢水可分為生產(chǎn)污水與生產(chǎn)廢水。生活污水一般不含有毒物質(zhì)，但是它有適合微生物繁殖的條件，含有大量的病原體，從衛(wèi)生角度來看有一定的危害性。生活污水是人類在日常生活中使用過的，并被生活廢料所污染的水。生產(chǎn)用水總量，其表示方法為一年中工業(yè)用水和農(nóng)業(yè)用水的總和。(1) 生產(chǎn)用水簡而言之，生產(chǎn)用水(Process Water)即生產(chǎn)過程所需用的水，通常包括工業(yè)生產(chǎn)用水和農(nóng)業(yè)生產(chǎn)用水。特別是當(dāng)多年平均降水量、多年平均徑流深、水資源總量、供水總量、生態(tài)環(huán)境用水的值愈小時，生態(tài)風(fēng)險等級就愈大，多年平均蒸發(fā)量的數(shù)值愈大時，生態(tài)風(fēng)險等級也愈大。進一步來說，維系生態(tài)系統(tǒng)最基本的生存條件及其最基本的環(huán)境服務(wù)功能所需求的水資源閾值，通常稱為最小生態(tài)環(huán)境需水。(6) 生態(tài)環(huán)境用水從廣義上講，維持全球生物化學(xué)平衡(諸如水熱平衡、水沙平衡、水鹽平衡等)所消耗的水分都是生態(tài)環(huán)境需水(劉昌明等，1996)。通常包括：河川徑流量、潛水蒸發(fā)量、地下潛流量之和。當(dāng)徑流深取多年平均值時，則為多年平均徑流深。(3) 多年平均徑流深徑流深是在某一時段內(nèi)，通過河流上指定斷面的徑流總量（W，以m3計）除以該斷面以上的流域面積（F，以km2計）所得的值。雨量稀少、地下水源及流入徑流水量不多的地區(qū)，如蒸發(fā)量很大，即易發(fā)生干旱。一般溫度越高、濕度越小、風(fēng)速越大、氣壓越低、則蒸發(fā)量就越大；反之蒸發(fā)量就越小。而蒸發(fā)量是指在一定時段內(nèi)，水分經(jīng)蒸發(fā)而散布到空氣中的量。通常是衡量一個地區(qū)在某段時間內(nèi)降水多少的數(shù)據(jù)，并取其多年平均值，稱為多年平均降水量。 自然因子自然因子包括研究區(qū)的多年平均降水量、多年平均蒸發(fā)量、多年平均徑流深、水資源總量、供水總量、生態(tài)環(huán)境用水等。因雅安處于長江上游，受水質(zhì)等因子影響較小，水量對其生態(tài)環(huán)境需水的影響比水質(zhì)要大得多，因而在進行水環(huán)境生態(tài)風(fēng)險評價時，可考慮能否為生態(tài)環(huán)境提供足夠的水量。通過已收集的資料和分析，初步選定各種影響雅安水環(huán)境生態(tài)風(fēng)險的各種風(fēng)險源。 雅安水環(huán)境生態(tài)風(fēng)險評價 問題分析與模型建立結(jié)合研究區(qū)域的實際情況，考慮各種自然因子和社會因子對水生態(tài)系統(tǒng)的影響范圍和程度、資料收集的難度和四川農(nóng)業(yè)大學(xué)倪福全的四川西緣山地農(nóng)村飲水安全的研究成果，在此，可按行政區(qū)進行劃分。特別是1996年以來，我國GIS技術(shù)在技術(shù)研究、成果應(yīng)用、人才培養(yǎng)、軟件開發(fā)等方面進展迅速，并力圖將GIS從初步發(fā)展時期的研究實驗、局部應(yīng)用推向?qū)嵱没?、集成化、工程化，為國民?jīng)濟發(fā)展提供輔助分析和決策依據(jù)。所謂空間數(shù)據(jù)，是指對出現(xiàn)在地球表面及其附近的地理現(xiàn)象的記錄，都與地形有關(guān)，都具有地理位置特征，均可以用地圖予以表示。它是以地理空間數(shù)據(jù)庫為基礎(chǔ)，在計算機硬件的支持下，對空間數(shù)據(jù)進行采集、管理、操作、分析、模擬和顯示，并采用地理模型分析方法，適時提供多種空間和動態(tài)地理信息的計算機系統(tǒng)技術(shù)。功能件工具箱主要用來擴充其符號計算功能、可視建模仿真功能及文字處理功能等；學(xué)科性工具箱專業(yè)性比較強，如控制系統(tǒng)工具箱(Control System Toolbox)、信號處理工具箱(Signal Processing Toolbox)、神經(jīng)網(wǎng)絡(luò)工具箱(Neural Network Toolbox)、最優(yōu)化工具箱(Optimization Toolbox)、金融工具箱(Financial Toolbox)、統(tǒng)計學(xué)工具箱(Statistics Toolbox)等，這些工具箱都是由該領(lǐng)域內(nèi)學(xué)術(shù)水平很高的專家編寫的，用戶可以直接利用這些工具箱進行相關(guān)領(lǐng)域的科學(xué)研究?；静糠謽?gòu)成了MATLAB的核心內(nèi)容，也是使用和構(gòu)造工具箱的基礎(chǔ)。MathWorks公司有心讓MATLAB成為新一代的通用軟件開發(fā)工具，并為此提供了將MATLAB源程序編譯為獨立于MATLAB集成環(huán)境運行的EXE文件以及將MATLAB程序轉(zhuǎn)化為C語言程序的編譯器。因此，對于從事數(shù)值計算、計算機輔助設(shè)計和系統(tǒng)仿真等領(lǐng)域的人員來說，用MATLAB編程的確是一個理想選擇。利用MATLAB圖形句柄操作，用戶可以更靈活地對圖形進行各種操作，為用戶在圖形表現(xiàn)方面開拓了一個沒有束縛的廣闊空間。MATLAB提供了兩個層次的繪圖操作：一種是對圖形句柄進行的低層繪圖操作，另一種是建立在低層繪圖操作之上的高層繪圖操作。MATLAB和著名的符號計算語言Maple相結(jié)合，使得MATLAB具有了符號計算功能。高質(zhì)量的數(shù)值計算功能為MATLAB贏得了聲譽。 概括來講，其主要的功能包含如下[49]： (1) 數(shù)值計算和符號計算功能 MATLAB以矩陣作為數(shù)據(jù)操作的基本單位，這使得矩陣運算變得非常簡捷、方便、高效。在歐美各高等院校MATLAB已經(jīng)成為線性代數(shù)、自動控制理論、數(shù)字信號處理、時間列序分析、動態(tài)系統(tǒng)仿真、圖像處理等諸多課程的基本教學(xué)工具，成為本科生，碩士生及博士生必須掌握的基本技能，在設(shè)計研究單位和工業(yè)部門，MATLAB已被廣泛的應(yīng)用與研究和解決各種工具的工程問題。 MATLAB MATLAB是英文MATrix LABoratory（矩陣實驗室）的縮寫。 計算機輔助工具 本次研究數(shù)據(jù)處理的工作量較大，特別是二維矩陣的分析和計算；同時，所得結(jié)果需要可視化的表達和顯示，因而都需要借助于一些計算機輔助工具。模糊綜合評價可表示為：式中為模糊綜合評價集；為權(quán)重集（權(quán)重向量），可視為一個行模糊矩陣。(5) 模糊綜合評判單因素模糊評價，僅反映了一個因素對評價對象的影響，顯然這是不全面的。表示和之間隸屬“合理”關(guān)系的程度，即按評價時，評價對象取的合理程度。顯然，為一模糊矩陣。(4) 單因素模糊評價單獨從一個因素出發(fā)進行評價，以確定評價對象對