求解振荡常微分方程的辛指数拟合Runge-Kutta(-Nystrom)型方法

来源 :南京农业大学 | 被引量 : 0次 | 上传用户:hongnanjing
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
物理学中的电路问题、弹性力学、天体力学、量子物理等领域中的许多问题都可以归结为解具有振荡性质的一阶或二阶的常微分方程(组)(ODEs)。本论文主要研究数值求解这类振荡常微分方程初值问题(ⅣPs)Runge-Kutta(RK)型及Runge-Kutta-Nystrom(RKN型)保结构方法。   RK方法从提出至今已有百年的历史,其主要用来求解一阶常微分方程(组)的初值问题(ⅣP)。二阶常微分方程(组)当然可以通过增加“速度”分量方程而化为一阶方程组,从而用RK方法求解。后来Nystrom提出了直接求解二阶方程组的RKN方法。另一方面,20世纪60年代J.Butcher建立了根数(rootedtree)与B级数的理论,使得推导高价方法成为现实。最近的科学研究与工程应用进一步要求数值方法能够保持系统原有的某些重要性质(如系统的某些不变量,特别是Hamilton系统的辛性、Hamilton函数)。我们的工作重点就是针对振荡微分方程的结构特点,设计有效的数值方法,保持精确解的定性行为。   本论文分为四章。   第一章概述了微分方程的基本概念,包括常微分方程(ODE)初值问题解的存在唯一性,讨论数值方法的相容性、收敛性、稳定性,介绍了求解Hamilton系统的辛Runge-Kutta方法,以及求解振动问题的保结构方法中最常用的技术之一-指数拟合。   第二章提出了一类新的RK型方法(我们将其记为RKNd),新方法的内级有比传统的RK方法更高一阶的精度。在适当的简化假设条件下利用阶条件导出了级数为二的三种显式和四种隐式的新方法,并对各个新方法的稳定区域,色散阶和耗散阶进行了分析。应用这些方法求解一些典型的一阶、二阶问题的数值结果显示:新方法更适合于求解二阶微分方程。   第三章研究了与第二章所给出的新方法相应的指数拟合型方法。我们推导出了方法的阶条件,进而给出了级数为二的三种指数拟合型方法,并对这三种方法的稳定区域,色散阶和耗散阶进行了分析。应用这些方法求解一些典型的二阶振荡问题所得到的数值结果表明,指数拟合的这些新方法的计算效率明显高于指数拟合的RK方法。   第四章主要研究了求解Hamilton问题的对称辛指数拟合RKN方法。我们给出了修正RKN方法的对称性条件、辛性条件、指数拟合条件及代数阶条件,并构造了二级二阶、三级二阶,四级四阶的新方法,对每个新方法的周期稳定区域,色散阶和耗散阶进行了分析。数值实验结果表明,本章的对称辛指数拟合RKN方法与非对称方法在计算效率和保持Hamilton函数值方面相比具有明显的优势。   概括起来,本论文主要贡献有:   ●给出了一类新的RK型方法,这类新方法的内级的精确性至少比传统RK方法高一阶;   ●基于这些新的RK型方法,构造了相应的指数拟合方法;   ●针对Hamilton问题,研究了一类对称辛指数拟合RKN方法。
其他文献
科学技术中大量的数学模型可以用偏微分方程来描述,但是大部分偏微分方程的解析表达式求解是极其困难,甚至不可能的,所以需要用近似方法来求解。差分方法是求解偏微分方程定解问
随着中国航空运输需求的急剧增加,运输需求与交通容量之间的矛盾日益突出。由于受恶劣天气,飞机故障,流量控制,地面保障不力等诸多不可控因素的影响,航班不正常情况时有发生,给航空
最优控制是现代控制理论的重要组成部分,它一直受到控制理论界的重视而得到不断深入的研究和发展。近几年来,有关Burgers方程、Kdv方程、Kdvb方程以及K-S方程的最优控制方面的
本文主要讨论了小波乘子的充分条件以及去噪过程中分解层数自适应确定的问题.  首先,给出了二进二元小波Fourier矩阵乘子的充分条件.通过对高维空间小波以及紧框架的研究,证
本文基于李(Lie)群理论研究了二阶、三阶和n阶自治系统的可积性。主要内容和结果包括以下几个方面:   (1)首先给出了二阶恰当自治系统所接受的单参数李群的生成元的一般形
学位
疲劳失效是机械和结构零件常见的失效形式,约占机械事故的50%以上。因此疲劳破坏问题是机械可靠性研究的一个重点领域。而疲劳失效的主要原因是存在漏检的宏观缺陷(如裂纹)。
幼儿园是一个民族未来的希望,也是社会发展中极为重要的一份子.钢琴课是学前教育专业学生必须掌握的一项艺术技能,在幼儿园教学工作中扮演着重要角色.但是当前我国学前教育专
超弹性杆方程是在研究弹性可压缩物质时得到的,广泛应用于弹性力学等领域,具有重要的研究价值。本文主要研究广义超弹性杆方程Cauchy问题解的爆破,以及广义超弹性杆方程初边值问
随着信息技术在教育领域中的广泛应用,为了进一步提高学前教育专业人才培养工作的实际效果,加强对学前教育专业教学资源库的构建成为必然发展趋势.本文从学前教育专业教学资
本文主要研究了带色散项Degasperis-Procesi方程的尖峰解的轨道稳定性,以及粘性色散Degasperis-Procesi方程的Cauchy问题的局部适定性理论。Degasperis-Procesi方程,(简称D-P方