时间分数阶微分方程数值求解的三次B样条方法

来源 :哈尔滨工业大学 | 被引量 : 0次 | 上传用户:calltt_stephy
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
时间分数阶微分方程是微分方程研究领域的一个重要分支,在各类扩散系统中有着广泛的应用。但是由于时间分数阶微分方程具有分数阶导数,不易求解,从现有的研究方法看,求解此类微分方程的主要方法是将其离散化进而求解数值解。由于B样条函数是对称单峰值函数,具有光滑性好、紧支集等特点,同时保留了分段插值多项式计算简单、稳定性好、收敛性有保证且易于在计算机上实现的优点。并且B样条配点法构造简单,数值精度高,易于处理复杂的边界问题,目前已成为求解偏微分方程的重要数值方法之一。所以本文提出使用三次B样条方法求解时间分数阶微分方程数值解以及参数反演,研究其求解此类问题的收敛性和稳定性。为了研究三次B样条方法能否较好地求解时间分数阶微分方程数值解,本文将利用Caputo分数阶微分定义,将时间分数阶微分转变为积分形式,同时对于时间维度的一阶偏导采用一阶差分格式,离散化时间分数阶微分方程。实验结果表明,三次B样条方法能够有效快速地求解时间分数阶微分方程,并且能够将误差控制在10-5数量级。随后,本文根据三次B样条插值方法的误差,通过理论证明来验证三次B样条方法求解时间分数阶微分方程的有效性。通过证明我们得到三次B样条方法求解时间分数阶微分方程的收敛阶为(Δt2-α+Δx2),进一步证明三次B样条方法求解时间分数阶微分方程的有效性。源项反演也是微分方程数值解的一个重要部分,在实际问题中通常需要通过观测值来推测源项、参数或者初边值,由于时间分数阶微分方程的特殊性,当前的分数阶微分与之前所有时间的函数值相关,所以很难使用整数阶微分方程的方法求解,因此本文提出了一种新的泛函,对方程进行平方计算并在时间空间上积分产生一个泛函,极小化该泛函所得到的源项即所求源项。但是仅仅对泛函偏微分等于零求源项具有不适定性,因此本文在这里引入正则化方法,然后通过最小二乘法的思想求解参数方程。实验结果表明,改进的三次B样条源项反演算法虽然因为三次B样条函数的分段特点增加了源项中的参数个数,但是算法不需要通过迭代,避免了复杂的计算,拟合效果整体上要优于多项式拟合算法。
其他文献
叠论是概型理论的推广。叠论与数论、几何Langlands等领域联系紧密,并为这些领域中的一些难题的解决提供了新方法。Higgs丛是极为重要的几何对象,近年来,很多重大问题的解决
在科技高速发展的21世纪,创新的意识和能力成为国家发展的最重要决定因素。瑞士虽作为一个资源匮乏,国土面积小的国家,其却以高效的科技创新体系闻名世界,其科学研究水平、自
市场竞争和全球化以及随着工业4.0、中国智造2025的不断推出,制造业面临新的挑战,快速并准时交货、需求品种的增加和特殊定制、不断改善质量及降低销售价格、产品的生命周期
2013年5月,场内质押业务正式开展,股权质押规模快速增长,2014年-2015年的牛市行情更是助推了股权质押业务的扩张。股权质押不仅具有便捷的特点,而且能够在不稀释大股东控制权
在市场化环境中,电价信号是市场参与者关注的焦点,它不仅是交易决策关键参考因素,而且是市场管理者的重要监管指标。因此,电价预测也是学术研究的热点。伴随着北京、广州两个
自2012年第一季度开始,中国乃至全球PC市场增长率持续下滑,甚至多个季度出现负增长。消费PC市场因受到手机等电子设备的影响,关注度持续走低。并且消费者也从对品牌的关注转
工程勘察作为工程建设领域的重要环节,自改革开放以来,伴随着我国不断深化的行业改革、不断完善的法规制度,以及市场经济的进一步发展,整个行业中的企业和从业人员数量和经营
地理学科的学习,实践性作业对于掌握学习内容、参与实践生活起到至关重要的作用,在完成实践作业的过程中,学生可以对理论知识进行实践检验,理解得更为深刻,形成长久记忆。在
无线通信在得到快速发展和广泛应用的同时,其自身所暴露出的安全问题日益趋于严峻。信息传输的媒介——无线信道,其开放式的传输环境,使得信息在传输过程中难免会受到潜在非法用户的窃听威胁,影响信息安全。以往保障信息安全的加密机制,会随着可用计算资源的改变而面临被破解的威胁。物理层安全技术的提出,则不同于依靠计算堆积的加密机制,而是创新性地利用无线信道的固有特性在物理层结合传输信号完成对数据的加密,保证信息
重卡驱动桥壳是重型卡车底盘的关键结构部件,在整车性能的提升中具有不可替代的作用。随着桥壳的加工质量要求不断提高,其对桥壳加工误差的测量精度要求也越来越高。桥壳测量