下载文档原格式
第四章有限差分方法4.1引言有限差分法：数值求解常微分方程或偏微分方程的方法。
物理学和其他学科领域的许多问题在被分析研究之后, 往往可以归结为常微分方程或偏微分方程的求解问题。
一般说来，处理一个特定的物理问题，除了需要知道它满足的数学方程外，还应当同时知道这个问题的定解条件，然后才能设计出行之有效的计算方法来求解。
有限差分法以变量离散取值后对应的函数值来近似微分方程中独立变量的连续取值。
在有限差分方法中，我们放弃了微分方程中独立变量可以取连续值的特征，而关注独立变量离散取值后对应的函数值。
但是从原则上说，这种方法仍然可以达到任意满意的计算精度。
因为方程的连续数值解可以通过减小独立变量离散取值的间格，或者通过离散点上的函数值插值计算来近似得到。
这种方法是随着计算机的诞生和应用而发展起来的。
其计算格式和程序的设计都比较直观和简单，因而，它的实际应用已经构成了计算数学和计算物理的重要组成部分。
有限差分法的具体操作分为两个部分：(1)用差分代替微分方程中的微分，将连续变化的变量离散化，从而得到差分方程组的数学形式； (2)求解差分方程组。
在第一步中，我们通过所谓的网络分割法，将函数定义域分成大量相邻而不重合的子区域。
通常采用的是规则的分割方式。
这样可以便于计算机自动实现和减少计算的复杂性。
网络线划分的交点称为节点。
若与某个节点P 相邻的节点都是定义在场域内的节点，则P 点称为正则节点；反之，若节点P 有处在定义域外的相邻节点，则P 点称为非正则节点。
在第二步中，数值求解的关键就是要应用适当的计算方法，求得特定问题在所有这些节点上的离散近似值。
有限差分法的差分格式：一个函数在x 点上的一阶和二阶微商，可以近似地用它所临近的两点上的函数值的差分来表示。
如对一个单变量函数f(x)，x 为定义在区间[a,b]的连续变量。
以步长h=Δx 将[a,b]区间离散化，我们得到一系列节点x = a , x = x + h , x = x + h = a + 212132Δx , ..., x = x + h = b , 然后求出 f(x)在这些点上的近似值。
有限差分法原理有限差分法（Finite Difference Method）是一种常见的数值计算方法，广泛应用于工程、物理、地质等领域的数值模拟和求解偏微分方程。
它的原理是将连续的微分方程转化为离散的差分方程，通过对网格节点上的数值进行逼近，从而求解微分方程的数值解。
在本文中，我们将介绍有限差分法的基本原理及其在实际问题中的应用。
首先，我们来看一维热传导方程的数值求解。
假设我们要求解一个长为L的均匀材料棒上的温度分布，其热传导方程可以写为：\[ \frac{\partial u}{\partial t} = \alpha\frac{\partial^2 u}{\partial x^2} \]其中，u(x, t)表示位置x上的温度分布，t表示时间，α为热扩散系数。
为了使用有限差分法求解这个方程，我们需要将空间和时间进行离散化。
假设我们在空间上取N个网格点，将材料棒分为N个小区间，每个小区间的长度为Δx。
在时间上也进行离散化，取时间步长为Δt。
这样，我们可以用u_i^n来表示位置为x_i的温度在时间t_n的值。
将热传导方程在离散点上进行近似，我们可以得到如下的差分格式：\[ \frac{u_i^{n+1} u_i^n}{\Delta t} = \alpha\frac{u_{i+1}^n 2u_i^n + u_{i-1}^n}{(\Delta x)^2} \]通过对时间和空间上的离散点进行迭代计算，我们可以逐步求解出温度在空间上的分布随时间的演化。
这就是有限差分法的基本原理。
除了一维热传导方程，有限差分法还可以应用于更加复杂的偏微分方程，比如二维热传导方程、波动方程、扩散方程等。
在这些情况下，我们需要在空间上取二维甚至三维的网格点，并相应地修改差分格式。
有限差分法的优点在于它简单易实现，而且可以直接应用于一般的偏微分方程，因此在实际工程和科学计算中得到了广泛的应用。
需要指出的是，有限差分法也有一些局限性。
