欢迎来到天天文库
浏览记录
ID:58452179
大小:16.50 KB
页数:2页
时间:2020-09-03
《幂法求矩阵主特征值.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、!程序说明:幂法求矩阵主特征值!日期:2010年11月30日PROGRAMMatrix_EigenValuePARAMETER(N=3)REALARR(N,N)CALLINPUT(ARR,N)CALLMATEV(ARR,N)ENDPROGRAMSUBROUTINEINPUT(ARR,N)REALARR(N,N)OPEN(1,FILE='MAT.TXT')READ(1,*)((ARR(I,J),J=1,N),I=1,N)ENDSUBROUTINESUBROUTINEMATEV(ARR,N)PARAMETER(EPS
2、=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)(主特征值)'DOWHILE(P/=N)WRITE(1,'(I6,A,F12.6,A,F12.6)')K,'(',X,')',MAXP=0MAX=0DOI=1,NX1(I)=0DOJ=1,NX1(I)=X1(I)+ARR(I,J)*X(J)!迭代过程ENDDOENDDODOI=1,NIF(A
3、BS(X1(I))>ABS(MAX))MAX=X1(I)!选取主特征值ENDDODOI=1,NIF(ABS(X(I)-X1(I)/MAX)4、2.8(0.0.1.)2.9(0.0.1.)2.10(0.0.1.)2.11(0.0.1.)2.12(0.0.1.)2.13(0.0.1.)2.14(0.0.1.)2.15(0.0.1.)2.16(0.0.1.)2.17(0.0.1.)2.18(0.0.1.)2.19(0.0.1.)2.20(0.0.1.)2.21(0.0.1.)2.22(0.0.1.)2.23(0.0.1.)2.24(0.0.1.)2.25(0.0.1.)2.26(0.0.1.)2.
4、2.8(0.0.1.)2.9(0.0.1.)2.10(0.0.1.)2.11(0.0.1.)2.12(0.0.1.)2.13(0.0.1.)2.14(0.0.1.)2.15(0.0.1.)2.16(0.0.1.)2.17(0.0.1.)2.18(0.0.1.)2.19(0.0.1.)2.20(0.0.1.)2.21(0.0.1.)2.22(0.0.1.)2.23(0.0.1.)2.24(0.0.1.)2.25(0.0.1.)2.26(0.0.1.)2.
此文档下载收益归作者所有