3.无光照时污染气体浓度分布图 Q=9.03;%输入源强 u=6.8;%输入风速 H=146;
d=10;%步长
x=0:d:10000;%下风向距离 y=-1000:d:1000;%横风向距离 [x,y]=meshgrid(x,y);
by0=0.08*x.*(1+0.0001*x).^(-1/2); bz0=0.06*x.*(1+0.0015*x).^(-1/2); fz=(-1/2).*((y./by0).^2+(H./bz0).^2); tempy1=2.718282.^fz;
c=Q/pi/u*((by0.*bz0).^(-1)).*(tempy1); Cs=20;%输入求解条数 contour(x,y,c,Cs); shading interp; colorbar; grid;
title('无光照时污染气体分布图') 4.低速风有光照时污染气体浓度分布图 Q=9.03;%输入源强 u=1.5;%输入风速 H=146;
d=10;%步长
x=0:d:3000;%下风向距离 y=-1000:d:1000;%横风向距离 [x,y]=meshgrid(x,y);
by0=0.22*x.*(1+0.0001*x).^(-1/2); bz0=0.2*x;
fz=(-1/2).*((y./by0).^2+(H./bz0).^2); tempy1=2.718282.^fz;
c=Q/pi/u*((by0.*bz0).^(-1)).*(tempy1); Cs=20;%输入求解条数 contour(x,y,c,Cs); shading interp; colorbar; grid;
title('低速风有光照时污染分布图')
二、逻辑斯蒂函数拟合源代码
t=0:1:500;
y=1./(1+100*exp(-0.02*t)); plot(t,y,'.')
title('空气污染指数与对人体危害的关系图')
19