含余项的常微分方程数值计算方法.pdfVIP

首都师范大学学报(自然科学版) Journalof Normal No．1 第36卷第1期 C印ital UniVersity Science Feb．．2015 2叭5年2月 (Natural Edition) 含余项的常微分方程数值计算方法 孟纯青 (山东省济南中学，山东济南250001) 摘 要 本文提出了一种新的微分方程计算方法的构建思路，并构建了一些公式，分析了算法特性，并给出了计算实 例．本文提出的算法具有稳定性好，精算精度高的特点，是一类具有使用价值的新方法． 关键词：含余项的常微分方程计算方法，常微分方程，余项． 中图分类号：0241．81 a；。，· 1 一阶偏差分的表达形式 2 含余项的常微分方程数值计算方法 以往的微分方程计算方法都会忽略余项，这无 疑会影响算法的精度，本文提出了包含余项的微分 设存在常微分方程初值问题 方程计算方法，并构建了一些公式．由于余项可以使 用一阶偏差分来表示，所以下面本文给出一阶偏差 』警2妣)，) (2．1) 分的表达式． 【)，(戈o)=yo 设存在函数八z，)，)，为了方便八石i，)，i)也写成 以算，y)为已知函数，对于八茁，)，)，通常假设， 八菇，y)在某个区域内连续，且关于)，三阶连续可微． ，，符号a以戈。，y。)=aZ，a以算i，yi)=aZ分别表 设^为步长．现在构建带余项的数值计算方法． 示八戈，y)在点(戈。，yi)处的一阶偏差分，设存在常数 r0，a以xi，y。)，a以石i，)，。)的具体表达式为： 包含余项的梯形公式的形式如下： a：+¨／石i，)，i)=／Ix“．，，)，i)一／l(石i，)，。)(1．1) y…=yi+要[以戈i，，，i)+以戈Ⅲ，)，…)]一 a：一¨／省i，y。)=／Ixi，)，i)一／xi一，，y。)(1．2) a：+¨／I戈i，)，i)=／xi，yi+，)一／石i，)，。)(1．3) 篙)，，『『(孝)孝∈(％％。) (2．2) a：一¨八菇i，y。)=八x。，y。)一“石i，yi一，)(1．4) (2．2)中最后一项称为梯形公式的余项，也就 (1．1)，(1．3)是向前偏差分，(1．2)，(1．4)是 是截断误差．关于(2．2)的推导可参阅相关的计算 向后偏差分．当r=1时称为一阶节点偏差分表达 方法的书籍，例如本文的参考文献[1]． 式，反之称为一阶非节点偏差分表达式．下面把 ay’八算i，)，i)，ay’以石i，yi)统一写成ay厂(菇i，yi)， 设y，，，(亭)：三坐i旦，则(2．2)变成： 凡 或者秽’，，a：+r’八戈。，y。)，a：_r’八髫i