气体高斯扩散模型模拟%气体高斯扩散模型模拟 Q=input(请输入泄漏源强(mg/s):Q=); u=input(请输入计算风速(m/s):u=); d=input(请输入计算精度(m):d=); Zo=input(请输入地面粗糙长度(m):Zo=); x,y=meshgrid(50:d:1000,-100:d:100); by0=0.08*x.*(1+0.0001*x).(-1/2); bz0=0.06 *x.*(1+0.0015*x).(-1/2); by=by0.*(1+0.38*Zo); fz=(2.53-0.13*log(x).*(0.55+0.042*log(x).(-1).*Zo.(0.35-0.03*log(x); bz=bz0.*fz; tempy1=-y.*y./by./by./2; tempy2=2.718282.(tempy1); c=Q/pi/u*(by.*bz).(-1).*tempy2; figure(1); Cs=input(请输入所有求解浓度(mg/m3):); contour(x,y,c,Cs); shading interp; colorbar; grid; xlabel(X轴向距离(m); ylabel(Y轴向距离(m); title(气体扩散下风向浓度分布图); set(gcf,color,white); figure(2); Cs=input(请输入所有求解浓度(mg/m3):); contour(x,y,c,Cs); shading interp; colorbar; grid; xlabel(X轴向距离(m); ylabel(Y轴向距离(m); title(气体扩散下风向浓度分布图); set(gcf,color,white); figure(3); Cs=input(请输入所有求解浓度(mg/m3):); contour(x,y,c,Cs); shading interp; colorbar; grid; xlabel(X轴向距离(m); ylabel(Y轴向距离(m); title(Cs气体扩散下风向浓度分布图); set(gcf,color,white); figure(4); Cs=input(请输入所有求解浓度(mg/m3):); contour(x,y,c,Cs); shading interp; colorbar; grid; xlabel(X轴向距离(m); ylabel(Y轴向距离(m); title(气体扩散下风向浓度分布图); set(gcf,color,white); figure(5); Cs=input(请输入所有求解浓度(mg/m3):); contour(x,y,c,Cs); shading interp; colorbar; grid; xlabel(X轴向距离(m); ylabel(Y轴向距离(m); title(气体扩散下风向浓度分布图); set(gcf,color,white); figure(6); Cs=input(请输入所有求解浓度(mg/m3):); contour(x,y,c,Cs); shading interp; colorbar; grid; xlabel(X轴向距离(m); ylabel(Y轴向距离(m); title(气体扩散下风向浓度分布图); set(gcf,color,white); figure(7); Cs=input(请输入所有求解浓度(mg/m3):); contour(x,y,c,Cs); shading interp; colorbar; grid; xlabel(X轴向距离(m); ylabel(Y轴向距离(m); title(气体扩散下风向浓度分布图); set(gcf,color,white);