高等数值分析结课论文--余帆.docVIP

《高 等 数 值 分 析》 结 课 论 文 学生姓名：余 帆 学 号：专业班级：数学研14-1班 2015年 1月7日问题重述 考虑如下双曲问题： 这里 求出精确解。 写出两种不同的数值格式并给出分析。 给定，计算结果。画出时刻的解并与精确解比较。 当时，计算时的最大误差，并画出log-log图像。 对你的结果进行解释。 问题分析 此问题即要求我们构造至少两种差分格式来计算此双曲问题，我们可以构造迎风格式、Lax-Friedrichs格式和Lax-Wendroff格式来进行计算比较。 问题的解答 利用特征线的方法我们可以计算出此问题的精确解为 我们首先进行网格剖分： 对于此问题的数值格式，我们可以构造如下三种不同的格式： 迎风格式 此格式的稳定性条件为。对于周期边值条件，我们可以利用来逼近。 Lax-Friedrichs格式 此格式的稳定性条件为。对于周期边值条件， 我们可以利用来逼近。 Lax-Wendroff格式 此格式的稳定性条件为。对于周期边值条件， 我们可以利用来逼近。 结果分析 当时，下面给出当时这三种格式的解和精确解。 当时，下面给出三种格式的解： 当时，下面给出三种格式的解： 当时，下面给出三种格式的解： 从上面的结果可以看出，精确解和数值解都是像波动一样在周期性的运动。 接下来我们分析这些格式的误差。 这里要说明的是当时网比，此时三种格式都不稳定，其数值结果会产生严重的数值振荡，在此我们不探究，只给出振荡的一个例子，如下 当时，对于我们给出了这三种格式的最大误差的log-log图及其对应的误差表格。 迎风格式： Maximum error Rate 0.10 18.2304 0.05 13.0428 0.4830925 0.025 8.5937 0.6019022 0.0125 3.8575 1.1556132 Lax-Friedrichs格式： Maximum error Rate 0.10 47.3863 0.05 26.9530 0.8140241 0.025 14.3051 0.9139162 0.0125 5.4916 1.3812311 Lax-Wendroff格式： Maximum error Rate 0.10 19.5787 0.05 11.8613 0.7230228 0.025 7.0696 0.7465616 0.0125 3.1685 1.1590584 从上面误差分析的结果可以看出，这三种格式的收敛阶几乎都在1左右，这与其格式是相符的。当空间步长和时间步长接近时，数值模拟的结果较好。 第5章 改进的分裂正定有限元格式的超收敛性分析