本文在前人对常微分方程初值问题的线性多步法公式研究的基础上,对于线性多步法公式中基于泰勒展开的构造方法进行了探究。我们尝试使用加权平均法构造得出了一个新的公式,随后对此给出了实例进行检验,并对构造生成新公式的原先两个公式和新公式进行了对照分析。经数值实验,该方法具有可行性,某些所得新公式在求解常微分方程初值问题中具有较高稳定性和较小误差。 Following the previous literature on multi-step formulae for initial value problems of ODEs (ordi-nary differential equations), we study the construction for linear multi-step methods based on Taylor expansion in this paper. We try the weighted average method and derive a new formula. Then we check this new method in an example, and compare the original two formulae and con-structed new formula. Via such numerical experiment, this method is reliable. Some constructed new formulae can have relatively high stability and small error while solving ODEs initial value problems.
黄志远1,胡志军2,王 承3
1广州大学,广东 广州
2广西师范大学,广西 桂林
3惠州学院,广东 惠州
收稿日期:2015年10月30日;录用日期:2015年11月13日;发布日期:2015年11月18日
本文在前人对常微分方程初值问题的线性多步法公式研究的基础上,对于线性多步法公式中基于泰勒展开的构造方法进行了探究。我们尝试使用加权平均法构造得出了一个新的公式,随后对此给出了实例进行检验,并对构造生成新公式的原先两个公式和新公式进行了对照分析。经数值实验,该方法具有可行性,某些所得新公式在求解常微分方程初值问题中具有较高稳定性和较小误差。
关键词 :初值问题,线性多步法,稳定性,加权平均法
解常微分方程在很多学科领域内都有着重要的应用,自动控制、各种电子学装置的设计、弹道的计算、飞机和导弹飞行的稳定性的研究、化学反应过程稳定性的研究等问题,都可以化为研究常微分解的性质的问题,或者化为求常微分方程的解。但大部分的常微分方程其真解一般难以通过解析的方法来获得,直到现在有许多类型的微分方程还不能给出解的解析表达式,通常只能用数值的方法进行计算。有关这一问题的研究早在十八世纪就已经开始了,现在计算机的发展更是为常微分方程的应用及理论研究提供了有力的工具,从而能使人们认识解的种种性质及其数值特征[
常微分方程的数值算法发展到今天已有了线性多步法、龙格–库塔法和在此基础上发展起来的单支方法、分块方法、循环方法、外推法、混合方法、二阶导数法以及各种常用的预估校正算法。其中比较经常用到的线性多步法公式有Euler公式、Heun公式、中点公式、Milne公式、Adams公式、Simpson公式、Hamming公式,Gear方法、Adams预估–校正法和Mile预估-Hamming校正法公式等 [
一般的,假设在第
这里
设
代入式(2.1.1),整理得
于是,欲使格式(2.1.1)成为
能符合到
如果格式(2.1.1)的系数满足这组条件,则将式(2.1.3)与式(2.1.2)相减,即得局部截断误差
其中
现研究基于泰勒展开的四步显示格式:
欲使这类格式成为四阶的,按格式(2.1.4),其系数应当满足条件
这里有七个待定系数
将数据代入式(2.1.5),解得截断误差为
根据文献[
其局部截断误差为
结合基于泰勒展开的四步显示格式(2.2.3)和Hamming格式(2.2.4)这两个公式,令
在式(2.2.5)中,取
该公式为一个隐式公式,对于其中的
预测
校正
进行求解。
下面通过求解一个例题,检验公式:
求解初值问题
下面我们依次使用基于泰勒展开的四步显示格式(2.2.3)、Milne-Hamming格式(2.2.7) (2.2.8)和新得到的公式(2.2.6)求解初值问题(2.3.1),画出它们在区间
其中,
其图像如下:
将精确值
2.3 | 2.366432 | 2.375243 | 2.373229 | 2.373033 |
2.31 | 2.370654 | 2.376292 | 2.377576 | 2.377376 |
2.32 | 2.374868 | 2.372462 | 2.381917 | 2.381714 |
2.33 | 2.379075 | 2.388074 | 2.386254 | 2.386047 |
2.34 | 2.383275 | 2.389802 | 2.390586 | 2.390375 |
2.35 | 2.387467 | 2.384688 | 2.394912 | 2.394698 |
2.36 | 2.391652 | 2.400788 | 2.399234 | 2.399016 |
2.37 | 2.39583 | 2.403332 | 2.403551 | 2.403329 |
2.38 | 2.4 | 2.396857 | 2.407863 | 2.407637 |
2.39 | 2.404163 | 2.413376 | 2.412171 | 2.41194 |
2.4 | 2.408319 | 2.416886 | 2.416474 | 2.41624 |
2.41 | 2.412468 | 2.408977 | 2.420773 | 2.420534 |
2.42 | 2.416609 | 2.425834 | 2.425068 | 2.424824 |
2.43 | 2.420744 | 2.430467 | 2.429358 | 2.42911 |
2.44 | 2.424871 | 2.421056 | 2.433645 | 2.433392 |
2.45 | 2.428992 | 2.438152 | 2.437927 | 2.437669 |
2.46 | 2.433105 | 2.44408 | 2.442205 | 2.441943 |
2.47 | 2.437212 | 2.433107 | 2.446479 | 2.446213 |
2.48 | 2.441311 | 2.450323 | 2.45075 | 2.450478 |
2.49 | 2.445404 | 2.457725 | 2.455017 | 2.45474 |
2.5 | 2.44949 | 2.44514 | 2.45928 | 2.458999 |
2.51 | 2.453569 | 2.462338 | 2.46354 | 2.463253 |
2.52 | 2.457641 | 2.471403 | 2.467797 | 2.467505 |
2.53 | 2.461707 | 2.457169 | 2.47205 | 2.471752 |
2.54 | 2.465766 | 2.474191 | 2.4763 | 2.475997 |
2.55 | 2.469818 | 2.485116 | 2.480547 | 2.480238 |
2.56 | 2.473863 | 2.469208 | 2.484791 | 2.484477 |
2.57 | 2.477902 | 2.485872 | 2.489032 | 2.488712 |
2.58 | 2.481935 | 2.49886 | 2.493271 | 2.492945 |
2.59 | 2.485961 | 2.481274 | 2.497507 | 2.497174 |
2.6 | 2.48998 | 2.497372 | 2.50174 | 2.501401 |
表1.
图1.
图2.
图3.
由上表可以发现在对例题方程求解过程中,使用新公式的解的误差小于使用基于泰勒展开的四步显示格式(2.2.3)和使用Hamming格式(2.2.4)的解的误差。因此,新公式比基于泰勒展开的四步显示格式(2.2.3)和Hamming格式(2.2.4)更为精确。
线性多步法是解决实际中所遇到的常微分方程初值问题的一种重要方法。本文在前人的基础上作了进一步的深入研究,对于线性多步法公式中基于泰勒展开的线性多步法构造法进行了探索,得出了一个新的公式。实验证明利用已有公式,进行加权平均所得的新公式确实可能出现某些具有较好性质,如局部截断误差较小、精确度较高的公式。
本工作受国家自然科学基金(项目号:11326096,11401247),广东省自然科学基金(项目号:2015A030313674),广东高校优秀青年创新人才培养计划项目(项目号:2013LYM_0089)和惠州学院博士启动项目(项目号:C511.0206)支持。
黄志远,胡志军,王承. 基于泰勒展开的线性多步法构造方法的研究A Study on Construction for Linear Multi-Step Methods Based on Taylor Expansion[J]. 应用数学进展, 2015, 04(04): 336-342. http://dx.doi.org/10.12677/AAM.2015.44042