偏微风方程数值解法.docVIP

偏微分方程数值解法作业 PAGE 3 偏微风方程数值解法 第2章 有限差分方法的基本概念 习题(p.44)第二题讨论对流方程,a0 的差风格式 的截断误差及稳定 解：(1)先计算截断误差：用Taylor级数在展开： ＝ = = 可见： 所以，截断误差为 （2）稳定性： 将代入上式得： = 得增长因子； = = 时对于任何值，所以此格式无条件稳定。 双曲型方程的差分方法 例：解方程在t=1.0时的解。 波动显式格式： （2）波动隐式格式 波动显式和波动隐式两个格式都是二阶精度的格式，但显示格式是条件稳定的，而隐式格式是无条件稳定的。而从计算结果来看，隐式格式更加精确些。 (p.81) 6. 讨论Crank-Nicolson型格式 的截断误差及稳定性。 解：(1)此格式的截断误差为： 可见： 此格式的截断误差为： （2）稳定性： 将代入上式得： = 所以，此格式无条件稳定。 数值计算例子(p.53)考虑初值问题 其中 取h=0.01, ，用迎风格式，Lax-Friedrichs格式和 Lax-Wendroff格式以及Beam-Warming格式计算时，方程的数值解。 对前两个格式（迎风格式和Lax-Friedrichs格式）把解抹平了。迎风格式为一阶精度，而Lax-Friedrichs格式的精度比一阶高，但低于二阶。而从解的结果来看，迎风格式对此题符合的更好些，它抹平的程度小些。也许对于其他题Lax-Friedrichs格式会好些。Lax-Wendroff格式以及Beam-Warming格式均为二阶精度，两者均出现了振荡，Lax-Wendroff格式在波前出现振荡，而Beam-Warming格式在低部出现振荡，这些现象是这些格式的正常现象。另外我用了蛙跳格式对此问题也进行了计算，蛙跳格式在波前的振荡非常厉害，蛙跳格式的计算结果与初值有很大的关系。另外，它的稳定性要求也高了些。 第4章 抛物型方程的有限差分法 10．对于初边值问题 用向前差分格式，向后差分格式，及Crank-Nicolson格式以及Du Fort-Frankel格式求解，取h=0.1,为0.1和0.5进行计算，并在t=0.1时与准确解比较。 解：(1)向前差分格式： t=0.100, lamda=0.1 t=0.100, lamda=0.1 0.100 0 0 0.200 0 0 0.300 0 0 0.400 0 0 0.500 0 0 0.600 0 0 0.700 0 0 0.800 0 0 0.900 0 0t=0.100,lamda=0.5 0.100 0 0t=0.100,lamda=0.5 0.100 0 0 0.200 0 0 0.300 0 0 0.400 0 0 0.500 0 0 0.600 0 0 0.700 0 0 0.800 0 0 0.900 0 0向前