构建自适应滤波法模型：
?yt?1?w1yt?w2yt?1???wNyt?N?1??wiyt?i?1
i?1Nyt?1为第t+1期的预测值，wi 为第t?i?1观测值权数，yt?i?1为第其中?t?i?1期的观测值，N为权数的个数。
权数调整模型为：
wi'?wi?2k?et?1yt?i?1
2、3?N、t?N、N+1、?n， 其中i?1、n 为序列个数，wi为调整前第i个观察值的权重，wi'为调整后的第i个观察值的权重。k为学习常数，ei?1为第t?1期的预测误差。
模型的优缺点：
1、对于单项污染指标模型，该模型的优点利用样本点的平均值作为整个区域，避免因某些取样点的值偏高而导致评价不合理，并且可以比较准确的评价各个地区各个重金属的污染程度，具有一定的实用性。缺点是背景的平均值作为标准背景值，并没有考虑平均背景值的偏差，还有不能对一个地区的重金属的污染进行综合评价，具有一定的片面性，不能改一个地区给一个确切的评价结果。
2、对于模型内梅罗综合指数法模型其优点是充分考虑了单项污染最严重的重金属元素对该地区的重金属污染情况影响。同时充分弥补了单项污染指标模型的缺点，对一个地区的污染进行综合评价。但在一定程度上有他的片面性，因为它只是用到背景值的平均值，没有把一些测量等误差屏蔽到最小化。
3、对于等指标污染负荷模型，其优点充分的考虑扩散原理，比较符合实际的找出浓度最高的地点为污染源，但是该模型对于判断污染源的C值没有一个恒定，因为这个C值指标需要根据具体的污染源的污染程度进行适应性变化，无法找到确切的真正的污染源，只能找到极为可能的污染源位置。
4、对于自适应滤波法模型，我们能够很好的进行时间方面的预测预估计，很有实用性。但我们只考虑时间的变化，没有深刻的考虑城市的外界因素，具有一定的局限性。
五、模型推广与评价
模型的推广：
1、对于城市污染情况进行综合评价，可以作为环保部门对环保政策实施的重要依据。
2、对于污染扩散模型，可以作为对污染防治提供参考。
3、对于建立的时间序列-自适应滤波模型，我们可以对地质环境进行预测，也可以作为对未来预测值的可靠预测方案。
模型的评价：
17
1、对于单项污染指标模型，该模型的优点利用样本点的平均值作为整个区域，避免因某些取样点的值偏高而导致评价不合理，并且可以比较准确的评价各个地区各个重金属的污染程度，具有一定的实用性。缺点是背景的平均值作为标准背景值，并没有考虑平均背景值的偏差，还有不能对一个地区的重金属的污染进行综合评价，具有一定的片面性，不能改一个地区给一个确切的评价结果。
2、对于模型内梅罗综合指数法模型其优点是充分考虑了单项污染最严重的重金属元素对该地区的重金属污染情况影响。同时充分弥补了单项污染指标模型的缺点，对一个地区的污染进行综合评价。但在一定程度上有他的片面性，因为它只是用到背景值的平均值，没有把一些测量等误差屏蔽到最小化。
3、对于等指标污染负荷模型，其优点充分的考虑扩散原理，比较符合实际的找出浓度最高的地点为污染源，但是该模型对于判断污染源的C值没有一个恒定，因为这个C值指标需要根据具体的污染源的污染程度进行适应性变化，无法找到确切的真正的污染源，只能找到极为可能的污染源位置。
4、对于自适应滤波法模型，我们能够很好的进行时间方面的预测预估计，很有实用性。但我们只考虑时间的变化，没有深刻的考虑城市的外界因素，具有一定的局限性。
六、参考文献
[1] 尹俊，上海市城郊梯度土壤重金属空间分布特征和评价研究，上海师范大学硕士学位论文，2010、05、01。
[2]林旭文，土壤重金属污染源定位模型，科技与生活：第7期，2012年。 [3] 赵瑾瑾，扭矩、转动惯量、能量大海 nbsp 2010年07月，曲阜师范大学。 [4] 石东伟、陈冬娜，高斯扩散模型在确定污染源位置的应用，河南科技学院报，第40卷：55-58，2012。 [5]饶运章、候运炳，尾矿库废水PH值对重金属污染的影响及治理技术研究，2003年。
[6] 卢铁光、杨广林、王立坤，基于相对土壤质量法的土壤质量变化评价与分析，东北农业大学报，第34期：56-59，2002年。
七、附录
附录1：
function drawcu(x,y,z,yuansu,q)
% x,y,z分别代表坐标，yuansu代表元素的浓度，q代表功能区类别 yuansu_XY= [x y yuansu];
[xx yy]= meshgrid(0:50:3e4,0:50:2e4);
yuansu1=griddata(yuansu_XY(:,1),Cu_XY(:,2),yuansu,xx,yy,'cubic'); [c h]=contour(xx,yy,yuansu,10,'black'); clabel(c,h); hold on
18
[c1 d1]=find(q==1); x1=x(c1,1); y1=y(c1,1);
plot(x1,y1,'r*'); [c2 d2]=find(q==2); x2=x(c2,1); y2=y(c2,1); hold on
plot(x2,y2,'black.') [c3 d3]=find(q==3); x3=x(c3,1); y3=y(c3,1);
plot(x3,y3,'yo') [c4 d4]=find(q==4); x4=x(c4,1); y4=y(c4,1);
plot(x4,y4,'g+') [c5 d5]=find(q==5); x5=x(c5,1); y5=y(c5,1);
plot(x5,y5,'black*') text(x5,y5,'△') hold on
pcolor(xx,yy,cucu) shading flat
legend('As浓度等值线','生活区','工业区','山区','交通区','绿地区'); title('As空间分布(颜色条表示浓度大小)') xlabel('X(m)','color','r'); ylabel('Y(m)','color','r'); colorbar
附录二：
等标负荷的代码：
function [aixu dxu]=dengbiao(aij,n) for i=1:319
b(i,1)=sum(aij(i,n)); end
c=sum(b(:,1)); for i=1:319
ai(i,1)=b(i,1)./c; ai(i,2)=i; end
aixu=sortrows(ai,1); for i=1:319
19
d(i,1)=sum(aixu(1:i,1)); d(i,2)=aixu(i,2); end
dxu=sortrows(d,-1);
20