资源描述:
《Fortran求逆矩阵子程序》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、!求逆矩阵subroutineinverse(A,IA)implicitnonereal::A(:,:),IA(:,:)real,allocatable::B(:,:)integer::i,j,NN=size(A,1)allocate(B(N,N))!先把IA设定成单位矩阵forall(i=1:N,j=1:N,i==j)IA(i,j)=1.0forall(i=1:N,j=1:N,i/=j)IA(i,j)=0.0!保存原先的矩阵A,使用B来计算B=A!把B化成对角线矩阵(除了对角线外,都为0)callUpper(B,IA,N)!先把B化成上三角矩阵callLo
2、wer(B,IA,N)!再把B化成下三角矩阵!求解forall(i=1:N)IA(i,:)=IA(i,:)/B(i,i)returnendsubroutine!求上三角矩阵的子程序subroutineUpper(M,S,N)implicitnoneinteger::Nreal::M(N,N)real::S(N,N)integer::I,Jreal::EdoI=1,N-1doJ=I+1,NE=M(J,I)/M(I,I)M(J,I:N)=M(J,I:N)-M(I,I:N)*ES(J,:)=S(J,:)-S(I,:)*Eenddoenddoreturnendsub
3、routineUpper!求下三角矩阵的子程序subroutineLower(M,S,N)implicitnoneinteger::Nreal::M(N,N)real::S(N,N)integer::I,Jreal::EdoI=N,2,-1doJ=I-1,1,-1E=M(J,I)/M(I,I)M(J,1:N)=M(J,1:N)-M(I,1:N)*ES(J,:)=S(J,:)-S(I,:)*EenddoenddoreturnendsubroutineLowerprogramqiuniimplicitnoneintegeri,jreal::a(3,3),aa(3
4、,3)print*,"请输入a矩阵"doi=1,3doj=1,3read*,a(i,j)enddoenddocallinverse(a,aa,3)print*,aaend!求逆矩阵subroutineinverse(A,IA,N)implicitnoneinteger::i,j,NrealA(N,N),IA(N,N),B(N,N)!先把IA设定成单位矩阵forall(i=1:N,j=1:N,i==j)IA(i,j)=1.0forall(i=1:N,j=1:N,i/=j)IA(i,j)=0.0!保存原先的矩阵A,使用B来计算B=A!把B化成对角线矩阵(除了对角
5、线外,都为0)callUpper(B,IA,N)!先把B化成上三角矩阵callLower(B,IA,N)!再把B化成下三角矩阵!求解forall(i=1:N)IA(i,:)=IA(i,:)/B(i,i)returnendsubroutine!求上三角矩阵的子程序subroutineUpper(M,S,N)implicitnoneinteger::Nreal::M(N,N)real::S(N,N)integer::I,Jreal::EdoI=1,N-1doJ=I+1,NE=M(J,I)/M(I,I)M(J,I:N)=M(J,I:N)-M(I,I:N)*ES(J
6、,:)=S(J,:)-S(I,:)*EenddoenddoreturnendsubroutineUpper!求下三角矩阵的子程序subroutineLower(M,S,N)implicitnoneinteger::Nreal::M(N,N)real::S(N,N)integer::I,Jreal::EdoI=N,2,-1doJ=I-1,1,-1E=M(J,I)/M(I,I)M(J,1:N)=M(J,1:N)-M(I,1:N)*ES(J,:)=S(J,:)-S(I,:)*EenddoenddoreturnendsubroutineLower