一类隐式的k+1阶线性k步法.pdf

2013年 11月 计 算 数 学 第 35卷第 4期 NOV．，2013 MATHEMATICA NUMERICA SINICA Vo1．35，No．4 一 类隐式的 +1阶线性 步法术1) 刘冬兵)马亮亮 (攀枝花学院数学与计算机学院，四川攀枝花617000) 摘 要 本文首先给出了一类比Adams—Moulton方法的绝对稳定区间大的隐式 k+1阶线性 k步 法基本公式．求出了3-9步新公式的分数形式的精确系数，阶数，局部截断误差主项系数和绝对 稳定区间，然后构造了由4阶隐式新公式和同阶显式 Nystr6m 公式组合而成的预估一校正方法， 比著名的 Adams—Bashforth-Moulton和 Nystr6m-Adams-Moulton预估校正方法的绝对稳定 区间大，最后用对比数值试验对结果进行了验证． 关键词：线性多步法；稳定性；预估一校正方法；数值试验 MR (2000)主题分类：65L05，65L20 1．引 言 线性多步法是求解常微分方程 (组)的初值问题： {【yx(o)Yo) (1．) 的重要经典算法之一 [1--23】，也可以用于求解偏微分方程 问题．线性多步法分为显式和隐式， 杨大地和刘冬兵在 1【3】中定义了线性多步法的基本公式，按照该定义基本公式只有有限个，属 于收敛的隐式线性多步法基本公式比显式的方法要多一些，实际应用中的线性多步法大多数 是基本公式，例如 Adams．Moulton公式、Gear(BDF)公式、YL公式 1【0]和 4阶的Hamming 公式等属于隐式方法．本文在此基础上进一步研究，给出了新的一类k+1阶隐式的线性 步 法基本公式，比同k步Adams—Moulton方法的绝对稳定区间都要大，新公式的发现丰富了具 有绝对稳定区间的隐式线性多步法基本公式的类型． 1．1．线性多步法的绝对稳定性 1．求解常微分方程初值问题的线性 k步法可用如下一般形式表示 (1．2) 其中0。，岛均为实常数，且OLk=1，I0I+I3／ol≠0 2013年 3月 15日收到． )基金项 目：国家 自然科学基金 资助项 目，四川省科技厅 自然科学基金项 目(2013JY0125)和攀 枝花学院自然科学基金项 目(2012YB21)资助． 。)通讯作者：刘冬兵 (1972一)，男，湖南宁乡人，研究方向：微分方程数值解．E—mail：liudongbingcqu@163．com． 394 计 算 数 学 2013住 定义 1．1[14]．一个数值方法以定步长 h求解试验方程 fY，： 【y(xo)、 Yo (1．3) 得到线性差分方程的解为Y．当咒_÷。。，若有Y ．÷0，则称该方法对步长h是绝对稳定 的 2．线性多步法的绝对稳定性可用其稳定多项式进行研究． 记 k k p(∈)=∑ i ()=∑ (1．4) i=0 i=0 定义 1．2_14]．分别称p(∈)，(∈)为线性多步法的第一特征多项式和第二特征多项式