数值分析幂法和反幂法_matlab程序.docVIP

下载本文档

45
0
约 5页
2017-08-17 发布于安徽
举报
版权申诉

数值分析幂法和反幂法_matlab程序.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文档。上传文档

数值分析幂法与反幂法 matlab程序随机产生一对称矩阵，对不同的原点位移和初值(至少取3个)分别使用幂法求计算矩阵的主特征值及主特征向量，用反幂法求计算矩阵的按模最小特征值及特征向量。要求 1）比较不同的原点位移和初值说明收敛性 2）给出迭代结果，生成DOC文件。 3）程序清单，生成M文件。wo 解答： A=rand(5) %随机产生5*5矩阵求随机矩阵 A = 0.7094 0.1626 0.5853 0.6991 0.1493 0.7547 0.1190 0.2238 0.8909 0.2575 0.2760 0.4984 0.7513 0.9593 0.8407 0.6797 0.9597 0.2551 0.5472 0.2543 0.6551 0.3404 0.5060 0.1386 0.8143 B=A+A %A矩阵和A的转置相加，得到随机对称矩阵B B = 1.4187 0.9173 0.8613 1.3788 0.8044 0.9173 0.2380 0.7222 1.8506 0.5979 0.8613 0.7222 1.5025 1.2144 1.3467 1.3788 1.8506 1.2144 1.0944 0.3929 0.8044 0.5979 1.3467 0.3929 1.6286 B= 编写幂法、反幂法程序： function [m,u,index,k]=pow(A,u,ep,it_max) % 求矩阵最大特征值的幂法，其中 % A为矩阵； % ep为精度要求，缺省为1e-5； % it_max为最大迭代次数，缺省为100； % m为绝对值最大的特征值； % u为对应最大特征值的特征向量； % index，当index=1时，迭代成功，当index=0时，迭代失败 if nargin4 it_max=100; end if nargin3 ep=1e-5; end n=length(A); index=0; k=0; m1=0; m0=0.01; % 修改移位参数，原点移位法加速收敛，为0时，即为幂法 I=eye(n) T=A-m0*I while k=it_max v=T*u; [vmax,i]=max(abs(v)); m=v(i); u=v/m; if abs(m-m1)ep; index=1; break; end m=m+m0; m1=m; k=k+1; end function[m,u,index,k]=pow_inv(A,u，ep,it_max) % 求矩阵最大特征值的反幂法，其中 % A为矩阵； % ep为精度要求，缺省为1e-5； % it_max为最大迭代次数，缺省为100； % m为绝对值最大的特征值； % u为对应最大特征值的特征向量； % index，当index=1时，迭代成功，当index=0时，迭代失败 if nargin4 it_max=100; end if nargin3 ep=1e-5; end n=length(A); index=0; k=0; m1=0; m0=0; % 修改移位参数，原点移位法加速收敛，为0时，即为反幂法 I=eye(n); T=A-m0*I; invT=inv(T); while k=it_max v=invT*u; [vmax,i]=max(abs(v)); m=v(i); u=v/m; if abs(m-m1)ep index=1; break; end m1=m; k=k+1; end m=1/m; m=m+m0; 修改输入的m0的值，所得结果：幂法： m u index k 0 5.2710 1 10 0.5 5.2710 1 13 0.3 5.2710 1 11 0.1 5.2710 1 10 0 5.2710 1 9 0.00001 5.2710 1 29 0.000005 5.2710 1 9 0.000001 5.2710 1 9 0 5.2710 1 10 0.00002 5.2710 0 1