【正文】 for k=1:nyc0(k+1)=ycl(k+1)ycl(k)。 ycl(1)=x0(1)。 y(i)=x0(i+1)。n=length(x0)。hold on。,linspace(0,18449),39。textdecoration:none}附錄7由于8種重金屬的空間分布圖的程序都差不多，下面以As的空間分布圖為例.x=data(:,1)。fontfamily:Monospaced。/SAVEITERATE(25)EIGENINITIALV7V3V8V4/VARIABLESNAME附件239。/TYPE=XLSfontweight:normal。/METHOD={color:0。ITERATE(25)/EXTRACTION/CRITERIAROTATION/PRINTV6V2LISTWISEV7V3/ASSUMEDSTRWIDTH=32767.DATASET/CELLRANGE=full39。D:\課件\畢業(yè)論文\A\39。 =sqrt((^2+^2)/2)公園路地區(qū)：= mean(a(5,:))。 =sqrt((^2+^2)/2)山區(qū)：= mean(a(3,:))。)附錄4矩陣a表示下表中各種重金屬在不同區(qū)域中的單項污染指數(shù)元素AsCdCrCuHgNiPbZn生活區(qū)工業(yè)區(qū)山區(qū)交通區(qū)公園綠地區(qū)生活區(qū)：= mean(a(1,:))。 y5=data5 (:,2)。 z4=data4(:,3)。scatter3(x3,y3,z3,39。ro39。),hold onx2= data2(:,1)。 y1=data1 (:,2)。v439。scatter3(x5,y5,z5,39。mp39。) ,hold onx4= data4(:,1)。 y3=data3 (:,2)。 z2=data2 (:,3)。contourf(X,Y,Z) ,hold onscatter3(x1,y1,z1,39。[X,Y,Z]=griddata(x,y,z,linspace(0,27823)39。)。 z=data (:,3)。三次多項式插值法分別對8種重金屬污染物的濃度進行數(shù)據(jù)插值，編程實現(xiàn)(附錄7).畫出圖像后，觀察圖像的極值點，根據(jù)附件2的數(shù)據(jù)，利用MATLAB的搜索函數(shù)find，先從重金屬的最高濃度開始查找，然后找到該點對應的x，y坐標，最后把找到的坐標跟圖像對比，符合的就是圖像的極值點，找到第一點后，再用次高濃度查找，直到找到所有的極值點為止. 2. 模型的結果圖5 As的濃度圖從圖5中，可以看出As的濃度圖中分別有2個極值點，可以大概確定極值點的位置為：(18134，10046)、(12696，3024).圖6 Cd的濃度圖從圖6中，可以看出Cd的濃度圖中分別有1個極值點，可以大概確定極值點的位置為：(21439，11383).圖7 Cr的濃度圖從圖7中，可以看出Cr的濃度圖中分別有2個極值點，可以大概確定極值點的位置為：(3299，6018)、(4592，4603).圖8 Cu的濃度圖從圖8中，可以看出Cu的濃度圖中分別有2個極值點，可以大概確定極值點的位置為：(2383，3692)、(3299，6018).圖9 Hg的濃度圖從圖9中，可以看出Hg的濃度圖中分別有4個極值點，可以大概確定極值點的位置為：(2383，3692)、(2708，2295)、(13694，2357)、(15248，29106).圖10 Ni的濃度圖從圖10中，可以看出Ni的濃度圖中分別有1個極值點，可以大概確定極值點的位置為： (3299，6018).圖11 Pb的濃度圖從圖11中，可以看出Pb的濃度圖中分別有2個極值點，可以大概確定極值點的位置為： (1647，2728)、(4777，4897).圖12 Zn的濃度圖從圖11中，可以看出Zn的濃度圖中分別有5個極值點，可以大概確定極值點的位置為： (9328，4311)、(13797，9621)、(12696，3024)、(4948，7293)、(2383，3692).根據(jù)各種重金屬濃度圖中的極值點，可以確定各種重金屬污染源的位置，如表9.表9 各種重金屬污染源的位置元素污染源位置/mAs(18134，10046)、(12696，3024)Cd(21439，11383)Cr(3299，6018)、(4592，4603)Cu(2383，3692)、(3299，6018)Hg(2383，3692)、(2708，2295)、(13694，2357)、(15248，29106)Ni(3299，6018)Pb(9328，4311)、(13797，9621)Zn(9328，4311)、(13797，9621)、(12696，3024)、(4948，7293)、(2383，3692)由上表可以看出，As和Zn有共同的污染源，來自交通區(qū)的(12696，3024)；Cr、Cu、Ni有共同的污染源，來自交通區(qū)的(3299，6018)；Cu、Hg、Zn有共同的污染源，來自工業(yè)區(qū)的(2383，3692)；Pb、Zn有共同的污染源，來自生活區(qū)的(9328，4311)和交通區(qū)的(13797，9621).[5]根據(jù)重金屬在土壤中的擴散特點，該問題的解決還可以選用高斯擴散模型進行建模. 污染源選為坐標原點，污染處重金屬元素的濃度記為，則時刻處的無窮空間中的任一點重金屬濃度為 ，單位距離內(nèi)通過單位法向面積的流量與濃度梯度成正比，則是擴散系數(shù)，grad表示梯度，負號表示由濃度高向濃度低的地方擴散. 考察空間域 ，其體積為，包圍 的曲面為， 為一規(guī)則的球面， 外法線向量為 . 則在內(nèi)通過的流量內(nèi)重金屬元素的增量為又由于污染物釋放的重金屬的總量為由質(zhì)量守恒定律，有根據(jù)曲面面積的Gauss積分公式有其中表示散度記號，由式至式并利用積分中值定理整理得這是無界區(qū)域的偏微分方程，根據(jù)假設(5)，初始條件為作用在坐標原點的污染源函數(shù)，可記作表示釋放的重金屬元素總量， 其中為時間內(nèi)傳播的速度. 問題4的分析與求解(1)模型一8種重金屬元素空間分布圖比較清晰、直觀，分析也比較詳細；對污染程度，采用內(nèi)梅羅單項指數(shù)評價法和綜合指數(shù)評價法，從不同的角度分析，得到相同的結論，更準確的判斷出各個區(qū)域的污染情況.(2)模型二在進行主成分分析前，先