正交变换与QR迭代矩阵特征值计算.ppt_第1页
正交变换与QR迭代矩阵特征值计算.ppt_第2页
正交变换与QR迭代矩阵特征值计算.ppt_第3页
正交变换与QR迭代矩阵特征值计算.ppt_第4页
正交变换与QR迭代矩阵特征值计算.ppt_第5页
已阅读5页,还剩10页未读 继续免费阅读

下载本文档

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

文档简介

1,第八章矩阵特征值计算,计算方法,正交变换与QR迭代,2,本讲内容,正交变换,QR迭代,Householder变换Givens变换QR分解Schur分解Hessenberg矩阵,3,Householder变换,性质,(1)对称:(2)正交:(3)对合:(4)保模:(5),定义:设且,称矩阵,为Householder变换,或初等反射矩阵。,4,Householder变换,定理:设x,yRn,xy且|x|2=|y|2,则存在n阶Householder变换H,使得y=Hx,5,Householder变换,定理:对任意的非零向量xRn,存在Householder变换H,使得Hx=e1其中=sgn(x1)|x|2,e1=(1,0,.,0)T,,的选取是为了防止在实际计算中与x1互相抵消若x1=0,则取=|x|2,6,Givens变换,定义:称矩阵,为Givens变换,或旋转变换。,i,j,7,Givens变换,性质,(1)只有四个元素与单位矩阵不同(2)正交:(3)用G左乘一个矩阵时,只改变该矩阵中两行的值(4)用G右乘一个矩阵时,只改变该矩阵中两列的值,8,Givens变换,定理:设x=(x1,.,xi,.,xj,.,xn)T,且xi,xj不全为零,则存在Givens变换G=G(i,j,),使得,9,QR分解,定理:(QR分解)设n阶实矩阵A非奇异,则存在正交分解A=QR其中Q是正交矩阵,R是非奇异上三角矩阵。若限定R的对角线元素为正数,则此分解唯一。,10,QR分解算法,设,(j=1,.,n),(1)构造H1使得H1a1=1e1,令,(2)构造使得,令,算法(QR分解),11,QR分解算法,以此类推,经过n-1步,可得Householder矩阵H1,H2,.,Hn-1,使得,令,即得,12,QR分解举例,例:用Householder变换计算的QR分解,解:(板书),13,Schur分解,定理:(Schur分解)设A为n阶实矩阵,则存在正交矩阵Q,使得其中Rii是一阶或二阶方阵。,若Rii是一阶方阵,则它就是A的特征值;若Rii是二阶方阵,则其特征值为A的两个共轭复特征值。,拟上三角矩阵,14,QR迭代,QR迭代算法,计算矩阵的所有特征值和特征向量计算过程,(1)令A1A(2)对k=1,

温馨提示

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

评论

0/150

提交评论