计算矩阵符号函数和不变子空间的Pad#233;切近亲近方法.pdfVIP

摘 要 本文主要研究矩阵特征值问题中的一个典型问题：不变子空间的计算。历 史上出现过的许多求解特征值问题的经典算法在求解此问题时，都碰到了很大 的困难，有时甚至不能采用。近年来不变子空间问题引起了许多学者的关注，并 产生了许多优秀的算法，矩阵符号函数技术就是其中的一种。 在本文中，我们首先介绍了矩阵符号函数的定义及性质，回顾了计算矩阵 符号函数的经典算法，然后分析了历史上的一些改进算法，在此基础上进一步 研究矩阵符号函数数值计算方法，及其在计算指定特征值范围的不变子空间 上的应用。首先，我们考虑传统的计算符号函数的迭代法的一种有理形式的改 进算法。这一改进形式并不增加计算量，但可以很好地改进算法的稳定性和计 算精度。同时，对于Hermite情形，这一有理方法适用于借助于矩阵QR分解来避 免算法中的求逆运算。其次，作为本论文的最主要部分，我们首次提出了利用 Pad6逼近计算矩阵符号函数的数值计算方法。一阶Pad6逼近被用作为基本的 迭代函数。与传统的计算符号函数的迭代法相比，一阶Pad6迭代方法具有明显 的优点：在基本不增加计算量的同时，它有更快的收敛性。我们分析并证明了 其三阶收敛速度。进而，我们也初步讨论了用计算量较少的高阶Pad6逼近，来 近似地计算不变子空间