附录
附录一： 烟尘下风向浓度分布图
q=input('请输入源强(kg/s)q=');
Q=q*exp((1.2*10^(-5))*(1935.8^0.5)); k=input('请输入风速(m/s)￡ou='); u=k*(8^(1/7));
h=81.514; %排烟口有效高度
d=input('请输入计算精度od='); %步长 x=0:d:10000; y=-500:d:500;
[x,y]=meshgrid(x,y);
thy=0.08*x.*(1+0.0001*x).^(-1/2); thz=0.06*x.*(1+0.0015*x).^(-1/2); %高斯扩散浓度模拟计算
t=exp(-0.5*((y.^2)./(thy.*thy)+(h^2)./(thz.*thz))); z=Q./(pi*u.*thy.*thz).*t; contour(x,y,z); shading interp; colorbar; grid;
title('烟尘下风向浓度分布图')
附录二： SO2下风向浓度分布图： q=input('请输入源强(kg/s)q=');
Q=q*exp((1.2*10^(-5))*(1935.8^0.5)); k=input('请输入风速(m/s)￡ou='); u=k*(8^(1/7));
h=81.514; %排烟口有效高度
d=input('请输入计算精度od='); %步长 x=0:d:10000; y=-500:d:500;
[x,y]=meshgrid(x,y);
thy=0.08*x.*(1+0.0001*x).^(-1/2); thz=0.06*x.*(1+0.0015*x).^(-1/2); %高斯扩散浓度模拟计算
t=exp(-0.5*((y.^2)./(thy.*thy)+(h^2)./(thz.*thz))); z=Q./(pi*u.*thy.*thz).*t; contour(x,y,z); shading interp; colorbar; grid;
19
title('SO2下风向浓度分布图')
附录三： NOx下风向浓度分布图 q=input('请输入源强(kg/s)q=');
Q=q*exp((1.2*10^(-5))*(1935.8^0.5)); k=input('请输入风速(m/s)￡ou='); u=k*(8^(1/7));
h=81.514; %排烟口有效高度
d=input('请输入计算精度od='); %步长 x=0:d:10000; y=-500:d:500;
[x,y]=meshgrid(x,y);
thy=0.08*x.*(1+0.0001*x).^(-1/2); thz=0.06*x.*(1+0.0015*x).^(-1/2); %高斯扩散浓度模拟计算
t=exp(-0.5*((y.^2)./(thy.*thy)+(h^2)./(thz.*thz))); z=Q./(pi*u.*thy.*thz).*t; contour(x,y,z); shading interp; colorbar; grid;
title('NOx下风向浓度分布图')
附录四： 高斯烟羽模型：
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)
20
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
Bi=[1 0.1 2.9 2.9;2.9 1 12.9 13.9 ;0.29 0.01 1 1;0.29 0.01 1 1] [V,D]=eig(Bi)
d=V(:,1)./sum(V(:,1)) %方案层对准测层Bi的权值%
21
附录六：层次分析法用图:
各监测点污染程度(A) 目标层A
准测层B 颗粒物SO2 NOx
C1 C2 C3 C4
方案层C
图1 层次结构图
表2 准则层对目标层的判断矩阵 A B1 B2 B3 准则层权值 B1 1 8/26 8/40 0.1081 B2 26/8 1 26/40 0.3514 B3 5 40/26 1 0.5405
表2方案层对准则层B1的判断矩阵
B1 方案C1 C2 C3 C4 层权值 1 0.1 2.9 2.9 0.1541 2.9 1 12.9 13.9 0.7567 0.29 0.01 1 1 0.0446 C1 C2 C3 22
k ?k C4 0.29 0.01 1 1 0.0446 表3方案层对准则层B2的判断矩阵 方案B2 C1 C2 C3 C4 层权值 C1 1 0.1 2.9 2.9 0.1541 C2 2.9 1 12.9 13.9 0.7567 C3 0.29 0.01 1 1 0.0446 C4 0.29 0.01 1 1 0.0446 表4方案层对准则层B3的判断矩阵 方案B3 C1 C2 C3 C4 层权值 C1 1 0.1 2.9 2.9 0.1541 C2 2.9 1 12.9 13.9 0.7567 C3 0.29 0.01 1 1 0.0446 C4 0.29 0.01 1 1 0.0446 表5特征值及一次性检验 0(n=3) 1(n=4) 2(n=4) 3.0000 0 0 3.1709 0.08545 0.09494 3.1709 0.08545 0.09494 3(n=4) 3.1709 0.08545 0.09494 CIk CRk 23