计算方法上机实习报告.docVIP

硕士研究生计算方法A班上机实习报告要求 实验一 一、实验题目 二、算法组织 1.算法分析： 分解法 矩阵A为对称正定矩阵，其各阶主子式皆大于0，进行高斯消去法和三角分解法时，不必选主元便可顺利、稳定的进行。对A的杜里特尔分解A=LR中，设R的对角元构成的矩阵为D，由于A对称正定，矩阵A可作如下分解： Cholesky分解法 由于A为对称正定矩阵，在分解的基础上进一步做如下分解： 2.程序代码： !题目二：LDL(T)分解法与CHOLESKY分解法 PROGRAM WORK02 IMPLICIT NONE INTEGER,PARAMETER :: NT=20 INTEGER :: I,J,K,ID DIMENSION :: A(NT,NT),L(NT,NT),LT(NT,NT),D(NT) REAL :: A,L,LT,D !**********构造正定对称矩阵A*********** DO I=1,NT DO J=1,NT IF (I==J)THEN A(I,J)=I ELSE A(I,J)=AMIN0(I,J) END IF END DO END DO PRINT *,请选择分解方法:（1-LDL(T)法分解，2-Cholesky法分解） READ (*,*) ID !**********分解正定对称矩阵A*********** SELECT CASE (ID) CASE (1) CALL DOOLITTLE(A,NT) CASE (2) CALL CHOLESKY(A,NT) CASE DEFAULT PRINT *, Unknown method STOP 1 END SELECT DO I=1,NT D(I)=A(I,I) L(I,I)=1 DO J=I+1,NT L(I,J)=A(I,J) LT(J,I)=A(J,I) END DO END DO !**********输出*********** SELECT CASE (ID) CASE (1) PRINT *,LDL(T)分解法结果： OPEN(10,FILE=PFG.DAT) CASE (2) PRINT *,Cholesky分解法结果： OPEN(10,FILE=CHOLESKY.DAT) CASE DEFAULT PRINT *, Unknown method,ID END SELECT DO I=1,NT WRITE(*,(20F4.1))(A(I,J),J=1,NT) WRITE(10,(20F4.1))(A(I,J),J=1,NT) END DO CLOSE (10) END PROGRAM WORK02 !*********LDL(T)法分解*********** SUBROUTINE DOOLITTLE(A,NT) IMPLICIT NONE INTEGER :: NT,I,J,K DIMENSION :: A(NT,NT) REAL :: A DO I=1,NT DO J=I,NT DO K=1,I-1 A(I,J)=A(I,J)-A(I,K)*A(K,J) END DO END DO DO J=I+1,NT DO K=1,I-1 A(J,I)=A(J,I)-A(J,K)*A(K,I) END DO A(J,I)=A(J,I)/A(I,I) END DO END DO DO I=1,NT DO J=I+1,NT A(I,J)=A(I,J)/A(I,I) END DO END DO END SUBROUTINE DOOLITTLE !*********Cholesky分解************ SUBROUTINE CHOLESKY(A,NT) IMPLICIT NONE INTEGER :: NT,I,J,K DIMENSION :: A(NT,NT) REAL :: A DO I=1,NT DO K=1,I-1 A(I,I)=A(I,I)-A(I,K)**2 END DO A(I,I)=A(I,I)**0.5 DO J=I+1,NT DO K=1,I-1 A(J,I)=A(J,I)-A(J,K)*A(I,K) END DO A(J,I)=A(J,I)/A(I,I) END DO END DO DO I=1,NT DO J=I+1,NT