rungekuttanystrm方法的若干新性质与阶结果-some new properties and order results of rungekkutta str m method.docxVIP

第一章引言二阶常微分方程初值问题常出现在科学与工程中的许多领域（如天体力学、量子力学、理论物理与化学等），研究其数值方法是重要而有意义的。数值求解二阶常微分方程有两条途径：一条是先将其化为一阶方程，然后再用相应数值方法求解，但这样会破坏二阶方程本身所具有的特征结构及方程所代表的物理意义；另一条途径是直接对二阶方程进行数值求解。对于有周期解或振荡解的二阶常微分方程初值问题的数值求解，人们有了大量的研究，提出了许多方法（如对称多步法、Runge?Kutta?Nystrom（简称RKN）方法、混合方法及相应的指数拟合型方法等），研究的主要方面包括精度、稳定性（如P-稳定性）、耗散性、相延迟性等（参见[3-6，11，14-16，33]）。一些可分离Hamilton系统（例如，多体问题）常可写为二阶常微分方程初值问题的形式。辛RKN方法是求解它们的有效方法，它是单步方法，具有自开始、可达高阶、保辛等优点。Suris[25]提出RKN方法保辛的代数充分条件，进一步，文[19]证明：在某些弱的条件下，这些代数条件也是必要的。自此以后，人们针对辛RKN方法的特征与构造及阶条件做了大量工作，例如文[18，19，20，23，24]，通过特殊的Nystrom带根树，系统研究了辛RKN方法的阶条件，并构造了一些有效的显式辛RKN方法。文[4]证明了一个显式RKN方法是辛的当且仅当它的伴随是显式的；文[37]构造了一族单参数2s-1阶扰动配置型隐式辛RKN方法；文[38]证明了Gauss-RKN方法是对称的且其对称性等价于辛性；文[1，28]获得了辛RKN方法的一些特征与阶性质，并通过变换技巧构造了2s、2s-1阶辛RKN方法及单隐、多隐、对角隐式RKN方法；文[15，39，40]进一步讨论了指数拟合辛RKN方法。其它进一步的工作还有文[26，30，31]等。文[23]构造了5阶、6阶显式辛RKN方法。我们要问：显式辛RKN方法能有更高的阶吗？文[26]构造了3级4阶单对角隐式辛RKN方法。同样，我们要问：单对角隐式辛RKN方法能有更高的阶吗？在本文第三章，我们证明了单对角隐式辛RKN方法和显式辛RKN方法的最高阶是6。NNN文[3，13，18，21]等表明：由一个Runge?Kutta（简称RK）方法Φ可以唯一生成一个RKN方法Φ，而RKN方法Φ也有伴随方法Φ*。我们在第二章中，考虑了由Φ生成Φ*的两条途径：Φ→Φ?→Φ*与Φ→Φ→Φ*，并证明这两条途径NN生成的Φ*是相同的。NNNNNNN文[7]等表明：若RK方法Φ是辛的，则其伴随方法Φ*也是辛的，且二者有相同的稳定函数。文[10]证明了：若RKN方法Φ是辛的，则它的伴随方法Φ*也是辛的。在第二章中，我们进一步讨论了具有辛性、对称性或P-稳定性的Φ、Φ、Φ*之间的一些关系，并证明：⑴若Φ是辛的，则由Φ生成的ΦN也是辛的。⑵若Φ是对称的，则由Φ生成的ΦN也是对称的。⑶由Φ的简化条件可以导出由Φ生成的ΦN的简化条件。NN⑷若Φ是辛的，则其伴随Φ*也是辛的（用不同于文[10]的证明方法）。NN⑸Φ和Φ*具有相同的P-稳定性。上述结论⑴、⑵提供了一条通过辛（或对称）RK方法构造辛（或对称）RKN方法的新途径。第二章Runge??Kutta??Nystrom方法的若干新性质2.1问题与方法为求解一阶常微分方程初值问题x′(t)??g(t,x)，t??t0，x(t0)??x0，(2.1.1)N1这里g:[t0,?∞)??R→?RN1，考虑s-级Runge??Kutta方法（以下简称为RK方法）：??X?xs??hag(t??ch,X)?in??∑j?1sijnjj(i??1,2,,s)，(2.1.2)??x?????n?1xn??h∑?dig(tn??cih,Xi)i?1其中h??0是积分步长，Xi≈?x(tn??cih)，xn≈?x(tn)，aij、di、ci（i,j??1,2,,s）为实常数，通常选取0≤?c1≤≤?cs()s?s≤?1。记TsTsA???aij∈R，c??(c1,c2,,cs)∈R，d??(d1,d2,,ds)∈?R，D??diag(d)，M??DA??ATD??ddT。下文中也记B??diag(b)，C??diag(c)，并设s∑?djj?1??1，di≠0s和∑aij?cij?1(i??1,,s)。此外，方法(2.1.2)可写成Butcher表的形式：，(2.1.2)′显然，当aij??0,j≥?i时，(2.1.2)是显式RK方法，否则它将是隐式的RK方法。下文中，方法(2.1.2)简记为Φ?或(c,A,d)。为求解二阶常微分方程的初值问题y?f(t,y)，y(t0)??y0，y(t0)??y0，t??t0，(2.1.3)N2这里f:[t0,?∞)??R→?RN2，考虑s-级Rung