【正文】 plot3(x,y,z,39。 z=data(45:80,10)。 zi=griddata(x,y,z,xi,yi,39。o39。 ti=0:10:22674。cubic39。) %2區(qū) Pb分布圖 x=data(45:80,2)。 tu=0:10:17980。)。 y=data(45:80,3)。 [xi,yi]=meshgrid(ti,tu)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 z=data(81:146,6)。 zi=griddata(x,y,z,xi,yi,39。o39。 ti=0:10:28654。cubic39。) %3區(qū) Cr分布圖 x=data(81:146,2)。 tu=0:10:18397。)。 y=data(81:146,3)。 [xi,yi]=meshgrid(ti,tu)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 z=data(81:146,10)。 zi=griddata(x,y,z,xi,yi,39。o39。 ti=0:10:28654。cubic39。) %3區(qū) Pb分布圖 x=data(81:146,2)。 35 tu=0:10:18397。)。 y=data(81:146,3)。 [xi,yi]=meshgrid(ti,tu)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 z=data(147:284,6)。 zi=griddata(x,y,z,xi,yi,39。o39。 ti=0:10:26424。cubic39。) 36 %4區(qū) Cr分布圖 x=data(147:284,2)。 tu=0:10:17949。)。 y=data(147:284,3)。 [xi,yi]=meshgrid(ti,tu)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 z=data(147:284,10)。 zi=griddata(x,y,z,xi,yi,39。o39。 37 ti=0:10:26424。cubic39。) %4區(qū) Pb分布圖 x=data(147:284,2)。 tu=0:10:17949。)。 y=data(147:284,3)。 [xi,yi]=meshgrid(ti,tu)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 z=data(285:319,6)。 zi=griddata(x,y,z,xi,yi,39。o39。 ti=0:10:20582。cubic39。) %5區(qū) Cr分布圖 x=data(285:319,2)。 tu=0:10:17133。)。 y=data(285:319,3)。 [xi,yi]=meshgrid(ti,tu)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 z=data(285:319,10)。 zi=griddata(x,y,z,xi,yi,39。o39。 ti=0:10:20582。cubic39。) %5區(qū) Pb分布圖 x=data(285:319,2)。 tu=0:10:17133。)。 y=data(285:319,3)。 40 [xi,yi]=meshgrid(ti,tu)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 z=data(1:319,4)。 zi=griddata(x,y,z,xi,yi,39。o39。)。 tu=0:10:18449。) %城區(qū)地形圖 x=data(1:319,2)。cubic39。 ti=0:10:20582。o39。 zi=griddata(x,y,z,xi,yi,39。 z=data(285:319,12)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 [xi,yi]=meshgrid(ti,tu)。 y=data(285:319,3)。)。 tu=0:10:17133。) %5區(qū) Hg分布圖 39 x=data(285:319,2)。cubic39。 ti=0:10:20582。o39。 zi=griddata(x,y,z,xi,yi,39。 z=data(285:319,8)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 [xi,yi]=meshgrid(ti,tu)。 y=data(285:319,3)。)。 tu=0:10:17133。) %5區(qū) As分布圖 x=data(285:319,2)。cubic39。 ti=0:10:26424。o39。 zi=griddata(x,y,z,xi,yi,39。 z=data(147:284,12)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 [xi,yi]=meshgrid(ti,tu)。 y=data(147:284,3)。)。 tu=0:10:17949。) %4區(qū) Hg分布圖 x=data(147:284,2)。cubic39。 ti=0:10:26424。o39。 zi=griddata(x,y,z,xi,yi,39。 z=data(147:284,8)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 [xi,yi]=meshgrid(ti,tu)。 y=data(147:284,3)。)。 tu=0:10:17949。) %4區(qū) As分布圖 x=data(147:284,2)。cubic39。 ti=0:10:28654。o39。 zi=griddata(x,y,z,xi,yi,39。 z=data(81:146,12)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 [xi,yi]=meshgrid(ti,tu)。 y=data(81:146,3)。)。 tu=0:10:18397。) %3區(qū) Hg分布圖 x=data(81:146,2)。cubic39。 ti=0:10:28654。o39。 zi=griddata(x,y,z,xi,yi,39。 z=data(81:146,8)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 [xi,yi]=meshgrid(ti,tu)。 y=data(81:146,3)。)。 tu=0:10:18397。) %3區(qū) As分布圖 x=data(81:146,2)。cubic39。 ti=0:10:22674。o39。 zi=griddata(x,y,z,xi,yi,39。 z=data(45:80,12)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 [xi,yi]=meshgrid(ti,tu)。 y=data(45:80,3)。)。 tu=0:10:17980。) 32 %2區(qū) Hg分布圖 x=data(45:80,2)。cubic39。 ti=0:10:22674。o39。 zi=griddata(x,y,z,xi,yi,39。 z=data(45:80,8)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 [xi,yi]=meshgrid(ti,tu)。 y=data(45:80,3)。)。 tu=0:10:17980。) %2區(qū) As分布圖 x=data(45:80,2)。cubic39。 ti=0:10:26453。o39。 zi=griddata(x,y,z,xi,yi,39。 z=data(1:44,12)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 [xi,yi]=meshgrid(ti,tu)。 y=data(1:44,3)。)。 tu=0:10:18449。) %1區(qū) Hg分布圖 x=data(1:44,2)。cubic39。 ti=0:10:26453。o39。 zi=griddata(x,y,z,xi,yi,39。 z=data(1:44,8)。 mesh(xi,yi,zi) 29 hold on plot3(x,y,z,39。 [xi,yi]=meshgrid(ti,tu)。 y=data(1:44,3)。)。 tu=0:10:18449。 八 、附錄 附錄一 ① 城區(qū)地形示意圖如下： 附錄二 ① 一區(qū)各元素的分布圖： 21 111 112 ② 二區(qū)各元素分布圖： 211 22 212 ③ 三區(qū)各元素分布圖： 311 312 ④ 四區(qū)各元素分布圖： 23 411 412 ⑤ 五區(qū)各元素分布圖： 511 24 512 附錄三 ① As 在城區(qū)情況分布圖： 121 ② Cd 在城區(qū)情況分布圖： 25 221 ③ Cr 在城區(qū)情況分布圖： 321 ④ Cu 在城區(qū)情況分布圖： 26 421 ⑤ Hg 在城區(qū)情況分布圖： 521 ⑥ Ni 在城區(qū)情況分布圖： 27 621 ⑦ Pb 在城區(qū)情況分布圖： 721 ⑧ Zn 在城區(qū)情況分布圖： 28 821 附錄四 %1區(qū) As分布圖 x=data(1:44,2)。 [7] 胡克林，張芬榮，呂貽忠，王茹，徐艷，北京市大興區(qū)土壤重金屬含量的空間分布特征，環(huán)境科學學報，第 24 卷第 3 期： 464～ 468,2020 年 4 月。