高振荡问题的高效数值方法研究

来源 :中南大学 | 被引量 : 6次 | 上传用户:truebyb
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
高振荡问题如高振荡数值积分、高振荡微分、积分方程数值解等广泛存在于量子化学、信号处理、医疗图像、分子动力学、流体力学等众多问题中,是公认的数学难题之一.本篇博士学位论文主要研究了高振荡问题的高效数值方法以及与高振荡相关的数值逼近问题.高振荡问题部分主要研究了高振荡函数数值积分如Fourier型和Bessel型积分,高振荡函数的Cauchy主值的计算以及具有高振荡Bessel核的Volterra第一型积分方程的数值解问题.逼近部分的研究主要包括Legendre系数的衰减估计和实现Legendre插值的快速稳定算法.全文由如下部分组成:   第一章简述了高振荡问题的应用背景和研究的意义.第二章介绍了近年来发展起来的几种高振荡积分的高效算法,如渐近方法、Filon型方法、Levin型方法、广义积分法则以及数值最速下降法等.   第三章讨论了高振荡积分的Levin迭代法.对于Fourier型积分,当被积函数在积分区间内没有驻点时,证明了Levin迭代法与渐近方法的等价性,同时将Levin迭代法推广到高振荡Bessel型积分.   第四章主要讨论了以下形式的高振荡有限Hilbert变换的计算,()提出了计算该型积分的三种新的算法.当函数f(x)在包含积分区间[-1,1]的一个充分大的复平面邻域内解析时,提出了一种可以利用经典的Gauss-Laguerre积分法则实现的数值最速下降法,此方法的精度随振荡频率ω的增加而快速提高.当f(x)在包含积分区间[-1,1]的一个小的邻域内解析时,提出了一种基于Chebyshev插值的快速实现算法,该算法可以由三项递推以及快速Fourier变换高效的实现.同时证明了所提出的Chebyshev插值型算法逼近振荡Hilbert变换是一致收敛的,也就是其误差界与τ无关.在前面两种方法的基础上,我们将Filon型方法推广到了当前积分的计算,Filon型方法不要求函数f(x)解析,它可以由递推关系快速的计算并且其精度随频率ω的增加而快速提高.   第五章主要研究了一个具有高振荡Bessel核的第一类Volterra卷积型积分方程的数值解问题.首先证明了对于积分区间含O的高振荡Bessel型积分,积分的渐近展开与Bessel函数的阶的关系,特别的,当Bessel函数的阶为整数时,提高了渐近方法以及Filon型方法的误差阶的估计.利用这个新的结果,研究了具有高振荡Bessel核的第一型Volterra积分方程的渐近展开和数值方法.首先利用Laplace变换技术将积分方程的解转化为一个高振荡的Bessel积分,然后利用提高的渐近展开给出了积分方程解的渐近展开.进一步还设计了一种计算积分方程解的Filon型方法,Filon型方法是一种对于固定频率收敛的方法,并且其精度可以用增加导数插值或增加插值节点而快速提高.   第六章主要研究了Legendre逼近问题.首先建立了Legendre系数的衰减估计,得到了Legendre的截断级数展开的误差估计.其次研究了Legendre插值公式的快速稳定算法.众所周知,Chebyshev插值公式可以通过其重心形式使用O(n)次运算量快速并且稳定的实现,然而对于Legendre插值公式的研究却没有Chebyshev插值那样受到重视.在本章中,我们使用Gauss-Legendre积分法则的节点和权显式的构造了Legendre重心权的显式形式,因此先使用Glaser-Liu-Rokhlin算法计算Gauss-Legendre积分法则的节点和权,这需要O(n)次运算量,然后再计算Legendre的重心权,从而Legendre插值公式可以用O(n)次运算量快速并且稳定的实现,这就使得Legendre插值公式和Chebyshev插值公式具有同阶的运算量.类似的,我们还显式的给出了重心Gauss-Legendre-Lobatto的权,从而Gauss-Legendre-Lobaao插值公式也可快速稳定的计算.
其他文献
学位
网络技术已广泛渗透于人类生活的各个领域,信息交流日渐频繁,人类社会步入了数字化时代。数据信息的交流使得它成为人类社会在信息利用方面重要的工具,所以通过网络安全传输
近年来,对随机种群扩散系统的研究引起了许多专家、学者的广泛关注,通常情况下,大多数随机种群扩散系统没有解析解,带poisson跳、Markov调制的随机种群扩散系统也是如此,因此,对随
分析以往针对教育系统建立的模型发现,以前所建立起来的模型,往往只是针对教学质量,在对教育创新扩散过程进行的研究方面还很欠缺。而微分方程模型在种群生态学和工程技术领
本文研究以下四阶强阻尼非线性波动方程的初边值问题utt-△u+△2u+μut-α△ut-△utt=f(u),μ,α>0,χ∈Ω,t>0,u|t=0=u0,ut|t=0=u1,χ∈Ω,u|aΩ=0,其中Ω是Rn中的有界开集,且
股票市场是一个高度复杂的动力学系统,基于股票间的相关性研究对掌握其内在结构和演化规律具有重要的价值。随机矩阵理论(Random matrix theory,RMT)在核物理、混沌系统以及无线通信等领域具有非常广泛的应用。通过将随机矩阵理论引入到股票市场,发现相关矩阵中确实蕴含着大量的噪声信息,同时学者们还提出可利用RMT对股票市场进行过滤,其去噪方法主要有LCPB法(Laloux L,Cizeau
经验似然作为一种非参数方法有很多优点,尤其是在用于估计和检验时,不需要假设总体分布.本文对于经验似然置信域体积,提出了重要抽样法计算置信域体积的方法;对于若干回归模型,利用
本文假设了渔业资源生长满足多代共存的Logistic模型,利用多个体博弈达到Nash均衡时的捕捞量得到总收获函数,建立了多个体博弈公共渔业资源的二维动力学模型。提出了把多个体博
我们主要研究的是微分动力系统中的一致双曲以外的一些动力系统的性质。   我们试图去理解那些在一致双曲系统中熟知的、非常理想的结论是否在在更大的一类系统中成立。我
由于事物的发展规律不仅依赖于当前时刻的状态,还依赖过去某时刻或某段时间内的状态,所以时滞微分方程对客观现象的描述比常微分方程更加准确和合理。现如今时滞微分方程被广泛