H第四次论文研读报告嫌犯围捕【DOC精选】.doc

城市表层土壤重金属污染分析（特等奖） 主要方法：第一问：地累积指数；克里金插值；ARCGIS软件。 图1各种区的分布图 代码：clear clear,clc A=xlsread(地点及功能区.xls); s1=A([1:44],:) s2=A([45:80],:) s3=A([81:146],:) s4=A([147:284],:) s5=A([285:319],:) plot(s1(:,1),s1(:,2),o,s2(:,1),s2(:,2),r*,s3(:,1),s3(:,2),gs,s4(:,1),s4(:,2),k+,s5(:,1),s5(:,2),m.) hold on legend(生活区,工业区,山区,交通区,公园区) hold on xi=linspace(min(A(:,1)),max(A(:,1)),500); yi=linspace(min(A(:,2)),max(A(:,2)),500); x=(A(:,1)); y=(A(:,2)); z=(A(:,3)); [X,Y,Z]=griddata(x,y,z,xi,yi,cubic) contour(X,Y,Z) 第一问： 4.3 用污染负荷指数的分析各功能区的污染程度 地积累指数分析了每种重金属的污染分布情况，为了评定不同功能区的重金属的综合污染程度，定义最高污染物指数： （4.2） 式中，为第i种金属浓度的实测值，第i种金属的背景值。 如果某点由多种重金属共同污染，该点的污染负荷指数为： （4.3） 为单个点的污染负荷指数，n为金属污染物种数。 由式（4.3），对有单点组成的区域来说，污染负荷指数为： （4.4） 为区域的污染负荷指数，m为采样点个数。 度如表4-3。 表4-3 各功能区污染程度 功能区 生活区 工业区 山区 主干道区 公园绿地区 污染负荷指数 1.83 2.35 1.06 1.92 1.58 污染程度 中等污染 强污染 中等污染 中等污染 中等污染 表格还原代码： %特等 表三数据的复现 clc; clear; A=load(金属浓度实测值.txt); B=load(地点及功能区.txt); B_GONG=B(:,5);%功能 BJ0=[3.6,130,31,13.2,35,12.3,31,69];%背景均值 SC=A(:,2:end);%浓度实测值 BJ=[]; for i=1:size(SC,1)%SC矩阵中行的长度 BJ=[BJ;BJ0];%将BJ变成一个和SC相同的矩阵 end %for i=1:319 % BJ=[BJ;BJ0];%每次在BJ0的基础上加一行相同得值 %end C=SC./BJ;%点除：矩阵对应的位置相除 PLI=prod(C,2).^(1/8);%prod连乘，2：列。(每一列的数据连乘再开8次方)%（C,1）：1表示行 PLI0=[PLI,B_GONG];%加入了功能区标识 flag=[]; for i=1:size(SC,1)%for i=1:319 if PLI0(i,2)==1 %查找功能区1 flag=[flag;i];%其中；表示分了两列 end end PLI1=PLI0(flag,1); m=size(PLI1,1);%计算功能区的样本点个数 PLI1=(prod(PLI1))^(1/m) flag=[]; for i=1:size(SC,1) if PLI0(i,2)==2 %查找功能区2 flag=[flag;i]; end end PLI2=PLI0(flag,1); m=size(PLI2,1);%计算功能区的样本点个数 PLI2=(prod(PLI2))^(1/m) flag=[]; for i=1:size(SC,1) if PLI0(i,2)==3 %查找功能区3 flag=[flag;i]; end end PLI3=PLI0(flag,1); m=size(PLI3,1);%计算功能区的样本点个数 PLI3=(prod(PLI3))^(1/m) %% flag=[]; for i=1:size(SC,1) if PLI0(i,2)==4 %查找功能区4 flag=[flag;i]; end end PLI4=PLI0(flag,1); m=size(PLI4,1);%计算功能区的样本点个数 PLI4=(prod(PLI4))^(1/m) %% flag=[]; for i=1:size(SC,1) if PLI0(i,2)==5 %查找功能