解大规模非对称矩阵特征问题Lanczos方式的两种精化版本.pdf

摘 要- 多科学与工程计算中，经常必须数值求解大规模矩阵特征问 题，理 析和大量的数值实验己经表明了求解此类问题的经典正交 投影、斜投影方法都存在着Ritz值收敛而Ritz向量可能不收敛的隐 患。为此，贾对经典的正交投影和斜投影方法进行了重大革新，使 用精化向量取代 Ritz向量作为待求特征向量的近似，由此导出的精 化投影方法不但克服了原有方法的隐患，使得算法的可靠性增强， 而且收敛的速度也大大加快。 进一步将精化策略和求解大规模矩阵问题的许多其它重要技术 或方法相结合，研究和开发出更多更高效的新算法，是一个十分重 要而且迫切需要解决的课叶本文吐是在这一思想的引导下夕灌下 述几个方面进行了研究: 1.将精化投影方法的思想和拟精化思想相结合，提出了拟精化 双正交Lanczos方法，建立了拟精化近似特征对和精化近似特 征对对应的残量范数之间的关系，给出了拟精化双正交 Lanczos算法 ，并进行了数值实验，具体算例表明，该算法比 双正交Lanczos算法优越; 2.根据精化投影方法的思想，提出了半精化双正交 Lanczos方 法，建立了半精化特征对和精化特征对对应的残量范数之间 的关系，给出了半精化双正交Lanczos算法，数值实验表明， 该算法比双正交Lanczos算法和拟精化双正交Lanczos算法 优越。 吮 关键词:非对称矩阵决征值问题，Ritz省Ritz向瑟正交 投影方法，斜投影方法，精化投影方法 Krylov子空间，双正交Lanczos 方法，拟精化双正交Lanczos方法，半精化双正交Lanczo。方洲 收敛性。 AMS分类:65F10,65F15,15A18. Abstract Largeunsymmetriceigenproblemsariseinmanyapplications inscientificandengineeringcomputing.Theoreticalanalysis andnumericalexperimentsshowthattheclassical.orthogonal andtheobliqueprojectionmethodshavepotentialdangerin findingeigenvectorsofthematrixinvolved.Thatis,Ritz vectorsmaynotconvergetothedesiredeigenvectorsevenif thecorrespondingRitzvaluesconverge.Inordertocircumvent theflaw,Jiasuggestedtousecertainrefinedvectorsto approximatethedesiredeigenvectors.Theresultingrefined projectionmethodsavoidtheabovedangerandconvergefaster. Itisimportanttocombinetherefinedstrategywithmany other techniques or methods.More robust and reliable algorithmscanbederivedfromthecombination.Thisthesis consistsofthefollowingparts 1.Thequasi-refinedbiorthog