第五章离散偏微分方程

1、 2013届毕业生外文文献翻译学 院 理学院 专 业 数学与应用数学 姓 名 孟文静 学 号 1 指导教师 张兆忠 第五章 离散偏微分方程5.1 介绍有多种方法求得的偏微分方程的数值解，然而本质上这些都涉及到构建一些近似值或获得求解的方程组的近似值，得到近似的离散系统的解决方案。这个过程通常有两个部分。第一个涉及发展近似值的一个未知的函数，它的导数以及其他涉及使用或应用这些算法计算出所需的数值的数量。虽然这过于简单化，但它提供了一个方便的研究技术的起点来求解偏微分方程。离散偏微分方程的最直观直接的方法是代替部分采用有限差分近似的导数，从而构造一个泛函微分方程（FDE）代替偏微分方程。本质上偏微

2、分方程的好处是被替换为一个泛函微分方程来产生一组代数方程组，原则上可以解决，以构建偏微分方程的一个近似解。完成离散化后，构建一个合适的计算方法来解决泛函微分方程的任务，具有挑战性，的确可能影响导数如何先被离散化。在开发一个解决偏微分方程的误差的数值方法时，需要一个数值解的质量评估。例如，确定提出的方法是否收敛，即，是否在有限的邻域内，近似解可以尽可能接近的偏微分方程的精确解，这是开发新方法的过程的一个重要部分。此外还有有效的计算方法和更多新的可平行性算法来检验和评估。本章介绍了这些思路，专注于开发取代有限差异的导数的技术，还展示了在开发一个简单的数值方案的背景下如何构建的一些偏微分方程的近似解

3、。其他的内容将在后面章节介绍。5.2 在实数域上构建有限差分近似对一个函数的导数的差分近似的建设可以直接找到多项式，使该多项式插值在指定的一列点上，。考虑使用二次多项式 差值在点和，即，使用局部坐标系，让和，然后 这个含三个未知数的三个方程可化为 解得，得到的差值多项式。为了得到导数的近似值，对求导 然后在（局部坐标系）计算该表达式，得 . 该方法构建方便，但无法说明近似误差，然而，而一阶导数的有限差分近似值可以使用第二章描述的的一些标准方法获得，在这些近似值中得到一个理想的误差是必要的。5.3 截断误差关于点，是足够光滑的，则的泰勒级数展开式可以表示为 把右侧的函数移到的左侧，然后除以，可得

4、标准差商（用来定义的导数），即， 令，大括号里的项消失，由介绍的导数定义得 考虑对于任意固定值，用差商代替产生的误差，即， 从的右边的项中提取，这个项可以视为代表一个线性算子作用于，即， 代表构造出了微分算子和近似线性算子之间的误差。误差或在作为截断误差被求出，并且一般来说当函数足够光滑时，通过一个离散的线性算子，截断误差被定义为代表线性算子中的误差。因此，定义5.1. （截断误差）让是一个到线性微分算子定义的有最大值的邻域上离散逼近；如果存在常数和那么 并且那么就有一个的截断误差。一个离散逼近到一个有的截断误差的线性算子，称为命令。此外，这个近似值可被称为阶精度。例如，在中对任意光滑的和充分

5、小的，最右侧的项被右侧绝对值中的首项所支配，即， （因为封闭，有界域和它的所有导数是连续有界的）因此，截断误差，当变量无限趋近于零时，可由给出的方法定义的收敛。 为得到一个有效的截断误差时，发现对任何固定值，值在中不线性依赖。这个行为只能在足够小时实现，即，当远远小于时。的确，这种行为在对计算近似值时不会数据性观察，若函数的二阶导数在有效范围内比其他阶导数远远小时，即，当 5.4 用特勒多项式构建导数在点附近再次考虑的泰勒级数，即， 已经表明，可以为重新安排以提供一个前向差分近似，即， 由命令的截断误差得到的一个前向差分近似公式 同样的，在方向上的点周围可扩展的泰勒级数 获得所谓的在点的向后差

6、分近似， 减得 求解得中心差分近似 只使用的两个点计算得到三个不同的离散表达式，像，和表示的一样通过扩展的泰勒多项式得的一阶导数。前两个式子，向前差分和向后差分，在计算模板中涉及使用两点，并且有截断误差。中心差分近似仅使用关于点的点和的中点，实际上使用了三个点，然而点因为它的系数为零而在式子中消失。表明中心差分近似在中为二阶精度，并且精确度的提高来自于向前和向后泰勒级数相加后偶数项的抵销。这个近似值可以看做向前和向后差分近似的平均值，即， 一般来说这个结果是正确的，即，向前和向后差分没有相应的中心差分精确。当然如果是这种情况，这个问题不必麻烦的使用片面的差分。答案更依赖于构建解决偏微分方程的问题，而不是只通过检验对一些方程的导数构建离散近似值这种狭窄问题来回答。对数值导数的构建不是唯一的，并且虽然中心差分近似的精度可能比相应的向前或向后差分近似值大，但是构建构建向前和向后差分近似可能需要一些截断误差，所以某种程度上，在这些近似值之间最根本的区别是定向的偏差，即，当评估导数时，向前或向后差分或没有差分被介绍。例5.1. 求解的向前近似差分，其中是二阶精度。解： 考虑关于点的在点，和的的泰勒级数展开式，那么,和结果，计算通过划分并重新安排，得即如果的导数在中使用的近似值的一个邻域内有界，例如如果，那么的所有导数可