线性高振荡常微分方程数值解法的若干研究

来源 :北京交通大学 | 被引量 : 0次 | 上传用户:miao4701730
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文以特殊的线性振荡方程y”+g(t)y=0(其中lim g(t)=+∞)为例讨论了高振荡常微分方程数值解问题。 高振荡微分方程是指其解含有高振荡函数的一类微分方程,它在分子动力学、天体力学、量子化学以及原子物理等方面有着广泛的应用。对于高振荡微分方程给出一种好的数值解法是一件非常困难的事情。例如,对于形如y”+g(t)y=0的线性高振荡方程,用经典的方法,如Runge-Kutta法、线形多步法等方法在处理该类问题时均会产生较大的误差。近来,Iserles利用Magnus展开方法详细研究了该类方程数值解法问题,给出了计算结果较好的数值算法。 在这篇论文中,我们首先介绍一些基本概念和基本知识,为后面的内容做准备工作。然后用梯形方法对线性振荡方程y”+g(t)y=0数值求解,理论分析及数值实验均显示,用梯形方法求解会产生较大的误差。我们就梯形格式作了几种修改,误差分析及数值结果均显示,修改后格式的数值解都要优于梯形格式的数值解。 我们还系统地介绍了Magnus展开方法及修正的Magnus展开方法,从修正的Magnus展开方法出发,我们考虑利用Cayley变换构造线性高振荡微分方程的数值解法。这样构造的解法涉及到高振荡函数的积分,我们采用Filon方法计算,数值结果显示,该数值解法具有好的长时间数值跟踪能力。另外,我们也利用数值实验比较了用Gauss方法及Filon方法计算高振荡函数积分所给出的数值解法,实验显示,Filon方法相应的数值解法优于Gauss方法所相应的数值解法。
其他文献
学位
基于身份的密码体制一直是公钥密码学中的研究热点。现存的大多数基于身份的加密与签名方案都是采用双线性对的,双线性对不仅运算非常耗时,大大降低了方案的效率,而且在如今量子
在全球经济一体化发展的市场环境中,企业的持续发展离不开高素质的人才队伍和科学高效的人力资源管理.而国有企业更应起带头的模范作用,革新管理理念,“科学发展观”为指导,
现代组合投资理论是关于在收益不确定条件下投资行为的理论,它是由美国经济学家马柯维茨在1952年首先提出的。此后人们做了不懈的研究,将该理论进行推广、改进、发展和完善。
1984年,Karmarkar提出求解线性规划的新方法称为内点法,并证明该方法不但具有多项式复杂性,而且实际计算对大规模线性规划问题的效果优于单纯形方法.该算法发表后,掀起研究内点法
本学位论文主要围绕这个猜想展开研究.  第一章,给出了一些基本概念以及图的(d,1)-全标号问题的研究背景和现状,并且介绍了本学位论文的主要结果. 第二章,给出了最大度为3的
我国地震界在过去三十多年的预报实践和几次地震预报攻关中,探索了许多新的地震学方法及预报指标。目前在地震中期和中短期预报中,常使用如反映介质所受应力的状态和介质的均匀
2013年11月,中共召开十八届三中全会,提出了鼓励发展混合所有制经济,推动国有企业深化改革的会议精神。在我国,混合所有制经济的发展和国有企业改革都经历了长期的历史发展过