差分方法的稳定性.docVIP

文档可能无法思考全面，请浏览后下载! 文档可能无法思考全面，请浏览后下载! PAGE / NUMPAGES 文档可能无法思考全面，请浏览后下载! 差分方法的稳定性 1.实验内容 对于一阶线性双曲线型方程： 其中初值 取空间长度h=0.01，对于不同的差分格式（迎风格式，Lax-Friedrichs格式，Lax-Wendroff格式，Beam-Warming格式以及蛙跳格式）及不同的网格比（时间长度与空间长度比）进行迭代计算。通过将计算结果与精确解进行比较，来讨论和分析差分格式的稳定性。 2.算法思想与步骤 2.1迎风格式 这种格式的基本思想是简单的，就是在双曲型方程中关于空间偏导数用在特征线方向一侧的单边差商来代替，格式如下： 运算格式： 2.2 Lax-Friedrichs格式 运算格式： 2.3 Lax-Wendroff格式 这种格式构造采用Taylor级数展开和微分方程本身得到 运算格式： 2.4 Bean-Warming格式（二阶迎风格式） 借助于双曲型方程的解在特征线上为常数这一事实，可以构造出多种差分格式。设在时间层上网格点A,B,C和D上的值已给定，要计算出在时间层上网格点P上的的值。假定C.F.L条件成立，过P点特征线与BC交于点Q，故微分方程解的性质知。 对于： 用B,C两点值进行线性插值，得到的是迎风格式; 用B,D两点值进行线性插值，得到的是Lax-Friedrichs格式； 用B,C和D三点值进行抛物型插值，得到的是Lax-Wendroff格式。 如果我们采用A,BC三点来进行抛物型插值，可以得到 这就是Beam-Warming格式。 2.5 蛙跳格式 运算格式： 由于它是个三层格式，需要先用一个二层格式计算出那一层的值。为了保持精度的阶数相同，一般我们用Lax-Wendroff格式或Beam-Warming格式。 2.6 目标点范围跟踪格式（迎风格式的改进） 其中是取整数部分，。下面的分析将会得到这是一个无条件稳定结构。 3.数据分析与作图 3.1迎风格式 稳定性分析： 记，则，得 即。 则在时，有，格式稳定。 3.2 Lax-Friedrichs格式 稳定性分析： 则在时稳定。 3.3 Lax-Wendroff格式 稳定性分析： 则在时稳定。 3.4 Beam-Warming格式 稳定性分析： 则在时稳定。 3.5 蛙跳格式 稳定性分析： 命 令 则 则时稳定。 3.6 目标点范围跟踪格式 稳定性分析： ，其中，的成立条件为。然而恒成立，故无条件稳定。