vb+access大气污染模型源代码论文

程序代码如下：i=input(Enteri (有风时（u101.5）气态污染物点源扩散模型输入 1;一次地面最大落地浓度及距离输入2;有风时（u10);ifi=1symsQUYKHeqyqzhnr1r2a1a2XQ=input(EnterQ);U=input(EnterU);Y=input(EnterY);r1=input(Enterr1);r2=input(Enterr2);He=input(EnterHe);h=input(Enterh);K=input(EnterK);a1=input(Entera1);a2=input(Entera2);X=input(EnterX);qy=r1*Xa1;qz=r2*Xa2;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)elseifi=2Q=input(EnterQ);r1=input(Enterr1);r2=input(Enterr2);a1=input(Entera1);a2=input(Entera2);U=input(EnterU);He=input(EnterHe);m1=-a1/a2;m2=(1+a1/a2)/2;m3=(1-a1/a2);m4=(1+a1/a2)m2*Hem3*exp(m3);P1=2*r1*r2m1/m4;m5=pi*exp(1)*U*He*He*P1;Cm=2*Q/m5m5=1/a1;m6=-1/2/a2;Xm=(He/r2)m5*(1+a1/a2)m6elseifi=3symsQUYKHeqyqzhnr1r2a1a2XHu10LpQ=input(EnterQ);u10=input(Enteru10);Y=input(EnterY);r1=input(Enterr1);r2=input(Enterr2);He=input(EnterHe);h=input(Enterh);K=input(EnterK);a1=input(Entera1);a2=input(Entera2);X=input(EnterX);p=input(Enterp);H=input(EnterH);L=input(EnterL);qy=r1*Xa1+L/4.3;qz=r2*Xa2+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)elseifi=4t=input(Entert);s=input(Enters);r01=input(Enterr01);r02=input(Enterr02);x=input(Enterx);y=input(Entery);He=input(EnterHe);Q=input(EnterQ);U=input(EnterU);symstmm=s;f=exp(-t2/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*Q*G/a6S=U*x/a7elseifi=5hf=input(Enterhf);He=input(EnterHe);qz=input(Enterqz);Q=input(EnterQ);U=input(EnterU);y=input(Entery);qy=input(Enterqy);symstPa1=sort(2*pi);P=(hf-He)/qz;f=exp(-t2/2);intf=int(f,t,-inf,P);intf=double(intf);B=intf/a1;qyf=qy+He/8;a2=Q/(a1*U*hf*qyf);a3=-y*y/(2*qyf*qyf);Cf=a2*exp(-a3)*BelseQc=input(EnterQc);A=input(EnterA);B=input(EnterB);C=input(EnterC);D=input(EnterD);L=input(EnterL);r=input(Enterr);Y=(B*LC+0.25*r2)2*L