特征值估计与表示_第1页
特征值估计与表示_第2页
特征值估计与表示_第3页
特征值估计与表示_第4页
特征值估计与表示_第5页
已阅读5页,还剩21页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

关于特征值估计与表示第1页,共26页,2022年,5月20日,7点57分,星期五§5.1特征值的估计

一、特征值的界

1.定理5.1:设A=(aij)Rn×n,若表示A的任一特征值,则其中。2.推论实对称矩阵的特征值都是实数。第2页,共26页,2022年,5月20日,7点57分,星期五3.引理1:设BCn×n,yCn为单位列向量,则证明:设B=(bij)

n×n,,则第3页,共26页,2022年,5月20日,7点57分,星期五4.定理5.2:设ACn×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=,因此由引理,于是有第4页,共26页,2022年,5月20日,7点57分,星期五例:估计矩阵特征值的上界。5.推论

Hermite矩阵的特征值都是实数,反Hermite矩阵的特征值为零或纯虚数.解:由定理5.2,对A特征值,有:|

|2,|Re()|2,|Im()|1.3,由定理5.1,知其虚部的另一逼近为:其特征值为:第5页,共26页,2022年,5月20日,7点57分,星期五6.定义5.1设A=(aij)Cn×n,记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设ACn×n,如果AT按行严格对角占优,则称A按列严格对角占优;如果AT按行(弱)对角占优、则称A按列(弱)对角占优。

第6页,共26页,2022年,5月20日,7点57分,星期五二、特征值的包含区域1.定义5.3

设A=(aij)Cn×n,称区域Gi:|z-aii|Ri为矩阵A的第i个盖尔圆,其中Ri=ji|aij|称为盖尔圆Gi的半径(i=l,…,n)。2.定理5.6矩阵A=(aij)Cn×n的一切特征值都在它的n个盖尔圆的并集之内。证明:设λ为其特征值,为对应特征向量,且为其绝对值最大者,则有即第7页,共26页,2022年,5月20日,7点57分,星期五3.定理5.7由矩阵A的所有盖尔圆组成的连通部分中任取一个,如果它是由k个盖尔圆构成的,则在这个连通部分中有且仅有A的k个特征值(盖尔圆相重时重复计数.特征值相同时也重复计数).证明思路:考虑由A的对角线元素构成的矩阵D=diag(a11,a22,…,ann),定义矩阵B(u)=(1-u)D+uA则其特征值变化连续依赖于参数u,D的盖尔圆连续变化成为A的盖尔圆。因此第8页,共26页,2022年,5月20日,7点57分,星期五例:讨论矩阵的特征值的分布。解:A的盖尔圆分别为|z-1|≤0.8和|z|≤0.5,这两个盖尔圆为连通的,因此包含两个特征值。其特征值为不在盖尔圆|z|≤0.5内。第9页,共26页,2022年,5月20日,7点57分,星期五考虑满秩对角阵则矩阵DAD-1与A具有同样的特征值,因此有若将Ri改作ri=ji(|aij|i/j)

,则两个盖尔定理仍然成立,其中i都是正数。第10页,共26页,2022年,5月20日,7点57分,星期五隔离矩阵特征值原则选取的一般方法是:观察A的n个盖尔圆,欲使第i个盖尔圆Gi的半径变大(或小)些,就取i>1(或i<1).而取其它正数=1。此时,B=DAD-1的第i个盖尔圆的半径变大(或小),而B的其余盖尔圆的半径相对变小(或变大).但是,这种隔离矩阵特征值的办法还不能用于任意的具有互异特征值的矩阵.比如主对角线上有相同元素的矩阵.如果矩阵A按行(列)严格对角占优,则detA0。第11页,共26页,2022年,5月20日,7点57分,星期五例:

隔离矩阵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的)一个特征值(见右图).第12页,共26页,2022年,5月20日,7点57分,星期五定理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第13页,共26页,2022年,5月20日,7点57分,星期五§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的普通特征值问题第14页,共26页,2022年,5月20日,7点57分,星期五因此,当B=GGT

正定时有正交矩阵P,使得令Q=G-TP,则有设A与B为正定对称阵,则A+B仍为正定对称阵,由以上结论,存在可逆矩阵Q,使得因此有第15页,共26页,2022年,5月20日,7点57分,星期五二、特征向量的共轭性1.在等价的普通特征值问题Sy=y中,特征向量系y1,y2,…,yn是完备的标准正交系。令xj=G-Tyj,j=1,2,…,n,则有xiTBxj=xiTGGTxj=(GTxi)T(GTxj)=yiTyj=ij,向量系x1,…,xn称为按B标准正交化向量系。2.

按B标准正交化向量系的性质:性质1xj0(j=1,2,…,n)(j=1,…,n);性质2x1,…,

xn线性无关。第16页,共26页,2022年,5月20日,7点57分,星期五§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)第17页,共26页,2022年,5月20日,7点57分,星期五性质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(x)=R(x2)

任取0yRn,令y0=y/||y||2,则y0S,根据性质3,有R(y)=R(y0),从而R(x1)R(y)R(x2)。第18页,共26页,2022年,5月20日,7点57分,星期五实对称矩阵A的特征值(都是实数)按其大小升序排列:12…n,对应的标准正交特征向量系设为P=[p1,…,pn],则有

定理:设A为实对称矩阵,则minxSR(x)=1,maxxSR(x)=n

证:任取xS,则x=Pc,||c||=1,Ax=APc=PcR(x)=xTAx=cTc1R(x)n,Api=ipiR(pi)=i。第19页,共26页,2022年,5月20日,7点57分,星期五推论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维子空间,1<k<n。第20页,共26页,2022年,5月20日,7点57分,星期五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)dim(VkWk)1

故存在x0=[pk,…,pn]VkWk,||||2=1满足||x0||2=1使得xTAx=Tk,即max{xTAx|xVk,||x||2=1}k

根据Vk的任意性,可得:第21页,共26页,2022年,5月20日,7点57分,星期五令Vk=L(p1,…,pk),取x=[pk,…,pn]Vk满足||x||2=l,则有xTAxk,即max{xTAx|xVk,||x||2=1}k

于是第22页,共26页,2022年,5月2

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论