特征值法

特征值法

对元素为实数或复数的n×n矩阵A，求数λ和n维非零向量x使A x=λx，这样的问题称为代数特征值问题，也称矩阵特征值问题，λ和x分别称为矩阵A的特征值和特征向量。代

数特征值问题的数值解法是计算数学的主要研究课题之一，它常出现于动力系统和结构系统的振动问题中。在常微分方程和偏微分方程的数值分析中确定连续问题的近似特征系，若用有限元方法或有限差分方法求解，最终也化成代数特征值问题。此外，其他数值方法的理论分析，例如确定某些迭代法的收敛性条件和初值问题差分法的稳定性条件，以及讨论计算过程对舍入误差的稳定性问题等都与特征值问题有密切联系。求解矩阵特征值问题已有不少有效而可靠的方法。

矩阵A的特征值是它的特征多项式P n(λ)det(λI-A)的根，其中I为单位矩阵。但阶数超过4的多项式一般不能用有限次运算求出根，因而特征值问题的计算方法本质上是迭代性质的，基本上可分为向量迭代法和变换方法两类。

向量迭代法是不破坏原矩阵A，而利用A对某些向量作运算产生迭代向量的求解方法，多

用来求矩阵的部分极端特征值和相应的特征向量，特别适用于高阶稀疏矩阵。乘幂法、反幂法都属此类，隆措什方法也常作为迭代法使用。

变换方法是利用一系列特殊的变换矩阵(初等下三角阵、豪斯霍尔德矩阵、平面旋转矩阵等)，从矩阵A出发逐次进行相似变换，使变换后的矩阵序列趋于容易求得特征值的特殊形式的矩阵(对角阵、三角阵、拟三角阵等)；多用于求解全部特征值问题，其优点是收敛速度快，计算结果可靠，但由于原矩阵A被破坏，当A是稀疏矩阵时，在计算过程中很难保持它的稀疏性，因而大多数变换方法只适于求解中小规模稠密矩阵的全部特征值问题。雅可比方法、吉文斯－豪斯霍尔德方法以及LR方法、QR方法等都属此类。

乘幂法计算矩阵的按模最大的特征值及对应特征向量的一种向量迭代法。设A为具有线性初等因子的矩阵，它的n个线性无关的特征向量是u i(i=1，2，…，n)，特征值排列

次序满足是一个n维非零向量，于是

若λ1>λ2，则当α1≠0，且k足够大时，A k z0除相差一个纯量因子外趋于λ1所对应的特征向量，这就是乘幂法的基本思想。实际计算中

最简单情况的乘幂法的迭代格式是：取初始向量z0，计算式

中指中绝对值最大的一个分量，这时

一般情形下，由于计算格式依赖于A的特征值的分布情况，实际使用很不方便，只是在求

阶数非常高的矩阵的个别特征值和相应的特征向量时偶尔使用。然而乘幂法的基本思想是重要的，由它可导出许多实用的计算方法，例如反幂法和子空间迭代法，它也是其他一些有效

计算方法(例如LR方法，QR方法)的理论基础。

子空间迭代法又称同时迭代法，乘幂法的直接推广，能同时求出模最大的一些特征值和相应的特征向量。它与乘幂法的差别主要在二个方面：①同时用p(p>1)个正交规范向量进行类似于乘幂法的迭代，将迭代向量看作一个p维子空间的正交规范基，每次迭代后得到一个新的子空间；②迭代过程中，在p维子空间上应用里茨原理进行加速。这个方法更便于使用计算机自动计算，而且加快了收敛速度，是大型稀疏矩阵特征值问题的有效解法。反幂法又称反迭代，其原理是:设矩阵A非奇异，为求A的模最小的特征值和相应的

特征向量而对A-1使用乘幂法。A-1的特征值次序是取z0为初始向

量，迭代格式为在一定条件下

反幂法的每次迭代要解一个线代数方程组。这种原始的反幂法在实际计算中很少应用，实际使用的反幂法总是带原点位移的，且位移常取为已求得的近似特征值，而用反幂法求其对应的特征向量。设i是A的某特征值λi的近似

的反幂法的迭代格式是：取初始向量z0，计算

值，带原点位移

当时，按

方向收敛于特征向量u i，收敛商是i越精确，收敛就越快，但

就越接近奇异矩阵，因而在迭代过程中需要求解非常病态的方程组。然而已经证明，这个病态性质对于反幂法的按方向收敛于特征向量是无关紧要的，而且当i相当精确时，一般只要一、二次反迭代就可求得相当精确的近似特征向量。反幂法是由已知近似特征值求对应近似特征向量的有效方法。

隆措什方法用相似变换将矩阵A约化为三对角矩阵的一种方法，其特点是不破坏原矩阵A。目前能实际使用的是A对称时的对称隆措什方法：取初始向量v1，‖v1‖=1，递推地计算式中αi和βi由使｛v i｝为正交规范化的条件确定，即当精确计算时v n+1=0，上述过程可写成矩阵形式其中T n是对称三

对角矩阵它的特征值可由二分法求得。由于舍

入误差的影响，隆措什方法所产生的隆措什向量很快失去正交性，因而这方法长期来被认为不稳定，很少用于实际计算。近年来对隆措什方法作了深入研究，进行了大量实际计算和细

致的误差分析，并且观察到如下的所谓隆措什现象：将m步递推关系写为

其中是m维行向量。对足够大的m，矩阵T m 的特征值包括A的所有相异特征值。注意，当出现舍入误差时，可能m>n。对高阶矩阵A，对不大的m(m<<<

雅可比方法对称矩阵可以通过正交相似变换化为对角阵，其对角元是原矩阵的全部特征值。雅可比方法就是通过一系列特殊的正交相似变换雅可比旋转，使对称矩阵近似对角化从而求得特征值和特征向量的方法。记A0=A，作正交相似序列

其中R k是(p，q)超平面的雅可比旋转矩阵，即

p、q的选取应使中

非对角元绝对值最大者。A k和A k-1仅在第p行(列)和第行(列)不同，它们之间的关系为

选

可使当k→时，A K趋于对角阵，这就是经典雅可比方法。此方法在第k次变换前要搜索A k-1的非对角元中绝对值最大者以确定雅可比旋转矩阵的旋转平面，但这很费时间。为避免这种消耗，可改用循环雅可比方法:在n(n-1)/2次的相继雅可比旋转(称为一次扫描)中，每个非对角元，不管按什么次序，恰好消去一次。其中最方便的是特殊循环雅可比方法，即每次扫描都按(1，2)，(1，3)，…，(1，n)；(2，3)，…，(2，n)；…；(n-1，n)的次序进行。实际计算中最广泛使用

的是特殊循环雅可比方法的一种变形，称为阈雅可比方法：确定一个正数作阈值，在特殊循环的一次扫描中只对绝对值超过阈值的非对角元所在的平面进行旋转变换，反复扫描，当所有非对角元的绝对值都不超过阈值时减小阈值，再按新的阈值进行扫描，如此继续下去，直至阈值充分小而达到近似对角化。雅可比方法的优点是：能在求特征值的同时求得相当精确的近似正交规范特征向量系；缺点是：与其他变换方法相比，收敛速度较慢。它适用于求解

低阶稠密矩阵的全部特征值问题，对近似对角(即非对角元素较小)的对称矩阵特别有效，常应用于子空间迭代法的里茨加速过程中。