附件
第一个程序
（城市地势图） 第一步：标准化
data_data=zeros (319,8);% 标准化后得到的数据 for i=1:8
data_data (:,i)=(data (:,i+5)-beijing (i,1)*ones (319,1))./beijing(i,2); end
第二步：引值 x=data (:,2); y=data (:,3); z=data (:,4); 第三步：插值
[X,Y,Z]=griddata (x,y,z,linspace (min (x),max (x),100)',linspace (min (y),max (y),100),'v4');%插值 surf (X,Y,Z)
第四步：命名x,y,z轴，
grid,xlabel ('x轴'),ylabel ('y轴'),zlabel ('海拔'),title (' 城市地势图') % axis ([min (x),max (x),min (y),max (y),min (z),max (z)]) 各种金属浓度分布图 第一步：标准化
data_data=zeros (319,8);% 标准化后得到的数据 for i=1:8
data_data (:,i)=(data (:,i+5)-beijing (i,1)*ones (319,1))./beijing(i,2); end
第二步：引值 x=data (:,2); y=data (:,3); z=data (:,6); 第三步：插值
[X,Y,Z]=griddata (x,y,z,linspace (min (x),max (x),100)',linspace (min (y),max (y),100),'v4');%插值 surf (X,Y,Z)
第四步：命名x,y,z轴，
grid,xlabel ('x轴'),ylabel ('y轴'),zlabel ('各种重金属的浓度分布图'),title (' 各种重金属的浓度分布图')
% axis ([min (x),max (x),min (y),max (y),min (z),max (z)]) 第二个程序
第一步：标准化处理后的数据
data_data=zeros(319,8);%标准化处理后的数据 for i=1:8
data_data(:,i)=(data(:,i+5)-beijing(i,1)*ones(319,1))./beijing(i,2);
26
end
第二步：定义X X=[data_data];
第三步：运用MATLAB中的corrcoef函数算出相关系数 r=corrcoef(X)
第四步：运行软件得到所求结果为：
r =
1.0000 0.2547 0.1890 0.1597 0.0644 0.3166 0.2899 0.2469 0.2547 1.0000 0.3524 0.3967 0.2647 0.3294 0.6603 0.4312 0.1890 0.3524 1.0000 0.5316 0.1032 0.7158 0.3828 0.4243 0.1597 0.3967 0.5316 1.0000 0.4167 0.4946 0.5200 0.3873 0.0644 0.2647 0.1032 0.4167 1.0000 0.1029 0.2981 0.1958 0.3166 0.3294 0.7158 0.4946 0.1029 1.0000 0.3068 0.4364 0.2899 0.6603 0.3828 0.5200 0.2981 0.3068 1.0000 0.4937 0.2469 0.4312 0.4243 0.3873 0.1958 0.4364 0.4937 1.0000
第三个程序
第一步：标准化
data_data=zeros (319,8); for i=1:8
data_data (:,i)=(data (:,i+3)-beijing (i,1)*ones (319,1))/beijing(i,2); end
第二步：定义相应的区域，并用插值法得出模型
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 55
stander_four _data=zeros (319,4);
stander_four _data (:,1)=data_data (:,2)+data_data (:,7)+data_data (:,8); % pb,cd ,zn stander_four _data (:,2)=data_data (:,3)+data_data (:,4)+data_data (:,6);% cr,cu,ni stander_four _data (:,3)=data_data (:,5);% hg stander_four _data (:,4)=data_data (:,1);% as
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
x=data (:,1); y=data (:,2);
z=stander_four _data (:,1);
[X,Y,Z]=griddata (x,y,z,linspace (min (x),max (x),200)',linspace (min (y),max (y),200),'v4'); surf (X,Y,Z)
27
标准化后的各个重金属元素的浓度值部分表 As (μg/g) Cd (ng/g) Cr (μg/g) Cu (μg/g) Hg (ng/g) Ni (μg/g) Pb (μg/g) Zn (μg/g) 4.7111 0.79333 1.4789 2.0444 28.875 1.5526 0.73 0.23929 2.5889 0.54 1.5611 2.5861 6.375 1.2895 0.86333 1.8279 1.4444 10.307 -0.21444 14.267 9.25 -0.44737 7.22 10.669 3.2889 3.13 1.0089 3.325 114.38 0.81579 0.21333 3.4536 3.0556 13.173 3.15 28.981 95.625 2.0789 23.16 46.93 11.644 32.097 4.1067 82.058 125.63 4.1842 67.3 64.124 5.9333 4.66 7.2033 8.7806 10.75 1.4474 5.3183 6.9807 6.6889 31.207 28.287 698.69 1683.1 7.7368 58.44 96.347 4.2333 33.13 6.3522 38.456 1995.6 3.5526 23.56 61.274 5.6889 4.57 3.84 4.5694 3.5 2.4737 0.99 2.2436 2.5889 2.38 1.5767 3.25 28 0.60526 0.81333 2.4036 6.1889 5.2333 1.4378 9.0472 16.625 1.9474 5.29 11.011 2.3556 2.1233 5.4833 3.7139 9.5 1.9737 4.44 1.4343 0.98889 7.65 25.239 30.575 5.25 0.15789 12.578 56.07 2.1222 12.88 6.7744 21.067 19.25 1.9737 15.12 30.414 8.7222 30.483 7.0867 34.381 20.875 2.6316 73.58 38.074 2.8222 10.513 5.7433 42.831 13.625 1.6053 13.373 22.914 4.7111 7.2633 2.9611 23.317 22.25 1.9211 6.6367 17.017 4.2333 7.19 14.272 16.064 6.25 1.5263 9.7233 22.28 5.4444 16.133 79.273 32.597 15.125 5.3947 32.94 67.462 2.1222 4.24 2.6267 4.3917 8.625 0.23684 9.4467 11.019 6.9333 36.117 98.871 375.46 10 34.263 25.08 124.96 6.4333 6.5267 15.699 25.469 5.875 5.0526 9.9833 25.735 -0.33333 2.7367 2.1256 7.0611 13 -0.44737 5.9967 8.3729 0.54444 -1.3167 0.44667 -0.38333 -2.375 -0.5 -0.31833 -1.5829 2.8222 15.113 7.1389 61.25 15 2.3158 11.078 17.287 1.9 7.88 1.26 14.292 19.125 1.3421 6.0183 8.1179 0.1 6.4633 0.46 5.9611 1.875 0.42105 5.7467 13.154 20.3 9.8167 4.7111 12.922 185.63 4.0789 8.7833 7.6221 16.422 16.667 7.2978 28.225 76.25 5.9211 16.56 111.22 7.7 16.843 3.67 24.486 19.375 4.2105 21.94 39.007 -0.11111 11.113 9.0211 16.592 10.375 0.47368 4.9333 8.8836 3.0556 13.4 2.9456 19.6 19.5 1.8947 7.0767 16.296 2.1222 21.623 4.8511 22.022 36.875 1.9474 13.2 20.188 0.98889 20.827 7.6533 23.533 26 3.3158 11.632 18.169 -0.11111 8.8767 11.93 12.714 16.875 3.1316 10.127 201.75 2.1222 18.593 6.0578 16.569 20.75 1.7632 12.033 23.876 0.76667 13.2 2.7011 18.953 7.25 2 11.608 21.486 0.76667 10.637 4.0244 10.678 35 0.89474 12.662 16.121 3.2889 24.09 4.6211 40.408 34.5 2.3421 15.54 22.01
28
14.422 10.967 7.0878 9.4361 233.13 2 4.2333 6.91 5.1411 65.458 6.875 2.0526 2.5889 14.603 4.9044 29.156 12.5 3.0263 1.2111 15.633 4.2278 30.272 10.75 1.9737 1.4444 16.85 4.1578 59.6 17.625 1.8947 1.9 15.69 1.5167 8.8611 2 0.84211 0.76667 14.587 3.25 9.8528 1.375 1 2.1222 3.2833 2.03 4.9028 3.375 2.7895 1.2111 14.62 30.558 15.892 108.13 1.1842 4 2.8233 2.1478 7.4889 15.125 2.2368 1.9 0.73 1.8044 3.1222 13.125 1.3158 1.4444 7.11 1.2233 12.669 5.625 0.39474 1.4444 5.4633 3.2544 10.508 2.25 0.078947 0.32222 6.0967 0.31 73.506 2.5 0.44737 0.1 6.1967 1.6033 5.7917 2.5 0.078947 -0.54444 9.5433 2.9867 9.5667 16.5 -0.10526 0.32222 8.1333 1.5744 10.275 0 0.71053 0.32222 7.1433 0.47667 9.6306 8.125 0.73684 -0.76667 4.0967 1.6644 16.206 -0.3575 0.55263 -0.33333 12.447 0.86 4.7944 0.17875 -1.3421 1.4444 5.7933 2.78 14.628 3.5263 7.3026 0.54444 -0.1 -0.38 2.9972 -0.625 -0.09737 -0.76667 4.5 0.51778 4.4972 -1.295 -0.80789 -0.97778 4.9633 1.3811 5.3917 3.6613 -0.63158 -0.54444 20.707 2.4567 11.278 -0.8925 -0.48421 -0.33333 7.71 1.8378 10.856 -1.1612 -0.83947 -0.33333 11.933 2.2422 5.9306 0.3125 -0.39474 2.8222 3.2333 1.2389 14.956 1.7862 1.0553 0.1 7.2467 0.75111 1.8806 -1.0263 -0.60263 0.98889 0.2 0.61778 2.7417 -1.6963 0.7 -0.54444 6.57 -0.55778 2.925 -1.1612 -0.66053 4.9556 -0.56333 2.3778 2.1139 3.795 1.9447 0.1 4.6833 0.23556 12.403 -1.1612 0.31579 0.1 1.01 0.81 3.5778 -1.1612 0.52368 -0.11111 5.85 0.94444 4.9056 -0.49125 0.64211 -0.97778 -1.97 -1.2833 -1.1333 -2.5 -1.3711 5.4444 4.3767 1.6633 3.2306 2.0538 1.5 -2.0333 -0.34 -1.7422 0.038889 -3.3038 -1.6079 -1.1889 11.293 0.67111 5.2306 1.25 -0.03947 0.1 0.69 3.1789 1.6111 -0.08875 3.3368 2.8222 -0.99 0.72111 1.9528 -0.08875 0.67105 8.2111 -0.67333 2.7856 15.517 2.9913 2.2105 3.0556 -1.2733 0.56889 1.0306 0.17875 0.05
29
6.72 11.43 13.423 11.953 10.888 5.8117 5.4567 2.4883 8.1117 2.7933 1.0817 8.0117 7.4883 6.2067 5.3067 14.198 4.5183 17.12 1.9983 3.7883 9.9483 4.4117 4.9233 3.7333 4.8783 13.743 3.9367 0.54667 3.9017 0.53167 5.4683 -0.24 5.1733 1.5217 5.1483 -1.3783 2.065 -0.78167 5.7083 1.245 -0.285 1.3117 -0.665 10.436 10.071 35.997 25.574 33.57 8.3807 9.9329 2.36 9.1236 23.927 1.6821 14.773 14.955 16.186 9.095 12.36 6.3107 5.1443 5.5157 2.4186 263.7 1.1864 3.8543 4.765 6.1429 10.662 4.0657 0.99714 10.873 0.71286 4.9329 -0.48286 3.5114 0.99 5.7057 -2.5814 1.1064 -1.5171 7.8557 1.4343 -0.50429 -0.40929 -1.5171 29.478 0.32222 -0.76667 -1.8222 3.7556 -0.76667 2.5889 2.5889 4.2333 0.76667 1.2111 2.3556 3.5222 0.98889 -0.54444 -0.76667 0.76667 4.4778 2.5889 -1.4 -0.76667 2.3556 -1.4 3.2889 1.2111 3.0556 1.6667 1.2111 -0.11111 1.2111 0.98889 -0.33333 -0.76667 0.54444 2.3556 1.4444 1.4444 0.54444 2.5889 -0.76667 -0.97778 20.463 9.56 7.9933 2.1333 -2.6633 2.2933 25.22 -0.03667 -0.52333 3.43 0.093333 49.66 5.0833 1.6867 8.5533 7.1667 -1.1467 -1.43 2.4333 7.4333 3.44 1.4767 -1.4133 3.8533 1.2633 -0.62 -1.18 -0.61333 -1.4833 1.3167 0.26667 0.046667 -2.9667 -0.01 0.6 -0.73 0.083333 1.0967 -1.3733 2.5567 -0.27333 2.0033 0.67111 0.59333 -1.4356 1.1133 -0.31 1.2989 1.8356 1.9267 -0.20333 0.56778 1.3867 1.2189 0.69222 0.54778 1.0511 -0.94556 1.6478 0.55222 -0.71889 -0.67556 0.31444 -1.3933 0.63667 0.23889 -0.24222 5.2133 -0.71444 -0.51889 0.90111 0.39111 0.55222 1.1967 -0.46333 0.63667 -0.92 5.3911 0.31 2.3789 -1.5044 -1.2244 3.6056 2.9389 10.911 0.825 0.29167 1.6778 4.1444 0.86389 2.2917 7.45 1.1333 0.63889 11 1.5278 3.6611 38.614 1.0972 0.45278 1.0222 -0.13889 2.3389 4.5139 -0.96389 13.361 1.6 12.769 1.9833 -0.30833 3.2444 2.5861 6.2056 -0.53056 -0.19444 -0.85556 0.28056 -1.0556 1.5194 -0.06667 0.58333 -0.32778 -1.525 -0.8925 -1.5625 -1.5625 -1.0263 -2.2325 -2.6338 1.1163 -2.7675 0.045 209.91 -1.83 -1.83 -0.8925 -0.8925 -1.295 -1.4288 6.3387 -0.625 -2.5 -2.9013 6.3387 -1.83 -2.6338 2.5888 3.125 21.474 -1.4288 -0.49125 0.17875 -1.1612 1.92 1.1163 -2.2325 0.71375 2.1875 0.045 0.9825 -0.625 -3.0362 -0.49125 -1.295 1.4105 -0.21579 -1.1342 -1.4 0.28684 -0.72105 2.3289 1.4105 1.1158 -0.66053 0.85 0.81842 1.5895 -0.09737 -0.01053 0.85 -0.57105 0.61316 0.58158 -0.86842 -1.0447 0.19737 -0.89737 0.49474 0.96842 -0.15789 4.3132 -0.98684 -0.24737 0.96842 0.078947 -0.33421 1.2632 -1.1947 0.078947 -0.89737 4.0474 0.16842 2.0026 -1.7868 -1.1053 6.96 3.9083 1.8367 1.545 -0.73833 2.8017 10.607 0.41833 0.805 1.275 1.005 -0.00167 3.0167 0.90833 4.8333 4.5083 -0.505 -0.29167 2.4717 4.6333 2.3667 1.85 -1.095 2.7067 0.41 -0.49667 -0.84667 -0.64667 -0.705 1.4217 -0.925 -0.15667 -1.8867 -0.36 -0.64167 0.23167 -0.45 -0.305 -0.39667 3.7983 -0.61667 8.0743 2.5714 1.0771 1.835 -2.0057 1.8429 6.7336 1.6236 -0.36571 1.8864 0.98286 -0.81786 2.5643 -0.35857 6.3029 7.265 -0.12571 -1.4443 1.4636 1.4343 2.2214 0.17357 -1.8307 3.235 -0.94929 -0.53357 -0.75214 -1.7943 -0.35143 0.99714 -1.1821 -0.50429 0.66929 -0.06 -0.79571 0.027857 1.66 0.42143 -0.36571 0.72714 -0.78857
30