分数阶扩散方程的预处理技术研究

来源 :电子科技大学 | 被引量 : 0次 | 上传用户:download1006
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
作为微分方程中很重要的一部分,分数阶扩散方程在力学,物理学上备受关注。除此之外,在化学,生物学,图像处理等众多领域也应用广泛。因为相比于整数阶微分模型,分数阶往往能够更加精确的描述其中的现象和结果,所以求解这类方程也变得十分重要。针对这类方程,一方面,由于可以作为封闭形式的解析解通常不可用,所以数值方法成了求解其近似解的主要方法。另一方面,由于分数阶算子具有非局部的性质,采用简单的离散化,即使是隐含的,也会导致无条件不稳定,并且大多数FDE的数值方法倾向于生成全系数矩阵,这导致采用直接计算法进行求解时,需要O(n~2)的存储和O(n~3)的计算复杂度,这里的n代表空间网格的数量。为了节约成本,加快计算速度,寻找一种高效快速的数值解法已经成为众多学者关注的重心。其中通过构造系数矩阵的预处理子,加速求解分数阶微分方程,成为求解这类方程的有力工具。并且通过众多研究者的努力,已经提出了很多经典且有效的预处理子。本文也是从寻找一种高效,快速的数值解法出发,通过引入线性方程组的矩阵分裂带代法,构造系数矩阵的预处理子,从而加速Krylov子空间迭代法求解(1)从空间分数阶扩散方程中离散出来的对角加Toeplitz线性系统;(2)从时间-空间分数阶扩散方程出发,在时空框架中通过有限差分进行离散,得到的完全耦合离散线性Kronecker乘积和的线性系统。针对(1)产生的线性系统,本文提出了双参数的双步分裂迭代法GDTS方法来处理;针对(2)产生的线性系统,提出了单参数的双步分裂迭代法TTS方法处理。根据这两种方法,分别求得其迭代矩阵,并通过迭代矩阵谱半径小于1,迭代法收敛这个条件,建立收敛理论,推导渐近收敛上界,求解迭代参数的最优值。最后引入Strang’s循环矩阵来代替迭代矩阵中的Toeplitz部分,构造出系数矩阵的预处理子。由于预处理子中将Toeplitz矩阵替换为循环矩阵,故矩阵向量的乘法计算量也降低到O(n log n)。通过理论推导,系数矩阵在经过文中所构造的预处理子处理后,可以表示为一个单位矩阵、一个低秩矩阵和一个小范数矩阵之和。数值实验也表明文中构造的预处理子可以有效集中原始系数矩阵的谱半径,并在采用Krylov子空间迭代法求解时,迭代步数明显减少,加快求解过程。
其他文献
邱华栋城市小说中的"新人物",写得最有深度的是"城市单面人",即被城市高度发达的物质文明所平面化了的人。这些艺术形象体现了邱华栋对现代文明与人的关系的深入思考,即物质挤压下城市人的异化。邱华栋着力塑造"城市单面人"的艺术典型,显然受到了马尔库塞关于"异化""单向度"等理论的深刻影响,触及城与人的哲学关联,因而在整个新时期小说发展中是一个重大的突破。
计算电磁学作为一门汇聚了电磁学、计算机科学和数值计算方法等学科知识的交叉学科,为研究当今时代无处不在的电磁现象提供了理论和技术支持。以Maxwell方程组为核心的电磁场正问题分析和反问题求解是计算电磁学研究中的重要部分,正问题可以归结为在不同的电磁环境参数下求解Maxwell方程组,反问题则可以看作在已知电磁响应的基础上反向寻找对应的电磁环境参数,针对它们的高效率、高精度的求解算法一直以来都是科学
主成分分析在计算机视觉和生物信息学等应用中起着重要作用。但传统的主成分分析容易受到严重损坏或异常的观测值的干扰,而这些现象在现实世界中无处不在,所以鲁棒性主成分分析的提出保证了一定程度的抗干扰能力。在现实世界中,许多数据,如彩色图像、视频和4-D fMRI数据是多维,简单的进行矩阵化会导致数据内部结构的破坏。张量分解是强大计算工具,它天然的可以刻画高维数据。自然地,通过建模分离数据中的低秩和稀疏成
图片理解的问题一直是计算机视觉领域的关键性问题,最近随着深度学习引领的计算机视觉的快速发展,相关的图片语义理解任务变得愈发关键。而边缘检测作为图片语义理解的基础任务,也得到了越来越多的关注。在传统的计算机视觉理论中,通过边缘的识别和检测,能够得到后续的图片分割等结果;而在深度学习时代,边缘检测和识别能得到后续的图片分割结果,并且在相关的图片分割、物体检测等领域,也能达到优化相关高层次任务效果的作用
地球物理反演作为获取地下岩石物理属性参数的重要手段,它的反演结果为后期油气勘探、储量预测、油气开发等奠定了坚实的基础。地质统计学反演作为一种重要的地球物理反演方法,在获得比常规方法更高分辨率反演结果的同时,还能获得与地震数据比较接近的横向连续性,因此得到了广泛的研究和应用。现有地质统计学反演方法主要是通过构建变差函数从已知的测井数据获取先验信息,然后再用于指导后续的反演过程。但是,该方法存在两个主
目前,随着科学技术的快速发展以及资源需求的日益增高,太空资源成为各国竞相开发的重要目标之一。航空航天飞行器作为开发太空资源的工具就显得尤为重要,这也意味这高空飞行器必然要面临通信中断的问题,也即是黑障问题,这是迄今为止航空航天飞行器遇到的最具有挑战性的难题之一。黑障问题主要体现为电磁波在色散介质等离子体中的传输问题,所以研究二者之间的作用关系对缓解黑障问题是非常必要的。因此本课题的研究对缓解黑障问
对称性在人们对世界的认识当中始终起着不可或缺的作用,而研究这一性质的数学工具群论也是近代数学中的一个重要分支。在上世纪80年代中期,量子群在应用反散射方法研究量子场论和统计物理中的精确可解模型中被人们提出,在此之后,各式各样的变形代数变成了物理学家们研究的热点之一。这些变形代数的不可分表示,不可约表示和振子实现等等问题有着重要的意义。本文讨论了多参数非线性变形so(3)代数Rq,r(v)以及多参数
近年来,随着计算机和电子信息技术的飞速发展,无人机领域也得到了飞跃式的发展。在无人机领域中,最基础也是最重要的研究方向就是无人机的导航定位。导航定位算法是无人机实现运动控制、自主导航等上层应用的基础。在实际工作和生产中,人们往往需要四旋翼无人机在山洞,隧道,野外等较为复杂的环境中完成作业,四旋翼无人机在此复杂的场景中工作时依靠单一的传感器无法完成可靠的定位。本文针对传统的导航定位方法鲁棒性、准确度
跌倒是危害老年人安全的主要原因之一,跌倒中头部撞击地面尤其危险,严重可导致死亡。近年来,意外跌倒对老年人肢体损伤的研究受到高度重视,对跌倒后严重性进行评估,可避免老年人因疼痛敏感性低而导致的延误就诊,在医生急救时给予有用信息,为后续跌倒保护研究提供相关参考。论文开展人体跌倒头部伤害评估模型研究,从跌倒过程不同人体姿态入手,面向视频和可穿戴传感的不同场景开展人体跌倒的头部伤害严重性评估,引入人体关节
属性基加密(Attribute-Based Encryption,ABE)作为一种将加密和访问控制的思想结合起来的非对称加密体制,在提供了一对多的加密方式的同时也提供了细粒度的访问控制。这些特性使得属性基加密在云计算、基于云的物联网以及无线体域网等领域得到了广泛的研究和应用。不幸的是,目前大多数的策略隐藏的属性基加密(Policy Hiding Attribute-Based Encryption