抛物型方程的一个实用有效的高阶差分算法.pdfVIP

第15卷 第1期 兰州工业高等专科学校学报 V01．15．No．1 2008年3月 J0哪 al of Lanzhou Polytechnic College M82．．2O08 文章编号：1009—2269(2008)01—0001—05 抛物型方程的一个实用有效的高阶差分算法 汪子莲 ，丁双平2 (1．兰州工业高等专科学校基础部，甘肃兰州730050；2．甘肃省礼县教育局，甘肃礼县 742200) 摘要：针对文献[1]的高阶差分算法，首先综合应用降维降阶法导出了抛物方程的一个不需要用 到u在边界上的二阶导数的值、计算量小、精度高的有限差分格式，并给出了该差分格式的截断误 差的表达式；然后通过引进过渡层给出了近似因式分解的交替方向隐差分格式的算法；最后是数 值例子，数值结果和理论结果是吻合的． 关 键 词：抛物型方程；高阶差分算法；差分格式 中图分类号：0 151 文献标识码：A 0 引言 1 差分格式的建立与截断误差的表 文献[1]建立了抛物型方程： 达式 a n 32n a2u ． ， 、 为了构造微分方程问题(1)一G)的有限差分 万 。 +6 +g( ，Y，‘)， 逼近，首先将求解区域西 ×(0，T)进行网格剖分 ( ，Y，t)∈n×[0，T]． (1) (为方便起见本文进行等距剖分)，选取正整数 u( ，Y，0)= ( ，Y)，( ，Y)∈n． (2) 和Ⅳ，并令 u( ，Y，t)= ( ，Y，t)， h= 1／M．r= ／N． ( ，Y)∈an，0t≤T． (3) 记 的紧交替方向格式，其中o，b是正常数，g( ，Y， i=ih，y{=jh，0≤ i，j≤M， )，声( ，Y， )， ( ，，，)是给定的函数，所得差分 t = nr，0≤ n≤ N， 格式具有0(r +h )阶精度．但求解所给的差分 n^={( ，y，)l 0≤i， ≤M}， 格式，需要用到u在边界上的二阶导数的值，计算 得到其近似解，降低了整体逼近精度，在每一时间 n。={t l 0≤n≤N}， 层上共需解4(M一1)个三对角方程组．为了克服 y={(0，．『)，(m， )l 0≤．『≤m}U{(i，0)， 上述缺陷，本文首先综合应用降维降阶法导出了 (i，m)l 1≤i≤m一1}． 该方程的一个不需要用到u在边界上的二阶导数 此外记 的值、计算量小、精度高的有限差分格式，并给出 {= 1( + +1)， + 了该差分格式的截断误差的表达式；然后通过引 进过渡层给出了近似因式分解的交替方向隐差分