8045 Cr 2383 4592 3299 10685 Cu 2383 3299 3052 3692 4603 6018 5528 3692 6018 2 4 Zn 4153 2383 3299 4948 9328 13797 4153 12696 2299 3692 6018 7293 4311 9621 2299 3024 7 附录5：土壤环境质量标准（GB15618-1995）
源程序：
A=xlsread('F:\\数模\\11A\\testdata.xls') for i=1:319; for j=1:8;
if A(i,j)>B(j,1); A(i,j)=B(j,1); end
if A(i,j) A=xlsread('F:\\数模\\11A\\testdata.xls') i=1; j1=1,j2=1,j3=1,j4=1,j5=1; while (i<320) if A(i,9)==1 for k=1:9 D1(j1,k)=A(i,k); end j1=j1+1; end if A(i,9)==2 for k=1:9 D2(j2,k)=A(i,k); end j2=j2+1; end if A(i,9)==3 for k=1:9 D3(j3,k)=A(i,k); end j3=j3+1; end if A(i,9)==4 for k=1:9 D4(j4,k)=A(i,k); end j4=j4+1; end if A(i,9)==5 for k=1:9 D5(j5,k)=A(i,k); end j5=j5+1; end i=i+1; end Variable x,y; Parameters dx,dy,ux,uy; Function z=(8*(pi*1)^(3/2))exp(-(x-ux)/(4*dx)-(y-uy)/(4*dy)); Data; //x y 1 15 2 13 3 11 4 9 5 8 6 6 7 3 8 1 A=xlsread('F:\\数模\\11A\\testdata.xls') j1=1,j2=1,j3=1,j4=1,j5=1,j6=1,j7=1,j8=1; for i=1:319; if A(i,1)>14.2 Z1(j1,1)=A(i,9); Z1(j1,2)=A(i,1); Z1(j1,3)=A(i,10); Z1(j1,4)=A(i,11); j1=j1+1; end if A(i,2)>600; Z2(j2,1)=A(i,9); Z2(j2,2)=A(i,2); Z2(j2,3)=A(i,10); Z2(j2,4)=A(i,11); j2=j2+1; end if A(i,3)>162; Z3(j3,1)=A(i,9); Z3(j3,2)=A(i,3); Z3(j3,3)=A(i,10); Z3(j3,4)=A(i,11); j3=j3+1; end if A(i,4)>57.6; Z4(j4,1)=A(i,9); Z4(j4,2)=A(i,4); Z4(j4,3)=A(i,10); Z4(j4,4)=A(i,11); j4=j4+1; end if A(i,5)>162; Z5(j5,1)=A(i,9); Z5(j5,2)=A(i,5); Z5(j5,3)=A(i,10); Z5(j5,4)=A(i,11); j5=j5+1; end if A(i,6)>51; Z6(j6,1)=A(i,9); Z6(j6,2)=A(i,6); Z6(j6,3)=A(i,10); Z6(j6,4)=A(i,11); j6=j6+1; end if A(i,7)>150; Z7(j7,1)=A(i,9); Z7(j7,2)=A(i,7); Z7(j7,3)=A(i,10); Z7(j7,4)=A(i,11); j7=j7+1; end if A(i,8)>330; Z8(j8,1)=A(i,9); Z8(j8,2)=A(i,8); Z8(j8,3)=A(i,10); Z8(j8,4)=A(i,11); j8=j8+1; end end k1=j1; for m=1:j1-1; for i=1:319; P=[Z1(m,3),Z1(m,4)]; Q=[A(i,10),A(i,11)] if norm(P-Q)<=2000&&norm(P-Q)>1; Z1(k1,1)=A(i,9); Z1(k1,2)=A(i,1); Z1(k1,3)=A(i,10); Z1(k1,4)=A(i,11); k1=k1+1; end end end