有限差分法的基本知识.pdfVIP

第二章有限差分法的基本知识 1、差分方程 1、差分方程 2、截断误差 2、截断误差 3、收敛性 3、收敛性 4、稳定性 4、稳定性 §1 差分方程 有限差分法和有限元法是解偏微分 方程的两种主要的数值方法。由于数字 电子计算机只能存储有限个数据和作有 限次运算，所以任何一种适用于计算机 解题的方法，都必须把连续问题离散化 ，最终化成有限形式的代数方程组。 以最简单一维对流方程 为例，引入用差分 方法求偏微分方程数值 解的一些概念，说明求 解过程和原理。 考虑对流方程的初值问 题 ⎧⎪∂u +c ∂u 0 x ∈R , t 0 ⎨∂t ∂x (1.1) ⎪ u x f x x ∈R ⎩ ( ,0) ( ) 1 区域的剖分(区域的离散化） 网格剖分可以采用两组 平行于x 轴和t轴 的直线形成的网覆盖区 域 ，它们的交点称 Ω 为网格结点（节点） t t nτ n 0,1,2, n x x j jh j 0,±1,±2, 节点(x j , tn )记为(j , n). 间距h 0称为空间步长，间距 τ 0称为时间步长。 t (x , t ) j n 0 x 1 微分方程离散(差分方程） 高等数学中，我们学习过Taylor公式： 设f (x ) 在x 的某个邻域 U(x ,δ) 内具有直 0 0 到n +1阶的导数，则 ∀x ∈U(x ,δ) 有 0 f (x ) f (x0 ) +f ′(x0 )(x −x0 ) ++ f (n) (x ) n 0 (x −x ) +R (x ) 0 n n! n R (x )是余项，且R (x ) o((x −x ) ) (x →x ). n n 0 0 设 是方程(1.1)的解，对于任何节点 ( , )， 的微商 u j n u 与差商之间的关系式 ( , ) ( , ) u x j tn+1 −u x j tn ∂ ( , ) (τ), (向前差商） (1.2) u x t +o