用矩阵符号函数和Pad#233;逼近求解稳定子空间.pdfVIP

摘 要 本文主要研究矩阵特征值中的不变子空间的计算问题。关于特征值的计算 问题，有很多经典的算法。然而这些经典算法用于计算不变子空间时，常会碰 到很大困难，有时甚至无法采用。产生于二十世纪七十年代早期的在控制论方 面得到广泛应用的矩阵符号函数，很快便被用于求解不变子空间。本文研究用 矩阵符号函数来求解不变子空间的～些迭代方法。 本文首先介绍矩阵符号函数的代数形式和几何形式的定义及其一些性质。 然后回顾计算矩阵符号函数的牛顿迭代和有理迭代方法并分析了牛顿迭代法的 收敛速度和加速收敛问题，有理迭代算法的收敛速度和Hermite矩阵情形的改 进问题。其次，作为本文的核心部分，主要研究用二阶Pad6逼近来计算矩阵符 号函数。与其它计算矩阵符号函数的迭代法相比较，二阶P甜6迭代方法具有更 快的收敛性，具有五次收敛速度，在此算法中，计算量稍有增加。Krylov子空间 方法是求解大型线性方程组和大型矩阵的特征值问题的一种很有效的方法。针 对大型矩阵求逆比较困难的问题，本文考虑如何用Krylov子空间来避免这一点。 文章的最后给出了数值实验，实验结果表明PadC-阶逼近确实可以用于解决不 变子空间问题，与本文前面所提到的几种方法相比，有更快的收敛速度。同时 对于高阶Pad6逼近的收敛特性，也给出了数值例子来验证。 关键词：矩阵符号函数，不变子空间，牛顿迭代，有理迭代，Pad6逼 近，Krylov子空间。 Abstract This dealswiththe ofmatrix the typical paper problem eigenvalue：compute invariant of matrixassociatedwithits subspacea舀Ven preindicatedeigenvalues． Thereare classical for the many algorithmssolvingeigenvalueproblem．when theseclassical to invariant issomuch using algorithmscompute subspace，there sometimes．Thematrixfunctionintroducedinthe difficulty,unfeasible sign early is usedincontrol canbeusedto invaxiant 1970s widely system，and compute studiesthe matrixfunction papermainly algorithmsusing sign subspace．The to invariant compute subspace． this introduce definitionofthe In the and paper，we algebraicgeometrical matrix withitscharacteristicsatfirst．We re- signfunction，together briefly Newtoniterationandtherationalite