高斯羽烟模型

已经是数学建模退役老油子了，偶尔发点模型计算方法出来vanvan 对高斯羽烟模型方程: 式中c为污染物浓度(单位：kg/m3) Q为源强(单位：kg/s) u为泄漏高度的平均风速(单位：m/s) y、z分别用浓度标准偏差表示的y轴及z轴上的扩散参数 H为泄漏有效高度(单位：m)

扩散系数系数

稳定度 a b c d A 0.527 0.865 0.28 0.90 B 0.371 0.866 0.23 0.85 C 0.209 0.897 0.22 0.80 D 0.123 0.905 0.20 0.76 E 0.098 0.902 0.15 0.73 F 0.065 0.902 0.12 0.67

Matlab 仿真代码:

clc;
clear all;

[x,y]=meshgrid(0:20:5000,-2000:20:2000);

Q=1000;%源强
%sigma_y=ax^b 水平扩散系数
a=0.123;
b=0.905;
sigy=a*x.^b;
%sigma_z=cx^d 垂直扩散系数
c=0.2;
d=0.76;
sigz=c*x.^d;
u=1;%风速
H=200;%泄漏有效高度
z=0;%高度

Qpi=2*Q./(2*pi*u.*sigy.*sigz+eps);
ex1=exp(-0.5*(y./(sigy+eps)).^2);
ex2=exp(-0.5*((z-H)./(sigz+eps)).^2)+exp(-0.5*((z+H)
C=Qpi.*ex1.*ex2;

contour(x,y,C,ShowText,on)

高斯羽烟模型 已经是数学建模退役老油子了，偶尔发点模型计算方法出来vanvan 对高斯羽烟模型方程: 式中c为污染物浓度(单位：kg/m3) Q为源强(单位：kg/s) u为泄漏高度的平均风速(单位：m/s) y、z分别用浓度标准偏差表示的y轴及z轴上的扩散参数 H为泄漏有效高度(单位：m) 扩散系数系数 稳定度 a b c d A 0.527 0.865 0.28 0.90 B 0.371 0.866 0.23 0.85 C 0.209 0.897 0.22 0.80 D 0.123 0.905 0.20 0.76 E 0.098 0.902 0.15 0.73 F 0.065 0.902 0.12 0.67 Matlab 仿真代码: clc; clear all; [x,y]=meshgrid(0:20:5000,-2000:20:2000); Q=1000;%源强 %sigma_y=ax^b 水平扩散系数 a=0.123; b=0.905; sigy=a*x.^b; %sigma_z=cx^d 垂直扩散系数 c=0.2; d=0.76; sigz=c*x.^d; u=1;%风速 H=200;%泄漏有效高度 z=0;%高度 Qpi=2*Q./(2*pi*u.*sigy.*sigz+eps); ex1=exp(-0.5*(y./(sigy+eps)).^2); ex2=exp(-0.5*((z-H)./(sigz+eps)).^2)+exp(-0.5*((z+H) C=Qpi.*ex1.*ex2; contour(x,y,C,ShowText,on)