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