大气污染模型源代码.docVIP

程序代码如下： i=input(Enter i (有风时（u101.5）气态污染物点源扩散模型输入 1 ;一次地面最大落地浓度及距离输入 2 ; 有风时（u101.5）气态污染物点源（面源）扩散模型输入 3 ; 小风、静风扩散模型输入 4 ; 逆温层破坏时熏烟模式输入 5 ; 卫生防护距离估算输入 6 ；)); if i==1 syms Q U Y K He qy qz h n r1 r2 a1 a2 X Q=input(Enter Q); U=input(Enter U); Y=input(Enter Y); r1=input(Enter r1); r2=input(Enter r2); He=input(Enter He); h=input(Enter h); K=input(Enter K); a1=input(Enter a1); a2=input(Enter a2); X=input(Enter X); qy=r1*X^a1; qz=r2*X^a2; sum1=Q/(2*pi*U*qy*qz); sum2=exp(-Y*Y/(2*qz*qz)); sum3=exp(-(2*h*n-He)*(2*h*n-He)/(2*qz*qz)); sum4=exp(-(2*h*n+He)*(2*h*n+He)/(2*qz*qz)); sum5=sum3+sum4; anst=symsum(sum5,n,-K,K); C=sum1*sum2*anst; C=double(C) elseif i==2 Q=input(Enter Q); r1=input(Enter r1); r2=input(Enter r2); a1=input(Enter a1); a2=input(Enter a2); U=input(Enter U); He=input(Enter He); m1=-a1/a2; m2=(1+a1/a2)/2; m3=(1-a1/a2); m4=(1+a1/a2)^m2*He^m3*exp(m3); P1=2*r1*r2^m1/m4; m5=pi*exp(1)*U*He*He*P1; Cm=2*Q/m5 m5=1/a1; m6=-1/2/a2; Xm=(He/r2)^m5*(1+a1/a2)^m6 elseif i==3 syms Q U Y K He qy qz h n r1 r2 a1 a2 X H u10 L p Q=input(Enter Q); u10=input(Enter u10); Y=input(Enter Y); r1=input(Enter r1); r2=input(Enter r2); He=input(Enter He); h=input(Enter h); K=input(Enter K); a1=input(Enter a1); a2=input(Enter a2); X=input(Enter X); p=input(Enter p); H=input(Enter H); L=input(Enter L); qy=r1*X^a1+L/4.3; qz=r2*X^a2+H/2.15; U=u10*(H/10)^p; sum1=Q/(2*pi*U*qy*qz); sum2=exp(-Y*Y/(2*qz*qz)); sum3=exp(-(2*h*n-He)*(2*h*n-He)/(2*qz*qz)); sum4=exp(-(2*h*n+He)*(2*h*n+He)/(2*qz*qz)); sum5=sum3+sum4; anst=symsum(sum5,n,-K,K); C=sum1*sum2*anst; C=double(C) elseif i==4 t=input(Enter t); s=input(Enter s); r01=input(Enter r01); r02=input(Enter r02); x=input(Enter x); y=input(Enter y); He=input(Enter He); Q=input(Enter Q); U=input(Enter U); syms t m m=s; f=exp(-t^2/2); intf=int(f,t,-inf,m); intf=double(intf); a1=r01*r01/(r02*r02); a2=sort(2*pi); a3=-U*U/(2*r01*r01); a4=exp(s*s/2); a5=(2*pi)^1.5; h2=x*x+y*y+a1*He*He; a6=a5*r02*h2; a7=r01*sort(h2); B=intf/a2; G=exp(a3)*(1+a2*s*a4*B); C=2