偏微分方程数值解法1105.pptVIP

偏微分方程数值解法（第四讲） 2014 年 11 月 05 日星期三 1. 一阶线性常系数双曲型方程组 考虑一阶线性常系数方程组 其中 为常 系 数矩阵。 ) 1 ( 0 0 ? ? ? ? ? ? ? ? t R x x U A t U p , )] , ( , ), , ( ), , ( [ ) , ( 2 1 T p t x u t x u t x u t x U U ? ? ? p p R A ? ? 定义： 如果 A 的特征值是实的，并存在非奇异矩阵 S 使得 其中 为 A 的特征值，称（ 1 ）为双曲型方程组， 如果 A 是对称阵，则称（ 1 ）为对称双曲型方程组，如果 A 的特征 值是实的且互不相同，则称（ 1 ）为严格双曲型方程组。 给定初始条件 与（ 1 ）构成初值问题。 , ] , , , [ 2 1 1 p diag AS S ? ? ? ? ? ? ? ? ] , , 2 , 1 ( p i i ? ? ? ) ( ) 0 , ( 0 x U x U ? 4.3 一阶双曲型方程的差分格式 常用求解该线性方程组的格式有 a. 后退的 Euler 方法 此格式不能用。 b.Lax — Friedrichs 格式 c. 蛙跳格式 ) ( 2 1 1 1 1 1 ? ? ? ? ? ? ? ? n j n j n j n j U U A h U U ? ) ( 2 ) ( 2 1 1 1 1 1 1 n j n j n j n j n j U U A h U U U ? ? ? ? ? ? ? ? ? ? ) 2 ( 2 ) ( 2 1 1 2 2 2 1 1 1 n j n j n j n j n j n j n j U U U A h U U A h U U ? ? ? ? ? ? ? ? ? ? ? ? ? 4.3.1 Lax-Friedrichs 格式 2. 二阶线性双曲型方程 为推导简单起见，令 令 因此我们有 ). ( ) , 1 ( ), ( ) , 0 ( : ), ( ) 0 , ( ), ( ) 0 , ( : 1 0 2 1 1 0 2 t g t u t g t u BC x u x u x u x u IC x u a u t xx tt ? ? ? ? ? ? ? ? ? ? ? ? x t u q u p 1 ? a x xx t tt q u p u ? ? ? ? ? ? ? ? ? ? ? x x t xt t x xx tt t p u u q q u u p ) ( 4.3.1 Lax-Friedrichs 格式 这里 A 的特征值为 -1,1 ，因此这个方程组为双曲型方程组。 边界条件 所以 所以 初始条件 x t q p q p ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 0 1 1 0 ) ( ) , 0 ( ) , 0 ( 1 t g t p t u t ? ? ? ? ? ? ? ? ? 0 1 1 0 A ) ( ) , 0 ( 1 t g t u ? ) ( ) , 1 ( 2 t g t u ? ) ( ) , 1 ( ) , 1 ( 2 t g t p t u t ? ? ) ( ) 0 , ( ) 0 , ( ) ( ) 0 , ( ) 0 , ( 0 1 0 x u x u x q x u x u x p x ? ? ? ? 4.3.1 Lax-Friedrichs 格式 也可直接求解： 截断误差为 稳定性条件为 ) ( 2 2 h O ? ? 1 ? ? a 2 1 1 2 2 1 1 2 2 h u u u a u u u n j n j n j n j n j n j ? ? ? ? ? ? ? ? ? ? 3 二维问题 求解的格式有： a. 隐格式 该截断误差为 ? ? ? ? ? ?? ? ? ? ? ? ? ?? ? ? ? ? ? ? ? ? ? ? ? ? ? y x y x g y x u t y x y u b x u a t u , , ) , ( ) 0 , , ( 0 , 0 ? ? ? ? 0 2 2 1 1 1 1 1 1 1 1 1 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? n jm n jm n m j n m j n jm n jm u u h b u u h a u u ? ) ( 2 h O ? ? 令 代入格式后可得， 格式是无条件稳定的。 b. Crank-Nicolson 格式 该截断误差为 格式为无条件稳定。 1 ) , ( 2 ? k G ? ) ( 2 1 mh k jh k i n n jm e v v ? ? h k ib h k