大规模矩阵特征值及线性系统的Krylov子空间算法研究

《大规模矩阵特征值及线性系统的Krylov子空间算法研究》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、中文摘要摘要大规模矩阵特征值及稀疏线性系统数值求解问题已成为许多科学、工程、工业模拟，以及金融等领域的核心问题，由于其求解时间在整个问题解决上占有相当大的比重．因此，高效地解决这些大规模矩阵计算问题能够很大程度上提高整个问题的求解效率．最近20多年来这些大规模矩阵计算问题的算法研究一直是计算数学的热点，国际、国内的研究十分活跃．其求解算法，特别是Krylov子空间类型的算法得到了长足的发展．然而，仍有许多问题尚待解决。基于问题的重要性及其长远意义，本文讨论这些大规模矩阵计算问题的数值求解算法．一方面，本文研究了大规模矩阵特征值问题的数值求解算法，算法的收敛性、稳定

2、性等问题；另一方面，还讨论了大规模稀疏线性系统的加速算法、预处理技术、以及相关的收敛性分析．全文共分八章．第一章分别介绍了大规模矩阵特征值和稀疏线性系统问题的来源、研究历史、发展现状以及解决这些问题的基本方法．第二章、第三章讨论了大规模非对称矩阵部分特征对的数值计算问题．其中，第二章提出了一种计算大规模矩阵部分内部特征对的Arnoldi类型的算法．分析了算法的收敛性，与调和Arnoldi算法之间的关系．最后给出数值算例．试验对比表明；当近似子空间空间维数较小时，本文的算法能够得到更快的收敛速度；当近似子空间维数逐渐增加时，两种算法的收敛速度都明显加快，子空间维数较

3、大时，两种算法的收敛效果相差不大．第三章讨论了求解大规模矩阵特征值问题的一类带压缩向量的块类型Arn01di算法．每次重新启动时，利用上一个循环得到的近似特征向量来构造初始块向量，在正交基向量的构造过程中采用‘非精确压缩’．一方面，在新的求解子空间中算法可以包含已经收敛的近似特征向量，另一方面，通过非精确压缩，初始‘块’的列数随着近似特征对的收敛而减小，从而使得新的近似子空间对于尚未收敛的特征对更有利．因此，这种方法能够克服单向量的KryloV子空间不能计算重特征值的缺陷，而且比传统的块Krylov子空间方法更稳定，更有效．近似子空间的性质分析表明这种方法渐进地具

4、备添加特征向量重新启动的Arnoldi算法的优势，随着近似特征对的收敛，每次重新启动生成的近似子空间对未收敛的特征对越来越有利．最后给出数值试验，结果表明本文的算法能够克服块II厦门大学理学硕学位论文的Krylov子空间不规则收敛的现象，具有非常光滑的收敛性质．第四至第七章讨论的是大规模稀疏线性系统的求解问题．其中，第四章首先通过实验发现添加近似特征向量重新启动的GMRES(GeneralizedMinimumRESidual)算法(GMRES—E)每隔一步生成的残向量角度往往很小．从而提出在GMRES．E的基础上添加修正向量的一种双重增广重新启动的GMRES算法

5、(LGMRES—E)．数值试验表明新的重新启动方式可以有效的纠正残向量经常出现的跳跃角很小的现象，从而可以使每次重新启动得到的近似子空间保持适当的正交性，一定程度上克服原来方法由于重新启动所造成的近似子空间整体维数下降问题，同时由于添加近似特征向量，新方法也能够有效地压缩掉影响收敛速度的小特征值．数值试验表明了这种双重增广方法对于系数矩阵具有少数几个相对较小特征值的线性系统具有非常好的效果．新方法可以加快GMRES—E的收敛速度．第五章我们提出在GCR0一DR每次重新启动时保留部分修正向量信息，并将其添加到新的求解子空间中．这种策略可以使前后近似求解子空间保持适当

6、的线性无关性，从而可以有效的缩短GCRo—DR经常出现的收敛较慢的现象．当只保留一个最新生成的修正向量时，分析表明，算法的改进是非常自然的，只需要极小的改动就能达到目的．当需要添加的修正向量的个数大于1时，本章重点讨论了一种非精确方法．它可以保证在近似子空间维数相同时，改进后的算法与原算法的运算量相差不大．然而改进后的方法具有明显的加速现象．我们对比了两种方法在求解单个线性系统以及序列线性系统方面的收敛速度，讨论了添加修正向量个数对于收敛性的影响．数值试验表明新方法能够有效加快GCR0一DR的收敛速度，在解决模拟机械疲劳断裂问题产生的系列线性系统等问题上的数值试验

7、效果表明，新方法的收敛速度提高接近10％．第六章讨论了切频率过滤分解与组合预处理技术．首先，我们观察到传统的右侧过滤分解同样可以在满足左侧过滤条件下来完成．在此基础上，本文提出了一种双侧切频率过滤分解预处理子．在过滤向量的选取上采用Dnes来作为过滤向量．如此选取有以下几个优点，一方面，可以节省其他类似算法的前处理过程来计算过滤向量．另一方面，对于使用。佗es作为左侧过滤向量所构造的过滤预处理子，如果选取适当的初始向量，那么预处理的Krylov子空间迭代算法能保证其残向量的和始终为零，即所谓的材料均衡误差为零的性质．将切频率过滤分解所构造的预处理子与传统的ILU(

8、0)预处理