波动方程的差分逼近.docVIP

第五章 双曲型方程的有限差分法 4.1 波动方程的差分逼近 1. 特征 针对波动方程 （1） 其初值条件为 其中是常数。其相应的特征方程为characteristic equation 即 得到两个特征方向：characteristic direction （3） 解（3），得到两族直线： 2. 显格式 取空间步长及时间步长，用两族平行直线two family of parallel lines 作矩形网格rectangle。在对方程（1）离散，得到 （5.1） 初始条件为 （5.2） （5.3） （5.1）式逼近的截断误差为。由于(5.3)式逼近截断误差为，因此对（5.3）的逼近可作适当改进。（5）可显示算出各网点的值。 （5.1）简化后可以写成 （6） 针对混合问题： 此时取空间步长及时间步长，同样建立离散格式（5），针对边值条件，可给出离散的边值条件 3. 稳定性分析 为了利用Fourier方法，令，将（1）化成一阶偏微分方程组： （7） 再令，则（7）变为 （8） 令及 则（8）变为 因此，差分方程（5）可写成 （10） 按照Fourier方法，设 ， 代入（10），消去公因子common factor和，得到 即 其中 为增长矩阵，其特征方程为 （14） 其根按模小于1的充要条件是absolute value of root （15） 即，此为必要条件。 另外，（14）的两个根为 因，故 设，从而 另一方面， 其F模为 根据上一章的结论，必须且只须存在常数使不等式成立： 而若，则当时，，不等式不能成立，因此应要求。因此格式（10）稳定的充要条件是。 稳定性条件的几何解释geometrical interpretation 从（6）可以看出，依赖于前两层的，，，，这四个值又依赖于，，，，，和，，，依次类推，可得到依赖于 因此称轴上位于区间的网格点成为差分解的依存域dependent region。它是轴上被过的两条直线 、 切下的区间所覆盖的网域。 而过的两条特征线为 因此，差分方程稳定的必要条件即差分解的依存域必须包含微分方程解的依存域。 利用依存域的概念，可以证明，当时，差分解不收敛。 如上图，意味着微分方程解的依存域大于差分解的依存域。固定，当步长变小，但不变，则依存域，不变，若改变，上的初值，但上的初值不变，则可取不同的值，而当时，是一串确定的数列，它不可能收敛到不同的。 当时，差分方程稳定，因而差分解收敛。 4. 隐格式 目的：得到绝对稳定的差分格式。 利用第层，层、层的中心差商的权平均去逼近，得到下列差分格式： （18） 其中是参数。当就是显格式。我们感兴趣的是，此时差分格式等价于 （19） 其增长矩阵为 可以证明的特征值按绝对值等于1,且是酉矩阵，因此，从而矩阵族一致有界，即（19）绝对稳定。 综合习题 1. 证明：方程