n_0=0.332;
n_1=0.6;
n_2=0.4;
T=35;
Pa=1010;
Ts=273.15;
Qv=96.30;
Hs=80;
Average=[1.905 1.848 1.509 1.614 2.808 2.704 2.533 2.637];
H_1=n_0.*(0.35.*Pa.*Qv.*T./Ts).^n_1.*Hs.^(n_2./Average);
H=Hs+H_1;
H=Hs+H_1;
n=0.25;
u_1=(10./H).^0.25.*Average;
sigma=[19.1 35.8 51.6 67.0 81.4 95.8 123 151 178 203 228 253 278 %%sigma_y
10.7 20.5 30.2 40.5 51.2 62.8 84.6 109 133 157 181 207 233];%%sigma_z
sigma_yz=sigma(1,:).*sigma(2,:);
sigma_zz=sigma(2,:).*sigma(2,:);
Q=[9027.78 14814.81 19340.28];
C=zeros(3,13,8);
for i=1:3
for j=1:13
for k=1:8
C(i,j,k)=Q(i)/(pi*u_1(k)*sigma_yz(j))*exp(-H(k)*H(k)/(2*sigma_zz(j)));
end
end
end
C(:,:,6);%%西南方向距离与污染浓度关系
plot_x=[0.1 0.2 0.3 0.4 0.5 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0]
hold on;
plot(plot_x,C(1,:,6),'Blue');
plot(plot_x,C(2,:,6),'Red');
plot(plot_x,C(3,:,6),'Black');
legend('烟尘','SO2','NOx')