3抛物型方程.pptVIP

§1最简差分格式 表1 r=1/2时Richardson格式的误差传播 §2 稳定性与收敛性 2.2 判别稳定性的直接法 §3 Fourier方法 4.1 ADI法 4.2 预-校法 ) / ( , ) 2 sin 2 (sin 4 1 ( 2 , 2 , 2 2 2 1 ) ( h a r v h h r v l p l p e v u n n y x i n n jk k j t b a p b p a b a = + - = = = = + + ，得 代到两端并消去公因子 取通项 解 的稳定性。 的向前差分方程： 物方程 考虑逼近带低阶项的抛 例 k j k j k j k j k j h u u u u u x u t u + - = - ? ? = ? ? - + + 2 1 1 1 2 2 2 3 t 。 一致有界的充要条件是 则知增长因子 代入 以 解 k p jh i k k j k j k j k j k j r G h r x G e v u u u r ru u 2 1 2 sin 4 1 ) , ( , ) 2 1 ( 2 1 1 1 ￡ - = = + - + = - + + a t a 得增长因子 令 ） ） （ ） ） （ ，则有 消去过渡阵 有逆 式左端的差分算子 都是正规阵，故以上二 和 注意 ） （ ） （ ） （ ） （ 我们将它改写成 的稳定性 为了检验格式 , ) 5 . 5 ( , 2 ( 2 2 ( 2 , ) 2 . 5 ( , 2 2 ) 2 . 5 ( , 2 2 , ) 2 . 5 ( ) ( 2 2 1 2 2 2 1 2 2 2 1 2 1 2 2 2 1 2 k j y x i n n jk n jk y x n jk y x n jk y x n jk x n jk y n jk y n jk x e v u u r I r I u r I r I u u r I u r I u r I u r I b a d d d d d d d d d d + + + + + + = + + = - - - - ￠ + = - ￠ + = - ) 7 . 5 ( 2 ( ) 6 . 5 ( , 1 2 ( 2 ) 5 . 5 ( 1 2 2 1 2 1 2 2 2 1 2 2 t d t d d t d d n jk n jk y n jk n jk n jk n jk y x n jk n jk y x u u r I u u u u h u u r I r I - - = - + = - - - + + + + ） ，使 引进过渡阵 ） （ ） ） （ 改写成 我们将 高维问题 2 ( 1 2 2 1 t d t n jk n jk y n jk n jk u u r I u u - - = - + + ） ） ） （ ） ） （ ) 5 . 5 ( , 2 ( 2 2 ( 2 2 2 1 2 2 n jk y x n jk y x u r I r I u r I r I d d d d + + + = - - * 第五章抛物型方程的有限差分法§1 最简差分格式 §2 稳定性与收敛性 §3 Fourier方法 §4 ADI法和预-校法 0 1 2 3 N-1 l 1 2 x y 考察Richardson格式的稳定性 ? -4 -3 -2 -1 0 ?1 2 3 ?4 0 ? ? ? ? ? ? ? ? ? 1 ? ? ? ? ? ? ? ? ? 2 ? ? ? ? ? ? ? ? ? 3 ? ? ? ? ? ? ? ? ? 4 ? ? ? ? ? ? ? ? ? 5 ? ? ? ? ? ? ? ? ? 6 ? ? ? ? ? ? ? ? ? ? ? ? -4 ? -3 ? -2 ? -1 ? 0 ? 1 ? 2 ? 3 ? 4 ? 0 ? ? ? ? ? ? ? ? ? ? 1 ? ? ? ? ? ? ? ? ? ? 2 ? ? ?