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

差分格式稳定性及数值效应比较实验一 实验目的：1.以一阶线性双曲线方程为例，使用Matlab工具分析4种差分格式的误差。2.了解4种差分格式的稳定性。二 实验问题：对于一阶线性双曲型方程：ut+aux=0,u0,x=fx=1 , x00 , x0取a=1,2,4, h=0.1, =0.08, 对不同的差分格式（迎风格式，Lax-Friedrichs格式，Lax-Wendroff格式，修正迎风格式）及不同的a值进行迭代计算。通过将计算结果与精确解来进行比较，来讨论分析差分格式的稳定性。三 实验原理：1. 迎风格式：这种格式的基本思想是简单的，就是在双曲型方程中关于空间偏导数用在特征线方向一侧的单边差商来代替，格式如下：ujn+1-ujn+aujn-uj-1nh=0,a0ujn+1-ujn+auj+1n-ujnh=0,a0ujn+1=1+aujn-auj+1n,a02. Lax-Friedrichs格式：ujn+1-12（uj+1n+uj-1n）+auj+1n-uj-1n2h=0运算格式：ujn+1=121-auj+1n+12(1+a)uj-1n3. Lax-Wendroff格式：这种格式构造是采用 Taylor 级数展开和微分方程本身得到运算格式：ujn+1=a2a-1uj+1n+1+a1-aujn+a2a+1uj-1n4. 修正迎风格式（目标点范围跟踪格式）：ujn+1=auj-a-1n+（1-a）uj-an其中a是a取整数部分，a=a-a。根据之后的理论分析可以得到这是一个无条件稳定结构。四 四种格式理论分析：通过求差分格式的增长因子G( , k)，来判定差分格式是否稳定。 1. 迎风格式：记ujn=vneijkh，则vn+1eijkh=vneijkh-a(vneijkh-vnei(j-1)kh)，得vn+1=vn1-a(1-e-ikh)，即G,k=1-a1-e-ikh=1-a1-coskh-asinkh。所以G,k2=1-4a(1-a)sinh2kh2。则在a1，满足von Neumann条件，格式稳定。以下格式用相同方法求解稳定性条件。2. Lax-Friedrichs格式：G,k=coskh-iasinkh，在a1时稳定。3. Lax-Wendroff格式：G,k=1-2a22sin2kh2-iasinkh，在a1时稳定。4. 修正迎风格式（目标点范围跟踪格式）：G,k=e-iakh1-a(1-e-ikh)，其中e-iakh=1，1-a(1-e-ikh)1的成立条件为a1。而a1恒成立，故格式无条件稳定。五 实验结果：a=1（a=0.8） 迎风格式 Lax-Friedrichs格式Lax-Wendroff格式 修正迎风格式a=2（a=1.6） 迎风格式 Lax-Friedrichs格式Lax-Wendroff格式 修正迎风格式a=4（a=3.2） 迎风格式 Lax-Friedrichs格式Lax-Wendroff格式 修正迎风格式六 总结：本次实验，通过4种差分格式求解T=4时的解并与解析解画图比较，可以看出：(1) a=1（a=0.81）时，迎风格式，Lax-Friedrichs格式，修正迎风格式的计算结果与解析解近似情况较好，而Lax-Wendroff格式则在间断点处出现了波前波，形成双波现象，这符合Lax-Wendroff格式为二阶迭代格式的性质。(2) a=2（a=1.61）时，迎风格式，Lax-Friedrichs格式，Lax-Wendroff格式都出现了比较强烈的震荡。这三种震荡中，Lax-Friedrichs格式震荡较小，迎风格式与Lax-Wendroff格式的震荡则较大。与之相对应的是修正迎风格式，保持着稳定的性质。(3) a=4（a=3.21）时，迎风格式，Lax-Friedrichs格式，Lax-Wendroff格式的震荡更加强烈。修正迎风格式则仍然保持着原有的稳定性不变。由上得出，稳定性对差分格式求解偏微分方程有重大意义。一个差分格式是否好，是否可用，首先要判定它是否稳定并找到稳定