- 1、本文档共43页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
* * 7.2 Jacobi方法 所以 再取i=2,j=3,aij(1)=a23(1)=2.020190,类似地可得 * * 7.2 Jacobi方法 以下依次有 * * 7.2 Jacobi方法 从而A的特征值可取为 ?1?2.125825, ?2?8.388761, ?3?4.485401 特征向量为 R1TR2T…RkT * * 7.3 QR方法 QR方法是目前求一般矩阵全部特征值的最有效并广泛应用的方法之一,在实际应用中,经常把一般的矩阵经过正交相似变换化成上Hessenberg矩阵,再应用QR方法求其特征值和特征向量. 6.3.1 基本QR方法 设A是n?n非奇异实矩阵。令 A1=A,对 A1作QR分解,得 A1= Q1R1 令 A2= R1Q1 又对A2作QR分解,得 A3= R2Q2 * * 7.3 QR方法 反复进行这种运算,就得到一个矩阵序列{Ak}, 由于 Rk= QTkAk 因而Ak+1与Ak(k=1,2,…)相似。所有Ak都与A 相似,它们都有共同的特征值。 在一定的条件下,序列{Ak}收敛于上三角(或分块上三角)矩阵,其主对角线元素(或子块)有确定的极限。如果收敛于上三角矩阵,则这个上三角矩阵的主对角线元素就是矩阵A的特征值;如果收敛于分块上三角矩阵,则主对角线上各个子块的特征值就是矩阵A的特征值 故有 Ak+1= QTkAk Qk, k = 1,2,… * * 7.3 QR方法 定理 8.2 设A 为n?n实矩阵,它的特征值λj(j=1,2,…,n)满足|λ1||λ2| … |λn| 0,则由基本QR方法产生的矩阵序列{Ak}收敛于一个上三角矩阵(其主对角线以上的元素可以不收敛)。若A是对称矩阵,则序列{Ak}收敛于一个对角矩阵. 定理 8.3 对于任何n?n 实矩阵A,由基本QR方法产生的矩阵序列{Ak}收敛于分块上三角矩阵(其主对角线子块以上的元素可以不收敛),并且每个主对角线子块有等模的特征值. * * 7.3 QR方法 要实现一步QR 方法的迭代,就需要做一次QR 分解,再做一次矩阵相乘,当A是一般矩阵时,计算量是很大的。 为了减少计算量,在实际计算时,总是先将原矩阵A经过相似变换,化为上Hessenberg矩阵 A1 =(a(1)ij)n?n, 其中,当i j+1 时,a(1)ij =0。然后对A1使用 QR 方法。容易看出,从上Hessenberg矩阵A1出发,经过迭代所得的Ak(k=2,3,…)全都是上Hessenberg矩阵。事实上,设Ak是上Hessenberg矩阵,则由于Rk?1是上Hessenberg矩阵,故 Qk= Ak Rk?1 是上Hessenberg矩阵,因而 Ak+1= RkQk 也是上Hessenberg矩阵。 * * 7.3 QR方法 6.3.3 带原点平移的QR 已经证明,基本QR方法的收敛速度一般是线性收敛的.为了加速收敛,引进带原点平移的QR 方法. 设A是n?n非奇异实矩阵,令A1= A,并取一适当的数u1,对矩阵(A1? u1I)作QR分解 A1? u1I = Q1R1 然后令 A2= R1Q1+u1I 又取一适当的数u2,对矩阵(A2? u2I)作QR分解 A2? u2I = Q2R2 反复进行这种运算,得到迭代公式 * * 7.3 QR方法 称为带原点平移的QR方法,称uk为第k次迭代中的原点平移量 产生的序列有下述性质 (1) Ak+1与Ak相似 (2)当Ak 为上 Hessenberg 阵时,Ak+1也是上 Hessenberg 阵,迭代若干次后,若 足够小,就可以认为 是特征值的近似值,然后将Ak+1的第n行第n列去掉 * * 7.3 QR方法 如果对(Ak?ukI)作QR分解时利用 Householder矩阵,则有 恰当的选用平移量uk可加速分解过程收敛,uk的选取一般有两种 (1) 令 uk= * * 7.3 QR方法 (2) 设A1已经是上 Hessenberg矩阵,则Ak的形式为 先求出Ak最右下角的二阶子矩阵 * * 7.3 QR方法 的两个特征值,然后取其中按模最小的
您可能关注的文档
最近下载
- 第一章 2.2 水量平衡.ppt
- 《GB/T 19326-2022锻制支管座》.pdf
- 2022年11月陕西省从优秀村社区干部中考试录用200名乡镇街道机关公务员上岸冲刺卷I含答案详解版(3套).docx VIP
- 2020年银行业从业人员职业操守和行为准则.pdf VIP
- 转预备党员思想汇报【银行】.pdf VIP
- 【新教材】人教版(2024)七年级上册英语Unit 4 My Favourite Subject教案.docx
- 米厂恒温仓库工程设计方案.docx
- 2024年党校入党积极分子培训考试必考重点知识汇编(共160题).doc VIP
- 《世界经典神话与传说故事》 测试题及答案.pdf
- 智能制造设备安装与调试职业技能等级标准(2021年).pdf
文档评论(0)