《矩阵特征特征向量计算》由会员分享,可在线阅读,更多相关《矩阵特征特征向量计算(17页珍藏版)》请在金锄头文库上搜索。
1、第三章 矩阵特征与特征向量的计算3.1 引言在科学技术的应用领域中,许多问题都归为求解一个特征系统。如动力学系统和结构系统中的振动问题,求系统的频率与振型;物理学中的某些临界值的确定等等。设A为n阶方阵,若,有数l使Ax= lx (5.1)则称l为A的特征值,x为相应于l的特征向量。因此,特征问题的求解包括两方面:1求特征值l,满足(5.2)2求特征向量,满足齐方程组(5.3)称j(l)为A的特征多项式,它是关于l的n次代数方程。关于矩阵的特征值,有下列代数理论,定义1 设矩阵A, BR nn,若有可逆阵P,使则称A与B相似。定理1 若矩阵A, BR nn且相似,则(1)A与B的特征值完全相同
2、;(2)若x是B的特征向量,则Px便为A的特征向量。定理2 设AR nn具有完全的特征向量系,即存在n个线性无关的特征向量构成Rn的一组基底,则经相似变换可化A为对角阵,即有可逆阵P,使其中li为A的特征值,P的各列为相应于li的特征向量。定理3 AR nn,l1, , ln为A的特征值,则(1)A的迹数等于特征值之积,即(2)A的行列式值等于全体特征值之积,即定理4 设AR nn为对称矩阵,其特征值l1l2ln,则(1)对任AR n,x0,(2)(3)定理5 (Gerschgorin圆盘定理) 设AR nn,则(1)A的每一个特征值必属于下述某个圆盘之中,(5.4)(5.4)式表示以aii为
3、中心,以半径为的复平面上的n个圆盘。(2)如果矩阵A的m个圆盘组成的并集S(连通的)与其余n m个圆盘不连接,则S内恰包含m个A的特征值。定理4及定理5给出了矩阵特征值的估计方法及界。例1 设有估计A的特征值的范围。解 由圆盘定理,A的3个圆盘为图5.1D1: D2: D3: 见图5.1。D2为弧立圆盘且包含A的一个实特征值l1(因为虚根成对出现的原理),则3l15。而l2,l3D1D2,则,即3.2 乘幂法与反幂法在实际工程应用中,如大型结构的振动系统中,往往要计算振动系统的最低频率(或前几个最低频率)及相应的振型,相应的数学问题便为求解矩阵的按模最大或前几个按模最大特征值及相应的特征向量问
4、题,或称为求主特征值问题。3.2.1 乘幂法乘幂法是用于求大型稀疏矩阵的主特征值的迭代方法,其特点是公式简单,易于上机实现。乘幂法的计算公式为:设AR nn,取初始向量x(0)R n,令x(1) = Ax(0),x(2) = Ax(1),一般有(5.5)形成迭代向量序列x(k)。由递推公式(5.5),有(5.6)这表明x(k)是用A的k次幂左乘x(0)得到的,因此称此方法为乘幂法,(5.5)或(5.6)式称为乘幂公式,x(k)称为迭代序列。下面分析乘幂过程,即讨论当k时,x(k)与矩阵A的主特征值及相应特征向量的关系。设A = (aij)nn有完全的特征向量系,且l1, l2, ln为A的n个
5、特征值,满足v1, v2, vn为相应的特征向量且线性无关,从而构成Rn上的一组基底。对任取初始向量x(0) Rn,可由这组基底展开表示为 (5.7)其中a1, a2, an为展开系数。将x(0)的展开式(5.7)代入乘幂公式(5.6)中,得 (5.8)利用(5.8)式为 (5.9)(1)如果A有唯一的主特征值,即,设l1 0,且由(5.9)式,有其中,由于,故当k充分大时,e k 0,此时(5.10)对i = 1, 2, , n,若(a1v1)i 0,考虑相邻迭代向量的对应分量比值,(5.11)即对i = 1, , n(5.12)这表明主特征值l1可由(5.11)或(5.12)式得到。由于迭
6、代序列x(k),当k充分大时,(5.10)式成立,x(k)与v1只相差一个常数因子,故可取x(k)作为相应于主特征值l1的特征向量的近似值。迭代序列x(k)的收敛速度取决于的大小。(2)如果A的主特征值不唯一,且可分三种情况讨论:a)l1 = l2;b)l1 = - l2;c)情况a)当l1 = l2时,A的主特征值为二重根,根据(5.9)式当k充分小时,由于,j = 3, n,e k 0,则对i = 1, 2, , n,如果,则(主特征值)且x(k)收敛到相应于l1 (=l2)的特征向量的近似值。这种重主特征值的情况,可推广到A的r重主特征值的情况,即当 且时,上述讨论的结论仍然成立。情况b
7、)当l1 = - l2时,A的主特征值为相反数,(5.9)式为当k充分大时,j = 3, 4, , n, e k 0,则(5.13)由于(5.13)式中出现因子(-1) k,则当k变化时,x(k)出现振荡、摆动现象,不收敛,利用(-1) k的特点,连续迭代两步,得从而,对i = 1, 2, , n,若,则(5.14)开方之后,便得到A的以上主特征值l1, l2 = - l1。为计算相应于l1, l2的特征向量,采取组合方式,(5.15)(5.16)可见分别为相应于l1与l2的特征向量。情况c)当时,A的主特征值为共轭复根。因A为实矩阵,于是由有即(v1与v2为互为共轭向量)。设,对任取x(0)
8、 R n,展开式(5.7)可为(5.17)将(5.17)式代入(5.9)式,同理,当k充分大时(5.18)对j = 1, 2, n,设复数表示则(5.18)式的复数表示可为连续迭代,得(5.19)利用三角函数运算性质及l1、l2的复数表示,不难验证。令(5.20)解方程(j = 1, 2, , n)(5.21)求出p, q后,再解出主特征值l1、l2,得 (5.22)同样,采取组合方式求相应于l1、l2的特征向量。由于(5.23)(5.24)则可分别取(5.23)、(5.24)左端的组合表达式作为相应于l1、l2的特征向量的近似值。通过上述分析,有定理6 设A Rnn有完全特征向量系,若l1,
9、 l2, ln为A的n个特征值且满足对任取初始向量x(0) Rn,对乘幂公式确定的迭代序列xk,有下述结论:(1)当时,对i = 1, 2, , n收敛速度取决于的程度,r 1收敛快,r 1收敛慢,且x(k)(当k充分大时)为相应于l1的特征向量的近似值。(2)当时a)若l1 = l2,则主特征值l1及相应特征向量的求法同(1);b)若l1 = -l2,对i = 1, 2, , n收敛速度取决于的程度。向量、分别为主特征值l1、l2相应的特征向量的近似值。c)若,则连续迭代两次,计算出x(k+1),x(k+2),然后对j = 1, 2, , n解方程求出、后,由公式解出主特征值l1、l2。此时
10、收敛速度取决于的程度。向量、分别为相应于l1,l2的特征向量的近似值。从分析乘幂过程可见,乘幂法可用于求矩阵按模最大的一个(或几个)特征值及相应的特征向量,当比值时,收敛速度快,r 1时,收敛速度慢,且计算公式简便,便于上机实现。分析中的假设、,在计算时可不用考虑,如果此条件不满足,则可通过迭代误差自行调整。在用乘幂法求矩阵的主特征值l1及对应的特征向量时,迭代向量的分量可能会出现绝对值非常大的现象,从而造成计算中溢出的可能。为此,需对迭代向量x(k)进行规范化。令max(x)表示向量x分量中绝对值最大者。即如果有某i0,使则max (x) = xi对任取初始向量x(0),记则一般地,若已知x
11、(k),称公式(5.25)为规范化的乘幂法公式或改进乘幂法公式,这里,乘幂迭代序列y(k)的分量绝对值最大者1。类似前面的分析乘幂过程,有定理7 设ARnn具有完全特征向量系,l1, l2, , ln为A的n个特征值,且满足则对任初始向量x(0),由规范化的乘幂法公式(5.25)确定的向量序列y(k),x(k)满足(1)(5.26)(2)y(k)为相应于主特征值l1的特征向量近似值(5.27)例2 用规范化乘幂法计算矩阵A的主特征值及相应特征向量解 A的特征值l1 = 6,l2 = 3,l3 = 2取初始值x(0) = (1, 1, 1)T,用规范化乘幂法公式(5.25)计算其它结果见表5.1
12、(表中的向量均为转置向量)。表5.1kmax(y(k)x(k) = y(k)/max(x(k)x(k+1) = Ay(k)01(1, 1, 1)(10, 8, 1)110(1, 0.8, 0.1)(7.2, 5.4, -0.8)27.2(1, 0.75, -0.111111)(6.5, 4.75, -1.222222)36.57(1, 0.730769, -0.203704)(6.230766, 4.499997, -1.407408)46.230766(1, 0.722222, -0.225880)(6.111108, 4.388886, -1.1451767)56.111108(1, 0.
13、718182, -0.237561)(6.054548, 4.336336, -1.475122)66.054548(1, 0.716216, -0.243639)(6.027024, 4.310808, -1.487278)76.027024(1, 0.715247, -0.246768)(6.013458, 4.298211, -1.483536)86.013458(1, 0.714765, -0.248366)(6.00671, 4.291945, -1.496732)96.00671(1, 0.714525, -0.249177)(6.00335, 4.28825, -1.496354)106.00335