特征值估计与表示.ppt

《特征值估计与表示.ppt》由会员分享，可在线阅读，更多相关《特征值估计与表示.ppt（26页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、关于特征值估计与表示现在学习的是第1页，共26页5.1特征值的估计特征值的估计 一、特征值的界一、特征值的界 1.定定理理5.1：设A=(aij)Rnn，若表示A的任一特征值，则其中 。2.推论推论 实对称矩阵的特征值都是实数。现在学习的是第2页，共26页3.引理引理1：设BCnn，yCn为单位列向量，则证明：设B=(bij)nn，则现在学习的是第3页，共26页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、是有现在学习的是第4页，共26页例：估计矩阵 特征值的上界。5.推论推论 Hermite矩阵的特征值都是实数，反Hermite矩阵的特征值为零或纯虚数解：由定理5.2，对A特征值，有：|2，|Re()|2，|Im()|1.3，由定理5.1，知其虚部的另一逼近为：其特征值为：现在学习的是第5页，共26页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

3、.2 设ACnn，如果AT按行严格对角占优，则称A按按列列严严格格对对角角占占优优；如果AT按行(弱)对角占优、则称A按按列列(弱弱)对角占优对角占优。现在学习的是第6页，共26页二、特征值的包含区域二、特征值的包含区域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个盖尔圆的并集之内。证明:设为其特征值，为对应特征向量，且 为其绝对值最大者，则有即现在学习的是第7页，共26页3.定定理理5.

4、7 由矩阵A的所有盖尔圆组成的连通部分中任取一个，如果它是由k个盖尔圆构成的，则在这个连通部分中有且仅有A的k个特征值(盖尔圆相重时重复计数特征值相同时也重复计数)证 明 思 路：考 虑 由A的 对 角 线 元 素 构 成 的 矩 阵D=diag(a11,a22,ann)，定义矩阵B(u)=(1-u)D+uA则其特征值变化连续依赖于参数u，D的盖尔圆连续变化成为A的盖尔圆。因此现在学习的是第8页，共26页例：讨论矩阵 的特征值的分布。解：解：A A的盖尔圆分别为的盖尔圆分别为|z-1|0.8|z-1|0.8和和|z|0.5|z|0.5，这两，这两个盖尔圆为连通的，因此包含两个特征值。其特征个盖

5、尔圆为连通的，因此包含两个特征值。其特征值为值为不在盖尔圆不在盖尔圆|z|0.5|z|0.5内。内。现在学习的是第9页，共26页考虑满秩对角阵则矩阵则矩阵DAD-1与与A具有同样的特征值，因此有具有同样的特征值，因此有 若将若将Ri改作改作ri=j i(|aij|i/j)，则两个则两个盖尔盖尔定理仍然定理仍然成立，其中成立，其中 i都是正数。都是正数。现在学习的是第10页，共26页隔离矩阵特征值原则隔离矩阵特征值原则选取的一般方法是：观察A的n个盖尔圆，欲使第i个盖尔圆Gi的半径变大(或小)些，就取i1(或i1)而取其它正数=1。此时，B=DAD-1的第i个盖尔圆的半径变大(或小)，而B的其余

6、盖尔圆的半径相对变小(或变大)但是，这种隔离矩阵特征值的办法还不能用于任意的具有互异特征值的矩阵比如主对角线上有相同元素的矩阵如果矩阵A按行(列)严格对角占优，则detA0。现在学习的是第11页，共26页例例:隔离矩阵隔离矩阵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的(也是

7、A的)一个特征值(见右图)现在学习的是第12页，共26页定理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现在学习的是第13页，共26页5.2 广义特征值问题广义特征值问题 定义定义:称 Ax=Bx 的特征值问题为(对称)矩阵矩阵A相对于相对于(对

8、称对称)矩阵矩阵B的广义特征值问题的广义特征值问题，称数为矩阵矩阵A相对于矩阵相对于矩阵B的特征值的特征值；而与相对应的非零解x称之为属于属于的特征向量的特征向量广义特征值由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 S 的普通特征值问题现在学习的是第14页，共26页因此，当 B=GGT 正定时有正交矩阵P，使得令 Q=G-TP,则有设A与B为正

9、定对称阵，则A+B仍为正定对称阵，由以上结论，存在可逆矩阵Q，使得因此有现在学习的是第15页，共26页二、特征向量的共轭性二、特征向量的共轭性1.在等价的普通特征值问题 Sy=y中，特征向量系y1,y2,yn是完备的标准正交系。令xj=GTyj，j=1,2,n，则有xiTBxj=xiTGGTxj=(GTxi)T(GTxj)=yiTyj=ij，向量系x1，xn称为按按B标准正交化向量系。标准正交化向量系。2.2.按按B标准正交化向量系标准正交化向量系的性质：性质1 xj0(j=1,2,n)(j=1，n)；性质2 x1，xn线性无关。现在学习的是第16页，共26页5.3 对称矩阵特征值的极性对称矩

10、阵特征值的极性一、实对称矩阵的一、实对称矩阵的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)现在学习的是第17页，共26页性质性质3 3 xL(x0)(x00)时，R(x)是一常数性性质质4 R(x)的最大值和最小值存在，且能够在单位球面S=x|xRn,|x|2=1上达到 证：S是闭集，在S上R(x)=xTAx连续，所以必有x1，x2S，使得 minxSR(x)=R(x1)maxxSR(

11、x)=R(x2)任取0yRn，令y0=y/|y|2，则y0S，根据性质3，有R(y)=R(y0)，从而R(x1)R(y)R(x2)。现在学习的是第18页，共26页实对称矩阵A的特征值(都是实数)按其大小升序排列：12n，对应的标准正交特征向量系设为 P P=p1，pn，则有 定理定理：设A为实对称矩阵，则minxSR(x)=1，maxxSR(x)=n 证：任取xS，则x=P Pc,|c|=1,Ax=AP Pc=P P c R(x)=xTAx=cT c1R(x)n，Api=ipiR(pi)=i。现在学习的是第19页，共26页推论1：在S上p1和pn分别是R(x)的一个极小点和极大点，即R(p1)

12、=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。现在学习的是第20页，共26页Courant-Fischer定理的证明的证明构造Rn的子空间Wk=L(pk,pn)，则dimWk=n-k+1由于Vk+Wk Rn，所以ndim(Vk+Wk)=dim(Vk)+dim(Wk)-dim(VkWk)=n+1-dim(VkWk)d

13、im(VkWk)1 故存在x0=pk,pn VkWk，|2=1满足|x0|2=1使得xTAx=T k，即maxxTAx|xVk,|x|2=1k 根据Vk的任意性，可得:现在学习的是第21页，共26页令Vk=L(p1,pk)，取x=pk,pn Vk满足|x|2=l，则有xTAxk，即maxxTAx|xVk,|x|2=1k 于是现在学习的是第22页，共26页二、广义特征值的极小极大原理二、广义特征值的极小极大原理1.定义定义：设A，B为n阶实对称矩阵，且B正定，xRn称R(x)=(xTAx)/(xTBx)，x0为矩阵矩阵A相对于矩阵相对于矩阵B的广的广义义Rayleigh商商 2.广义Raylei

14、gh商可以只在椭球面SB=x|xRn,xTBx=1上讨论。3.定定理理：非零向量x0是R(x)的驻点的充要条件是x0为Ax=Bx的属于特征值的特征向量。4.推推论论 若x是Ax=Bx的特征向量，则R(x)是与之对应的特征值。.现在学习的是第23页，共26页5.定定理理：设Vk是Rn的任意一个k维子空间，则广义特征值问题Ax=Bx的第k个特征值和第n-k+1个特征值具有下列的极小极大性质6.推推论论1 设Vk是Rn的任意一个k维子空间，则实对称矩阵A的第k个特征值和第n-k+1个特征值具有极性质现在学习的是第24页，共26页7.推推论论2 设Vn-k+1是Rn的任意一个n-k+1维子空间，则现在学习的是第25页，共26页感感谢谢大大家家观观看看现在学习的是第26页，共26页