论文部分内容阅读
科学与工程的很多领域如流体力学,高阶微分方程求解,计算电磁学,最优化问题和油藏模拟等都涉及到大规模稀疏线性代数系统的求解.大规模稀疏线性代数系统求解方法的研究甚至是大规模科学与工程计算的核心问题之一,具有重要的理论意义和实际应用价值.本文对与大规模稀疏线性代数系统迭代求解有关的特殊矩阵数值特征和求解方法进行了深入的研究.特别地,研究了矩阵分裂迭代法的收敛性和比较理论及迭代求解预处理技术.全文共六章,分四个部分:
第一部分研究了严格对角占优M—矩阵逆的无穷大范数的上界估计.利用严格对角占优M—矩阵的特殊结构,逆M—矩阵和M—矩阵元素之间的关系,得到逆M—矩阵无穷大范数上界估计.进一步,得到对其谱半径的估计.
第二部分研究了鞍点问题迭代求解预处理技术.首先,提出了松弛不精确Uzawa算法和预条件的Uzawa算法,这两种算法扩充原有算法,并且讨论了算法的收敛性,数值实验验证了这两种算法的有效性.其次,建立了两类求解(1,1)块矩阵为高奇异对称和非对称鞍点问题的预条件子,深入研究这两类预条件子的谱性质,通过数值例子说明所建立的预条件能很有效的解决此类病态鞍点问题.最后,对混合型时谐Maxwell方程离散得到的线性方程组,根据其系数矩阵特殊性质提出了带参数的免增广和免Schur余块三角预处理技术,理论分析说明其构造及应用代价和已有的免增广和免Schur余块对角预条件子相当,但有更好的特征值聚集性质,特别是在给出的相对最优参数的时候更聚集.数值实验说明其性能大大优于免增广和免Schur余块对角预处理技术,而且也验证了在理论上给出最优参数情形下性能最佳.
第三部分,首先研究了矩阵多分裂方法的收敛性和比较理论,提出了κ+1参数的非定常多分裂算法和κ+1参数的非定常二级多分裂算法.详细地讨论了在系数矩阵为H—矩阵时算法的收敛性,也研究了通过不同零模式得到的不完全LU分解形成的多分裂矩阵对上述算法的收敛性,并且利用迭代矩阵构造出预条件子.其次,对块三对角H—矩阵,根据其特殊结构和性质构造块LU预条件子,数值实验显示此预条件子是非常有效的.
第四部分,首先分析了两类修正预条件子结合Gauss—Seidel方法和SOR方法对系数矩阵为L—矩阵情况的收敛性,并给出了比较结果,进而得到这两类预条件子的最优结构.其次,研究了Mixed—type分裂方法对系数矩阵为Z—矩阵的线性方程组的求解方法.最后,提出了松弛的交替迭代方法,并对系数矩阵为单调矩阵,Hermitian正定矩阵的情况下研究了其收敛性.