关于稳定性分析的附注.docVIP

第 第 PAGE 9 页 共 NUMPAGES 9 页 20** 年 秋 季学期研究生课程考核 （读书报告、研究报告） 考核科目 : 偏微分方程数值解法 学生所在院（系） : 理学院数学系 学生所在学科 : 数学 学 生 姓 名 : 学 号 : 学 生 类 别 : 考核结果 阅卷人 关于稳定性分析的附注 摘要 本篇论文从运用Fourier方法研究有限差分格式稳定性的出发，具体介绍了如何运用Fourier方法得到差分格式的稳定性情况，并给出了一些判别有限差分格式稳定性的准则，如von Neumann条件，并列举了一个具体的应用Fourier方法判断有限差分格式稳定性的例子。然后，又介绍了两个研究有限差分格式稳定性的方法，Saul’ev算法和分组显式方法。 关键字：Fourier方法，稳定性分析，差分格式 Abstract This paper from the use of Fourier method to study the stability of finite difference scheme, specifically describes how to use the Fourier method to get the stability of the difference scheme, and gives some criteria for the stability of the finite difference scheme, such as Neumann von conditions, and cited a specific application of Fourier method to determine the stability of finite difference scheme example. Then, it introduces two methods to study the stability of the finite difference scheme, ev Saul algorithm and grouping explicit method. Keywords: Fourier method, stability analysis, difference scheme 1 前言 在《偏微分方程数值解法》的课程中，我们已经跟着老师学习过如何分析差分格式的稳定性问题。特别是在第二章“有限差分格式的基本概念”这一章中，我们学过用矩阵方法对初值问题的一些差分格式的稳定性进行过分析。此外也可以用能量不等式方法来进行差分格式的稳定性分析。一般说来，这样做较为麻烦。在实际应用中，仍可以用Fourier方法来进行分析。但必须注意到，Fourier方法仅适用于线性常系数初值问题。对于一般的初值问题，我们假定边界处理是准确的，再用Fourier方法来作近似分析。这种近似性方法是简单的，但对于许多实际问题来说，可对计算提供很大的帮助。因此，在本篇论文中，将先回顾一下有限差分格式稳定性的Fourier方法，然后介绍两种分析差分格式稳定性的算法，分别是Saul’ev算法和分组显式方法。 2 运用Fourier方法研究有限差分格式的稳定性 2.1 Fourier方法 以对流方程的初值问题[1] 的差分格式 （2.1） 为例进行讨论。注意到（2.1）式可以写成 （2.2） 其中的解以及初值只是在网格点上有意义。 为了应有Fourier方法进行讨论，必须扩充这些函数的定义域，是的它们在整个实轴上都有定义。令 这样，对任意都有定义了。这里使用大写字母仅为区别微分方程初值问题中使用的小写字母。由此，（2.2）式可以写为 （2.3） 显然，上式对有意义。而且对任意，上式也是有意义的。对（2.3）式两边用Fourier积分来表示，可以得到 由此得出 （2.4） 上面推导方法可以推广到一半形式的差分格式，可以得到 （2.5） 上式中的称为增长因子。显然，差分格式（2.1）的增长因子为 上式等号右边中可以通过和来表示。 由于增长因子不依赖与时间层，因此有（2.5）式可以得出 如果增长因子的任意次幂是一致有界的，并设其界为，则应用Parseval等式有 再次应用Parseval等式有 由的定义可知有 由此得到，常系数的差分格式是稳定的。同样地应用Parseval等式，可以证明，如果差分格式为常系数，那么差分格式的稳定性可以推出其增长因子的任意次幂是一致有界的。这样就得到了下面的重要结论。 常系数