抛物型方程的一族精高度恒稳格式.pdfVIP

维普资讯 第 21卷 第 2期 华 侨 大 学 学 报 (自 然 科 学 版 ) V0}．21 No 2 2000年 4月 JournalofHuaqiaoUniversity (NaturalScience) A Dr 2000 文章编号 1000-5013(2000)02—0124—05 抛物型方程的一族高精度恒稳格式 f4 f2 黄浪扬 曾文平 (华侨太学管理信息科学系 ．泉州 36201l 摘要 对 二 阶抛物 型方程 = ．构造 丁一旗新 的三层 隐式差分 格式 (在特殊情况下是两层) 它 们含有 非 负参数 ， 和 -其截 断误差至少可达 D((At)4-( ))．对三层格式 ，在条件 ≥ ≥ o·≤音及 口-4 十 一1之下绝对稳定 ．特别地．在蕞件 =0．-二 或 口1一 ，-二0之下成为 两层 不舍参效 的臆 式格 式 ．且 也是绝对稳 定 的 这 些格 式均 可用追赶法求解 ．在谖 格 式 中．选 取 适 当的参数 ．可得抛物 型方程初曲值 差分格 式 中的高精度格式 关键词 三堕地塑型左显r鲁堕堡室，高精度 壁蔓董坌整 中田分类 号 O 241．82 文献标识码 A ，7 本文考虑 如下 的抛物 型方程初边 值 问题 ，即 害=筹， 0， ￡0， “(，0)一 _，()， 0 -r ． f “(O，)一 l(f)， “( ，)= ：(f)， 0．J 解上述 问题 ，人们 已建立 了各种各样的差分格式 。‘．本文构造 了一族三层 (在特殊情况下是 两层 )含参数 、高精度 、绝对稳定 、三对角 型 的隐式差 分格式 ，其截断误差 阶至少为 O((At)。 ( ))(特殊情 况下还可 以提高)，它包 含 了文 (1～4]中的所有 高精度恒 稳格式 ． 1 差分格式的构造 现设 问题 (1)的解 “ )充分光滑 ，时 间步长为 ，空间步长为 △ ，并设 r一 为网格 比，在网点 (靠 ，)处 的网格 函数 “( ，f)记为 “：．利用 Taylor展开不难验证下列数值微分公 式为 ‘ 一 一 + 十 + 一．一， ：，= c 一 一 +一-·， ㈣ ：～= c～ ～ 一 一一 。+-．一． ㈨ 收稿 日期 1999一l1—18 作者简介 黄 浪扬 (1974一)，男．助教 维普资讯 第 2期 黄浪扬等 ：抛物型方程 的一族高精度恒稳格式 l25 引入记 号 △^“ 一 l+ 10￡ + ， ：一 “+}一 2 + “ 】， (5) 则 有 A(J：一1 + (雾 十 蕞fJ+0cc ． e 将 式 (6)用于方程 (1)，便 得到 下面 3个差 分格式为 = —— —一 一丽忐 m’