微分方程数值解第五章答案微分方程数值解第五章答案.pdf

第五章 双曲型方程的数值解 - 1 - 第五章 ⎧ 1, x 0, u u ⎪ ∂ ∂ 1. 对初值问题 + 0, u( x,0) ⎨1/ 2, x 0, 试分别用左偏心格式、 ∂ ∂ t x ⎪ ⎩0, x 0. LW 格式计算其数值解 uk , k=1,2,3,4, 取τ/ h 1/ 2 . 解: 矩形网格剖分区域. 取空间步长 , 时间步长 的矩形网格剖分 h τ 区 域 , 用 节 点 (j ,k ) 表 示 坐 标 点 (x , t ) ( jh, kτ) , j 0, ±1, ±2,...; j k k 0,1,2,3,4. k k ∂u ∂u ⎡ ⎤ ⎡ ⎤ （1）左偏心格式： + 0 ，在 t 上用向前差商，x 上用向后 ⎢ ⎥ ⎢ ⎥ ∂t ∂x ⎣ ⎦ ⎣ ⎦ j j k +1 k k k u u u u − − 差商，得 j j + j j −1 0 / 1τ/ 2h ，因为 ，整理得到 τ h k +1 1 k 1 k u u u + j 2 j −1 2 j 1, ⎧ 0 j 0 ⎪ 1/ 2, 0 j 把已知条件离散成uj ⎨ ，则可以根据下一层求上一层的值得