气体高斯扩散模型模拟

1、%气体高斯扩散模型模拟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).A(-1/2);bz0=0.06 *x.*(1+0.0015*x).a(-1/2); by=by0.*(1+0.38*Zo);fz=(2.53-0.13*

2、log(x).*(0.55+0.042*log(x)A(-1).*Zo.A(0.35-0.03*log(x); bz=bz0.*fz;tempy1=-y.*y./by./by./2; tempy2=2.718282.A(tempy1);c=Q/pi/u*(by.*bz).A(-1).*tempy2;figure(1);Cs=i nput('请输入所有求解浓度(mg/mA3):'); contour(x,y,c,Cs);shading interp; colorbar;grid;xlabel('X 轴向距离 (m)');ylabel('Y 轴向距离 (m)

3、');title(' 气体扩散下风向浓度分布图 ');set(gcf,'color','white');figure(2);Cs=input(' 请输入所有求解浓度 (mg/mA3):'); contour(x,y,c,Cs);shading interp; colorbar;grid;xlabel('X 轴向距离 (m)');ylabel('Y 轴向距离 (m)');title(' 气体扩散下风向浓度分布图 ');set(gcf,'color','w

4、hite');figure(3);Cs=input(' 请输入所有求解浓度 (mg/mA3):'); 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(' 请输入所有求解浓度contour(x,y,c,Cs);sh

5、ading interp;colorbar;grid;xlabel('X 轴向距离 (m)');ylabel('Y 轴向距离 (m)');title(' 气体扩散下风向浓度分布图 set(gcf,'color','white');figure(5);Cs=input(' 请输入所有求解浓度contour(x,y,c,Cs);shading interp;colorbar;grid;xlabel('X 轴向距离 (m)');ylabel('Y 轴向距离 (m)');title('

6、; 气体扩散下风向浓度分布图 set(gcf,'color','white');figure(6);Cs=input(' 请输入所有求解浓度contour(x,y,c,Cs);shading interp;colorbar;grid;xlabel('X 轴向距离 (m)');ylabel('Y 轴向距离 (m)');title(' 气体扩散下风向浓度分布图 set(gcf,'color','white');figure(7);Cs=input(' 请输入所有求解浓度contour(x,y,c,Cs);shading interp;colorbar;grid;xlabel('X 轴向距离 (m)');ylabel('Y 轴向距离 (m)');(mg/mA3):');');(mg/mA3):');');(mg/mA3):&#