幂法求矩阵主特征值.doc

《幂法求矩阵主特征值.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.