求解大规模非对称矩阵特征值问题的加权Arnoldi方法

来源 :南京航空航天大学 | 被引量 : 2次 | 上传用户:zuiainivip
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
求解大规模矩阵特征值问题是当今科学与工程计算的热点之一。最近一、二十年来,在大规模非对称矩阵特征值问题的数值求解方面已经取得了较大的进展。但是,仍然有一些问题没有得到彻底的解决。正因为如此,本文在如下方面做了一些工作:首先,本文介绍了大规模矩阵特征值问题的来源和解决这类问题的一些经典数值方法,并对问题的研究现状作了简要的回顾。Arnoldi方法非常适合计算矩阵的端部特征对,对于实际问题中产生的大规模矩阵内部特征问题,这种方法却常常失效。Shift-and-invert Arnoldi方法对于大规模矩阵内部特征问题的求解非常有效,且能达到快速收敛。但是由于在计算过程中需要对大规模矩阵进行分解,一方面会破坏原矩阵的稀疏结构,增加存储量;另一方面,Arnoldi方法在迭代过程中对矩阵的扰动非常敏感,而求逆过程很难保证达到可靠的精度,因此这种方法被认为是不实用的。由Morgan提出的调和Arnoldi方法目前被认为是求解大规模矩阵内部特征问题的最有效方法之一。为了进一步减少该方法的存储量,加快收敛速度,改善精度,本文第二章我们针对调和Arnoldi方法进行加权,得到加权调和Arnoldi方法,并且考虑权矩阵的选择,对它进行理论分析和数值试验,表明了新算法的优越性。由Arnoldi算法的过程知,形成Krylov子空间的一组标准正交基时需要将新产生的向量与已产生的所有基向量正交,同时还要计算一个m阶Hessenberg阵的全部特征对,其存储量很大甚至无法实现。因此,在实际计算时,由内存、收敛速度和计算速度等方面考虑,Krylov子空间的维数m不应该太大。而对于选定的m ,用Arnoldi方法计算的Ritz值和Ritz向量又未必满足精度要求。克服这一困难的办法是用一个Ritz向量或几个Ritz向量的线性组合代替初始向量v1重复使用Arnoldi方法,即采用重新开始技术。本文第三章讨论几种重新开始Arnoldi方法,并运用加权的思想提出一种新算法,通过数值试验验证其收敛性和精度。
其他文献
“西厢调”是流传于甘肃白银、榆中一带的一种民间曲子戏的特定称谓。“西厢调”集歌唱、表演与器乐演奏于一体,以表演短小精悍的故事为主要内容。从其曲牌、唱腔来看,属于曲
很多朋友反映,有时移动光电鼠标,光标的反应不像往常那般迅捷,感觉不听指挥了,这是因为鼠标灵敏度下降,最常见的可能就是透镜通路脏了,光线不能顺利到达。光电鼠标的工作时间
网易邮箱(mail.163.com)推出极速3.0版有一段时间了,笔者在使用过程中,发现极速3.0版有一些快速操作技巧,可以为我们通过邮件与好友联系提供方便。打开163邮箱登录页面,单击
一、资产取得时的计价资产取得时的计价可分为以公允价值入账和以账面价值入账两种方法,对于支付现金、接受投资、接受捐赠、融资租赁和盘盈等5种取得方式,以资产的公允价值
2016年7月22日,翠屏区经信科技局党委副书记赵启宽率局相关人员一行赴宋家镇大地村、洋坪村开展帮扶慰问活动,为21户建档贫困户发放了第二批鸡苗共420只,专业技术人员将不定期对
镍多金属矿的地形测量工作包括矿区地形测量、野外地形地貌测量两项,根据地形测量结果,可以为找矿工作提供丰富的资料与参考依据。本文以某地区镍多金属矿为例,首先分析了矿区地