三类结构化离散系统的高效迭代法与预处理技术

来源 :兰州大学 | 被引量 : 0次 | 上传用户:sukoo777
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
在大规模科学计算中,许多问题如PDE约束优化问题,Navier-Stokes方程,最小二乘解问题以及分数阶微分方程等经离散得到一些具有特殊结构的线性系统或者矩阵方程.本文针对三类具有特殊结构的离散系统:包括鞍点结构的线性系统,Toeplitz类结构的线性系统以及Sylvester结构的矩阵方程,提出了一系列新的迭代方法与预处理技术.本文第一章详细介绍了三类结构化离散系统的背景和研究意义,研究现状,以及本文的主要工作和创新点.本文第二章的前半部分对PDE约束的优化问题的离散系统,根据其特殊结构讨论了有效解法及预处理技术:对以热方程为约束条件的优化问题离散得到的特殊结构的线性系统,采用了自然降阶的方法,并对降阶后的线性系统提出了一类新的加性块对角预处理技术;对以非稳态Burgers方程为约束条件的优化问题离散得到的鞍点线性系统,提出了一种非标准内积意义下对应Schur补近似的新的预处理技术,得到一种有效的预处理算法;对以Stokes方程为约束条件的速度追踪优化问题的离散线性系统提出了一种松弛分裂迭代法,建立了双参数的松弛分裂预条件子.详细分析了上述三类预条件子对应的预处理矩阵的谱性质,在本章的后半部分,我们针对2 × 2块结构的线性系统,提出了三种预处理技术:即广义的加性块对角预条件子,广义的位移分裂预条件子和一类加参旋转块预条件子,分析了相应迭代法的收敛性条件和预处理矩阵的谱性质.数值例子验证了本章的新算法及新预条件子的有效性.本文第三章针对空间分数阶对流扩散方程离散得到的具有Toeplitz类结构的线性系统提出了基于不完全的循环与反循环矩阵分裂的三步交错迭代方法.讨论了新迭代法的收敛性条件,数值实验验证了新算法的有效性.该方法的主要优点在于,每一步迭代过程只需进行两次快速Fourier变换和一次对角阵与向量的乘积,计算复杂度小.本文第四章首先针对大型稀疏的Sylvester方程提出了预处理的不对称的Her-mitian与反Hermitian分裂(PAHSS)迭代法和不精确的PAHSS(IPAHSS)迭代法,给出了收敛性条件和最优迭代参数的选取.接着针对时间-空间分数阶对流扩散方程离散所得的特殊Sylvester方程,提出了带状预处理的向后代入迭代方法,对预条件子进行了相关的理论分析.最后针对时间周期的二维分数阶扩散方程离散得到的低秩的 Sylvester 方程,基于单步的 HSS(SHSS)迭代法和 Krylov-Plus-Inverted-Krylov(KPIK)子空间迭代方法,提出了 SHSS-KPIK迭代方法,讨论了新算法的相关理论结果,数值例子验证了新算法及新预条件子是有效的.
其他文献
学位
学位
学位
学位
学位
表观遗传沉默机制在植物基因表达调控,生长发育以及胁迫应答中具有重要的作用。与此同时,植物中也存在一些抗沉默因子,它们通过DNA去甲基化,组蛋白修饰和染色质重塑等途径来抑制基因的表观遗传沉默。为研究植物抗基因沉默的机制,本研究以模式植物拟南芥为材料,采用一套基于荧光素酶(LUC)的报告基因系统,筛选拟南芥抑制基因沉默的重要因子。在这套筛选系统中,报告基因LUC的启动子为双35S启动子(d35S),在
学位
学位
学位
丛代数与量子群,泊松几何,整系统等领域都有紧密的联系,特别是利用丛代数研究量子群的典范基和利用丛结构与泊松结构的相容性研究李群的丛结构.本文主要研究丛代数在量子仿射代数和泊松几何中的一些应用.全文共分为四章.第一章介绍了研究背景和预备知识.第二章主要介绍了型An和Bn的M-系统和对偶M-系统,并且介绍了 M-系统(对偶M-系统)与丛代数的关系.证明Hernandez-Leclerc猜想对于型An和