用QR算法求矩阵的特征值.docVIP

下载本文档

95
0
约3.37千字
约 7页
2017-02-13 发布于重庆
举报
版权申诉

用QR算法求矩阵的特征值.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

用QR算法求矩阵的特征值

一、实验名称：用QR算法求矩阵的特征值二、实验目的：1、通过实验进一步熟悉掌握求矩阵特征值的QR方法及原理。 2、理解QR方法的计算流程。 3、能够编程实现QR方法。三、实验内容：给定矩阵 , ，采用QR方法计算A和H矩阵的全部特征值。四、实验要求：（1）根据QR算法原理编写程序求矩阵及矩阵的全部特征值（要求误差＜10）。（2）直接用MATLAB的内部函数eig求矩阵及矩阵的全部特征值，并与（1）的结果比较。五、QR方法计算矩阵特征值的程序： function [namda,time,data_na]=qr_tz(A,tol) if nargin==1; tol=1e-5; end wucha=1; time=0; while (wuchatol)(time500) [q,r]=qr(A); A1=r*q; tz0=diag(A1); tz1=diag(A); wucha=norm(tz0-tz1); A=A1; time=time+1; data_na(time,:)=tz1; end namda=tz1; disp(‘特征值为’) namda disp(‘第一个特征在值’) time n1=length(data_na); n2=(1:n1)’; temp1=[n2,data_na]; subplot(2,2,1:2) plot(date_na(:,1)) title(‘迭代次数为’) grid subplot(2,2,3) plot(data-na(:,2)) title(‘第二个特征值’) grid subplot(2,2,4) plot(data-na(:,3)) title(‘第三个特征值’) grid 六、实验结果： A=[6,2,1;2,3,1;1,1,1];[namda,time,data_na]=qr_tz(A,1e-5); 特征值为 namda = 7.287984301265196 2.133082001665311 0.578933697069494 迭代次数为 time = 6 图 1 A=[6,2,1;2,3,1;1,1,1];[V,D]=eig(A,nobalance), V = -0.043168204294791 -0.497425032351127 -0.866432249704755 -0.350731446032482 0.819589100011505 -0.453057567982585 0.935480603167125 0.284327354176152 -0.209842790596346 D = 0.578933385691052 0 0 0 2.133074475348525 0 0 0 7.287992138960417 A=[2,3,4,5,6;4,4,5,6,7;0,3,6,7,8;0,0,2,8,9;0,0,0,1,0];[namda,time,data_na]=qr_tz(A,1e-5); 特征值为 namda = 13.172341736418078 6.551888013598945 1.595655961778729 -0.929097677537164 -0.390788034258589 迭代次数为 time = 22 图 2 A=[2,3,4,5,6;4,4,5,6,7;0,3,6,7,8;0,0,2,8,9;0,0,0,1,0];[V,D]=eig(A,nobalance), V = 0.724952325211240 0.639133139255961 0.336044072268703 0.541826077153744 -0.303591564596599 1.000000000000000 1.000000000000000 0.803166939605544 -1.000000000000000 0.314250026625090 0.792999044338331 0.470313263418813 -1.000000000000000 0.659286335558985 -0.352004942324224 0.353299625949037 -0.333345861484