差分方法的稳定性.docVIP

差分方法的稳固性 差分方法的稳固性 差分方法的稳固性 差分方法的稳固性 实验内容 对于一阶线性双曲线型方程： u u t 0,x0,1,t0,T x ux,0 u0x 此中初值u0 1,x 0 x 0 0,x 取空间长度h=0.01,对于不一样的差分格式（迎风格式，Lax—Friedrichs格式,Lax—Wendroff格式，Beam—Warming格式以及蛙跳格式)及不一样的网格比（时间长度 与空间长度比）进行迭代计算。经过将计算结果与精准解进行比较，来讨 h 论和分析差分格式的稳固性。 算法思想与步骤 2.1迎风格式 这类格式的基本思想是简单的，就是在双曲型方程中对于空间偏导数用在特色线方向一侧的单边差商来取代，格式以下： unj 1 unj a unj unj1 0,a 0 h unj 1 unj a unj 1 unj 0,a 0 h unj 1 (1 a)unj aunj1,a0 运算格式： unj 1 1 aunj aunj1,a0 2.2Lax—Friedrichs 格式 unj 1 1 unj1 unj1 n n 1 2 a uj1 uj 0 2h 运算格式： unj 1 1 1 aunj1 1 1 aunj1 2 2 2。3Lax—Wendroff 格式 这类格式结构采纳Taylor 级数睁开和微分方程自己获取 运 算 格 式 ： unj 1 a a1unj11a1aunj a a1unj1 2 2 2.4Bean—Warming格式（二阶迎风格式） 借助于双曲型方程的解在特色线上为常数这一事实,能够结构出多种差分格式。 设在ttn时间层上网格点A,B,C和D上u的值已给定，要计算出在ttn1时间层上网格点P上的u的值。假设C。F。L条件建立,过P点特色线与BC交于点Q，故微分方程解的性质知uPuQ. 对于uQ： ①用B,C两点值进行线性插值，获取的是迎风格式; ②用B，D两点值进行线性插值,获取的是Lax-Friedrichs格式; ③用B，C和D三点值进行抛物型插值，获取的是Lax-Wendroff格式. 假如我们采纳A,BC三点来进行抛物型插值，能够获取 unj11a1aunja2aunj1a1aunj1 22 这就是Beam—Warming格式。 2。5蛙跳格式 unj 1 unj 1 a unj1 unj1 0 2 2h 运算格式： unj 1 unj 1 a unj1 unj 1 因为它是个三层格式，需要先用一个二层格式计算出 t 那一层的值u1j。为了 保持精度的阶数同样，一般我们用 Lax—Wendroff格式或Beam—Warming格式。 2.6目标点范围追踪格式(迎风格式的改良） unj1aunja11aunja 此中a是a取整数部分，aaa.下边的分析将会获取这是一个 无条件稳固结构。 数据分析与作图 3.1迎风格式 1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 00.10.20.30.40.50.60.70.80.91 0.5 1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 00.10.20.30.40.50.60.70.80.91 1 400 300 200 100 0 -100 -200 -300 -400 -500 00.10.20.30.40.50.60.70.80.91 1.1 10 x10 8 6 4 2 0 -2 -4 -6 -8 00.10.20.30.40.50.60.70.80.91 2 稳固性分析： 记unjvneijkh，则vn1eijkhvneijkhavneijkhvneij1kh，得 vn1vn1a1eikh 即 G,k 1a 1eikh 1a1coskhaisinkh 。 则在a 1时，有G ,k1,格式稳固。 3.2Lax—Friedrichs 格式 1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 00.10.20.30.40.50.60.70.80.91 0.5 1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 00.10.20.30.40.50.60.70.80.91 1 10 5 0 -5 -10 -15 00.10.20.30.40.50.60.70.80.91 1.1 6 x10 4 3 2 1 0 -1 -2 -3 -4 00.10.20.30.40.50.60.70.80.91 2 稳固性分析： G,kcoskhihasinkh 则在a1时稳固。 3。3Lax—Wendroff格式 1.4 1.2 1 0.8 0.6 0.4 0.2 0 00.10.20.30.40.50.60.70.80.91 0.5 1