43多项式方法求特征值问题.docVIP

4．3多项式方法求特征值问题 4.3.1 F-L方法求多项式系数 我们知道，求n阶方阵A的特征值就是求代数方程 （4.3.1） 的根。称为A的特征多项式。上式展开为 （4.3.2） 其中为多项式的系数。 从理论上讲，求A的特征值可分为两步： 第一步 直接展开行列式||求出多项式； 第二步 求代数方程的根，即特征值。 对于低阶矩阵，这种方法是可行的。但对于高阶矩阵，计算量则很大，这种方法是不适用的。这里我们介绍用F-L（Faddeev-Leverrier）方法求特征方程（4.3.2）中多项式的系数。由于代数方程求根问题在第2章中已经介绍，所以本节中解决特征值问题的关键是确定矩阵A的特征多项式，所以称这种方法为多项式方法求特征值问题。 记矩阵A=的对角线元素之和为 （4.3.3） 利用递归的概念定义以下n个矩阵 （4.3.4） 可以证明,(4.3.4)式中即是所求A的特征多项式的各系数。用（4.3.4）式求矩阵的特征多项式系数的方法称为F-L方法。相应特征方程为： （4.3.5） 而且可证矩阵A的逆矩阵可表示为 (4.3.6) 例1 求矩阵 的特征值与. 解 用F-L方法求得 所以A的特征方程为 此方程的根,即特征值为 从例1中的计算结果可知Faddeev曾经证明: 对n阶矩阵A,按(4.3.4)式计算出的总有 (4.3.7) 4.3.2 特征向量求法 当矩阵A的特征向量确定以后,将这些特征值逐个代入齐次线性程组()x=0中,由于系数矩阵的秩小于矩阵的阶数n,因此虽然有n个方程n个未知数,但实际上是解有n个未知数的相互独立的r个方程(rn). 当矩阵A的所有特征值互不相同时,这样的问题中要解的齐次方程组中有n-1个独立方程,其中含有n个特征向量分量,因此特征向量分量中至少有一个需要任意假设其值,才能求出其他特征分量. 在计算机中解这样的齐次线性程组,可用高斯-若当消去法,以便把一组n个方程简化为等价的一组n-1个方程的方程组.然而,用高斯-若当消去法简化一个齐次线性程组时,方程之间不都是独立的,在消去过程中系数为零的情况较多.必需交换方程中未知数的次序,以避免主元素位置上为零的情况.因此,为了提高精度和避免零元素的可能性,我们总是用主元素措施把绝对值最大的系数放于主元素位置. 例如,假设矩阵A为 其特征方程为 =0 展开后为 故特征值分别为 下面求特征向量，将代入方程组中，得 （4.3.8） 以-5为主元素,交换上式第一与第二个方程得 (4.3.9) 用高斯-若当消去法消去-5所在列中的,并把主元素所在行调到最后,得 (4.3.10) 再以16/5为主元素,消去它所在列中的,并把主元素所在的行调到最后,得 (4.3.11) 这就是用高斯若当消去法实现把一组三个方程简化为等价的一组两个独立方程的情形.因为这个等价的方程组包含两个独立的方程,而有三个未知数，所以只要假定其中一个值,则其它两个值就可以通过两个独立方程解出.比如,令,则得到矩阵A的对应于的一个特征向量为 对另外两个特征值的对应特征向量求法与上述对的推导过程相同. 计算机中实现求解这样的齐次线性方程组的消去步骤是,用第3章讨论过的高斯-若当消去法的公式,方程组(4.3.9)的系数矩阵经过第一次消去后的矩阵B为 (4.3.12) 以矩阵为方程组(4.