发布时间 : 星期三 文章数学建模模板 - 图文更新完毕开始阅读2fb08b7376232f60ddccda38376baf1ffc4fe3aa
八、参考文献
[1] 曲向荣.土壤环境学.北京：清华大学出版社,2010
[2] 陶玲，任珺，乔国栋等.兰州市蔬菜基地土壤重金属含量及评价.兰州交通大学环境与市政工程学院环境生态研究所，2010年06期
[3] 国家环境保护局.土壤环境质量标准(GB15618－1995).www.fert.cn/1017/2005/9/26/20059261164099229.shtml,2011.9.11
[4] 百度百科.重金属污染.baike.http://m.china-audit.com//view/282970.htm，2011.9.10 [5] 杨希东.实验数据异常值的剔除方法.唐山师专学报，1998.5
[6] 钱天伟，刘春国.饱和-非饱和土壤污染物运移.北京：中国环境科学出版社，2007 [7] 黄铭洪等.环境污染与生态恢复.北京：科学出版社，2003
[8] 彭应登.区域开发环境影响评价.北京：中国环境科学出版社，1999
[9] 何峰.重庆市农田土壤-粮食作物重金属关联特征与污染评价.西南农业大学，中国学术期刊网络出版总库，2004年03期
[10] 褚卓栋.土壤-中草药重金属含量及中药中砷汞生物可给性研究.河北农业大学，中国学术期刊网络出版总库，2008
[11] 李胤，毛义伟，周立晨等.上海世博会规划区域城市土壤重金属调查及评价研究.上海市城市化生态过程与生态恢复重点实验室，2009年04期
[12] 周爱国，周建伟，梁合诚等.地质环境评价.武汉：中国地质大学出版社，2008 [13] 张丹.贵州主要药材基地土壤及中药材重金属污染状况调查研究.中国科学院上海冶金研究所，2006
[14] 毛文永.生态环境影响评价概论.北京：中国环境科学出版社，1998 [15] 张志涌等.精通MATLAB 6.5版.北京：北京航空航天大学出版社，2002
[16] 欧阳志云，郑华，高吉喜等.区域生态环境质量评价与生态功能区划.北京：中国环境科学出版社，2009
[17] 邹乔，王瑶，杜显元等.内梅罗指数法在石油开采区土壤PAHs污染分级评价中的应用.安徽农业科学，2011
23
九、附录
9.1 程序
9.1.1各元素在该城区内的浓度平面分布图的程序（以As为例）
%各元素在该城区内的浓度平面分布 clc close all clear
den_As=xlsread('x-y-z-w.xls'); x=den_As(:,1); y=den_As(:,2);
Ndu=xlsread('wuranwu.xls','A1:A319'); nx=500:1000:29500; ny=500:1000:19500; Zmean=nan(29,19); Zcount=nan(29,19); for ii=nx for jj=ny
tem=Ndu( ii-500<=x & x pcolor(nx,ny,Zmean') xlabel('x坐标(m)') ylabel('y坐标(m)') colorbar colormap(cool) 9.1.2各重金属元素的三维分布图的程序（以As为例） %各元素的浓度分布图---三维图 clc close all clear a=xlsread('x-y-z-w.xls'); x=a(1:319,1); y=a(1:319,2); z=a(1:319,3); b=xlsread('wuranwu.xls'); z_as=b(1:319,1); nx=linspace(min(x),max(x),100); ny=linspace(min(y),max(y),100); [xx,yy]=meshgrid(nx,ny); zz=griddata(x,y,z,xx,yy); zz_as=griddata(x,y,z_as,xx,yy); 24 surf(xx,yy,zz_as); colorbar('horiz'); shading interp; colormap(cool); title('As含量'); 9.1.3各重金属元素的空间四维浓度分布图（以As为例） clc close all clear a=xlsread('x-y-z-w.xls'); x=a(1:319,1); y=a(1:319,2); z=a(1:319,3); b=xlsread('wuranwu.xls'); z_as=b(1:319,1); nx=linspace(min(x),max(x),100); ny=linspace(min(y),max(y),100); [xx,yy]=meshgrid(nx,ny); zz=griddata(x,y,z,xx,yy); zz_as=griddata(x,y,z_as,xx,yy); surf(xx,yy,zz,zz_as); colorbar('horiz'); shading interp; title('As含量'); 9.1.4 计算梯度的程序 9.1.4.1 海拔梯度 %输入海拔矩阵a，求梯度（八个方向的平均值） clc close all clear c=xlsread('x-y-z-w.xls'); x=c(1:319,1); y=c(1:319,2); aa=c(1:319,3); nx=linspace(min(x),max(x),100); ny=linspace(min(y),max(y),100); [xx,yy]=meshgrid(nx,ny); a=griddata(x,y,aa,xx,yy); b=nan(100,100); for i=2:99 for j=2:99 sur=[a(i-1,j-1) a(i-1,j+1) a(i+1,j-1) a(i+1,j+1) a(i-1,j) a(i,j-1) a(i,j+1) a(i+1,j)]; blo=nan(1,8); for k=1:4 25 blo(k)=abs((sur(k)-a(i,j)))/100/sqrt(13); end for k=[5 8] blo(k)=abs((sur(k)-a(i,j)))/200; end for k=[6 7] blo(k)=abs((sur(k)-a(i,j)))/300; end if sum(~isnan(blo))==0 b(i,j)=nan; else b(i,j)=max(blo); end end end surf(xx,yy,b); shading interp; title('海拔梯度'); 9.1.4.2 浓度梯度（以As为例） %输入浓度矩阵a，求梯度（八个方向的平均值） clc close all clear c=xlsread('x-y-z-w.xls'); x=c(1:319,1); y=c(1:319,2); z=c(1:319,3); d=xlsread('wuranwu.xls'); aa=d(1:319,1); nx=linspace(min(x),max(x),100); ny=linspace(min(y),max(y),100); [xx,yy]=meshgrid(nx,ny); a=griddata(x,y,aa,xx,yy); b=nan(100,100); for i=2:99 for j=2:99 sur=[a(i-1,j-1) a(i-1,j+1) a(i+1,j-1) a(i+1,j+1) a(i-1,j) a(i,j-1) a(i,j+1) a(i+1,j)]; blo=nan(1,8); for k=1:4 blo(k)=abs((sur(k)-a(i,j)))/100/sqrt(13); end for k=[5 8] blo(k)=abs((sur(k)-a(i,j)))/200; end 26