有限差分法的Crank-Nicolson方法

有限差分法是数值计算中一种常用的方法，它将偏微分方程中的连续问题离散化为有限个离散问题，用数值方法求解。Crank-Nicolson方法是有限差分法的一种，它是一种隐式的求解方法。

在使用Crank-Nicolson方法求解偏微分方程时，先将空间离散化，将时间离散化。对于时间上的每一个离散点，都可以使用数值方法求解空间上的离散点。

假设求解的偏微分方程为：

$$\frac{\partial u}{\partial t}=D\frac{\partial^2 u}{\partial x^2}$$

其中，$D$为常数，$u(x,t)$是未知函数。将空间离散化，假设$x$方向上的网格点数为$N$，则可以将空间上的离散点表示为$u_j^n$，其中$j=1,2,...,N$，$n$为时间离散点。假设网格点的间距为$\Delta x$，时间步长为$\Delta t$，则可以将偏微分方程离散化为：

$$\frac{u_j^{n+1}-u_j^n}{\Delta t}=\frac{D}{2\Delta x^2}(u_{j+1}^{n+1}-2u_j^{n+1}+u_{j-1}^{n+1}+u_{j+1}^n-2u_j^n+u_{j-1}^n)$$

这个式子就是Crank-Nicolson方法的离散化表达式。由于该表达式中既包含了$n+1$时刻的值，也包含了$n$时刻的值，因此它是一种隐式的求解方法。可以将该表达式写成矩阵形式：

$$(I-\frac{D\Delta t}{2\Delta x^2}L)u^{n+1}=(I+\frac{D\Delta t}{2\Delta x^2}L)u^n$$

其中，$L$是一个$N\times N$的三对角矩阵，表示离散化后的偏微分方程中的系数矩阵，$I$是一个$N\times N$的单位矩阵，$u^n$表示$n$时刻的离散点向量，$u^{n+1}$表示$n+1$时刻的离散点向量。

Crank-Nicolson方法是一种比较稳定的方法，它可以避免一些数值解法中可能出现的数值震荡问题。但是，由于它是一种隐式的求解方法，每次求解需要解线性方程组，因此计算量比较大