%气体高斯扩散模型模拟 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/m^3):'); 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/m^3):'); 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/m^3):'); 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/m^3):'); 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/m^3):'); 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/m^3):'); 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/m^3):'); contour(x,y,c,Cs); shading interp; colorbar; grid; xlabel('X 轴向距离(m)'); ylabel('Y 轴向距离(m)'); 

title('气体扩散下风向浓度分布图'); set(gcf,'color','white');