二阶常微分方程一类配置方法-计算数学专业论文.docxVIP

1 绪论 1.1 选题背景 二阶常微分方程被大范围地运用于力学，电学，工程技术领域，如何精确地求 解二阶微分方程是人们努力追求的目标。然而只有极少部分的微分方程能够得到精 确的解析解，大部分方程是不能通过积分求解的。有些方程即使有解析解，但表达 式过于繁琐，这也将在实际运用中失去意义。于是，寻求方程的达到精度要求的数 值解是十分重要的。在求解大型微分方程组的过程中，发现了所谓刚性现象。运用 具有有限稳定区域的一些数值方法求解这类方程时，只有步长取得非常小，计算才 会成功。而步长取小后，计算步数增加了，以至于增加到方程本身的刚性比的程度。 于是我们所寻找的数值方法应该具有尽可能高的阶以及尽可能大的稳定区域。 求解二阶常微分方程初值问题，一种可行的做法是先将其先化成一个一阶微分 方程组，然后再用一阶方程中已有的各种数值算法进行求解。这样的说法叫做间接 方法，关于其收敛性，稳定性，误差的分析完全平行于一阶方程的结果。不过，在 1925 年，Nystr?m[10]提出将二阶方程独立于一阶方程进行专门研究，他构造了一种 如下所示的新方法,现在称之为 Runge-Kutta-Nystr?m 方法： i n i n  s ∑ ij n j j Y = y + c hy + h2 a j =1 f (t + c h,Y ), n+1 n s ∑ j n j j y = y + h d j =1 f (t + c h,Y ), n+1 n n s ∑ j n j j y = y + hy + h2 b j =1 f (t + c h,Y ). 1978 年，Hairer 和 Wanner[11]给出了 RKN 方法的阶条件。同年，Battin 在[14] 中独立地给出了 8 阶 RKN 方法的阶条件。1982 年，Hairer[16]给出了一个具体的 10 阶方法。后面的章节中，我们将看到方法的稳定函数是一个二阶矩阵，这样我们确 定稳定区域就比较困难。1980 年，L.Kramarz[9]研究了基于直接配置的 RKN 方法的 稳定性,他提出了一种相对简单的方法来确定方法的稳定区域，他指出凡使用对称配 置点的配置方法，至少存在一个有限的稳定区域。1986 年，M M Chawla[18]研究了 两步四级的 P-稳定的 RKN 方法。1990 年，Sharp,Fine 和 Burrage 在文[40]中研究了 2 级和 3 级对角隐式 RKN 方法，并且分析了这两类方法的阶条件和 R 稳定性。1991 年，P.J.Van Der Houwen[15]等人研究了直接配置方法与间接配置方法，发现对于采 用相同的配置点所分别构成的两类配置方法，只要满足一定条件，直接配置比间接 配置级阶升高一阶。1999 年，Kazufumi Ozawa 在文[34]中研究了一类解决二阶微分 方程周期初值问题的变系数四级隐式 RKN 方法，在他所考虑的方法中，系数是关于 频率与步长的函数。2001 年，Beatrice Paternoster[15]中研究了一般的两步 RKN 方法 及其 P-稳定性。 2009 年，Gonzalez 等人在文[4]中提出了一类新的解一阶刚性微分方程的龙格 库塔配置方法，称之为 SAFERK 方法。这种方法的配置点取自如下方程的根： m(x) =  2s + 1(P * (x) ? P  * (x)) + α  2s ?1(P  * (x) ? P  * (x)) = 0, s s?2 s?1 s?3 其中 P * (x) = 1 P (x) ， P (x) 是位移勒让德多项式，即 n 2n + 1 n k P (x) = 2k + 1 d k  (xk (x ?1)k ). k k! dxk 当α = 0 时，m(x) 的零点就是 Lobatto IIIA 方法的节点。当 s ≥ 4 时，多项式 m(x) 在[0,1]上与所有次数不超过 s ? 4 的多项式关于权函数 ω(x) = 1正交，从而由 m(x) 的零 点 s s {ci }i=1 确定的插值求积公式{bi ,ci }i=1 至少具有 2s ? 3 阶代数精度，当然前提是这些零 点都是[0,1]上互不相同的实数。又由于 c1 = 0 ，cs = 1 都是 m(x) = 0 的根，所以由此生 成的方法具有首级显式，刚性精确的特点。对于 s 级的该方法与 s ?1级的 Radau IIA 方法，计算量是相当的，两种方法的阶都是 2s ? 3 ，但是前者的级阶比后者高一阶。 作者指出，这种方法还是强 A 稳定的。 1.2 本文的主要工作 本文主要就是把 Gonzalez 等人在 2009 年提出的新方法运用到二阶微分方程