高振荡常微分方程的高效数值方法

来源 :南京农业大学 | 被引量 : 0次 | 上传用户:q43372958
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
高振荡微分方程大量存在于天文学、量子力学、生物化学、分子动力学、生物学等领域.由于大部分振荡微分方程的解析解难以得到,所以对其数值解法的研究显得尤为重要.虽然对于微分方程数值解的研究已经有了丰富的成果,如Runge-Kutta方法,Nystr(o)m方法和线性多步法等,但是这些基于Taylor级数分析的传统方法并不能有效地求解高振荡问题.最近10多年来数值分析学家们已经对高振荡数值积分以及高振荡微分方程数值解做了大量的研究,形成了一系列的理论和方法.已经被证明为有效的方法主要有:渐近方法,Filon-型方法,Levin-型方法,最速下降法等.本文前人工作的基础上,发展了若干新方法和新技巧.   作为本文工作的起点,第一章概述了现有的理论结果以及数值方法.主要包括:求解高振荡积分的渐近方法(asymptotic method),Filon-型方法(Filon-type method)以及Levin-型方法(Levin-type method);求解线性高振荡微分方程的渐近方法(asymptoticmethod)和Filon-型方法;求解非线性高振荡微分方程的波形松弛Filon-型方法(wav-form relaxation Filon-type method);指数拟合方法的误差分析.   第二章建立了求解一阶振荡微分方程y=f(x,y)初值问题的三角配置法.基本思想是用三角多项式来逼近振荡微分方程的解.首先证明了三角配置多项式等价于一类Runge-Kutta方法,并导出系数表达式,其形式与传统的配置法类似.然后导出了三点三角配置法的系数,分析了相应的误差.实验结果表明新的三角配置法优于传统的方法,包括Matlab软件包中求解微分方程的方法ode113.   第三章对线性高振荡微分方程组提出了基于绝热变换的复节点Filon-型方法.首先通过变换将高振荡微分方程y(x)=Aωy(x)+f(x)转化成一个非高振荡的一阶微分方程组η(x)=exp(-i(Φ)(x))Q-1f(x)并在区间[xn,xn+1]上积分.对Q-1f(x)的不同分量,在区间端点附近选取不同的节点,在选取的节点上分别作各个分量的Hermite插值多项式,在积分式中以插值多项式代替原函数精确积分.二阶高振荡微分方程组y"+ω2y=f(x)的初值问题等价地转换成一阶初值问题求解.实验结果表明我们的方法优于传统的Filon-型方法.   第四章对二阶高振荡微分方程组考虑两步Levin-型方法.通过常数变易法给出方程的解在连续三个时间点处值的关系,将问题转化为一个高振荡积分.首先用Levin-型方法求解这个高振荡积分,继而导出了求解微分方程的数值格式,并给出了相应的误差分析.   最后简要地总结了对本文的主要工作,并提出了未来研究的课题.
其他文献
如果说,姚文元于1965年11月10日发表的《评新编历史剧〈海瑞罢官〉》,是“文化大革命”的导火索的话,那么,以聂元梓为首的7个人于1966年5月25日在北京大学贴出的后来被称之为
混沌是自然界中广泛存在的一种复杂的非线性运动形式。由于混沌信号具有初值和参数的极端敏感性、非周期性、连续宽带频谱、类噪声以及长期不可预测性等优良特性,因而特别适用
Ni-W-P composite coatings reinforced by CeO2 and SiO2 nano-particles on the surface of common carbon steels, were prepared by double pulse electrodeposition. Th
均衡是一个研究许多实际生活现象中某些系统的一个核心概念,这包含了从 经济、网络到力学等许多领域,均衡在现实中的应用研究也促进了不动点和最优 化理论的发展.抽象变分不等
随着信息的快速发展,图像正在以惊人的速度增长,如何在海量的图像中检索出需要的图像是一个关注的问题。近年来,图像的检索技术正在不断的发展,由最初的基于文本的图像检索技
偏微分方程在科学和工程技术中有着广泛的应用,许多实际问题的数学模型都可以用偏微分方程来描述,但很多偏微分方程无法求出解析解,只能用各种方法求出其数值解。格子Boltzma
期刊
Anodized composite films containing SiC nanoparticles were synthesized on Ti6Al4V alloy by anodic oxidation procedure in C4O6H4Na2 electrolyte. Scanning electro
不确定性问题作为人工智能最核心的研究任务,将不确定性问题的求解方法大致分为两类:一类是基于概率的方法,一类是基于非概率的方法。因果图推理是一种概率的方法,因果图以图
由于非自治动力系统以及自治动力系统与其对应的超空间动力系统一直以来是广大学者研究的热点课题,而对超空间非自治动力系统的研究就成为一个比较新的课题.本文主要研究Li-Yo