关于不完全正交三角分解的算法构造与数值研究

来源 :中国科学院数学与系统科学研究院 | 被引量 : 0次 | 上传用户:wjz5201
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
快速、准确且稳定地求解大型稀疏、非奇异非对称线性代数方程组是科学与工程计算研究领域中的最基本问题之一。研究关于这类问题的预处理迭代算法,特别是预处理子的构造和性质,是一项具有重要的理论意义和很高的实用价值的,且十分有挑战性的课题。   本文对基于Givens变换的不完全正交三角分解方法的理论和性质进行了进一步发展和深入的研究。结合矩阵的具体结构和性质,本文作者提出了比较系统、高效且实用的不完全正交三角分解算法,并进行了深入的理论分析和大量的数值实验。   理论结果表明,基于Givens变换的不完全正交三角分解方法数值上稳定,且总能得到一个正交矩阵Q和一个上三角矩阵R,使得A=QR—E,其中E为满足一定约束的误差矩阵。特别,该正交矩阵Q可以存储为若干Giveils变换;这给相应的迭代算法的预处理过程的具体实施带来了很大的方便。对于正规方程组的预处理共轭梯度法而言,由于Q的正交性,在构造预处理子的过程中,不需要存储矩阵Q,只需要在迭代的每步求解两个稀疏的三角方程组,再完成相应的广义残量方程的求解过程。   为了节省存储空间,首先提出了一类不引入非零元素的不完全正交三角分解算法,称之为IGO(0)方法。构造了基于该方法的预处理迭代算法,并对离散的偏微分方程和矩阵市场中的例子做了大量的数值实验。结果表明,IGO(0)方法是有效的,并且由它所诱导出的预处理迭代算法比相应的不完全三角分解方法具有更快的收敛速度。   其次,改进了不完全正交三角分解方法;通过将丢弃的非零元素以一定的策略补偿到对角元上,得到了MIGO方法。本文还给出了在不完全正交三角分解过程中关于对角补偿元素的数值选取的一个一般性公式。特别,对于求解离散的常系数偏微分方程,MIGO方法相当于给不完全正交三角分解方法增加了一个额外的强制约束。对于某些实际问题,这种预处理方法可以更加准确而快速地得到它们的数值解。数值实验也进一步佐证了这一点。   再次,利用矩阵的稀疏结构,提出了一个预先指定稀疏结构的、一般的、不完全正交三角分解(GIGO)方法。特别,构造了一个以矩阵的半带宽l为参数的不完全正交三角分解方法。该方法允许引入主对角线附近的宽为2l+1的带状区域内的非零元素,称之为IGO(l)方法。如果参数l足够大,则可能会得到一个完全的正交三角分解。对于一些实际问题,本文做了大量的数值实验,并对IGO(l)方法中参数l的选取进行了详细的讨论。   之后,提出了一个在不完全正交三角分解过程中考虑非零元素的大小,且允许引入数值较大的非零元素的不完全正交三角分解(IGO(τ))方法。考虑到计算机的实际存储和计算能力,本文通过引入另外一个参数p来限制非零元素的个数的增长。这种新的更为实用的方法被称为IGO(τ,p)方法。理论分析表明,IGO(τ)方法和IGO(τ,p)方法都是数值稳定的。大量的数值实验表明,IGO(τ)方法和IGO(τ,p)方法是可行且有效的。值得指出的是,在实际计算中,可以通过选择参数τ和p使总的迭代时间达到最优。   最后,对于最小二乘问题,将上述不完全正交三角分解方法推广到了长方矩阵的情形。理论分析表明,这类新方法同样是数值稳定的,且能够保持原不完全正交三角分解方法的优点。本文构造了用新方法做预处理子的,求解最小二乘问题的共轭梯度法和广义极小残量法。大量的数值实验表明,这种预处理方法是可行且十分有效的。
其他文献
磁流体动力学(MHD:magnetohydrodynamics)在科学研究和工程应用中都有着很重要的意义,对理想MHD方程和浅水磁流体方程(SMHD:shallowwaterMHD)进行数值模拟有助于了解导电流体的
随着科技化时代的快速发展,在不断探索自然界物种的生理结构和活动行为过程中,人们受到很多启发.大脑是人体十分重要的器官,它能够进行信息处理和信息加工,同时也能完成复杂的学
时频分析作为分析时变非平稳信号的有力工具,成为现代信号处理研究的一个热点。它的主要任务是表述信号的频率成分随时间变化的规律,并进一步建立一种时频分布,其能够在时间和频
随着科学的进步,信息技术在当今教学过程中占据着举足轻重的地位.信息技术发展飞速,教学方式灵活多样.作为网络时代下的信息技术教师,我们在教学过程中要立足课堂,创造有利条
本文主要讨论的是多元多项式的近似因式分解和有关单变元和多变元多项式的近似最大公因子的一些问题。所得到的主要结果包括三个方面:   第一,计算多变元多项式的近似因式分
凭借最新的PMI(Pilz人机界面)6 Control产品系列,Pilz自动化公司推出了配备PLC软件(符合IEC 61131-3标准)的首款操作终端。此终端不仅可以提供专业诊断和可视化,还可以在一个
在一些实际问题中,群体的数目是已知的.但在有些科学研究的统计问题中,群体数目和相关参数是未知的,这就需要进行估计.本文研究了一个带异质性的捕获移出frailty模型,在个体间
本文研究的是二维调和映照热流方程与四维Yang-Mills热流方程解的长时间行为,主要探讨内容如下所示:  第一章为绪论,简单介绍了二维调和映照热流方程与四维Yang-Mills热流方
实代数曲线曲面的拓扑结构确定和表示不仅是一个有趣的数学问题,也是计算机图形学和计算机辅助几何设计的一个核心问题。论文的内容由三部分组成:三角列多项式系统的实根隔离算
为了使Hilbert变换在信号分析的应用中具有坚实的数学基础,本文研究了乘积函数的Hlibert变换问题。在前人研究结果的基础上,给出了L2(R)上的Bedrosian恒等式成立的新的充要条件