一类显式k阶线性k步法基本公式.pdfVIP

2013年 2月 计 算 数 学 第 35卷第 1期 Feb．．2013 MATHEMATICA NUM ERICA SINICA Vo1．35．No．1 一 类显式的k阶线性 步法基本公式 ) 刘冬兵 )马亮亮 (攀枝花学院数学与计算机学院，四川攀枝花617000) 摘 要 本文给出了一类比Adams--Bashforth方法的局部截断误差主项系数小和绝对稳定区间大的 显式k阶线性 k步法基本公式．作者求出了公式的分数形式的精确系数，阶数和局部截断误差主 项系数，给出了3-9步公式的绝对稳定区间，构造了由新公式的4阶显式公式和一个同阶隐式基 本公式组合而成的特殊预估 一校正方法，它的绝对稳定区间大于预估公式而且等于校正公式，比 著名的Adams-Bashforth-M·oulton预估校正方法的绝对稳定区间大，最后用数值试验对结果进 行了验证，适合于求解常微分方程初值问题． 关键词：线性多步法；稳定性；预估一校正方法；数值试验 MR (2000)主题分类：65L05，65L20 1．引 言 线性多步法是求解常微分方程 (组)的初值问题 【：(0)：YoYn 6 、(1．1) 的十分重要的方法，用显式方法可以直接求解，如著名的Adams-．Bashforth公式，由于误差较大 且稳定性不好 [1-8l，一般不常用．用具有绝对稳定区间或 ()稳定的隐式线性多步法求解时， 误差和稳定性都有很大地改进，如文献 9【]中的Adams-．Moulton公式、Gear(BDF)公式、[10】 中的YL公式和 [11]中新的隐式线性多步法公式等，但一般不能直接求解，而需要解一个非线 性代数方程．一般采用迭代法求解这个代数方程，这时需要一个迭代的初值，称为预估值，通 常用显式公式作计算提供，然后用隐式公式生成的迭代公式进行迭代，组成预估校正方法．文 献 1【2—18]中研究了选取显式公式作预估，隐式公式作校正，刘冬兵和杨大地在 [19]中研究并 给出了4阶的 Adams--Bashforth--Moulton预估 一校正公式，其绝对稳定区间为 (一1．2848，0)； 4阶的Adams--Bashforth-Hamming预估 一校正公式，其绝对稳定区间为 (一1．5644，0)；4阶的 Milne-．Hamming预估 一校正公式，其绝对稳定区间为 (一0．5，0)；4阶的Adams．Bashforth—Gear 预估 一校正公式，其绝对稳定区间为 (一1，0)等著名的预估 一校正方法的绝对稳定区间．杨大 地和刘冬兵在 2【0]中给出了线性多步法基本公式的定义，按照该定义基本公式只有有限个，实 际应用中的线性多步法公式大多数是基本公式，在一般文献中报道的数十个常用的重要算法， 包括上面提到的方法大部分都是属于线性多步法基本公式，而属于收敛的显式线性多步法基 本公式非常少，只有 Adams．．Bashforth公式、NystrSm公式和 4阶的Milne公式等少数几类 (个)，本文在此基础上进一步研究，给出了新的一类 k步 k阶显式的线性多步法基本公式． 2012年 4月23日收到． )基金项目：国家 自然科学基金项 目、攀枝花学院自然科学基金项目(2012YB21)和攀枝花学院 数学与计算机学院自然科学预研基金项目(Y2011--05)资助． 。)通讯作者：刘冬兵 (1972一)，男，研究方向：微分方程数值解．E-mail：liudongbingcqu@163．com． 32 计 算 数 学 2013危 1．1．线性多步法的绝对稳定性 1．求解常微分方程初值问题的线件 k步法可用如下一般形式表示