二阶非线性微分方程两种混合数值格式.pdfVIP

二阶非线性微分方程的两种混合数值格式1 吴奇，化存才 云南师范大学数学学院，昆明（650092 ） E-mail ：Cuncai-hua@ 摘 要：在构造算法时，通过同时运用数值积分和差商近似导数的方法，给出解二阶非线性 微分方程的两种混合数值格式，分析得到了它们的局部截断误差和整体误差；数值算例表明： 这两种方法与Euler 方法和Runge-Kutta 法相比，都具有一定的优势。 关键词：二阶非线性微分方程；数值解法；误差分析 1. 引言 对于形如y′′ f (y , y, t ) 的二阶非线性微分方程，不论是求解其初值问题，还是边值问′ 题都是非常困难的，仅有少量特殊二阶微分方程可直接求出解析解，而更多的方程很难求 得解的明确表达式，甚至根本无法求解。这就促使学者们寻找某种“逼近”方法去近似的得 到微分方程的解，如从渐近分析的角度去构造和研究它们的渐近解，这可参看文献[1]；又 如从数值的角度将方程的初边值问题转化为在给定的有限个离散节点上的近似差分或有限 元方程求解，这方面的方法和论述在文献[2]的第一章中有详细介绍。此外，对于一些特殊 的二阶微分方程，还有一些更为深入的数值方法研究，如：对于方程形如y′′ f (x, y)，P. 亨利西在[3]中对线性多步法进行了研究，尔后，Lambert 和Watson 在[4]中针对S-C 方法提 出了对称线性多步法。在本文中，我们在传统的将二阶方程化为一阶方程组的基础上，对 两个方程的一阶数分别用数值积分和差商进行近似处理，给出求解二阶非线性微分方程的 两种混合计算格式，分析得到了它们的局部截断误差和整体误差。 2. 整体误差为一阶的混合数值格式 我们考虑如下的二阶非线性微分方程的初值问题： ′′ ( ′ ) y f y , y, t , ⎧ ⎨ ′ ′ (1) ( ) ( ) y t 0 y0 , y t 0 y0 , ⎩ ′ ′ 其中f (y , y, t )对t 充分光滑。作变量替换令y z ，将(1)化为如下的一阶微分方程组： ′ y z, ⎧ ⎪ ′ ( ) ⎨z f z, y, t , (2) ⎪ ( ) ( ) y t 0 y 0 ,z t 0 z0 , ⎩ 传统的数值方法是用Euler 法或Runge-Kutta 法对(2)进行求解。 在这里，为得到新的计算格式，我们先对(2)中第一个方程的两边在区间[t m-1，t m ]求 积分，并利用数值积分中的梯形公式得： h