matlab衰减指数的程序是,Matlab数值计算最简单的一个例子——指数衰减

放射性衰变是指数衰减的典型例子。另外还有化学反应某反应物的减少，RC电路电流的减小，大气压随海拔高度的减小等。

指数衰减的方程：

\begin{equation}

\frac{dN(t)}{dt}=-\frac{N(t)}{\tau}

\label{eq1}

\end{equation}

其中，$N(t)$为$t$时刻的物理量$N$，对于放射性衰变，$N$就是未衰变的原子核数目。$\tau$为时间常数。

方程\eqref{eq1}有解析解：

$$N(t)=N(0)\exp(-t/\tau)$$

这个解可以通过Matlab符号计算求得：

dsolve('DN=-N/tau')

ans =

C3*exp(-t/tau)

数值求解方程\eqref{eq1}，可用欧拉格式将方程离散化。

$$t_i=(i-1) \Delta t,\quad i=1,2,\dots,\mathrm{npoints}$$

$$\frac{dN(t)}{dt}\approx\frac{N(t)-N(t-\Delta t)}{\Delta t}$$

将以上两式带入方程\eqref{eq1}，得离散之后的方程：

$$N(t_{i+1})=N(t_i)-N(t_i)\frac{\Delta t}{\tau}$$

代入初始条件，即可得解。

下面写个Matlab 脚本文件，重复出Computational Physics_Giordano 2nd Edition的图1.1，pp11

%

% Exp