线性方程组和鞍点问题的迭代法与预处理技术研究

来源 :电子科技大学 | 被引量 : 4次 | 上传用户:linxinrudo
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
大规模科学计算和工程技术中许多问题的解决,最终归结为大型稀疏线性方程组的求解,其求解时间在整个问题求解时间中占有很大的比重,有的甚至达到80%.由于现今科学研究和大型项目中各种复杂的课题对计算精度和计算速度的要求越来越高.因此,作为大规模科学计算基础的线性代数方程组的高效数值求解引起了人们的普遍关注.这种方程组的求解一般采用迭代法,所以,迭代法的收敛性和收敛速度就成为人们关注的焦点,为许多专家和学者所研究.本文对与大型稀疏线性方程组迭代求解有关的特殊矩阵迭代法进行了深入和系统的研究,特别研究了松弛型矩阵多分裂迭代法的收敛性,两种Krylov子空间方法的性能和鞍点问题的预处理技术.全文共六章,分四个部分: 研究了H-矩阵的松弛型矩阵多分裂法,并给出详细的理论分析和敛散速度的比较.一方面,给出了松弛型矩阵多分裂TOR法,研究了方法的收敛性,比较了他们的敛散速度,分别进行了串行和并行试验,验证了所提方法的优越性。另一方面,给出了松弛型矩阵多分裂USAOR法,研究了方法的收敛性,给出了实现算法的例子,并用数值试验与存在的方法进行了比较. 进一步研究了一些H-矩阵松弛型矩阵多分裂法新的收敛性结果.分别为非线性方程组的非定常矩阵多分裂法,线性互补问题的矩阵多分裂法,松弛型矩阵多分裂SSOR法和松弛型矩阵多分裂TOR法,构建了相应方法的收敛性理论,得到了新的更弱的收敛性条件,进行了数值试验的比较.基于多分裂法的并行性,矩阵多分裂的研究和理论分析对于多分裂预处理子的构造有一定的理论和应用价值.我们的方法选取参数的余地更大,当选取近似最优参数时,能实现更快的收敛速度和构造出更有效的预处理子. 基于BiCR算法设计了求解非对称线性方程组Krylov子空间平方共轭残差(CRS)算法和适合分布式并行计算的改进的平方共轭残差(ICRS)算法,并对两种算法进行了理论分析和算法比较,串行和并行数值试验表明所提方法具有较好的收敛速度和并行性能. 研究了鞍点问题特殊矩阵迭代求解预处理技术.首先。对内点优化问题产生的一类鞍点问题给出了一种预处理技术,进行了相应的理论分析和数值试验.接着,基于离散化混合型时谐Maxwell方程的块三角鞍点问题和特殊矩阵的结构,提出了带多个参数的预处理技术,并进行了理论分析和参数的理论选取.理论分析表明提出的预处理子有更好的特征值聚集性。数值试验也表明本章提出的预处理子性能大大优于免增广和免Schur余块对角预处理子.
其他文献
圆填充是常曲率曲面上具有特定相切模式的一种圆格局。其理论在复分析和离散微分几何的交叉学科中是一个快速发展的研究领域。网填充理论起源于费尔兹(Fields)奖获得者W.Thur
框架理论发展至今已广泛应用于光学、信号处理、图像处理、数据压缩、采样理论等领域。由于再生核空间具有很多良好的性质,再生核函数对研究再生核空间中的性质起着至关重要的
代数表示论高维理论是Iyama等人推广经典Auslander-Reiten理论,引入n-Auslander代数[63],n-Auslander-Reiten平移函子[64]等建立发展起来的。作为平移代数的推广,郭引入了n-平移
非方常数表示空间的非方状态,它们的取值与一致正规结构和空间的一些其他几何性质密切相关。空间几何常数的表示与计算能够更精确的描述空间的性质。本文主要研究了非方常数的
本文主要研究了一类系数为非常数的时空白噪声驱动的双边反射随机偏微分方程和一类系数为常数的、关于时间是分式的分式噪声驱动的双边反射随机偏微分方程,而且双边反射壁为光
学位
强自吸收的C*—代数是一类具有特殊性质的C*—代数,因此对其性质的归纳、总结、优化对于研究其它代数有重要意义。 本文主要研究了强自吸收C*—代数的一些性质和C*—代数强
Loewner微分方程是单叶函数中的一个重要内容,它被证明是解决单叶函数中极值问题最有用的工具之一。为了研究统计物理中一些模型的Scaling极限,1999年,Schramm建立了一个随机版