geta.m

污染物连续衰变浓度变化 基于matlab的编程

function aa=GetA(v,D,num,a0,k,x,t)
% %global v,D,num,c0,k,a0,t
% [v,D,num,c0,k]=Readfile;
% a0=Geta0(num,c0,k);

global aa
aa=zeros(num,1);
% kk=1;
% for x=0:0.1:100
  
for i=1:num
beta=sqrt(v^2/(4*D^2)+k(i)/D);
aa(i) =a0(i)*0.5*exp(v*x/2/D)*(exp(-beta*x)*erfc((x-sqrt(v^2+4*k(i)*D)*t)/2/sqrt(D*t))+exp(beta*x)*erfc((x+sqrt(v^2+4*k(i)*D)*t/2/sqrt(D*t))));
end
% kk=kk+1;
% end