对称正定线性方程组的多级分裂预处理

来源 :长沙理工大学 | 被引量 : 0次 | 上传用户:zhanghongyingyxl
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
对称正定线性方程组在通信工程等工科领域中扮演着非常重要的角色.事实上,工科许多问题的计算最终会转化为大规模的对称正定线性方程组的求解问题.因此,设计求解这类方程组的高效算法是数值计算领域的一个重要的研究方向.   直接求解和迭代求解是求解线性方程组的主要两种方法.直接求解法一般是通过求解系数矩阵的逆来求解方程组的.当系数矩阵维数较小、逆易求时,直接法比较占优势.然而,当系数矩阵阶数较大时,迭代法一般是首选.但迭代法的收敛速度依赖于系数矩阵的条件数.因此对于系数矩阵条件数较大时,有效的预条件子构造就是很有必要的了.   本文主要研究系数矩阵为对称正定及块对称正定线性方程组基于多级分裂基础上,构造相应预条处理子的问题.首先,对(块)对称正定矩阵进行多级分裂,在分裂过程中通过设置不同的嵌套数目,推导出最终的迭代分裂公式A=MTp-NTp.然后,选取M-1Tp作为方程组共轭梯度法的预处理子进行加速,并证明了M-1Tp作为预处理子的合理性.对于块对称正定线性方程组,可类似地构造相应的块预处理子.最后我们做了一些数值试验,并对比了多级迭代法,共轭梯度法,循环预处理共轭梯度法,和多级分裂预处理共轭梯度法四种方法的计算结果,试验显示了多级分裂预处理子在计算这类方程组上的优势.本文一共六章,具体如下:   第一章介绍了多级分裂预处理的研究背景、研究现状及本文的创新之处,并介绍了文章的主要研究内容.   第二章介绍了相关预备知识,如多级分裂法、块多级分裂法以及共轭梯度预处理技术.其次,介绍了一些基本定理、引理和推论.   第三章通过系数矩阵的多级分裂构造出预处理子M-1Tp,并证明了该预处理子的合理性.   第四章主要是将预处理子的构造法延伸到块对称正定线性方程组,并证明了该预处理子的合理性.   第五章主要是进行数值试验,对比了级迭代法,共轭梯度法,循环预处理共轭梯度法,和多级分裂预处理共轭梯度法四种方法的数值算例.且数值试验显示了多级分裂预处理子的有效性.   第六章总结了本文的主要工作,并对该算法在Toeplitz矩阵上的应用进行了探讨.
其他文献
匹配理论是图论的主要研究专题之一,并且与其他理论课题具有密切联系.鉴于n-可扩图、导出匹配可扩图、PM-紧邻图的研究工作,我们提出两个新的概念:圈唯一可扩图和导出圈可扩图.
本文介绍了左对称代数的概念和两个特殊的左对称代数结构。一个是A Dzhumal dild daev和C.Lof wall给出的在根树上的一个生成元的自由左对称代数结构,另一个是由A Connes和D.K
约束矩阵方程问题是指在满足一定约束条件的矩阵集合中求矩阵方程(组)的解的问题.作为当代数值代数领域中的前沿方向,约束矩阵方程问题以及其延伸出来的新问题普遍应用于结构设
作为整数阶阻尼系统的推广,分数阶阻尼系统具有更加广泛的实用价值.近年来,分数阶微分方程的可控性已被学者大量研究,并且得到了许多结果,研究成果被广泛地应用于生物、工程和化
随着数学在实际应用中的深入发展,数学物理反问题越来越成为一门重要的学科。近些年来,反问题的应用渗透在工业设计、地球物理、材料科学、生物医学等各个领域中,但是由于其
主要研究了具有不连续系数的三阶奇摄动边值问题.首先研究了具有不连续系数的三阶拟线性奇摄动边值问题(此处公式省略)  式中式B,D是给定的常数,(此处公式省略)  而函数f1,f
本文研究了混合KN与AKNS方程族的超Bi-Hamilton结构及无穷守恒律。文中首先从3×3矩阵谱问题出发,借助零曲率方程导出与这个谱问题相关联的一组非线性微分方程,其中第一个非平
随着土壤盐碱化和土壤污染的日益加剧,有关土壤溶质问题的研究引发了众多学者的广泛关注。寻找高效的科学方法来控制和调节土壤溶质浓度对土壤养分管理及环境污染监测问题具