用矩阵符号函数和Padé逼近求解不变子空间
作者单位:浙江大学
学位级别:硕士
导师姓名:张振跃
授予年度:2007年
学科分类:07[理学] 070102[理学-计算数学] 0701[理学-数学]
主 题:矩阵符号函数 不变子空间 牛顿迭代 有理迭代 Padé逼近 KryloV子空间
摘 要:本文主要研究矩阵特征值中的不变子空间的计算问题。关于特征值的计算问题,有很多经典的算法。然而这些经典算法用于计算不变子空间时,常会碰到很大困难,有时甚至无法采用。产生于二十世纪七十年代早期的在控制论方面得到广泛应用的矩阵符号函数,很快便被用于求解不变子空间。本文研究用矩阵符号函数来求解不变子空间的一些迭代方法。 本文首先介绍矩阵符号函数的代数形式和几何形式的定义及其一些性质。然后回顾计算矩阵符号函数的牛顿迭代和有理迭代方法并分析了牛顿迭代法的收敛速度和加速收敛问题,有理迭代算法的收敛速度和Hermite矩阵情形的改进问题。其次,作为本文的核心部分,主要研究用二阶Padé逼近来计算矩阵符号函数。与其它计算矩阵符号函数的迭代法相比较,二阶Padé迭代方法具有更快的收敛性,具有五次收敛速度,在此算法中,计算量稍有增加。Krylov子空间方法是求解大型线性方程组和大型矩阵的特征值问题的一种很有效的方法。针对大型矩阵求逆比较困难的问题,本文考虑如何用Krylov子空间来避免这一点。文章的最后给出了数值实验,实验结果表明Padé二阶逼近确实可以用于解决不变子空间问题,与本文前面所提到的几种方法相比,有更快的收敛速度。同时对于高阶Padé逼近的收敛特性,也给出了数值例子来验证。