点污染源空气污染扩散模型

8 点、中午12 点、晚上9 点都没有排放气体，该怎么算，是不是需要找到一个关于时间t的函数，来计算多长时间之后污染还剩下多少

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

这个函数对吗？该调用什么函数？

问题：

建立单污染源空气污染扩散模型，描述其对周围空气污染的动态影响规律。

现有河北境内某一工厂废气排放烟囱高50m，主要排放物为氮氧化物。早上9 点至下午

3 点期间的排放浓度为406.92mg/m3，排放速度为1200m3

/h；晚上10 点-凌晨4 点期间

的排放浓度为1160mg/m3，排放速度为5700m3

/h；通过你的扩散模型求解该工厂方圆51

公里分别在早上浓度8 点、中午12 点、晚上9 点空气污染分布和空气质量等级。

源代码

clear all

clc

[x,y]=meshgrid(0:20:5100,0:20:5100);

Q=135.64; z=1.5; H=50; u=1.94;

sigy=0.3914238*x.^0.865014;

sigz=0.0757182*x.^1.00770;

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

mesh(x,y,c);

xlabel('X'),ylabel('Y'),zlabel('C'),

clear all

clc

[x,y]=meshgrid(-5100:20:5100,-5100:20:5100);

Q=1836.7; z=1.5; H=50; u=1.7;

sigy=0.3914238*x.^0.865014;

sigz=0.0757182*x.^1.00770;

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

mesh(x,y,c);

xlabel('X'),ylabel('Y'),zlabel('C'),

分享到：

2015-05-29 16:32

提问者采纳

clear all

[x,y]=meshgrid(-51000:100:51000,-51000:100:51000);

Q=135.64; z=1.5; H=50; u=1.94;

sigy=0.3914238*x.^0.865014;

sigz=0.0757182*x.^1.00770;

%c=@(x,y)Q./(2*pi*sigy.*sigz*u+eps).*exp(-0.5*(y.^2)./((sigy+eps).^2)).*(exp(-0.5*(z-H).^2./((sigz+eps).^2))+exp(-0.5*(z+H).^ 2./((sigz+eps).^2)));

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

%g=dblquad(c,-51000,51000,-51000,51000,100),

%g=abs(g/51000/51000*1000),

mesh(x,y,c);

xlabel('X'),ylabel('Y'),zlabel('C'),

clear all

clc

[x,y]=meshgrid(-51000:100:51000,-51000:100:51000);

Q=1836.7; z=1.5; H=50; u=1.7;

sigy=0.3914238*x.^0.865014;

sigz=0.0757182*x.^1.00770;

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

mesh(x,y,c);

xlabel('X'),ylabel('Y'),zlabel('C'),

clear all

[x,y]=meshgrid(-51000:100:51000,-51000:100:51000);

Q=1836.7; z=1.5; H=50; u=1.7;

sigy=0.3914238*x.^0.865014;

sigz=0.0757182*x.^1.00770;

c=@(x,y)Q./(2*pi*sigy.*sigz*u+eps).*exp(-0.5*(y.^2)./((sigy+eps).^2)).*(exp(-0.5*(z-H).^2./((sigz+eps).^2))+exp(-0.5*(z+H).^2./ ((sigz+eps).^2)));

%mesh(x,y,c);

%xlabel('X'),ylabel('Y'),zlabel('C'),

g=dblquad(c,-51000,51000,-51000,51000,100),

g=abs(g/51000/51000*1000),

其具体标准如下:当空气污染指数达0—50时为1级；51—100时为2级；101—200时为3级；201—300时为4级；300以上时为5级。其中3级属于轻度污染，4级属于中度污染，5级则属于重度污染

按照国家统一规定，空气质量达到优良标准即达到国家质量二级标准是指空气污染指数小于等于100，如果空气污染指数小于等于50，说明空气质量为优。空气污染指数大于50且小于等于100时，说明空气质量为良好。具体如何判断，下面就给出空气质量的评定标准。

空气质量的评定标准：