求解约束Hamilton系统的指数拟合方法

来源 :南京农业大学 | 被引量 : 1次 | 上传用户:honeykaka
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
随着科学的发展和工程技术的进步,单纯的微分方程不足以对某些复杂系统进行正确的描述,这就促使人们去研究由微分方程组和代数方程组混合而成的系统,即微分-代数方程组(简记为DAEs)。另一方面很多应用研究中迫切要求数值方法能够长期保持问题精确解内在的物理性质或几何性质,这就导致了保结构算法的研究。   一切真实的、耗散可以忽略不计的物理过程都可以表示成Hamilton系统。理论与实践证明,辛方法(symplectic methods)在保持Hamilton系统解的定性性质方面要比非辛方法优越得多,这种优势在长时间积分中尤其明显。约束Hamilton系统(constrained Hamiltonian systems)是一类重要的微分-代数方程组,其中微分方程的解的轨迹始终保持在由约束方程定义的流形上。对于这一类约束问题的数值求解,前人已经作了一些较为系统的工作,如Ryckaert,Ciccotti和Berendsen给出了SHAKE方法,Jay考虑辛分块Runge-Kutta方法(PRK)。另一方面,近年来出现了一种称为指数拟合的方法,这种方法能精确积分某些基本的振荡函数,应用于振荡微分方程时数值解具有相当好的相性质,计算效率高于传统方法。   本文主要研究具有振荡解的约束Hamilton系统的指数拟合方法。   第一章是预备知识,概述了微分方程的基本概念,包括常微分方程(ODEs)初值问题解的存在唯一性,数值方法的相容性、收敛性、稳定性,介绍了求解Hamilton系统的辛Runge—Kutta方法以及辛Partitioned Runge-Kutta方法,以及求解振动问题的保结构方法中最常用的技术之一——指数拟合方法。   第二章讨论求解微分-代数方程组的RK方法。在简要地介绍有关微分-代数方程组的基本概念之后,分别阐述了求解指数-1,指数-2和指数-3的微分-代数方程组的RK方法,以及在隐式方法实现过程中需要使用的简化Newton法.对指数-1微分-代数方程组及其数值解系统地建立了双色有根树(bi-colored rooted tree)理论与相关的B-级数理论,并导出了求解指数-1 DAEs的RK方法的阶条件。   第三章讨论求解约束Hamilton系统的辛方法。这章介绍Hamilton系统精确流的一些性质以及满足这些性质的数值方法,着重分析了LabattoⅢA-ⅢB相延迟阶。由于是辛的,这个方法是零耗散的。   第四章主要研究求解约束Hamilton系统的指数拟合方法.利用辛条件以及指数拟合条件,我们构造一个极大化相延迟阶的三级辛指数拟合PRK方法。数值结果验证,与传统的PRK方法相比,新的方法能更精确地保持Hamilton能量,计算效率也更高。   最后总结本文的主要内容贡献,指出约束振荡Hamilton系统的指数拟合方法方面有待进一步研究的问题。
其他文献
论文主要研究了基于多尺度的非采样Contourlet变换(NonSubsampled Contourlet Transform,简称NSCT)与小波变换(Wavelet Transform)相结合的图像去噪方法。小波变换具有各向同
近几十年,随着计算机科学技术的飞速发展,大维数据分析在现代科学研究中越来越突显其重要性,比如在生物学的微阵列数据,金融学的股票市场分析,无线通讯网络等新兴领域中,都出
学位
在大量的自然和社会现象中不可避免地存在时滞现象,亦即事物的发展趋势不仅依赖于当前的状态,而且还依赖于事物过去的情况。时滞系统的控制是控制理论应用的一个重要领域。时
关联规则挖掘是数据挖掘领域中一个重要的研究方向,揭示数据集中不同领域或属性间的有价值联系,具有重要的理论价值和广泛的应用前景。本文系统地讨论了关联规则挖掘的相关理
基于身份的密码体制可弥补传统的基于证书密码体制的不足,因此近年来发展迅速。基于身份密码体制的基本工具是双线性对,其有效实现取决于双线性对的计算效率。但就目前而言,双线
我国是世界上煤炭储量最丰富的国家之一,煤炭能源在国民经济发展中起着重要的基础作用。随着煤炭资源开采强度的不断加大,煤层开采效率越来越受到关注。水溶液对煤岩损伤弱化
近年来,随着消费者环保意识的增强,非政府环保组织的监督,政府不断出台措施鼓励厂商生产绿色环保的产品。再制造作为绿色制造的典型形式,已经受到越来越多制造商的关注。再制造产
通过人们对社交系统,信息系统等的整体研究,发现它们背后都可以用复杂网络来表示,复杂网络除了具有无标度特性和小世界特性外,社区结构也是其中一种重要的结构特性.社区结构中的