附录附录一：程序1．城区各功能区的地形二维等高线污染浓度散点图分布图程序clcclearcloseallD=load('坐标数据．dat');M=load('浓度含量shuju．dat');x=D(:,2);y=D(:,3);z=D(:,4);c=D(:,5);xi=linspace(min(x),max(x),300);yi=linspace(min(y),max(y),300);[xi,yi]=meshgrid(xi,yi);zi=griddata(x,y,z,xi,yi);ci=griddata(x,y,c,xi,yi);marker={'d','^','s','o','V'};color={'b','y','c','r','m'};mat={'As','Cd','Cr','Cu','Hg','Ni','Pb','Zn'};str={'等高线','生活区','工业区','山林区','交通区','绿地区'};forj=1:8%等高线图形figurecontourf(xi,yi,zi,0:10:500);%set(h,'ShowText','on','TextStep',get(h,'LevelStep')*2)%clabel(C,h,[0:10:50,50:50:300,300:100:500]);legend('生活区','工业区','山林区','交通区','绿地区');%legend({'k','r','y','c','b'});title(['重金属',mat{j},'空间二维等高线分布图'])xlabel('X---轴')ylabel('Y---轴')colormapsummercolorbargridonholdonfori=1:5loc=c==i;plot(x(loc),y(loc),marker{i},'markerfacecolor',color{i},'MarkerEdgeColor',color{i});end%legend(str,'location','best')fork=1:length(x)text(x(k)-200,y(k)+200,num2str(M(k,j+1)),'fontsize',8);endlegend('功能一区——生活区','功能二区——工业区','功能三区——山林区','功能四区——交通区','功能五区——绿地区');end2．重金属‘**’含量程序clc;clear;loadshuju1．txtA=shuju1;x=A(:,1);y=A(:,2);z=A(:,3);scatter(x,y,5,z)%散点图figure(1)[X,Y,Z]=griddata(x,y,z,linspace(min(x),max(x),200)',linspace(min(y),max(y),200),'v4');%插值pcolor(X,Y,Z);shadinginterp%伪彩色图figure(2),contourf(X,Y,Z)%等高线图legend('污染程度分布');title('重金属As含量空间分布平面图')xlabel('x—轴'),ylabel('y—轴');%zlabel('z—轴')figure(3),surf(X,Y,Z)%三维曲面legend('污染程度分布');title('重金属As含量空间分布平面图')xlabel('x—轴'),ylabel('y—轴');zlabel('z—轴')3重金属综合污染指数程序%Cd综合污染指数Grandvalue=190;maxvalue=0;sum=0;C=[352．60556．20113．70240．50120．30514．10246．30158．00211．30236．50];%可以输入任意长度的污染浓度数据C=C'n=length(C);fori=1:n;ifC(i)>maxvaluemaxvalue=C(i);elsesum=sum+C(i)/Grandvalue;endendp=sqrt(((sum/n)^2+(maxvalue/Grandvalue)^2)/2)附录二：污染状况空间分布图
