关于两类线性鞍点问题的预处理迭代法和两类线性矩阵方程的低秩并行化方法的研究

来源 :兰州大学 | 被引量 : 0次 | 上传用户:monishibi
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文主要针对两类具有鞍点结构的3×3块线性系统以及两类线性矩阵方程的数值求解展开研究.由于其广泛的应用背景,构造快速高效的数值求解方法解决这些问题具有重要的理论价值和实际意义.第二章中对于一类大型稀疏的3×3块鞍点线性系统,通过分析松弛维数分解(Relaxed dimensional factorization,简记为RDF)预处理子的性质,我们发展了一类参数化松弛维数分解(Parameterized RDF,简记为PRDF)预处理子去进一步改善RDF预处理子求解这类3×3块鞍点线性系统的预处理效果.类似于RDF预处理子,PRDF预处理子的预处理效果也密切地依赖于所含参数的取值.因此,为了使得PRDF预处理子的计算效率达到最佳,我们提出了一类快速而有效的方法用于计算PRDF预处理子中所含参数的最优值.最后,结合这些最优参数值并通过使用两个不同的数值算例,即Navier-Stokes方程和偏微分方程的约束优化问题,我们验证了PRDF预处理子求解大型稀疏3×3块鞍点线性系统的有效性和可行性.第三章中对于一类由液晶定向特性的连续介质模型离散而得的双鞍点线性系统,通过分析块三角(Block triangular,简记为BT)预处理子的性质,我们构造了一个两参数块三角(Two-parameter BT,简记为TPBT)预处理子,并证明了TPBT预处理矩阵的谱都是实的且完全落在了区间(0,2)上,与矩阵-1-1的谱半径取值无关.同时,还估计了TPBT预处理矩阵极小多项式次数的上界及其相应Krylov子空间的维数.此外,由于TPBT预处理子的预处理效果依赖于两个待定参数的取值,所以我们进一步推导出了一类快速而又实用的公式去估计这两个参数的拟最优值.最后,通过数值算例验证了含有拟最优参数值的TPBT预处理子对加速Krylov子空间方法的收敛速率是非常有效的.第四章中对于两类具有带状Topelitz结构的广义线性矩阵方程,通过利用Topelitz矩阵与循环矩阵之间的联系,我们提出了一种低秩并行化方法去求解这两类具有特殊结构的线性矩阵方程.同时,提出了一种有理扩展的Krylov子空间方法(Rational variant of extended Krylov subspace method,简记为REKSM)去求解相应的低秩校正方程,并且提供了一种非常实用的方法去计算REKSM中所含位移参数的最优值.此外,对于一类出现在2D线性波方程中的广义Sylvester方程来说,通过将单边扩展的Krylov子空间方法(EKSM)与Sherman-Morrison-Woodbury(SMW)公式结合起来,从而提出了一种新的方法(记为EKSM-SMW)求解相应的低秩校正方程及其压缩方程.最后,通过数值算例验证了所提出的低秩并行化方法用于求解这两类线性方程时的有效性.第五章我们考虑的是一类出现在时间依赖偏微分方程的时空离散化中的线性矩阵方程,例如,在并行时间积分的背景下,一类Para Diag算法已经被提出用于求解这类矩阵方程.我们发展了两类新的方法求解此类方程.第一种方法主要依赖于这样的一个观察,即Para Diag算法为了并行求解这些方程而对其系数矩阵所做的校正具有较低的秩.基于之前对矩阵方程低秩更新所做的工作,这使得我们可以利用张量化Krylov子空间方法来求解相应的低秩更新方程.第二种方法主要利用几个校正方程的解的线性组合来逼近初始线性矩阵方程的解.这两种方法都避免使用Para Diag及其相关的时空方法为了获得较高的计算精度而产生的迭代加细.反过来,我们的新方法有可能显著地优于现有的方法.最后,通过几个不同的数值算例验证了所提出的方法总是非常有效的.
其他文献
分汊河槽是自然河流地貌中的主要地貌单元,河口分汊河槽是水、沉积物和营养盐从上游流向下游直至外海的关键节点,在河流地貌塑造理论和河流治理工程应用中发挥着核心作用。河流通过河口分汊系统控制水沙再分配,分汊节点处的不均衡水沙输运导致节点上游大部分或全部流量流向下游中的一个汊道,致使河流改道甚至发生河道的废弃。因此,河口分汊河槽极不稳定,易发洪水灾害,并威胁河槽沿程水资源供应和航行安全。而且,径流下泄与潮
学位
在自然界中,信息处理能力对于生物体的生存与进化至关重要,如大脑中神经网络的感知、模式识别及决策生成能力。在大脑进化之前,单细胞有机体基于复杂的生物分子电路也能够处理信息以产生生存所需的“智能”行为。基于分子“思考”方式研究已经产生了一系列计算模型及人工化学系统相关应用,但开发功能性分子电路以提升其信息处理能力依然是纳米技术领域内的重要科研目标之一。核酸分子具有可用于工程化人工分子电路的诸多理想特性
学位
时间周期反应扩散传染病模型是理解具有季节性特征传染病空间分布演化的有力工具,其空-时动力学性质也是发展方程中的重要前沿问题之一.然而,除了系统的非单调性耦合结构和非自治性共同导致的基础困难外,其它现实因素的引入也对问题的分析带来了更多新的挑战.本文研究时间周期环境下S-I型反应扩散传染病模型的地理传播性质,以及有界域上的阈值动力学.具体研究如下:第二章研究了一类考虑由潜伏期内宿主随机移动引起的非局
学位
疾病的发生和发展与细胞命运密切相关,越来越多的治疗手段开始聚焦于通过改变细胞命运实现治疗目的。然而,细胞命运受细胞内多种因素影响,导致现有治疗手段难以通过调节体内单一参数实现细胞命运的定向改变,无法做到将不同疾病的治疗方法系统化。细胞命运改变的本质,是细胞内信号传输系统将外界刺激转化为内在的输入性信号,用于改变细胞的生命活动并最终决定细胞的生死存亡。钙离子(Ca2+)是细胞内第二信使之一,其细胞内
学位
本学位论文主要研究以下两类带磁场的二阶椭圆方程:磁Schr(?)dinger方程和临界磁Choquard方程,并利用变分方法得到了解的存在性,多解性和集中现象在第一章,介绍带磁场的椭圆方程的物理背景及国内外研究现状,给出本文所需的预备知识以及主要结果介绍.在第二章,研究在Neumann边界条件下磁Schr(?)dinger方程#12解的存在性,其中Ω是RN中有界的C1区域,v是(?)Ω上的单位外法
学位
本文从带算子代数角度研究了左余单位双代数,广义罗巴代数,DN-双代数和微分李代数的泛包络微分代数,全文共分为六章.第一章介绍了本文相关的研究背景,研究动机.同时为了本文内容的完整性,本章还回顾了所用到的基本概念和事实.第二章首先介绍了左余单位双代数的基本概念与性质.通过扩展的1-闭链条件,在双装饰平面根森林上定义了余乘,从而构造了双装饰平面根森林上的左余单位双代数结构.然后结合带算子代数和左余单位
学位
电子转移现象广泛存在于物理、化学、生命医学和材料科学等领域;其中,(类)电子转移反应在生命活动中起到十分重要的作用。一般而言,高度有序的(类)电子转移反应是机体稳定实现特定生物学功能的基础,一旦特定电子转移反应被破坏,极易影响机体生物学功能,引发一系列疾病。由此可知,精准调控关键电子转移反应,继而干扰或重塑机体特定生物学功能,有望为研发疾病治疗新策略提供借鉴。值得一提的是,纳米材料具有独特的物理化
学位
可比较图是表示偏序关系的无向图.余可比较图是可比较图的补图.这两类图与偏序集之间存在自然的一一对应关系.字典序广度优先搜索(LBFS)最早由Rose等人提出,用于解决弦图的识别问题.LBFS+是LBFS使用了“+”这一特定的破局机制后的一个变体.设σ0是图G的任意一个LBFS顶点排序,{σi}i≥1是G的一个顶点排序序列,满足σi=LBFS+(G,σi-1).字典序广度优先搜索循环数则定义为通过扫
学位
塑料垃圾与微塑料污染遍及各类环境介质当中,并且在高原、深海等人迹罕至之处均有报道。由于塑料自身的环境持久性以及潜在的生态风险,已经成为全球性的环境污染问题,引发科学家、政府、企业、公民和环保组织等相关各方的共同关注。舟山渔场作为我国重要的渔场,受人类活动影响显著,但是目前对该塑料垃圾与微塑料污染热点区域的深入研究十分缺乏,制约了研究者对该区域塑料垃圾与微塑料污染的来源、归趋以及评估生态风险的理解。
学位
帕金森病(Parkinson’s disease,PD)是一种运动障碍性疾病,常见于中老年人群,其典型临床表现为静止性震颤、僵直和运动迟缓。脑深部电刺激(deep brain stimulation,DBS)手术是治疗PD的重要手段。但鉴于PD个体化诊疗具有主观性与异质性,患者手术靶点的精准定位和手术受益率的评估仍然充满挑战。定量磁化率成像(quantitative susceptibility
学位