【数值分析--有限差分法】

引入–问题

(1)一维有限差分

推导：
对于函数$T\left(t\right)，T\left(t\right)$连续。
使用$T(t)+∆t$的泰勒展开式来逼近$T\left(t\right)$：
移项得：
在将∆t移项到左边：

因为差分方程是将连续变为离散，而题干中只需要对$\frac{dT\left(t\right)}{dt}$做近似，因此考虑只保留等式右边得第一项，将等式右边得其余项视为截断误差：

如果截断误差为0，则有：

这个公式称为欧拉前向公式已知题干条件：

$\frac{dT\left(t\right)}{dt}=-T\left(t\right)$

将其带入差分方程：

整理得到：

已知初始条件：
T(0)=1
假设将t的整个区间分成n份，步长为∆t需要自行设定，理论上∆t的值越小近似效果越好，则分别对应$t_1=0,t_2=∆t,\dots ,t_n=n-1∆t,t_{n+1}=n∆t$。当$t=t_1$时：

即：

此时根据初始条件T(0)=1，计算出$T\left(t_2\right)$时刻的近似值$T_{t_2}$然后以此类推，根据$T_{t_2}$求出$T_{t_3}$，…,直到求出$T_{n+1}$。这样就使用有限差分方法将连续的$T\left(t\right)$近似成了离散的$T\left(t_1\right)，T_{t_2},\dots ,T_{n+1}$。

(2)二维有限差分

从题干中可以发现：t是时间变量，x,y是空间变量，不仅需要对t差分做近似，还需要对x,y差分做近似。，将x分为m份，y分为u份，t分为n份，则在x方向上$x_i\left(i=1,2,\dots ,m+1\right)$，在y方向上有$y_j\left(j=1,2,\dots ,u+1\right)$，在t方向上有$t_k\left(k=1,2,\dots ,n+1\right)$，因此为了方便表达将$T\left(x,y,t\right)$离散后的值记作$T_{i,j}^k$。∆x，∆y，∆t为x,y,t方向的步长，$x_i=i∆x，y_j=j∆y，t_k=k∆t$。
使用有限差分近似$\frac{{\partial }^{2}T}{\partial x^2}$：