【正文】 end disp(uint16(yc0(2:1:n+1)))； 。 end yc0(1)=x0(1)。 for k=1:n c=x0(1)au(2)/au(1)。 au=b\y。 end b(:,2)=1。 end for i=1:n1 b(i,1)=z(i+1)。 for i=1:n x1(i)=sum(x0(1:i))。) 附錄 8 由于 8 種重金屬的高斯擴(kuò)散模型的程序都差不多，下面以 As 的高斯擴(kuò)散模型為例 . x0= DATA(:,1)。plot3(x,y,a,39。 mesh(X,Y,A)。v439。 [X,Y,A]=griddata(x,y,a,linspace(0,27823)39。y= data(:,2)。fontweight:normal。fontsize:13pt。 /CELLRANGE=full /READNAMES=on /ASSUMEDSTRWIDTH=32767. DATASET NAME DataSet1 WINDOW=FRONT. FACTOR /VARIABLES V1 V2 V3 V4 V5 V6 V7 V8 /MISSING LISTWISE /ANALYSIS V1 V2 V3 V4 V5 V6 V7 V8 /PRINT INITIAL EXTRACTION /PLOT EIGEN /CRITERIA FACTORS(3) ITERATE(25) /EXTRACTION PC /ROTATION NOROTATE /SAVE REG(ALL) /METHOD={color:0。 /SHEET=name 39。textdecoration:none} 附錄 6 GET DATA /TYPE=XLS /FILE=39。fontstyle:normal。fontfamily:Monospaced。附件 239。D:\課件 \畢業(yè)論文 \A\濃度表 .xls39。 maxiP =。 maxiP =。 maxiP =。 maxiP =。 maxiP =。v39。 z5=data5(:,3)。) ,hold on x5= data5(:,1)。 scatter3(x4,y4,z4,39。 y4=data4 (:,2)。ks39。 z3=data3(:,3)。),hold on x3= data3(:,1)。 scatter3(x2,y2,z2,39。 y2=data2 (:,2)。*39。 z1=data1 (:,3)。 hold on x1= data1(:,1)。)。,linspace(0,18449),39。v39。 z5=data5(:,3)。) ,hold on x5= data5(:,1)。 scatter3(x4,y4,z4,39。 y4=data4 (:,2)。ks39。 z3=data3(:,3)。),hold on x3= data3(:,1)。 scatter3(x2,y2,z2,39。 y2=data2 (:,2)。*39。)。,linspace(0,18449),39。 z1=data1 (:,3)。 figure,surf(X,Y,Z) 附錄 2 利用 excel 把采樣點(diǎn)位置數(shù)據(jù)按不同區(qū)域分開，分別存為 data1， data2， data3， data4， x,y,z 與附錄 1 中表示的數(shù)據(jù)相同 . x1= data1(:,1)。v439。 [X,Y,Z]=griddata(x,y,z,linspace(0,27823)39。 y=data (:,2)。v439。v439。v439。 3Pi? :嚴(yán)重污染 . 下面用 Excel 統(tǒng)計(jì)得出不同功能區(qū)中重金屬元素的含量平均值，并將得出的數(shù)據(jù)作為 Ci ， 8 種重金屬的背景值分別為： ， ， 31， ， ， ， 31， 69 (mg/kg)，把它們作為 Si .得到各種金屬在不同區(qū)域中的污染程度，如表 2. 表 1 各種金屬不同區(qū)域中的平均值 (mg/kg) 元素 As Cd Cr Cu Hg Ni Pb Zn 生活區(qū) 工業(yè)區(qū) 山區(qū) 交通區(qū) 公園綠地區(qū) 表 2 各種 金屬不同區(qū)域中的污染指數(shù) 元素 As Cd Cr Cu Hg Ni Pb Zn 生活區(qū) 工業(yè)區(qū) 山區(qū) 交通區(qū) 公園綠地區(qū) 從表 1 中，可以得出，在生活區(qū)， Cu、 Zn 屬于嚴(yán)重污染， Hg 、 Cd、 Cr、 Pb 屬于中度污染， As、 Ni 屬于輕度污染；在工業(yè)區(qū)， Hg、 Cu、 Zn、 Cd、 Pb 屬于嚴(yán)重污染， As 屬于中度污染， Cr、 Ni 屬于輕度污染；在山區(qū)，污染較輕， 8 種重金屬都屬于輕度污染；在交通區(qū)， Hg、 Cu、 Zn 屬于嚴(yán)重污染， Cd、 Pb 屬于中度污染， Cr、 As、 Ni 屬于輕度污 9 染；在公園綠地區(qū)， Hg 屬于嚴(yán)重污染， Cu、 Zn、 Cd 屬于中度污染， Pb、 As、 Cr、 Ni 屬于輕度污染 . 模型 [3] 由于單因子指數(shù)評價(jià)法只能估計(jì)每種重金屬在某一區(qū)域的污染程度，但如果要評價(jià)某一區(qū)域的綜合污染程度，就要把各種重金屬在同一區(qū)域的污染指數(shù)結(jié)合起來，下面就用內(nèi)梅羅指數(shù)評價(jià)法來評價(jià)五個(gè)功能區(qū)的污染程度 . 綜合污染評價(jià)不僅兼顧了單因子污染指數(shù)平均值和最高值，還可以突出污染較重的重金屬污染物的作用 .其計(jì)算公式為：,2 2m a x2 iPPP ?? ）（綜 其中 P綜 是采樣點(diǎn)的綜合污染指數(shù)； maxiP 為 i 采樣點(diǎn)重金屬污染物單項(xiàng)污染指數(shù)中的最大值； ???ni iPnP 11 為單因子污染指數(shù)平均值 .采用《土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)》中國家一級標(biāo)準(zhǔn)作為環(huán)境質(zhì)量的標(biāo)準(zhǔn)， 8 種重金屬的標(biāo)準(zhǔn)值分別為： 15， ， 90，35， ， 40， 35， 100(mg/kg).用 Excel 表格統(tǒng)計(jì)得到各種金屬在不同區(qū)域中的單項(xiàng)污染程度，如表 3. 表 3 各種金屬不同區(qū)域中的單項(xiàng)污染指數(shù) 元素 As Cd Cr Cu Hg Ni Pb Zn 生活區(qū) 工業(yè)區(qū) 山區(qū) 交通區(qū) 公園綠地區(qū) 得出上表后，利用 Excel 表格統(tǒng)計(jì)得出單因子污染指數(shù)平均值 P ，再通過內(nèi)羅梅綜合評價(jià)公式計(jì)算得到各種金屬不同區(qū)域中的綜合污染指數(shù)，如表 4，編程實(shí)現(xiàn) (附錄 4). 內(nèi)梅羅指數(shù)以 5 為界，劃分為清潔、輕污染、污染、重污染、嚴(yán)重污染 5 個(gè)等級 . 10 表 4 各種金屬不同區(qū)域中的綜合污染指數(shù) 功能區(qū) 污染指數(shù) 污染等級 生活區(qū) 輕污染 工業(yè)區(qū) 重污染 山區(qū) 清潔 交通區(qū) 污染 公園路地區(qū) 輕污染 從上表可以看出，工業(yè)區(qū)的污染最嚴(yán)重，污染等級為重污染，其次，交通區(qū)，污染等級為污染，生活區(qū)和公園綠地區(qū)的污染較輕，為輕污染，山區(qū)的污染最輕，清潔 . 問題 2 的分析與求解 [4] 由問題的求解中，知道山區(qū)的人為污染比較輕，而其他四個(gè)不同功能區(qū)都受到不同程度的人類活動(dòng)污染 .其中，生活區(qū)的重金屬污染可能來源于污水的排放、煤煙的排放和垃圾的污染，工業(yè)區(qū)的重金屬污染可能來源冶煉金屬的廢氣、廢水排 放和廢棄物中的金屬污染，公園路地區(qū)的重金屬污染可能來源于農(nóng)藥的使用 . 由于不同的重金屬可能來源于相同的污染源，下面先進(jìn)性 KMO檢驗(yàn) .該檢驗(yàn)從比較觀測變量之間的簡單相關(guān)系數(shù)和偏相關(guān)系數(shù)的相對大少出發(fā)，其值的變化從 0到 1，一般 KMO大于 ，如果 KMO大于 ，說明還可以做因子分析，若 KMO過小，表明不適合做因子分析 .相關(guān)性檢驗(yàn)結(jié)果如表 5， 編程實(shí)現(xiàn) (附錄 5). 表 5 KMO檢驗(yàn)和 and Bartlett檢驗(yàn) KaiserMeyerOlkin Measure of Sampling Adequacy. .778 Bartlett39。1 2Pi? :輕度污染 。 Gaussian diffusion model。 Correlation analysis。