求高斯气体扩散模型 matlab编程代码，求每个得解释

meshgrid(-3000:20:3000,0:20:5000);
Q=14831;
sigy=0.314278*x.^0.876;
sigz=0.057025*x.^1.0935;
H=80;
u=1.72;
c=Q./(2*pi*sigy.*sigz*u).*exp(-0.5*(y.^2)/sigy^2).*exp(-0.5*(x.^2)/sigz^2);
mesh(x,y,c);
xlabel('X'),ylabel('Y'),zlabel('C')

[x,y]=meshgrid(0:20:5000,-2000:20:2000); %设定网格Q=0.015079644737231;%源强
sigy=0.16*x.*(1+0.0001*x).^(-0.5);
sigz=0.12*x;%
H=[209.981728 245.982304 215.960608 385.389424 162.662320 172.409248 224.093152 162.138016];%泄漏源有效高度
u=1./[0.763762 1.013766 0.805282 1.981871 0.435155 0.502842 0.861758 0.431514];
f=[0.067692 0.033846 0.070769 0.043076 0.080000 0.233846 0.083076 0.387692];
s=size(u);%计算u的元素个数，以便设置循环次数for i=1:s(2)
    Qpi=Q./(pi*u(i)*sigy.*sigz+eps);%计算公式第一部分
    ex1=exp(-0.5*(y./(sigy+eps)).^2);%计算公式第二部分
    ex2=exp(-0.5*(H(i)./(sigz+eps)).^2);%计算公式第三部分
    X1(:,:,i)=Qpi.*ex1.*ex2*f(i);%计算扩散气体浓度高维数组
end
X=X1;
for j=2:s(2)
    X(:,:,j)=X(:,:,(j-1))+X(:,:,j);%计算扩散气体浓度高维数组各二维数组之和
endX2=1000000*X(:,:,s(2));%计算求和后最终扩散气体浓度mesh(x,y,X2);