论文部分内容阅读
科学与工程计算的很多领域,诸如计算流体力学、约束优化、计算电磁学、PDEs的混合有限元近似、非线性规划、中子输运理论等问题的求解,最终都可归结为大型稀疏代数系统的求解.因此,对大型稀疏代数系统的求解研究就具有非常重要的理论意义和实际应用价值.由于许多实际问题产生的大型稀疏代数系统往往具有某种特殊结构,对具有这些特殊结构的大型稀疏代数系统的数值求解研究引起了国内外众多专家和学者的关注.本文对几类特殊的大型稀疏代数系统的数值求解方法进行了深入系统地研究.特别研究了求解线性鞍点问题的迭代法和预处理技术.求解非线性鞍点问题的迭代法以及求解代数黎卡提方程的迭代法.本文共六章,分四个部分:
研究求解线性鞍点问题的Uzawa类迭代法.首先,提出了一个修正的非线性Uzawa迭代法,讨论了算法的收敛性,并给出了理论和数值比较,数值实验也验证了修正方法的有效性.其次,给出了求解(2,2)块不为零的非对称广义鞍点问题的GMLHSS迭代方法,并探讨了算法的收敛条件.最后,给出了非精确Uzawa方法、GSSOR方法和MLHSS方法求解奇异鞍点系统时的半收敛性分析.
研究求解非对称鞍点问题的预处理技术.首先,给出了含参数的广义非精确块三角预条件子,对预处理矩阵的特征对性质给出了分析,并给出了预处理矩阵的特征值扰动分析,其次,基于系数矩阵的部分HS分裂和PS分裂,提出了PHSS预条件子和PPSS预条件子,并详细研究了预处理矩阵的谱性质,指出对于一个充分小的正参数,预处理矩阵特征值聚集在两个点附近:一个是(0,0)点,另一个是(2,0)点,并且通过大量的数值实验验证了理论分析和此两类预条件子的有效性.最后,给出了SIMPLE预条件子,利用特征值理论,研究了预处理矩阵两种不同表达形式的谱之间的联系.
研究求解非线性鞍点问题的迭代法.在求解线性鞍点问题的迭代解法基础上,给出了几个求解非线性鞍点问题的迭代法,并对方法进行了收敛性分析,数值实验验证了所提算法的有效性.
研究求解代数黎卡提方程的迭代法.事实上,输运理论中的代数黎卡提方程可以写成与之等价的向量方程,首先,基于松弛思想和牛顿方法,提出了求解向量方程的松弛Newton-like方法,并给出了算法的收敛性分析.其次,利用拟牛顿思想,结合已有的牛顿类型方法,给出了两个修正牛顿方法求解向量方程,并得到了算法的收敛结果,数值实验也表明所给出的三个算法能有效改进和提高已有算法的收敛性.