Gossip-gas 高斯烟羽扩散模型的matlab实现

高斯烟羽扩散模型的matlab实现，可以用来模拟气体扩散的过程。

Q=input（‘请输入泄漏源强（mg／s）：Q=‘）;
u=input（‘请输入计算风速（m／s）：u=‘）;
d=input（‘请输入计算精度（m）：d=‘）;
Z0=input（‘请输入地面粗糙长度（m）：Z0=‘）;
[xy]=meshgrid（0:d:1000-100:d:100）;%定义解空间和计算精度
by0=0.08*x.*（1+0.0001*x）.^（-1/2）;%计算Y轴向的基本扩散参数
bz0=0.06*x.*（1+0.0015*x）.^（-1/2）;%计算z轴向的扩散参数
by=by0.*（1+0.38*Z0）;%对Y轴向的扩散参数按地面粗糙长度进行修正。
fz=（2.53-0.13*log（x））.*（0.55+0.042*log（x））.^（-1）.*Z0.^（0.35-0.03*log（x））;%按地面粗糙长度计算=轴向的扩散参数修正系数
bz=bz0.*fz;%对z轴向的扩散参数按地面粗糙长度进行修正。
%高斯扩散浓度模拟计算％
tempy1=-y.*y./by./by./2;
tempy2=2.718282.^（tempy1）;
c=Q/pi/u*（（by.*bz）.^（-1））.*tempy2

Csa=input（‘请输入初始求解浓度（mg／m3）：‘）;
NS=20;
Cs=[Csa Csa-NS];
Csi=Csa-NS;

for i=2:（Csa/NS-1）
    Cs（i）=Csi;
    Csi=Csi-NS;
contour（xy

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1587  2015-04-16 19:57  高斯烟羽模型\gsyy.m

     目录          0  2015-11-09 22:13  高斯烟羽模型

----------- ---------  ---------- -----  ----

                 1587                    2