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

点污染源空气污染扩散模型.doc8点、中午12点、晚上9点都没有排放气体，该怎么算，是不是需要找到一个关于时间 t的函数，来计算多长时间之后污染
还剩下多少
c=Q./(2*pi*sigy.*sigz*u+eps)*exp(-*()./((sigy+eps).A2)).*(exp(-*(z-H).A2./((sigz+eps).A2))+exp(-*(z+H).A2./((sigz+
eps)A2)));
这个函数对吗？该调用什么函数？
问题：
建立单污染源空气污染扩散模型，描述其对周围空气污染的动态影响规律。
现有河北境内某一工厂废气排放烟囱高 50m，主要排放物为氮氧化物。早上 9点至下午
3点期间的排放浓度为 ，排放速度为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=; z=; H=50; u=;
sigy=*;
sigz=*;
c=Q./(2*pi*sigy.*sigz*u+eps).*exp(-*()./((sigy+eps).A2)).*(exp(-*(z-H).A2./((sigz+eps).A2))+exp(-*(z+H).A2./((sigz+
eps)A2)));
mesh(x,y,c);
xlabel('X'),ylabel('Y'),zlabel(C),
clear all
clc
[x,y]=meshgrid(-5100:20:5100,-5100:20:5100);
Q=; z=; H=50; u=;
sigy=*;
sigz=*;
c=Q./(2*pi*sigy.*sigz*u+eps).*exp(-*()./((sigy+eps).A2)).*(exp(-*(z-H).A2./((sigz+eps).A2))+exp(-*(z+H).A2./((sigz+
eps)A2)));
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=; z=; H=50; u=;
sigy=*;
sigz=*x".00770;
%c=@(x,y)Q./(2*pi*sigy*sigz*u+eps).*exp(-*()./((sigy+eps).A2)).*(exp(-0.