main.m

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

[v,D,num,c0,k,t,xmin,xmax,interval]=Readfile;

a0=Geta0(num,c0,k);
C=zeros(num,1);
%tmin=0;tmax=400;interval=1;x=10;
i=1
%for t=tmin:interval:tmax
 for x=xmin:interval:xmax
    aa=GetA(v,D,num,a0,k,x,t);
    Con=GetC(num,aa,k);
    if i>1
    C=[C Con];
    end
    plot(x,Con);
    %plot(t,Con);
    i=i+1;
hold on;
end
% plot(x,aa);
% hold on



% for xd=0:1:100
%     AV=zeros(num,1);
%     AV =Avariable(ve,Di,n,cn0,kr,xd,ti);
% end

% for x=0:0.1:100
%     
%     for i=1:num
%     plot(x,Con(i,kk));
%     hold on
%     end
%     kk=kk+1;
% end