有限差分方法基础-研究生课程.pptVIP

第一节差分原理及逼近误差/逼近误差（4/4）将与的Taylor展开式相加可得（1-18）第一节差分原理及逼近误差/非均匀步长（1/3）在有些情况下要求自变量的增量本身是变化的，如图2-1中的，是不相等的，相应的差分和差商就是不等距的。Ox

图1-1非均匀步长差分3．非均匀步长一阶向后差商 一阶中心差商 （1-23）第一节差分原理及逼近误差/非均匀步长（2/3）图1-2均匀和非均匀网格实例1第一节差分原理及逼近误差/非均匀步长（3/3）图1-3均匀和非均匀网格实例2（2-1）第二节差分方程、截断误差和相容性/差分方程（1/3）从上节所述可知，差分相应于微分，差商相应于导数。只不过差分和差商是用有限形式表示的，而微分和导数则是以极限形式表示的。如果将微分方程中的导数用相应的差商近似代替，就可得到有限形式的差分方程。现以对流方程为例，列出对应的差分方程。0102图2-1差分网格第二节差分方程、截断误差和相容性/差分方程（2/3）若时间导数用一阶向前差商近似代替，即01空间导数用一阶中心差商近似代替，即02则在03点的对流方程就可近似地写作04（2-2）05（2-3）06（2-4）07第二节差分方程、截断误差和相容性/差分方程（3/3）08第二节差分方程、截断误差和相容性/截断误差（1/6）按照前面关于逼近误差的分析知道，用时间向前差商代替时间导数时的误差为,用空间中心差商代替空间导数时的误差为，因而对流方程与对应的差分方程之间也存在一个误差，它是 这也可由Taylor展开得到。因为（2-5）（2-6）第二节差分方程、截断误差和相容性/截断误差（2/6）一个与时间相关的物理问题，应用微分方程表示时，还必须给定初始条件，从而形成一个完整的初值问题。对流方程的初值问题为 这里为某已知函数。同样，差分方程也必须有初始条件： 初始条件是一种定解条件。如果是初边值问题，定解条件中还应有适当的边界条件。差分方程和其定解条件一起，称为相应微分方程定解问题的差分格式。（2-7）（2-8）第二节差分方程、截断误差和相容性/截断误差（3/6）FTCS格式（2-9）FTFS格式（2-10）（2-11）FTBS格式FTCS(b)FTFS(c)FTBS图2-2差分格式第二节差分方程、截断误差和相容性/截断误差（5/6）第二节差分方程、截断误差和相容性/截断误差（6/6）1FTCS格式的截断误差为2FTFS和FTBS格式的截断误差为3（2-12）4（2-13）53种格式对6都有一阶精度。7第二节差分方程、截断误差和相容性/相容性（1/3）一般说来，若微分方程为 其中D是微分算子，f是已知函数，而对应的差分方程为其中是差分算子，则截断误差为这里（2-14） （2-15）（2-16）如果当、时，差分方程的截断误差的某种范数也趋近于零，即则表明从截断误差的角度来看，此差分方程是能用来逼近微分方程的，通常称这样的差分方程和相应的微分方程相容（一致）。如果当、（2-17）第二节差分方程、截断误差和相容性/相容性（2/3）1若微分问题的定解条件为2其中B是微分算子，g是已知函数，而对应的差分问题的定解条件为3其中4是差分算子，则截断误差为5（2-18）6（2-19）7（2-20）8 第二节差分方程、截断误差和相容性/相容性（3/3）1只有方程相容，定解条件也相容，即2和3整个问题才相容。4（2-21）5无条件相容条件相容6以上3种格式都属于一阶精度、二层、相容、显式格式。7第三节收敛性与稳定性/收敛性（1/6）所谓相容性，是指当自变量的步长趋于零时，差分格式与微分问题的截断误差的范数是否趋于零，从而可看出是否能用此差分格式来逼近微分问题。然而，方程（无论是微分方程或是差分方程）是物理问题的数学表达形式，其目的是为了借助数学的手段来求问题的解。因此，除了必须要求差分格式能逼近微分方程和定解条件（表明这两种数学表达方法在形式上是一致的）外，还进一步要求差分格式的解（精确解）与微分方程定解问题的解（精确解）是一致的（表明这两种数学表达方法的最终结果是一致的）。即当步长趋于零时，要求差分格式的解