【正文】 o39。)。 zi=griddata(x,y,z,xi,yi,39。 tu=0:10:18449。 z=data(1:319,4)。) %城區(qū)地形圖 x=data(1:319,2)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。cubic39。 40 [xi,yi]=meshgrid(ti,tu)。 ti=0:10:20582。 y=data(285:319,3)。o39。)。 zi=griddata(x,y,z,xi,yi,39。 tu=0:10:17133。 z=data(285:319,12)。) %5區(qū) Pb分布圖 x=data(285:319,2)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。cubic39。 [xi,yi]=meshgrid(ti,tu)。 ti=0:10:20582。 y=data(285:319,3)。o39。)。 zi=griddata(x,y,z,xi,yi,39。 tu=0:10:17133。 z=data(285:319,10)。) %5區(qū) Hg分布圖 39 x=data(285:319,2)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。cubic39。 [xi,yi]=meshgrid(ti,tu)。 ti=0:10:20582。 y=data(285:319,3)。o39。)。 zi=griddata(x,y,z,xi,yi,39。 tu=0:10:17133。 z=data(285:319,8)。) %5區(qū) Cr分布圖 x=data(285:319,2)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。cubic39。 [xi,yi]=meshgrid(ti,tu)。 ti=0:10:20582。 y=data(285:319,3)。o39。)。 zi=griddata(x,y,z,xi,yi,39。 tu=0:10:17133。 z=data(285:319,6)。) %5區(qū) As分布圖 x=data(285:319,2)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。cubic39。 [xi,yi]=meshgrid(ti,tu)。 ti=0:10:26424。 y=data(147:284,3)。o39。)。 zi=griddata(x,y,z,xi,yi,39。 tu=0:10:17949。 z=data(147:284,12)。) %4區(qū) Pb分布圖 x=data(147:284,2)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。cubic39。 [xi,yi]=meshgrid(ti,tu)。 37 ti=0:10:26424。 y=data(147:284,3)。o39。)。 zi=griddata(x,y,z,xi,yi,39。 tu=0:10:17949。 z=data(147:284,10)。) %4區(qū) Hg分布圖 x=data(147:284,2)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。cubic39。 [xi,yi]=meshgrid(ti,tu)。 ti=0:10:26424。 y=data(147:284,3)。o39。)。 zi=griddata(x,y,z,xi,yi,39。 tu=0:10:17949。 z=data(147:284,8)。) 36 %4區(qū) Cr分布圖 x=data(147:284,2)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。cubic39。 [xi,yi]=meshgrid(ti,tu)。 ti=0:10:26424。 y=data(147:284,3)。o39。)。 zi=griddata(x,y,z,xi,yi,39。 tu=0:10:17949。 z=data(147:284,6)。) %4區(qū) As分布圖 x=data(147:284,2)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。cubic39。 [xi,yi]=meshgrid(ti,tu)。 ti=0:10:28654。 y=data(81:146,3)。o39。)。 zi=griddata(x,y,z,xi,yi,39。 35 tu=0:10:18397。 z=data(81:146,12)。) %3區(qū) Pb分布圖 x=data(81:146,2)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。cubic39。 [xi,yi]=meshgrid(ti,tu)。 ti=0:10:28654。 y=data(81:146,3)。o39。)。 zi=griddata(x,y,z,xi,yi,39。 tu=0:10:18397。 z=data(81:146,10)。) %3區(qū) Hg分布圖 x=data(81:146,2)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。cubic39。 [xi,yi]=meshgrid(ti,tu)。 ti=0:10:28654。 y=data(81:146,3)。o39。)。 zi=griddata(x,y,z,xi,yi,39。 tu=0:10:18397。 z=data(81:146,8)。) %3區(qū) Cr分布圖 x=data(81:146,2)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。cubic39。 [xi,yi]=meshgrid(ti,tu)。 ti=0:10:28654。 y=data(81:146,3)。o39。)。 zi=griddata(x,y,z,xi,yi,39。 tu=0:10:18397。 z=data(81:146,6)。) %3區(qū) As分布圖 x=data(81:146,2)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。cubic39。 [xi,yi]=meshgrid(ti,tu)。 ti=0:10:22674。 y=data(45:80,3)。o39。)。 zi=griddata(x,y,z,xi,yi,39。 tu=0:10:17980。 z=data(45:80,12)。) %2區(qū) Pb分布圖 x=data(45:80,2)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。cubic39。 [xi,yi]=meshgrid(ti,tu)。 ti=0:10:22674。 y=data(45:80,3)。o39。)。 zi=griddata(x,y,z,xi,yi,39。 tu=0:10:17980。 z=data(45:80,10)。) 32 %2區(qū) Hg分布圖 x=data(45:80,2)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。cubic39。 [xi,yi]=meshgrid(ti,tu)。 ti=0:10:22674。 y=data(45:80,3)。o39。)。 zi=griddata(x,y,z,xi,yi,39。 tu=0:10:17980。 z=data(45:80,8)。) %2區(qū) Cr分布圖 x=data(45:80,2)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。cubic39。 [xi,yi]=meshgrid(ti,tu)。 ti=0:10:22674。 y=data(45:80,3)。o39。)。 zi=griddata(x,y,z,xi,yi,39。 tu=0:10:17980。 z=data(45:80,6)。) %2區(qū) As分布圖 x=data(45:80,2)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。cubic39。 [xi,yi]=meshgrid(ti,tu)。 ti=0:10:26453。 y=data(1:44,3)。o39。)。 zi=griddata(x,y,z,xi,yi,39。 tu=0:10:18449。 z=data(1:44,12)。) %1區(qū) Pb分布圖 x=data(1:44,2)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。cubic39。 [xi,yi]=meshgrid(ti,tu)。 ti=0:10:26453。 y=data(1:44,3)。o39。)。 zi=griddata(x,y,z,xi,yi,39。 tu=0:10:18449。 z=data(1:44,10)。) %1區(qū) Hg分布圖 x=data(1:44,2)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。cubic39。 [xi,yi]=meshgrid(ti,tu)。 ti=0:10:26453。 y=data(1:44,3)。o39。)。 zi=griddata(x,y,z,xi,yi,39。 tu=0:10:18449。 z=data(1:44,8)。) %1區(qū) Cr分布圖 x=data(1:44,2)。 mesh(xi,yi,zi) 29 hold on plot3(x,y,z,39。cubic39。 [xi,yi]=meshgrid(ti,tu)。 ti=0:10:26453。 y=data(1:44,3)。o39。)。 zi=griddata(x,y,z,xi,yi,39。 tu=0:10:18449。 z=data(1:44,6)。 八 、附錄 附錄一 ① 城區(qū)地形示意圖如下： 附錄二 ① 一區(qū)各元素的分布圖： 21 111