[工学]chapter 5 -20440926 有限差分法.ppt

1 有限差分法 (Finite Difference Method) 1 简介 2 差分原理及逼近误差 3 差分方程 4 截断误差和相容性 5 收敛性与稳定性 6 Lax等价定理 7 基于差分的分析系统的组成 8 有限差分分析软件 1 有限差分法 (Finite Difference Method) 1 简介 2 差分原理及逼近误差 3 差分方程 4 截断误差和相容性 5 收敛性与稳定性 6 Lax等价定理 7 基于差分的分析系统的组成 8 有限差分分析软件 1 有限差分法 (Finite Difference Method) 1 简介 2 差分原理及逼近误差 3 差分方程 4 截断误差和相容性 5 收敛性与稳定性 6 Lax等价定理 7 基于差分的分析系统的组成 8 有限差分分析软件 2 1 简介 有限差分法：用差分代替微分，是一种微分方程 和积分微分方程数值解的方法。 有限差分法以变量离散取值后对应的函数值来近 似微分方程中独立变量的连续取值。 3 1 简介 基本思想： 把连续的定解区域用有限个离散点构成的网格来代替，这些 离散点称作网格的节点； 把连续定解区域上的连续变量的函数用在网格上定义的离散 变量函数来近似； 以Taylor 级数展开等方法把控制方程和定解条件中的微商 用网格节点上的差商代替进行离散，于是原微分方程和定解 条件就近似地代之以代数方程组，即有限差分方程组，解此 方程组就可以得到原问题在离散点上的近似解。 然后再利用插值方法便可以从离散解得到定解问题在整个区 域上的近似解。 4 1 简介 主要内容： 如何根据问题的特点将定界区域作网格划分； 如何把原微分方程离散化为差分方程组； 如何解此方程组； 为了保证计算过程的可行性和计算结果的准确性，还需 从理论上分析差分方程组的性态，包括解的唯一性、存 在性和差分格式的相容性、收敛性和稳定性等问题。 5 有限差分法的具体操作分为两个部分： （1）用差分代替微分方程中的微分，将连续变化的变量离散化， 从而得到差分方程组的数学形式； （2）求解差分方程组。 在第一步中，我们通过所谓的网格分割法，将函数定义域分成 大量相邻而不重合的子区域。通常采用的是规则的分割方式。 这样可以便于计算机自动实现和减少计算的复杂性。 在第二步中，数值求解的关键就是要应用适当的计算方法，求 得特定问题在所有这些节点上的离散近似值。 1 简介 = lim dy dx f (x + Δx) − f (x) Δx = lim Δx→0 Δy Δx→0 Δx ) 2 差分原理及逼近误差 差分原理 设有 x 的解析函数y = f (x，从微分学的角度函数 y 对 x 的导数为： 导数 差商 其中，dy、dx 分别是函数及自变量的微分， 相应的Δy、 Δx 分别称为函数及自变量的差分。 6 Δx)− f (x − Δx) 7 一阶差分 差分原理 在导数的定义中，Δx 是以任意方式趋近于零的，因而 Δx可正可负。在差分方法中，Δx总是取一小的正数。 这样一来，与微分对应的差分可以有3种形式： 向前差分：Δy = f (x + Δx) − f (x) 向后差分：Δy = f (x) − f (x − Δx) 1 2 1 2 中心差分：Δy = f (x + 8 Δ2y = Δ(Δy) = Δ[f (x + Δx)− f (x)] = Δf (x + Δx) − Δf (x) =[f (x + 2Δx)− f (x + Δx)]−[f (x + Δx)− f (x)] = f (x + 2Δx) − 2 f (x + Δx) + f (x) 以向前差分为例可得到二阶差分： 差分原理 9 ] = Δ[Δ(Δn−2y) } = Δ{Δ⋅⋅⋅[Δ(Δy)] } = Δ{Δ⋅⋅⋅[Δ( f (x+Δx)− f (x)] 差分原理 依此类推，任何阶差分都可由其低一阶的差分得 到，例如 n 阶向前差分为： Δny = Δ(Δn−1y) … f (x + Δx) − f (x − Δx) 10 差分原理 函数的差分与自变量的差分之比，即为函数对自 变量的差商，如： = = = Δy Δx Δy Δx Δy Δx f (x + Δx) − f (x) Δx f (x) − f (x − Δx) Δx 1 1 2 2 Δx = Δy