大规模稀疏线性方程组的预条件迭代法的研究

来源 :南昌大学 | 被引量 : 0次 | 上传用户:jonelove0000
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
在科学计算与工程应用领域,如核能工业、石油工业、电路计算机辅助设计和分析、偏微分方程数值解、图像处理等,许多问题的计算最后往往归结为大规模线性方程组的求解,而这也恰恰是计算中最耗时的部份.因此,设计求解线性方程组的有效算法是大规模科学计算领域一个非常重要的课题.   众所周知,求解线性方程组的数值方法主要有两类,即直接法和迭代法.直接法主要是基于将系数矩阵分解成容易求逆的矩阵.当系数矩阵为小型稠密的矩阵时,直接法比较受欢迎.然而,当系数矩阵为大规模稀疏矩阵时,迭代法往往被使用.迭代法在求解大规模线性方程组时较直接法有许多优势,但是迭代法也存在一个普遍的缺点,即收敛速度较慢.因此,结合各种预条件技术的迭代法得到了广泛的研究和应用.   本文主要研究大规模线性方程组的预条件迭代求解算法.首先,利用离散余弦变换给出了对称Toeplitz线性方程组的一个预处理子;其次,结合Krylov子空间方法和预处理技术,给出了变预处理子SOR-双共轭残量算法,理论和数值实验验证了算法的有效性;由于HSS迭代在求解非埃米特正定线性方程组时非常有效,因此在变预处理子SOR-双共轭残量算法的基础上,研究并给出了求解非埃米特正定线性方程组的变预处理子HSS-广义共轭残量法;最后,研究了下三角Toeplitz矩阵的数值求逆问题,给出了一个基于离散正弦变换的快速求逆算法.本文共分六章,组织如下:   第一章介绍了求解大规模稀疏线性方程组的预条件迭代法的研究背景、研究现状及相关预备知识,同时介绍了本文的主要研究内容.   第二章考虑了对称Toeplitz线性方程组的求解问题,研究给出了一个基于离散余弦变换的预处理矩阵,同时给出了该预处理矩阵的快速构造算法.数值实验显示了该预处理矩阵的有效性.   第三章结合Krylov子空间方法和预处理技术,给出了变预处理子SOR-双共轭残量算法,并且给出了算法的收敛性分析.通过数值实验,验证了算法的有效性.   第四章在变预处理子SOP-双共轭残量算法的基础上,研究并给出了求解非埃米特正定线性方程组的变预处理子HSS-广义共轭残量法.   第五章研究了下三角Toeplitz矩阵的数值求逆问题,给出了一个快速求逆算法.   第六章对全文的工作进行了总结,并对今后的研究方向作了一些展望.  
其他文献
安全协议是实现信息安全的基础,是网络安全通信的核心技术,它的正确性对网络的安全起着非常重要的作用,因此其自身的安全性问题已成为安全研究的重要内容。然而怎样保证安全
谣言传播是网络上的一个经典问题,文章共分为三部分。首先介绍了复杂网络的概念以及研究动态。然后在经典谣言传播模型的基础上,研究了社会网络上的谣言传播行为,并建立了数学模
本文研究几类不连续动力系统-脉冲微分系统,时滞脉冲系统以及时滞脉冲切换系统的稳定性。   第一章介绍了非线性系统以及脉冲系统的稳定性理论,阐述了研究该类问题的理论
当CT系统在数据采集过程中受到严重的噪声干扰,或者投影数据采集不完全时,用解析重建算法重建出的图像有伪迹。而统计重建算法具有物理模型准确、对噪声不敏感、易于加入约束
<正>日前,央视新闻频道推出八集系列主题报道《大国工匠》,聚焦八位行业顶级的普通技工,讲述了他们立足岗位、追求卓越的感人故事,从中挖掘我国火箭、高铁、国产大飞机、LNG