关于波动方程的差分格式.pdfVIP

维普资讯 反 程 差 移 ，差方军子／ 钱 期 删蹦扎。 矗 G业 )F撮一一 7 关于波动方程的差分格式 一 张天德 孙传灼 (=)z z - - — — — — 一 (山东I业大学数理系 济南 250014) 1 建立差分格式 对波动方程 L—一等j一一笔一一0 (疋定解胖条寐件1千略J) L(l1．。1) 取矩形网格，r为时间步长，h为空间步长 ，设 r与 h同阶，网比，一r，^为常数．如对方程 L 一 一4。 一0，则 ，一4r／h 4o所得结论同样成立． 用 表示 u(x，f)关于 的 阶偏导数 ， 表示 “关于 t的一阶偏导数关于 z的n 一 1阶偏导数．以下的 及其各阶导数如无特别指明均表示在点 ( ，tD的取值． 构造如下形式的三层差分格式 L^U一南军(r5j+~‘一6，u+，一qUj+I‘1)一0 (1．2) 一 告 (Ⅱ， )≠o， ∈n={～N’．．．，0，…，Ⅳ} 将 (1．1)的充分光滑的解代入 (1_2)并在 ( ，)用Taylor级数展开 (1．3) 其中 D=r3+lh D一 rO+lh 并注意到u )_出m 。Ⅲ+1)!u墨赫”，则 收稿 日期 ；1993一n一03 维普资讯 284 山 东 工 业 大 学 学 报 1995年 9月 在( 有 l 4 ) ( )一 + 手 ) 中 当 + + 手((，一 一 “+r(口+c1) + (口一 一c，)U+ (口+c1) ] 三 + 兰手 [()(一一c) + T 一()一2(n，一日)U 二Ⅲ)兀罢二一；而 一-1(拍)一 (，+)U ] (1．4) Uo一 三( 一 一 )一 0 一 三(口l— bi— q) 一 0 一 六手[(口一一)+r(口J一)]一0