第五章特征值估计与表示.ppt

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

5. 推论推论 Hermite矩阵的特征值都是实数，反Hermite矩阵的特征值为零或纯虚数．解：由定理5.2，对A特征值 ，有：|  | 2，|Re()|2， |Im()|1.3，由定理5.1，知其虚部的另一逼近为：其特征值为：6. 定义定义5.1 设A=(aij)Cn×n,记Rr=sr|ars|,r=1,…,n，如果|arr|＞Rr (r=1, 2, … , n)，则称矩阵A按行严格按行严格对角占优对角占优；如果|arr|Rr (r=1，…，n)，且有lron，使得|aroro|>Rro成立，则称矩阵A按行按行(弱弱)对角占优对角占优7. 定定义义5.2 设ACn×n，如果AT按行严格对角占优，则称A按按列列严严格格对对角角占占优优；如果AT按行(弱)对角占优、则称A按列按列(弱弱)对角占优对角占优 二、特征值的包含区域二、特征值的包含区域1. 定定义义5.3 设A=(aij)Cn×n，称区域 Gi: |z-aii|Ri 为矩阵A的第第i个个盖盖尔尔圆圆，其中 Ri=ji|aij| 称为盖尔圆盖尔圆Gi的半径的半径(i=l，…，n)2. 定定理理5.6 矩阵A=(aij)Cn×n的一切特征值都在它的n个盖尔圆的并集之内。

证明:设λ为其特征值， 为对应特征向量，且 为其绝对值最大者，则有即3.定定理理5.7 由矩阵A的所有盖尔圆组成的连通部分中任取一个，如果它是由k个盖尔圆构成的，则在这个连通部分中有且仅有A的k个特征值(盖尔圆相重时重复计数．特征值相同时也重复计数)．证 明 思 路 ： 考 虑 由 A的 对 角 线 元 素 构 成 的 矩 阵D=diag(a11,a22,…,ann)，定义矩阵B(u)=(1-u)D+uA则其特征值变化连续依赖于参数u，D的盖尔圆连续变化成为A的盖尔圆因此例：讨论矩阵 的特征值的分布解：解：A A的盖尔圆分别为的盖尔圆分别为|z-1||z-1|≤0.8≤0.8和和|z||z|≤0.5≤0.5，这两，这两个盖尔圆为连通的，因此包含两个特征值其特征个盖尔圆为连通的，因此包含两个特征值其特征值为值为不在盖尔圆不在盖尔圆|z||z|≤0.5≤0.5内考虑满秩对角阵则矩阵则矩阵DAD-1与与A具有同样的特征值，因此有具有同样的特征值，因此有 •若将若将Ri改作改作ri= j i(|aij| i/ j) ，，则两个则两个盖尔盖尔定理定理仍然成立，其中仍然成立，其中 i都是正数。

都是正数隔离矩阵特征值原则隔离矩阵特征值原则•选取的一般方法是：观察A的n个盖尔圆，欲使第i个盖尔圆Gi的半径变大(或小)些，就取i＞1(或i＜1)．而取其它正数=1•此时，B=DAD-1的第i个盖尔圆的半径变大(或小)，而B的其余盖尔圆的半径相对变小(或变大)．•但是，这种隔离矩阵特征值的办法还不能用于任意的具有互异特征值的矩阵．比如主对角线上有相同元素的矩阵．•如果矩阵A按行(列)严格对角占优，则detA0例例: 隔隔离离矩矩阵阵A= 的的特特征征值值．． •A的3个盖尔圆为G1: |z-20|5.8，G2: |z-10|5，G3: |z-10j|3G1与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的)一个特征值(见右图)．定理5.11：设矩阵A=(aij)Cn×n的， 0α1，λ是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/2[R1(AT)]1/2=0.41/2=0.6324|λ |  [R2(A)]1/2[R2(AT)]1/2=0.41/2=0.6324§§5.2 广义特征值问题广义特征值问题 •定义定义: 称 Ax=Bx 的特征值问题为(对称)矩阵矩阵A相对于相对于(对对称称)矩阵矩阵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 的普通特征值问题因此，当 B=GGT 正定时有正交矩阵P，使得令 Q = G-TP , 则有设A与B为正定对称阵，则A+B仍为正定对称阵，由以上结论，存在可逆矩阵Q，使得因此有二、特征向量的共轭性二、特征向量的共轭性1. 在等价的普通特征值问题 Sy=y中，特征向量系y1, y2,…,yn是完备的标准正交系。

令xj=G－Tyj，j=1,2,…, n，则有xiTBxj=xiTGGTxj=(GTxi)T(GTxj)=yiTyj=ij，向量系x1，…，xn称为按按B标准正交化向量系标准正交化向量系2.2. 按按B标准正交化向量系标准正交化向量系的性质：•性质1 xj0 (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 3 xL(x0)( x00)时，R(x)是一常数．•性性质质4 R(x)的最大值和最小值存在，且能够在单位球面S={x|xRn,||x||2=1}上达到． 证：S是闭集，在S上R(x)=xTAx连续，所以必有x1，x2S，使得 minxSR(x)=R(x1)maxxSR(x)=R(x2) 任取0yRn，令y0=y/||y||2，则y0S，根据性质3，有R(y)=R(y0)，从而R(x1)R(y)R(x2)。

实对称矩阵A的特征值(都是实数)按其大小升序排列：12…n，对应的标准正交特征向量系设为 P P = [p1，…，pn]，则有 •定理定理：设A为实对称矩阵，则minxSR(x)=1，maxxSR(x)=n 证：任取xS，则x=P Pc, ||c||=1, Ax=AP Pc=P P c R(x)=xTAx=cT c1R(x)n，Api=ipiR(pi)=i•推论1：在S上p1和pn分别是R(x)的一个极小点和极大点，即R(p1)=1，R(pn)=n •推论2 若1=…=k (1kn)．则在||x||2=l上R(x)的所有极小点为[p1，…，pk] ，|| ||2=1•定理：设xL(pr，…，ps) ， 1rsn ，则有minxR(x)=r，maxxR(x)=s •Courant-Fischer定理定理: 设实对称矩阵A的特征值按升序排列，则A的第k个特征值其中Vk是Rn的任意—个k维子空间，1＜k＜nCourant-Fischer定理的证明的证明•构造Rn的子空间Wk=L(pk,…,pn) ，则dimWk=n-k+1．由于Vk+Wk Rn，所以ndim(Vk+Wk)=dim(Vk)+dim(Wk)-dim(VkWk)=n+1-dim(VkWk)  dim(VkWk)1 故存在x0=[pk,…,pn] VkWk，|| ||2=1满足||x0||2=1使得xTAx= T k，即max{xTAx |xVk, || x||2=1}k 根据Vk的任意性，可得:•令Vk=L(p1,…,pk)，取x =[pk,…,pn] Vk满足||x||2=l，则有xTAxk，即max{xTAx|xVk, || x||2=1}k 于是二、广义特征值的极小极大原理二、广义特征值的极小极大原理1.定义定义 ：设A，B为n阶实对称矩阵，且B正定，xRn．称R(x)=(xTAx)/(xTBx)， x0为矩阵矩阵A相相对于矩阵对于矩阵B的广义的广义Rayleigh商商． ． 2.广义Rayleigh商可以只在椭球面SB={x|xRn, xTBx =1}上讨论。

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