差分格式稳定性及数值效应比较实验.docVIP

差分格式稳定性及数值效应比较实验 5090719044 张 赟 F0907102 一 实验目的： 1.以一阶线性双曲线方程为例，使用Matlab工具分析4种差分格式的误差。 2.了解4种差分格式的稳定性 二 实验问题： 对于一阶线性双曲型方程： 取a=1,2,4, h=0.1, τ=0.08, 对不同的差分格式（迎风格式，Lax-Friedrichs格式，Lax-Wendroff格式，修正迎风格式）及不同的a值进行迭代计算。通过将计算结果与精确解来进行比较，来讨论分析差分格式的稳定性。 三 实验原理： 迎风格式： 这种格式的基本思想是简单的，就是在双曲型方程中关于空间偏导数用在特征线方向一侧的单边差商来代替，格式如下： 运算格式： Lax-Friedrichs格式： 运算格式： Lax-Wendroff格式： 这种格式构造是采用 Taylor 级数展开和微分方程本身得到，运算格式： 修正迎风格式（目标点范围跟踪格式）： 其中是取整数部分，=。根据之后的理论分析可以得到这是一个无条件稳定结构。 四 四种格式理论分析： 通过求差分格式的增长因子G(τ , k)，来判定差分格式是否稳定。 迎风格式： 记，则， 得, 即。 所以。 则在，满足von Neumann条件，格式稳定。 以下格式用相同方法