幂法求矩阵主特征值.docxVIP

下载本文档

12
0
约2.59千字
约 2页
2017-06-09 发布于北京
举报
版权申诉

幂法求矩阵主特征值.docx

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

!程序说明：幂法求矩阵主特征值!日期：2010年11月30日PROGRAM Matrix_EigenValuePARAMETER(N=3)REAL ARR(N,N)CALL INPUT(ARR,N)CALL MATEV(ARR,N)END PROGRAMSUBROUTINE INPUT(ARR,N)REAL ARR(N,N)OPEN(1,FILE=MAT.TXT)READ(1,*)((ARR(I,J),J=1,N),I=1,N)END SUBROUTINESUBROUTINE MATEV(ARR,N)PARAMETER(EPS=1E-7)REAL :: ARR(N,N),X(N),X1(N),MAX=0INTEGER :: K=0,P=0X=RESHAPE((/1,1,1/),(/3/))WRITE(1,*) 迭代次数 U（规范化向量） MAX(V)（主特征值）DO WHILE(P/=N)WRITE(1,(I6,A,NF12.6,A,F12.6)) K, (,X, ),MAXP=0MAX=0DO I=1,NX1(I)=0DO J=1,NX1(I)=X1(I)+ARR(I,J)*X(J)!迭代过程ENDDOENDDODO I=1,NIF(ABS(X1(I))ABS(MAX)) MAX=X1(I)!选取主特征值ENDDODO I=1,NIF(ABS(X(I)-X1(I)/MAX)EPS) P=P+1X(I)=X1(I)/MAX!规范化ENDDOK=K+1ENDDOEND SUBROUTINE输出结果：1 1 0.51 1 0.250.5 0.25 2 迭代次数 U（规范化向量） MAX(V)（主特征值） 0 ( 1.000000 1.000000 1.000000 ) 0.000000 1 ( 0.909091 0.818182 1.000000 ) 2.750000 2 ( 0.837607 0.743590 1.000000 ) 2.659091 3 ( 0.799016 0.703035 1.000000 ) 2.604701 4 ( 0.777415 0.680338 1.000000 ) 2.575267 5 ( 0.765108 0.667406 1.000000 ) 2.558792 6 ( 0.758025 0.659963 1.000000 ) 2.549406 7 ( 0.753925 0.655655 1.000000 ) 2.544003 8 ( 0.751544 0.653153 1.000000 ) 2.540876 9 ( 0.750158 0.651697 1.000000 ) 2.539060 10 ( 0.749351 0.650848 1.000000 ) 2.538003 11 ( 0.748880 0.650354 1.000000 ) 2.537387 12 ( 0.748606 0.650065 1.000000 ) 2.537028 13 ( 0.748445 0.649897 1.000000 ) 2.536819 14 ( 0.748352 0.649799 1.000000 ) 2.536697 15 ( 0.748298 0.649741 1.000000 ) 2.536626 16 ( 0.748266 0.649708 1.000000 ) 2.536584 17 ( 0.748247 0.649688 1.000000 ) 2.536560 18 (