稀疏矩阵方程组预处理迭代技术研究

来源 :华南理工大学 | 被引量 : 0次 | 上传用户:houwenjin
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
大型稀疏线性方程组的求解是对自然科学和社会科学中许多问题进行数值模拟的关键技术之一。而GMRES算法是目前求解大型稀疏非对称线性方程组最为有效的迭代算法之一。在执行整体的GMRES算法时,所需的计算量和存储量会随着迭代步数的增加而变得不可接受。为了克服这一困难,可以使用重新开始策略或混合迭代策略。最近,重新开始GMRES算法在迭代过程中表现出的补足收敛性质引起了人们的兴趣。特别地,基于这一性质所提出的积混合GMRES算法能够显著改善混合迭代策略求解方程组的效率。积混合GMRES算法在执行过程中,需要首先计算出多次迭代GMRES迭代循环的残量多项式,然后重复使用这些多项式的乘积进行Richardson迭代。然而,当迭代循环的步长较大时,计算出的残量多项式可能是不稳定的,从而导致Richardson迭代的发散。为了提高积混合GMRES算法的稳定性以及收敛速率,本文提出了多项式预处理积混合广义极小剩余算法。首先,本文介绍了求解大型稀疏矩阵的预处理Krylov子空间方法的原理以及在此基础上发展起来的各种迭代法,包括共轭梯度法,广义极小剩余法。其次,本文重点介绍了在重新开始GMRES循环的残量多项式在矩阵的谱上收敛的互补性以及在此基础了得到的积混合广义极小剩余算法,构造出了多项式预处理矩阵,并将该矩阵作为积混合广义极小剩余算法的预处理矩阵,改善其系数矩阵谱的性质,提高了该算法的收敛速率和稳定性。最后,本文对预处理后的新算法做了数值实验模拟与分析,将新算法与经典的成熟算法进行了对比,结果均表明,新算法更适合大型稀疏矩阵问题的求解,在计算量和存储量方面都有相应的改进。求解大型稀疏矩阵的积混合广义极小剩余算法得到了进一步的改善。
其他文献
<正> 目的本试验研究肾脏移植术后应用免疫抑制剂治疗,比较应用普乐可复和环孢素的临床疗效。方法48例患者均为首次移植,按服用免疫抑制剂的不同分为普乐可复组和环孢素组。
会议
在阐述煤化工废水来源、特征与危害的基础上,对目前治理手段进行了综述,介绍了煤化工废水“零排放”技术,并提出了煤化工废水治理存在的问题与展望,从而希望为煤化工废水治理
在切削加工过程中,切削温度是影响被加工工件的表面质量和刀具磨损的重要因素,因此对切削温度的研究有着重要意义。以人工热电偶为温度传感器,设计了热电偶的安装方法和以P87LPC
从国家统计局广东调查总队了解到,广东省城乡居民收入在去年实现新突破:农村居民人均纯收入首次过万元,城镇居民人均可支配收入首破3万元大关。
以葡萄霜霉病生防菌株BCJB01和BMJBN02为供试菌,分别在山东省济南和烟台地区的葡萄园对不同葡萄品种的霜霉病进行了小区防治试验,以期明确BCJB01和BMJBN02对不同葡萄品种霜霉
通过密度秤的制作使用 ,引导学生钻研一些物理问题 ,加深对物理知识的理解。1 密度秤的构造  密度秤的构造如图 1 ,它与普通的杆秤类似 ,但比杆秤更简单。它有秤杆、秤钩
汉语口语课堂教学的任务是训练学生汉语口头表达能力,目的是培养学生用汉语做事的能力,研究设计行之有效的汉语课堂教学评价标准尤为关键。构建汉语口语课堂教学评价应具有导
作为全球碳排放总量最多的国家,中国面临碳减排的国际压力与日俱增,而中西部欠发达地区所处的工业化和城市化的阶段性特征又使得中国碳减排更加艰难,因此,中西部欠发达地区碳
何谓精密司法?大凡了解日本刑事司法者,大多知道“精密司法”一词。第二次世界大战之后,日本刑事司法在很大程度上受到美国的影响。尽管日本司法 What is fine justice? Alm