收缩和扩张Krylov子空间方法

来源 :复旦大学 | 被引量 : 0次 | 上传用户:nayinian2688
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文研究解大规模稀疏线性方程组的收缩和扩张Krylov子空间方法。在科学计算中,尤其是在解大规模稀疏线性方程组时,Krylov子空间方法显示出与众不同的有效性。当矩阵是对称正定时,常用的方法是具有短递推的共轭梯度方法(CG)。但是在许多情况下,系数矩阵不是对称的,这时常用的方法中有完全正交化方法(FOM)和广义最小残量方法(GMRES)。矩阵的非对称性导致这两种方法不具有短递推的性质。由于存储量和计算量的限制,这两种方法通常需要重开始。研究表明,如果系数矩阵具有模很小的特征值,那么Krylov子空间方法一般会收敛得比较慢。对重开始方法来说,Krylov子空间维数比较小,有时并不含有跟模很小的特征值对应的特征向量,或者不含有相应的好的近似向量。因此,重开始方法收敛得更慢,甚至会停滞。收缩和扩张的Krylov子空间方法正是因为这个原因而被研究者提出来。其基本思想是用跟模最小的特征值对应的近似特征向量扩张Krylov子空间,以达到收缩小特征值,从而加快收敛速度的目的。本文对收缩和扩张Krylov子空间方法作了全面的介绍,并研究了它们的收敛性。本文根据前人的思想提出了解广义Sylvester方程的完全正交化方法和最小残量方法。在此基础上把重点放在应用收缩和扩张Krylov子空间技术于Sylvester方程和广义Sylvester方程。近年来,许多人对如何快速求解这两个方程作了深入的研究,提出了不同的方法。但是据作者所知,本文提出的方法应该是解Sylvester方程和广义Sylvester方程的第一个加速方法。本文所使用的近似解空间是由两个扩张的Krylov子空间作Kronecker积得到的子空间。解空间的基表示为这两个扩张的Krylov子空间的基的Kronecker积,称为Kronecker乘积基。这种方法在具有加速收敛的同时,也比应用于线性方程组的通常的扩张Krylov子空间方法需要少很多的存储量。非常适合大规模Sylvester方程和广义Sylvestcr方程的求解。
其他文献
20世纪的钢琴音乐反应着所有创作思想、技术、美学观念的变革与尝试。普罗科菲耶夫即将钢琴处理为强烈的打击性,又具有在高度不协和背景上的抒情性。练习曲是我们学习钢琴的
旗帜鲜明讲政治是我们党作为马克思主义政党的根本要求,是我们党克服艰难险阻、从胜利走向胜利的重要保证,也是我们党在革命性锻造中走向成熟的突出优势。党的十九大把党的政
<正>党的十九大作出了实施乡村振兴战略的重大决策部署,绘就了新时代"三农"事业新征程的宏伟蓝图。实施乡村振兴战略,是当前和今后一个时期"三农"工作的总抓手。做好乡村振兴