论文部分内容阅读
许多工程和科学计算的数值求解均可转化为大型线性和非线性方程组的求解.如不可压液态流问题的数值求解,具有广泛应用的PDE约束问题以及病态的反问题数值求解等.通过适当的数值离散,它们最终均可转化为具有某些特殊结构的大型线性方程组.这些线性方程组一般具有某种特殊结构,如块结构:2×2块或3×3块结构,大型、稀疏病态,有些甚至是奇异的.高效快速求解这类线性系统一直是科学计算的核心问题之一,具有非常重要的意义.面对实际问题中涉及的形态各异的大型系数结构线性系统,如何高度利用问题自身的结构和特性来构造稳定、高效的算法是现代数值计算的重要课题和研究热点之一.本论文分共有七章,主要研究离散PDE约束优化问题、几类PDE及一些不适定问题中涉及的大型线性和非线性方程组的求解,重点针对具有块2×2结构的线性系统,提出了一系列具有针对性的快速、稳定的迭代算法及预处理技术.第一章介绍了问题的研究背景、研究意义以及研究现状,并简单介绍了本文的主要研究内容和创新点.第二章主要研究由几类特殊的偏微分方程(如Helmholtz方程)离散得到的复对称线性方程组的数值求解问题.针对复线性方程组的实等价块2×2结构的线性系统,得到了两类具有更好特性的块预处理子.所得到的块预处理子在某些特定条件下有更聚集的特征值分布,理论证明了预处理矩阵特征值的界与问题参数无关,体现了它的鲁棒性和有效性.数值例子进一步验证了理论的正确性以及新预处理子的有效和鲁棒性.第三章主要研究了由PDE约束优化问题导出的离散结构线性系统的快速算法.针对离散得到的特殊的块2×2结构复线性方程组进一步研究了目前最有效的三个块预处理子,发现并导出了三个预处理子间的关系,即相应的预处理矩阵的特征值都是受同一矩阵的影响.通过讨论最优参数的选取得到了三个预处理的最有效的形式.此外,基于多步预处理的思想,又构造了两个两步预处理子.理论分析及数值实验都证明了新的两步预处理子比现有的预处理子有更聚集的谱分布,且与问题参数无关,说明了所提出的预处理子的有效性和鲁棒性.第四章针对图像复原问题所得到的结构化线性系统我们提出了有效的数值算法.基于Pan等人(SIAM,2014)提出的一类新的HSS(NHSS)预处理子,通过引入两个迭代参数,提出了一个广义的两半步交替方向(GNHSS)的迭代方法,理论分析了该迭代方法的收敛性质以及预处理矩阵的谱性质,讨论了双参数的选择方式.数值试验表明GNHSS迭代法及相应预处理子都是非常有效的.第五章主要针对不可压缩Navier-Stokes方程导出的鞍点结构的线性系统,设计了一种有效的预处理子,研究了预处理矩阵的谱性质,给出了新的预处理子与已有预处理子的算法复杂性比较,表明给出的预处理子在计算方面具有很大的优势.此外,我们还讨论了最优参数的选择并给出了一种比较实用的参数选择方法.第六章针对一般广义鞍点问题,考虑了一种块乘积型高效预处理子及迭代方法.研究了迭代方法的收敛性条件以及预处理矩阵的相关性质,讨论了参数的选择办法.数值实验结果表明新预处理子及参数选择办法是非常有效的.第七章对本文做了简要总结并对未来的工作进行了展望.