4．3 多项式方法求特征值问题4.3.1 F-L 方法求多项式系数我们知道，求 n 阶方阵 A 的特征值就是求代数方程( ) | A I | 0（ 4.3.1）的根。

( )称为 A 的特征多项式。

上式展开为( )np 1n 1p2 n2 .....p n（ 4.3.2）其中 p 1, p 2 ,...p n 为多项式 ( )的系数。

从理论上讲，求 A 的特征值可分为两步：第一步 直接展开行列式 |AI|求出多项式( ) ；第二步 求代数方程 (x) 0的根，即特征值。

对于低阶矩阵， 这种方法是可行的。

但对于高阶矩阵，计算量则很大，这种方法是不适用的。

这里我们介绍用 F-L （ Faddeev-Leverrier ）方法求特征方程（4.3.2）中多项式 () 的 系数。

由于代数方程求根问题在第 2 章中已经介绍， 所以本节中解决特征值问题的关键是确 定矩阵 A 的特征多项式()，所以称这种方法为多项式方法求特征值问题。

记矩阵 A=(a ij)n n的对角线元素之和为trAa11a22... a nn（4.3.3 ）利用递归的概念定义以下n 个矩阵 B k ( k 1,2,...., n) :p 1 trB 1B 1A,p 21B 2A(B 1p 1I ),trB22p 31trB 3B 3A(B 2 p 2 I ),3...............p k1trB k B k A(B kp k 1 I ),k1................p n1B nA(B n 1p n 1 I ),trB nn（ 4.3.4）可以证明 ,(4.3.4)式中p k, k 1,2,..., n,即是所求 A 的特征多项式( )的各系数。

用（ 4.3.4）式求矩阵的特征多项式系数的方法称为 F-L 方法。

相应特征方程为：( 1)n (np 1n 1p 2n 2..... p n ) 0（ 4.3.5）A 的逆矩阵可表示为而且可证矩阵 A 11( B n 1p n 1I )p n(4.3.6)例1 求矩阵3 2 4A2 0 24 2 3的特征值与 A 1 .解 用 F-L 方法求得3 24 B 1A2 0 24 2 3p 1trB 1 611 2 4 B 2A( B 1 p 1 I )2 8 242 11p 21trB 21528 0 0 B 3A( B 2 p 2 I )0 8 00 08p 31trB 38所以 A 的特征方程为3( 1)3(36215 8) 0此方程的根 ,即特征值为18,2 1,311 1 12 4 2 A 11 ( B2 p 2 I )1 7 1p 34 8 4 1 1 1242从例 1 中的计算结果可知 B 3p 3I .Faddeev 曾经证明 : 对 n 阶矩阵 A, 按 (4.3.4)式计算出的Bn 总有B np n I(4.3.7)4.3.2 特征向量求法当矩阵 A 的特征向量确定以后 ,将这些特征值逐个代入齐次线性程组(AI)x=0 中,由于系数矩阵 AI 的秩小于矩阵 A I 的阶数 n,因此虽然有 n 个方程 n 个未知数 ,但实际上是解有 n 个未知数的相互独立的 r 个方程 (r<n). 当矩阵 A 的所有特征值互不相同时 ,这样的问题中要解的齐次方程组中有n-1 个独立方程 ,其中含有 n 个特征向量分量 ,因此特征向量分 量中至少有一个需要任意假设其值 ,才能求出其他特征分量 . 在计算机中解这样的齐次线性程组,可用高斯 -若当消去法 ,以便把一组 n 个方程简化为 等价的一组 n-1 个方程的方程组 .然而 ,用高斯 -若当消去法简化一个齐次线性程组时 ,方程之间不都是独立的 ,在消去过程中系数为零的情况较多 .必需交换方程中未知数的次序,以避免 主元素位置上为零的情况 .因此 ,为了提高精度和避免零元素的可能性 ,我们总是用主元素措 施把绝对值最大的系数放于主元素位置 .例如 ,假设矩阵 A 为4 2 2 A5 3 2 2 41其特征方程为42 2 53 2241=0展开后为(1)(2)(5)0 故特征值分别为11, 22,35下面求特征向量，将1 代入方程组( AI )x 0 中，得3x 12 x 2 2x 35x 1 2x 22x 32x 1 4x 2 0x 3（4.3.8 ）以-5 为主元素 ,交换上式第一与第二个方程得5x 1 2x 2 2x 3 0 3x 1 2 x 2 2x 32x 1 4x 2 0x 3(4.3.9)用高斯 -若当消去法消去 -5 所在列中的x1,并把主元素所在行调到最后,得0x 116x 2 4x 3 05 50x 116x 24 x 3 055x 12x 2 2x 3 05 5(4.3.10)再以 16/5 为主元素 ,消去它所在列中的x 2,并把主元素所在的行调到最后,得0x 1 0x 2 0 x 30 x 1 0x 21 x 320x 1 x 21x 34(4.3.11)这就是用高斯若当消去法实现把一组三个方程简化为等价的一组两个独立方程的情形这个等价的方程组包含两个独立的方程 ,而有三个未知数，所以只要假定其中一个值 两个值就可以通过两个独立方程解出 .比如 ,令x 31,则得到矩阵 A 的对应于 1特征向量为.因为 ,则其它1的一个1 2 1 4 1对另外两个特征值的对应特征向量求法与上述对11的推导过程相同.计算机中实现求解这样的齐次线性方程组的消去步骤是,用第 3 章讨论过的高斯 -若当消去法的公式 ,方程组 (4.3.9)的系数矩阵经过第一次消去后的矩阵 B 为16 45 516 4B552 25 5 (4.3.12)以矩阵为方程组 (4.3.10) 的系数矩阵 ,其中省略了有0 和 1 元素的第一列 .在进行第二次消元之前,要应用完全主元素措施对前两行进行最大主元素选择,然后再进行必要的行或列交换 .每完成一次消元过程 ,总省略只有 0 和 1 元素的第一列 ,并且计算机仅寻找矩阵的前 n-k 行中的最大主元素 ,其中 k 是消元过程应用的次数.对 (4.3.12)式再进行一次消元过程 ,则得到列矩阵B1 1214 (4.3.13)此矩阵是对应于方程组(4.3.11)的系数矩阵 ,不过省略了含 0 和 1 元素的前两列 .一般来说 ,最后矩阵列的数目等于矩阵 A I 的阶数和秩的差值.由于方程组 (4.3.8) 有三个未知数 ,两个独立方程 ,所以计算机必须任意给定一个未知数的值,以便可以从其他两个独立方程中解出另外两个未知数.为方便，在计算机决定特征向量时 ,要恰当地设定任意选取的未知数的值.例如 ,令x31, 由方程组 (4.3.11)知道 ,其他两个分量的值正好能从含x3的非零系数项得出.为此 ,从计算机所存储的最终矩阵中,令B1 最上面的 0( 1 ,1, 1)T元素为 -1,并把它顺次调到最下面第三行的位置上,就得到所求的特征向量 2 4 .在工程问题中 ,从特征方程所求出的特征值,少数情形也有相同的 .一般地 ,当一个特征方程有 k 重根时 ,矩阵A I 的秩可能比其阶数少1,或 2,或 3, ,或 k,当然对应于的线性无关的特征向量的个数也就是1,或 2,或 3, ,或 k,下面通过一个特征值对应两个线性无关特征向量的例子进一步说明计算机求特征向量的方法.设矩阵A为3 2 4A 2 0 24 2 3其特征方程为3 24 22 0423展开后得(1)2( 8)0 所以特征值为121,38为了决定1的特征向量 ,将1 代入方程组 ( AI )x=0, 得4 2 4 x 12 1 2 x 2 042 4 x 3(4.3.14)应用一次高斯 -若当消去法 ,得0 0 0 x 1 00 x 21 1 /2 1 x 3(4.3.15)写成矩阵形式 ,(4.3.15) 式的系数矩阵为0 0 B1/2 1(4.3.16)1有两个线性无 因为方程组 (4.3.15)的系数矩阵的秩为 1,它比矩阵阶数少 2,因此对应于关的特征向量 ,必须给两个未知数任意规定值 ,才能确定这两个线性无关的特征向量 , 由（ 4.3.15）式可看出 ,一般总是选择 x 21, x 3 0求一个特征向量;选择x 20, x31 求另一个特征向量 ;这样有两个线性无关的特征向量1/ 2 1 10 01,,在 (4.3.16) 式的 B 中 ,把第一列中第一个 0计算机中求两个线性无关的特征向量的办法是 元素用 -1 代替 ,第二列中第二个 0 元素也用 -1 代替 ,然后把第一、 第二行顺次调到最下面一行 的位置上， 第三行自然就成了第一行， 如此调换后矩阵的第一列和第二列就是所求的两个线性无关的特征向量。

对应于1 的全部特征向量为1/ 2 1 k 11 k2 01其中k1 与k2 是任意常数，且不同时为零。

为了说明列交换的必要性，避免主元素为零，再举一个例子，设矩阵A 为2 8 12A1 4 4 01其特征方程为( 2) ( 1) 0特征值为1 2,2 0,3 1对应于 2 的特征向量可由解下列方程组而求得4 8 12 x11 2 4 x2 00 0 1 x3(4.3.17)用一次高斯 -若当消去法，得0 0 1 x10 0 1 x2 01 2 3 x3(4.3.18)若不进行列交换，则下一个消元过程只能在第一行的第二个元素与第二行的第二个元素中找最大主元素，而它们都是零，我们不得不对 (4.3.17) 式进行列交换，即交换未知数之间的次序，之后再进行消去过程 .对 (4.3.17)式进行列交换，即把绝对值最大系数放在主元素位置，显然是第一列与第三列的交换，交换后成为1284x3421x20100x1其中未知数列矩阵中x1 与x3 也进行了交换，这样才能保证(4.3.19) 式进行一次高斯 -若当消去法，得0 2 / 3 1/ 3 x30 2 / 3 1 / 3 x2 01 2 / 3 1 / 3 x1再进行一次消去过程，得(4.3.19)(4.3.17) 式与 (4.3.19) 式等价，对(4.3.20)0 0 0 x31 0 0 x2 00 1 1/ 2 x1(4.3.21)在计算机中计算，剩下一个最终的列矩阵B 01/ 2(4.3.22)将(4.3.22) 式中的列矩阵 B 中第一个0 元素用 -1 代替，并随即调到最下面一行，便得到1/ 21（ 4.3.23）这就是对应于方程组(4.3.19)的解，在计算机程序中应把原来进行列交换的列号次序记住，重新把 (4.3.23)式中各分量排列一下，即交换第一行和第三行的元素，就得到对应于 2 的特征向量11/ 2对应于的全部的特征向量为11/ 2k其中 k 为不等于零的任意常数.。