计算最小奇异三元组的隐式重启精化Lanczos双对角化方法

来源 :厦门大学 | 被引量 : 0次 | 上传用户:tonykaier
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
随着科学和技术的发展,越来越多的应用和计算问题需要求解大规模矩阵的少数几个最小奇异值及其相应的左(右)异向量.众所周知,求矩阵的奇异值分解通常用Lonczos Bidiagonatization(LBD)先把矩阵转化为双对角阵,然后用Golub-KahanSVD求出双对角阵的奇异值分解,从而求出矩阵的奇异值分解.最早LBD算法是Golub等人在文献中提出的。后来Soremen在文献(25]提出了-种有效的隐式重新启动策略并成功应用到Arnoldi方法和Lanczos三对角化方法中,这种方法可看作是有一种截断QR算法(curtailed QR).再后来Bjock等人把这种策略应用到LBD并得到相应的递推公式,接着贾仲孝通过添加“精化向量”提出了精化算法.最近Kokiopoulou等人应用LBD的递推公式并添加“精化向量”来求解-些奇异值及相应的奇异向量.本文就是在Kokiopoulou等人的工作上进行改进的. 本文的主要工作是给出-种计算少数几个最小奇异三元组的隐式重新启动精化Lallcozs双对角方法,与上述方法区别在于:我们求的是小的奇异值,因此我们用调和Ritz值作为位移,能有效地逼近大规模矩阵的最小奇异值的奇异三元组.在算法中我们还用精化残量,精化奇异向量和精化Rayleigh商,同时采取压缩技术压缩掉已经求出的小的奇异三元组.最后,我们用数值实验表明,我们的算法可以成功地求解大规模矩阵的小的奇异三元组,收敛速度也较快. 本文的结构和内容如下:第一章,我们在第一节对奇异值的背景及数值解法的研究历史给出一个比较详细的概述,第二节介绍了怎样用低维的奇异三元组来逼近大规模矩阵的奇异三元组.在第二章,我们首先介绍LBD,然看给出隐式重启LBD及讨论位移的选择问题.在第三章中我们证明正交压缩变换(ODT)的几个结论,给出精化向量、精化残量及精化Rayleigh商的求解,第四章,我们给出完整的算法,最后-章是数值实验.
其他文献
期刊
近几十年来,函数逼近在理论研究和实际应用中均获得重大进展,它不仅是数值分析的基础,同时在微分方程数值解等方面起着重要作用.具体说来,函数逼近所涉及的问题是:在某一区间
本文介绍了锚杆支护在回采巷道中的应用情况,从支护设计、施工工艺、设备配套及质量管理监测等方面做了全面论述。 This paper introduces the application of bolt support
该论文对属性值分别以区间灰数和模糊数这两种常用形式的多属性决策问题进行了深入系统的研究: (1)对于灰色多属性决策问题,给出了区间灰数夹角余弦公式,并给出灰色多属性决
本论文是对线性混合效应模型中参数的谱分解估计方法的深入讨论.我们知道,由谱分解方法得到的参数的估计有很多优良的性质.其中,对于观测向量协方差阵的谱分解估计,我们很容易得
研究型教学作为一种新的教学理念、教学方法、教学模式,是针对中国传统教学所存在的系列问题提出的创新性解决方案,是当下培养高素质创新型人才、实现创新型国家的必然要求,是高
Bi-Cayley图是一类新定义的图,它的连通性已被深入的研究.本文主要研究了一些Cayley有向图的邻接矩阵和Bi-Cayley图的一些代数性质:特征值和生成树数. 第一章介绍了背景和一些
随着Web服务技术的研究、发展和应用,众多学者试图寻找Web服务组合验证的有效方法,开始运用形式化方法描述Web服务组合并对其进行分析。目前国内外有使用进程代数、petri网、sp