[9]孙志宽.环境与可持续发展2013.(5).107-109 [10]孙志宽.高斯烟羽扩散模型再研究.2013
[11]胡世明.气体释放源的三维瞬态重气扩散研究.劳动保护科学技术,2002,3（20）：28-30 [12]我国垃圾处理补贴费标准过低.http://news.solidwaste.com.cn/view/id_21314.2014-5-18 [13]关于完善核电上网电价机制有关问题的通知.2013.7.8 [14]《数学建模与实验MATLAB教程》.第11讲，回归分析
[15]设备管理与设备故障率.http://www.gztaiyou.com/html/2012102694259.html.2014-5-18
20
附录
附录一：高斯烟羽模型代码：
x=input('请输入目标监测点相对横坐标(m)：x='); y=input('请输入目标监测点相对纵坐标(m)：y=0'); H=input('请输入污染源有效高度(m)：H=');
d=input('请输入大气稳定度（提示：以0~5分别代表A~F六类）：d=3'); Q=input('请输入污染物释放速率(kg/s)：Q='); u=input('请输入平均风速(m/s)：u='); %扩散系数的计算% if(d==0)
m=0.22*x/(1+0.0001*x)^0.5; n=0.2*x; end
if(d==1)
m=0.16*x/(1+0.0001*x)^0.5; n=0.12*x; end
if(d==2)
m=0.11*x/(1+0.0001*x)^0.5; n=0.08*x/(1+0.0002*x)^0.5; end
if(d==3)
m=0.08*x/(1+0.0001*x)^0.5; n=0.06*x/(1+0.0015*x)^0.5; end
if(d==4)
m=0.06*x/(1+0.0001*x)^0.5; n=0.03*x/(1+0.0003*x); end
if(d==5)
m=0.04*x/(1+0.0001*x)^0.5; n=0.016*x/(1+0.0003*x); end
c=(Q/pi*u*m*n)*exp((-1/2)*(y/m)^2)*exp(H^2/(2*(n^2)))
附录二：层次分析法对各监测点污染程度进行综合排名的代码
%层次分析法计算污染程度%
A=[1 8/26 8/40;26/8 1 26/40;5 40/26 1] [V,D]=eig(A)
d=V(:,1)./sum(V(:,1)) %准则层对目标层的权值% clear
B1=[1 3 15 2.5 27.5 12.7 23.6 9.1 19.4;1/3 1 5 0.85 9.2 4.2 7.9 3.1 6.5;1/15 0.2 1 0.17 1.8 0.85 1.6
21
0.16 1.3;0.4 1/0.85 5.9 1 10.8 5 9.3 3.6 7.6;0.036 0.1 0.54 0.09 1 0.46 0.86 0.33 0.71;0.079 0.24 1.2 0.2 2.2 1 1.8 0.72 1.5;0.042 0.13 0.64 0.11 1.2 0.54 1 0.39 0.82;0.109 0.33 1.6 0.28 3 1.4 2.6 1 2.1;0.51 0.15 0.77 0.13 1.4 0.65 1.2 0.47 1] [V,D]=eig(B)
d=V(:,1)./sum(V(:,1)) %方案层对准测层B1的权值% clear
B2=[1 30 3 1.9 28.5 13 25 9.3 20;0.033 1 0.1 0.063 0.95 0.43 0.83 0.31 0.67;0.33 10 1 0.63 9.5 4.3 8.3 3.1 6.7;0.53 16 1.6 1 15 6.8 13 4.9 10;0.035 1.1 0.1 0.067 1 0.45 0.87 0.34 0.69;0.077 2.3 0.23 0.15 2.2 1 1.9 0.72 1.5;0.042 1.2 0.12 0.77 1.1 0.52 1 0.38 0.79;0.1 3.2 0.32 0.2 2.9 1.4 2.6 1 2.1;0.05 1.5 0.15 0.1 1.45 0.67 1.3 0.48 1] [V,D]=eig(B2)
d=V(:,1)./sum(V(:,1)) %方案层对准测层B2的权值% clear
B3=[1 29 15 2.5 27 12.5 23 8.9 19;0.034 1 0.5 0.086 0.92 0.43 0.8 0.3 0.65;0.067 2 1 0.17 1.8 0.86 1.6 0.61 1.3;0.4 1.2 5.9 1 10.7 5 9.3 3.5 7.6;0.037 1.09 0.56 0.093 1 0.46 0.87 0.33 0.7;0.08 2.3 1.16 0.2 2.2 1 1.9 0.71 1.5;0.043 1.3 0.625 0.11 1.2 0.53 1 0.38 0.81;0.11 3.3 1.6 0.29 3.03 1.4 2.6 1 2.1;0.053 1.54 0.77 0.13 1.4 0.67 1.2 0.48 1] [V,D]=eig(B3)
d=V(:,1)./sum(V(:,1)) %方案层对准测层B3的权值% %据此进行一致性检验，一致性可接受%
附录三： 层次分析法求解的判断矩阵
表1 方案层对准则层B2的判断矩阵
B2 C1 C2 C3 C4 C5 C6 C7 C8 C9 C1 1 0.033 0.33 0.53 0.035 0.077 0.042 0.1 0.05 C2 30 1 10 16 1.1 2.3 1.2 3.2 1.5 C3 3 0.1 1 1.6 0.1 0.23 0.12 0.32 0.15 C4 1.9 C5 28.5 C6 13 0.43 4.3 6.8 0.45 1 0.52 1.4 0.67 C7 25 0.83 8.3 13 0.87 1.9 1 2.6 1.3 C8 9.3 C9 20 方案层权值 0.4476 0.063 0.95 0.63 1 0.067 0.15 0.77 0.2 0.1 9.5 15 1 2.2 1.1 2.9 1.45 0.31 0.67 0.0149 3.1 4.9 6.7 10 0.1488 0.2345 0.34 0.69 0.0157 0.72 1.5 0.0344 0.38 0.79 0.0344 1 0.48 2.1 1 0.0469 0.0229
22
表2 方案层对准则层B3的判断矩阵
B3 C1 C2 C3 C4 C5 C6 C7 C8 C9 C1 1 0.034 0.067 0.4 0.037 0.08 0.043 0.11 0.053 C2 C3 C4 2.5 0.086 0.17 1 0.093 0.2 0.11 0.29 0.13 C5 27 0.92 1.8 10.7 1 2.2 1.2 3.03 1.4 C6 12.5 0.43 0.86 5 C7 23 0.8 1.6 9.3 C8 8.9 0.3 0.61 3.5 C9 19 0.65 1.3 7.6 0.7 1.5 方案层权值 0.5600 0.0191 0.0381 0.2009 0.0208 0.0448 0.0242 0.0626 0.0295 29 15 1 2 1.2 0.5 1 5.9 1.09 0.56 2.3 1.16 0.46 0.87 0.33 1 0.53 1.4 0.67 1.9 1 2.6 1.2 0.71 1.3 0.625 3.3 1.6 0.38 0.81 1 0.48 2.1 1 1.54 0.77
23
24