《43多项式方法求特征值问题》由会员分享，可在线阅读，更多相关《43多项式方法求特征值问题（7页珍藏版）》请在金锄头文库上搜索。
1、43多项式方法求特征值问题4.3.1 F-L方法求多项式系数 我们知道，求n阶方阵A的特征值就是求代数方程 （4.3.1）的根。称为A的特征多项式。上式展开为 （4.3.2）其中为多项式的系数。 从理论上讲，求A的特征值可分为两步： 第一步 直接展开行列式|求出多项式； 第二步 求代数方程的根，即特征值。 对于低阶矩阵，这种方法是可行的。但对于高阶矩阵，计算量则很大，这种方法是不适用的。这里我们介绍用F-L（Faddeev-Leverrier）方法求特征方程（4.3.2）中多项式的系数。由于代数方程求根问题在第2章中已经介绍，所以本节中解决特征值问题的关键是确定矩阵A的特征多项式，所以称这种方
2、法为多项式方法求特征值问题。 记矩阵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的特征向量确定以后,将这些特征值逐个代入齐次
3、线性程组()x=0中,由于系数矩阵的秩小于矩阵的阶数n,因此虽然有n个方程n个未知数,但实际上是解有n个未知数的相互独立的r个方程(rn). 当矩阵A的所有特征值互不相同时,这样的问题中要解的齐次方程组中有n-1个独立方程,其中含有n个特征向量分量,因此特征向量分量中至少有一个需要任意假设其值,才能求出其他特征分量. 在计算机中解这样的齐次线性程组,可用高斯-若当消去法,以便把一组n个方程简化为等价的一组n-1个方程的方程组.然而,用高斯-若当消去法简化一个齐次线性程组时,方程之间不都是独立的,在消去过程中系数为零的情况较多.必需交换方程中未知数的次序,以避免主元素位置上为零的情况.因此,为了
4、提高精度和避免零元素的可能性,我们总是用主元素措施把绝对值最大的系数放于主元素位置. 例如,假设矩阵A为 其特征方程为 =0展开后为 故特征值分别为 下面求特征向量，将代入方程组中，得 （4.3.8）以-5为主元素,交换上式第一与第二个方程得 (4.3.9)用高斯-若当消去法消去-5所在列中的,并把主元素所在行调到最后,得 (4.3.10)再以16/5为主元素,消去它所在列中的,并把主元素所在的行调到最后,得 (4.3.11)这就是用高斯若当消去法实现把一组三个方程简化为等价的一组两个独立方程的情形.因为这个等价的方程组包含两个独立的方程,而有三个未知数，所以只要假定其中一个值,则其它两个值就
5、可以通过两个独立方程解出.比如,令,则得到矩阵A的对应于的一个特征向量为 对另外两个特征值的对应特征向量求法与上述对的推导过程相同. 计算机中实现求解这样的齐次线性方程组的消去步骤是,用第3章讨论过的高斯-若当消去法的公式,方程组(4.3.9)的系数矩阵经过第一次消去后的矩阵B为 (4.3.12)以矩阵为方程组(4.3.10)的系数矩阵,其中省略了有0和1元素的第一列. 在进行第二次消元之前,要应用完全主元素措施对前两行进行最大主元素选择,然后再进行必要的行或列交换.每完成一次消元过程,总省略只有0和1元素的第一列,并且计算机仅寻找矩阵的前n-k行中的最大主元素,其中k是消元过程应用的次数.对
6、(4.3.12)式再进行一次消元过程,则得到列矩阵 (4.3.13)此矩阵是对应于方程组(4.3.11)的系数矩阵,不过省略了含0和1元素的前两列.一般来说,最后矩阵列的数目等于矩阵的阶数和秩的差值. 由于方程组(4.3.8)有三个未知数,两个独立方程,所以计算机必须任意给定一个未知数的值,以便可以从其他两个独立方程中解出另外两个未知数.为方便，在计算机决定特征向量时,要恰当地设定任意选取的未知数的值.例如,令,由方程组(4.3.11)知道,其他两个分量的值正好能从含的非零系数项得出.为此,从计算机所存储的最终矩阵中,令最上面的0元素为-1,并把它顺次调到最下面第三行的位置上,就得到所求的特征
7、向量. 在工程问题中,从特征方程所求出的特征值,少数情形也有相同的.一般地,当一个特征方程有k重根时,矩阵的秩可能比其阶数少1,或2,或3,或k,当然对应于的线性无关的特征向量的个数也就是1,或2,或3,或k,下面通过一个特征值对应两个线性无关特征向量的例子进一步说明计算机求特征向量的方法. 设矩阵A为 其特征方程为 展开后得 所以特征值为 为了决定的特征向量,将代入方程组()x=0,得 (4.3.14)应用一次高斯-若当消去法,得 (4.3.15)写成矩阵形式,(4.3.15)式的系数矩阵为 (4.3.16)因为方程组(4.3.15)的系数矩阵的秩为1,它比矩阵阶数少2,因此对应于有两个线性
8、无关的特征向量,必须给两个未知数任意规定值,才能确定这两个线性无关的特征向量,由（4.3.15）式可看出,一般总是选择求一个特征向量;选择求另一个特征向量;这样有两个线性无关的特征向量 , 计算机中求两个线性无关的特征向量的办法是,在(4.3.16)式的B中,把第一列中第一个0元素用-1代替,第二列中第二个0元素也用-1代替,然后把第一、第二行顺次调到最下面一行的位置上，第三行自然就成了第一行，如此调换后矩阵的第一列和第二列就是所求的两个线性无关的特征向量。对应于的全部特征向量为 其中与是任意常数，且不同时为零。为了说明列交换的必要性，避免主元素为零，再举一个例子，设矩阵A为 其特征方程为 特
9、征值为 对应于的特征向量可由解下列方程组而求得 (4.3.17)用一次高斯-若当消去法，得 (4.3.18)若不进行列交换，则下一个消元过程只能在第一行的第二个元素与第二行的第二个元素中找最大主元素，而它们都是零，我们不得不对(4.3.17)式进行列交换，即交换未知数之间的次序，之后再进行消去过程.对(4.3.17)式进行列交换，即把绝对值最大系数放在主元素位置，显然是第一列与第三列的交换，交换后成为 (4.3.19)其中未知数列矩阵中与也进行了交换，这样才能保证(4.3.17)式与(4.3.19)式等价，对(4.3.19)式进行一次高斯-若当消去法，得 (4.3.20)再进行一次消去过程，得 (4.3.21)在计算机中计算，剩下一个最终的列矩阵 (4.3.22)将(4.3.22)式中的列矩阵B中第一个0元素用-1代替，并随即调到最下面一行，便得到 （4.3.23）这就是对应于方程组(4.3.19)的解，在计算机程序中应把原来进行列交换的列号次序记住，重新把(4.3.23)式中各分量排列一下，即交换第一行和第三行的元素，就得到对应于的特征向量 对应于的全部的特征向量为 k 其中k为不等于零的任意常数.
