特征值估计与表示.ppt

《特征值估计与表示.ppt》由会员分享，可在线阅读，更多相关《特征值估计与表示.ppt（26页珍藏版）》请在三一办公上搜索。

1、第五章 特征值估计及极性,知识要点：特征值的估计；广义特征值问题；特征值的极小极大原理;特征值和奇异值的扰动；广义特征值分析的应用。,5.1特征值的估计,一、特征值的界 1.定理5.1：设A=(aij)Rnn，若表示A的任一特征值，则其中。,2.推论 实对称矩阵的特征值都是实数。,3.引理1：设BCnn，yCn为单位列向量，则,证明：设B=(bij)nn，则,4.定理5.2：设ACnn，则A的任一特征值 满足(1)|A|m(2)|Re()|0.5|A+AH|m(3)|Im()|0.5|AAH|m。证明：设A属于的单位特征向量为y，则有Ay=y，即 yHAy=yHy=，因此,由引理，于是有,例：

2、估计矩阵 特征值的上界。,5.推论 Hermite矩阵的特征值都是实数，反Hermite矩阵的特征值为零或纯虚数,解：由定理5.2，对A特征值，有：|2，|Re()|2，|Im()|1.3，由定理5.1，知其虚部的另一逼近为：,其特征值为：,6.定义5.1 设A=(aij)Cnn,记Rr=sr|ars|,r=1,n，如果|arr|Rr(r=1,2,n)，则称矩阵A按行严格对角占优；如果|arr|Rr(r=1，n)，且有lron，使得|aroro|Rro成立，则称矩阵A按行(弱)对角占优。7.定义5.2 设ACnn，如果AT按行严格对角占优，则称A按列严格对角占优；如果AT按行(弱)对角占优、则

3、称A按列(弱)对角占优。,二、特征值的包含区域,1.定义5.3 设A=(aij)Cnn，称区域 Gi:|z-aii|Ri 为矩阵A的第i个盖尔圆，其中 Ri=ji|aij|称为盖尔圆Gi的半径(i=l，n)。2.定理5.6 矩阵A=(aij)Cnn的一切特征值都在它的n个盖尔圆的并集之内。证明:设为其特征值，为对应特征向量，且 为其绝对值最大者，则有即,3.定理5.7 由矩阵A的所有盖尔圆组成的连通部分中任取一个，如果它是由k个盖尔圆构成的，则在这个连通部分中有且仅有A的k个特征值(盖尔圆相重时重复计数特征值相同时也重复计数)证明思路：考虑由A的对角线元素构成的矩阵D=diag(a11,a22

4、,ann)，定义矩阵B(u)=(1-u)D+uA则其特征值变化连续依赖于参数u，D的盖尔圆连续变化成为A的盖尔圆。,因此,例：讨论矩阵 的特征值的分布。,解：A的盖尔圆分别为|z-1|0.8和|z|0.5，这两个盖尔圆为连通的，因此包含两个特征值。其特征值为不在盖尔圆|z|0.5内。,考虑满秩对角阵,则矩阵DAD-1与A具有同样的特征值，因此有 若将Ri改作ri=ji(|aij|i/j)，则两个盖尔定理仍然成立，其中i都是正数。,隔离矩阵特征值原则,选取的一般方法是：观察A的n个盖尔圆，欲使第i个盖尔圆Gi的半径变大(或小)些，就取i1(或i1)而取其它正数=1。此时，B=DAD-1的第i个盖

5、尔圆的半径变大(或小)，而B的其余盖尔圆的半径相对变小(或变大)但是，这种隔离矩阵特征值的办法还不能用于任意的具有互异特征值的矩阵比如主对角线上有相同元素的矩阵如果矩阵A按行(列)严格对角占优，则detA0。,例:隔离矩阵A=的特征值,A的3个盖尔圆为G1:|z-20|5.8，G2:|z-10|5，G3:|z-10j|3。G1与G2相交；而G3孤立，其中恰好有A的一个特征值，记作3(见左图)选取D=diag(1，1，2)，则B=DAD-1的三个盖尔圆为G1:|z-20|5.4，G2:|z-10|4.5，G3:|z-10j|6。易见，这是3个孤立的盖尔圆，每个盖尔圆中恰好有B的(也是A的)一个特

6、征值(见右图),定理5.11：设矩阵A=(aij)Cnn的，01，是A的任一个特征值，则存在i使得,|aii|Ri(A)Ri(AT)1-例：讨论矩阵 的特征值的分布。解：R1(A)=0.8,R2(A)=0.5;R1(AT)=0.5,R2(AT)=0.8.取=0.5,则A的特征值满足不等式|1|R1(A)1/2R1(AT)1/2=0.41/2=0.6324|R2(A)1/2R2(AT)1/2=0.41/2=0.6324,5.2 广义特征值问题,定义:称 Ax=Bx 的特征值问题为(对称)矩阵A相对于(对称)矩阵B的广义特征值问题，称数为矩阵A相对于矩阵B的特征值；而与相对应的非零解x称之为属于的

7、特征向量广义特征值由det(A-B)=0的根给出。一、广义特征值问题的等价形式1.等价形式1：B可逆时B-1Ax=x，等价地化为非对称阵B-1A的普通特征值问题。2.等价形式2：B正定时 B=GGT 使得 Sy=y，其中 y=GTx，对称阵 S=G-1AG-T。等价地转化为对称矩阵 S 的普通特征值问题,因此，当 B=GGT 正定时有正交矩阵P，使得,令 Q=G-TP,则有,设A与B为正定对称阵，则A+B仍为正定对称阵，由以上结论，存在可逆矩阵Q，使得,因此有,二、特征向量的共轭性,1.在等价的普通特征值问题 Sy=y中，特征向量系y1,y2,yn是完备的标准正交系。令xj=GTyj，j=1,

8、2,n，则有xiTBxj=xiTGGTxj=(GTxi)T(GTxj)=yiTyj=ij，向量系x1，xn称为按B标准正交化向量系。2.按B标准正交化向量系的性质：性质1 xj0(j=1,2,n)(j=1，n)；性质2 x1，xn线性无关。,5.3 对称矩阵特征值的极性,一、实对称矩阵的Rayleigh商的极性1.定义：设A是n阶实对称矩阵，x Rn称为矩阵A的Rayleigh商2.Rayleigh商的性质：性质1 R(x)是x的连续函数性质2 R(x)是x的零次齐次函数即，对任意的实数0，有R(x)=R(x)=0R(x),性质3 xL(x0)(x00)时，R(x)是一常数性质4 R(x)的最

9、大值和最小值存在，且能够在单位球面S=x|xRn,|x|2=1上达到 证：S是闭集，在S上R(x)=xTAx连续，所以必有x1，x2S，使得 minxSR(x)=R(x1)maxxSR(x)=R(x2)任取0yRn，令y0=y/|y|2，则y0S，根据性质3，有R(y)=R(y0)，从而R(x1)R(y)R(x2)。,实对称矩阵A的特征值(都是实数)按其大小升序排列：12n，对应的标准正交特征向量系设为 P=p1，pn，则有 定理：设A为实对称矩阵，则minxSR(x)=1，maxxSR(x)=n 证：任取xS，则x=Pc,|c|=1,Ax=APc=Pc R(x)=xTAx=cTc1R(x)n

10、，Api=ipiR(pi)=i。,推论1：在S上p1和pn分别是R(x)的一个极小点和极大点，即R(p1)=1，R(pn)=n 推论2 若1=k(1kn)则在|x|2=l上R(x)的所有极小点为p1，pk，|2=1。定理：设xL(pr，ps)，1rsn，则有minxR(x)=r，maxxR(x)=s Courant-Fischer定理:设实对称矩阵A的特征值按升序排列，则A的第k个特征值其中Vk是Rn的任意个k维子空间，1kn。,Courant-Fischer定理的证明,构造Rn的子空间Wk=L(pk,pn)，则dimWk=n-k+1由于Vk+Wk Rn，所以ndim(Vk+Wk)=dim(V

11、k)+dim(Wk)-dim(VkWk)=n+1-dim(VkWk)dim(VkWk)1 故存在x0=pk,pnVkWk，|2=1满足|x0|2=1使得xTAx=T k，即maxxTAx|xVk,|x|2=1k 根据Vk的任意性，可得:,令Vk=L(p1,pk)，取x=pk,pnVk满足|x|2=l，则有xTAxk，即maxxTAx|xVk,|x|2=1k 于是,二、广义特征值的极小极大原理,定义：设A，B为n阶实对称矩阵，且B正定，xRn称R(x)=(xTAx)/(xTBx)，x0为矩阵A相对于矩阵B的广义Rayleigh商 广义Rayleigh商可以只在椭球面SB=x|xRn,xTBx=1

12、上讨论。定理：非零向量x0是R(x)的驻点的充要条件是x0为Ax=Bx的属于特征值的特征向量。推论 若x是Ax=Bx的特征向量，则R(x)是与之对应的特征值。.,定理：设Vk是Rn的任意一个k维子空间，则广义特征值问题Ax=Bx的第k个特征值和第n-k+1个特征值具有下列的极小极大性质推论1 设Vk是Rn的任意一个k维子空间，则实对称矩阵A的第k个特征值和第n-k+1个特征值具有极性质,推论2 设Vn-k+1是Rn的任意一个n-k+1维子空间，则,三.矩阵奇异值的极小极大性质,实矩阵A的奇异值(A)和实对称半正定矩阵ATA的特征值(ATA)有关系定理：设A的奇异值排列为0=1=n-r n-r+1n，则A的第k个奇异值和第n-k+1个奇异值具有极性质:,