【正文】 *39。 cc(i,2)=c(kkk(i),2)。 n=size(kkk,1) for i=1:n rr(i)=pp(kkk(i))。 kkk=find(kk500)。 for i=1:n kk(i)=k(ii(i),jj(i))。 c=[ii,jj]。 end。 end。(z(i,j)z(i+1,j+1))。(z(i,j)z(i+1,j1))amp。(z(i,j)z(i1,j+1))amp。(z(i,j)z(i1,j1))amp。(z(i,j)z(i,j1))amp。(z(i,j)z(i,j+1))amp。(z(i,j)z(i+1,j))amp。 for i=2:99 for j=2:99 20 if (z(i,j)z(i1,j))amp。%插值 z=Z39。v439。 [X,Y,Z]=griddata(x,y,z,linspace(0,30000)39。y=A(:,3)。}。,39。,39。,39。,39。,39。,39。,39。 ss={39。B4:I32239。F:\A\cumcm2020A 附件 _數(shù)據(jù) .xls39。)。,1,39。),title(ss{k}) end 附錄五： 篩選 污染源 的 MATLAB程序 （以砷為例，其它元素的程序只需調(diào)整參數(shù)即可）： A=xlsread(39。 subplot(2,2,k4),scatter(ii,jj,39。disp(jj39。 disp(ii39。 end。 end。(z(i,j)z(i+1,j+1))。(z(i,j)z(i+1,j1))amp。(z(i,j)z(i1,j+1))amp。(z(i,j)z(i1,j1))amp。(z(i,j)z(i,j1))amp。(z(i,j)z(i,j+1))amp。(z(i,j)z(i+1,j))amp。 for i=2:99 for j=2:99 if (z(i,j)z(i1,j))amp。)。,linspace(0,20200),39。 for k=5:8 z=B(:,k)。 x=A(:,2)。Zn39。Pb39。Ni39。Hg39。Cu39。Cr 39。Cd39。As 39。)。,2,39。 B=xlsread(39。A4:E32239。F:\A\cumcm2020A 附件 _數(shù)據(jù) .xls39。*39。)。)。 [ii,jj]=find(z==1000)。 end。 z(i,j)=1000。amp。amp。amp。amp。amp。amp。amp。%插值 z=Z39。v439。 [X,Y,Z]=griddata(x,y,z,linspace(0,30000)39。y=A(:,3)。}。,39。,39。,39。,39。,39。,39。,39。 ss={39。B4:I32239。F:\A\cumcm2020A 附件 _數(shù)據(jù) .xls39。)。,1,39。 c(i,j)=a(i,j)/b(j) end end 附錄 三 ： 由 各區(qū)的平均單項污染指數(shù) Pj,ave 和最大單項污染指數(shù) Pj,max 求各區(qū)的綜合污染指 數(shù)的 MATLAB 程序： Pj,ave Pj,max function f=fun4(x) a=((x(1)^2+x(2)^2)/2)^(1/2) 附錄 四 ： 重金屬元素砷、鎘、鉻、銅的污染源分布圖的 MATLAB 程序： A=xlsread(39。 b=b39。 n=size(a,2)。功能區(qū) 39。功能區(qū) 39。%插值 pcolor(X,Y,Z)。v439。 scatter(x,y,5,z)%散點圖 figure [X,Y,Z]=griddata(x,y,z,linspace(0,30000)39。y=A(:,3)。)。,2,39。 B=xlsread(39。A4:E32239。F:\A\cumcm2020A 附件 _數(shù)據(jù) .xls39。公園綠地區(qū) 39。主干道區(qū) 39。山區(qū) 39。工業(yè)區(qū) 39。生活區(qū) 39。功能區(qū)分布 39。x39。 y=y(x5(:))。 x5=find(z==5)。 x=A(:,2)。p39。 y=y(x4(:))。 x4=find(z==4)。 x=A(:,2)。s39。 y=y(x3(:))。 x3=find(z==3)。 x=A(:,2)。h39。 y=y(x2(:))。 x2=find(z==2)。 x=A(:,2)。d39。 y=y(x1(:))。 x1=find(z==1)。y=A(:,3)。)。,1,39。) figure,surf(X,Y,Z)%三維曲面 功能區(qū)分布散點圖的 MATLAB 程序： A=xlsread(39。) figure,contourf(X,Y,Z) %等高線圖 figure,contour(X,Y,Z) title(39。shading interp%偽彩色圖 title(39。)。,linspace(0,20200),39。z=A(:,4)。 x=A(:,2)。A4:E32239。F:\A\cumcm2020A附件 _數(shù)據(jù) .xls39。 [5] 周濤發(fā) 岳書倉 柏林，《城市地質(zhì)環(huán)境及其評價與保護(hù)》 [J].安徽：合肥工業(yè)大學(xué)學(xué)報，第 20卷第 3期， 2227， 1997。 [3] GB156181995， 《 土壤環(huán)境質(zhì)量標(biāo)準(zhǔn) 》 [S].北京：中國標(biāo)準(zhǔn)出版社 ， 1995。 七、 參考文獻(xiàn) [1] 張志涌 ， 《精通 MATLAB 版》 [M].北京：北京航天航空大學(xué)出版社， 234302，2020。 對于問題三所建模型，我們在研究城市地質(zhì)環(huán)境的演變模式時，我們僅對城市海拔進(jìn)行了分析。 對于問題二所建模型，我們在分析重金屬污染的主要原因時，僅考慮了城市內(nèi)各區(qū)造成的污染，忽略了該城市周邊農(nóng)田中農(nóng)藥的使用等因素造成的污染。 六、模型的改進(jìn)與推廣 對于問題一所建模型，我們在求各區(qū)污染程度的時候，僅僅拿各區(qū)重金屬的平均濃度進(jìn)行分析，得出的結(jié)果只能反映各區(qū)的平均污染等級，不能反映各區(qū)在不同位置的污染等級。城市地質(zhì)環(huán)境是一個涉及到地球巖石圈表層的巖石、土壤、大氣、水和生物的復(fù)雜系統(tǒng)。 附表所給數(shù)據(jù)是靜態(tài)的 ， 無法根據(jù)所給數(shù)據(jù)建立城區(qū)污染的動態(tài)演化過程。 缺點：解決 問題一 第二小問時，我們把各區(qū)內(nèi)采樣點重金屬濃度實測值的平均值用作各區(qū)重金屬濃度的實測值 ，經(jīng)過內(nèi)梅羅綜合污染指數(shù)評價法進(jìn)行求解得出的各區(qū)污染等級只能反映各區(qū)的平均污染等級，不能反映各個采樣點各自的污染等級。然后結(jié)合國家土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)對污染源進(jìn)行篩選 ， 能方便的求出各種重金屬元素的主要污染源。模型三中依據(jù)模 型一中建立的濃度分布矩陣建立了遍歷搜索模型。通過綜合比較之后選定三維 surf曲面建立重金屬元素的空間分布模型 ， 直觀明了。由于所給數(shù)據(jù)采樣點的不規(guī)則性 ， 首先使用 griddata函數(shù)對所給數(shù)據(jù)進(jìn)行插值規(guī)整得出一個X,Y分別等步長的某種元素的濃度分布矩陣。 得到各種重金屬元素 空間分布極大值點的散點圖 如下： 11 圖 9：砷、鎘、鉻、銅的空間分布極大值散點圖 圖 10：汞、鎳、鉛、鋅的空間分布極大值散點圖 結(jié)合國家土壤環(huán)境 三個等級的 質(zhì)量標(biāo)準(zhǔn)通過 MATLAB軟件對極值點進(jìn)行 分級 篩選 ：首先用國家 土壤環(huán)境 一級 質(zhì)量標(biāo)準(zhǔn) 進(jìn)行篩選，得出篩選結(jié)果，再用國家 土壤環(huán)境 二級 質(zhì)12 量標(biāo)準(zhǔn) 對一級指標(biāo)得出的點進(jìn)行篩選 ， 依次類推，最終得到篩選結(jié)果 如表 7 表 7：不同國標(biāo)等級下的極大值個數(shù)