衰减方程如下:
matlab程序:
C0=1;
k=0.36;
t=linspace(0,15,90);
dt=t(2)-t(1);
C=C0*exp(-k*t);%解析解
%下面求解数值解
C1=zeros(1,length(t));
C1(1)=C0;
%下面用显式求解
for i=2:length(t)
C1(i)=C1(i-1)-k*dt*C1(i-1);
end
%下面用隐式求解
C2=C1;
for i=2:length(t)
C2(i)=C2(i-1)/(1+k*dt);
end
%四阶龙格库塔
lambda=-k;
C3=C1;
for i=2:length(t)
C3(i)=C3(i-1)*(1+dt*lambda+0.5*(dt*lambda)^2+1/6*(dt*lambda)^3+1/24*(dt*lambda)^4);
end
%半隐半显
C4=C1;
for i=2:length(t)
C4(i)=C4(i-1)*(1-0.5*k*dt)/(1+