附录二
中国空气质量等级分类 空气质量指数值 中国级别 中国分类 25
0-50 51-100 101-150 151-200 201-250 一级 二级 三级 四级 五级 优 良 轻度污染 中度污染 重度污染 251-300 301-350 351-400 401-450 451-500 500+ 附录三
问题一对应MATLAB代码
%利用层次分析加权法计算各污染物的权重 A=[1,6,2,5,1/2,1/2 1/6,1,1/3,5/6,1/12,1/12 1/2,3,1,2/5,1/4,1/4 1/5,6/5,2/5,1,1/10,1/10 2,12,4,10,1,1 2,12,4,10,1,1]; [D,X]=eig(A); ZD=max(max(X)); [m,m]=size(D); p=max(X); for i=1:m if p(i)==ZD q=i end end
六级 严重污染 26
z=sum(D); for i=1:m for j=1:m
D(i,j)=D(i,j)/z(j); end end B=(D(:,q))'
RI=[0,0,0.58,0.9,1.12,1.24,1.32,1.41,1.45,1.49,1.51]; CI=(ZD-m)/(m-1) CR=CI/RI(m) Q=zeros(11,1); for i=1:11 Q(i,1)=0; for j=1:6
Q(i,1)=Q(i,1)+sum(i,j)*B(1,j); end End
问题三对应matlab码 %中午12点图形 x=0:100:20000;
oy=0.433384.*(x).^0.565188; oz=0.0570251.*(x).^1.09356; H=60.8; u=2.5; q=135.64;
C=q./(2*pi*u*oy.*oz).*exp(-H^2./(oz.^2)); plot(x,C,'bd',x,C,'b-') %晚上9点图形 x=0:100:20000;
oy=0.433384.*(x).^0.565188;
27
oz=0.0570251.*(x).^1.09356; H=60.8; u=2.5; q=1836.7;
C_=q*3600*6./((4*pi*10*3600)^(3/2).*((oy.*oy.*oz).^(1/2))).*exp(-x.^2./(4*oy*10*3600))*2.*exp(-H^2./(4*oz*10*3600)) plot(x,C_,'bd',x,C_,'b-') %晚上9点图形 x=0:100:20000;
oy=0.433384.*(x).^0.565188; oz=0.0570251.*(x).^1.09356; H=60.8; u=2.5; q=135.64;
C_=q*3600*6./((4*pi*12*3600)^(3/2).*((oy.*oy.*oz).^(1/2))).*exp(-x.^2./(4*oy*12*3600))*2.*exp(-H^2./(4*oz*12*3600)) plot(x,C_,'bd',x,C_,'b-')
问题四对应MATLAB代码 %二环早上8点污染物浓度曲线 x=0:10:1000; oz=0.24*x; H=1; u=2.5;
q=8000*40.05/3600;
C=(2/pi)^(1/2)*q/u./oz.*exp(-H^2/2./((oz).^2)); plot(x,C,'bd',x,C,'b-') hold on
%二环中午12点污染物浓度曲线
28
x=0:10:1000;
oz=0.24*(x)./((1+0.001*x).^(-1/2)); H=1; u=2.5;
q=4000*40.05/3600;
C=(2/pi)^(1/2)*q/u./oz.*exp(-H^2/2./((oz).^2)); plot(x,C,'r*',x,C,'r-') %二环晚上9点污染物浓度曲线 x=0:10:1000;
oz=0.08*(x)./((1+0.0015*x).^(-1/2)); H=1; u=2.5;
q=5600*40.05/3600;;
C=(2/pi)^(1/2)*q/u./oz.*exp(-H^2/2./((oz).^2)); plot(x,C,'y+',x,C,'y-') xlabel('距线源距离x') ylabel('CO浓度C')
title('二环CO污染物浓度曲线')
legend('蓝色8点','红色12点','黄色21点') x=0:100:25000;
oy=0.04.*(x)./((1+0.0001*x).^(-1/2)); oz=0.016.*(x)./((1+0.0003*x).^(-1)); H=60.8; u=2.5; q2=135.64
C_=q2*3600*6./((4*pi*12*3600)^(3/2).*((oy.*oy.*oz).^(1/2))).*exp(-x.^2./(4*oy*12*3600))*2.*exp(-H^2./(4*oz*12*3600)) plot(x,C_,'bd',x,C_,'b-') title('晚上9点污染物浓度分布'); xlabel('距工厂距离X');
29
ylabel('污染物浓度C')
30