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

差分格式稳定性及数值效应比较实验 一实验目的: 1?以一阶线性双曲线方程为例，使用 2?了解4种差分格式的稳定性。 Matlab工具分析4种差分格式的误差 二实验问题: 对于一阶线性双曲型方程： ut + aux = 0, 如=心）=唸; 取a=1,2,4, h=0.1, t =0.08,对不同的差分格式（迎风格式， Lax-Friedrichs格 式，Lax-Wendroff格式，修正迎风格式）及不同的 a值进行迭代计算。通过将计 算结果与精确解来进行比较，来讨论分析差分格式的稳定性。 三实验原理: 迎风格式： 这种格式的基本思想是简单的，就是在双曲型方程中关于空间偏导数用在特 征线方向一侧的单边差商来代替，格式如下： TOC \o "1-5" \h \z ..n+1 1tn ..n ..n U -11； U 一片?[ J——+a丄■宀= I h ..n+l .11 ..n .,11 I h I h 运算格式： 二（1 +训曙為加＜0 Lax-Friedrich格式： 护-扌婢1 +略）I嚅-％ “ i +a^i^=0 运算格式： 申=（1 ■就用土 （1 + a喀1 Lax-Wendroff格式： 这种格式构造是采用Taylor级数展开和微分方程本身得到 运算格式： = —(aX -l)u^1 + (l + aX)(l-aX)l^ + —(aX + l)u^x 小 是d 小 是d取整数部分，fa^ 其中。根据之后的理论分析可以 其中 得到这是一个无条件稳定结构。 四四种格式理论分析： 通过求差分格式的增长因子G(t , k),来判定差分格式是否稳定 迎风格式： 记:F，一乙、；二：「.卩「三爲；■ ：：：「-,：； 得「-「. ：「「, G(T,k) = l-a^(l-e_lkh)=l-a^(l-眺 kh 卜小 sinldi 所以 |訓!:;;:卜二 1 f : I - .S '、，::二好二 则在以下格式用相同方法求解稳定性条件。^<1 ,满足von Neumann 则在 以下格式用相同方法求解稳定性条件。 Lax-Friedrich格式： G（;k） = cosli-iahinkli，在 J时稳定。 Lax-Wendroff格式: G（t*）二1 — 2讥sin2学—i讥血kh，在？ X < 1时稳定 2 * ■亠 修正迎风格式（目标点范围跟踪格式）： G（调，咖卜岡（1-严）], 其中「皿恤卜1,卜訓1一严壮1的成立条件为< ] 而\《1恒成立，故格式无条件稳定。 五实验结果: 迎风格式 Lax-Friedrich格式 Lax-Wendroff格 式 Lax-Wendroff格 式 修正迎风格式 a=22 Lax-Wendroff格式 修正迎风格式 迎风格式 Lax-Friedrich格式 Lax-Wendroff格式 修正迎风格式 六总结： 本次实验，通过4种差分格式求解T=4时的解并与解析解画图比较，可以看 出： ⑴a=1 (a入=0.8v 1 )时，迎风格式，Lax-Friedrichs格式，修正迎风格式 的计算 结果与解析解近似情况较好，而 Lax-We ndroff格式则在间断点处出现了波前波， 形成双波现象，这符合Lax-We ndroff格式为二阶迭代格式的性质。 ⑵ a=2 (a 入=1.6> 1)时，迎风格式，Lax-Friedrichs格式，Lax-Wendroff 格式都 出现了比较强烈的震荡。这三种震荡中， Lax-Friedrichs格式震荡较小，迎风 格式与Lax-Wendroff格式的震荡则较大。与之相对应的是修正迎风格式，保 持着稳定的性质。 ⑶ a=4 (a 入=3.2> 1)时，迎风格式，Lax-Friedrichs格式，Lax-Wendroff 格式的 震荡更加强烈。修正迎风格式则仍然保持着原有的稳定性不变。 由上得出，稳定性对差分格式求解偏微分方程有重大意义。一个差分格式是 否好，是否可用，首先要判定它是否稳定并找到稳定性条件。修正迎风格式强大 的稳定性在解决一阶线性双曲线方程中有着很强的实用价值。 七程序: 迎风格式： 2 3 4 5 6 function yingfeng (a? h』tj mi rut j jm= ^jnaxx-ninx) /h; T=4; p-t/h; n-I/t; ul=ones(m+n+lJ 1); ul fn+l: Jn+n+ll=O; g io n 12 13 14 15 16 17 u2=ul; for i=l:1:n for j=i+l: 1 ：Jit+n+l u2(j)=a?p*ul (j-l>+(l-a*p)*u2(j); -ertd ul=u2; -end yl=u2(n+l：R+n+l)； xl=ninx:h:ikaKX; pl