第三讲 消去法求逆矩阵.doc

第三讲 消去法求逆矩阵 并行变换法 AX=B 经n次初等变换 (k) 表示第k步消元 k=1……..m I=1,2,….,k-1,k+1…….m (i≠k) J=1,2,….,n,n+1,……..,2m+1 SUB INVERSE1(M as Integer) for K=1 To M for I=1 to M IF IK then c = A(I, K) / A(K, K) For J=1 to 2*M+1 A(I,J)=A(I,J)-C* A(K,J) Next j Endif Next I Next k For k=1 to M C=a(K,K) For J=1 to 2*M+1 A(K,J)=A(K,J)/C Next j Next k END sub 紧凑变换法 AX=B 经n次初等变换 I,J≠K J≠K I≠K (k) 表示第k步消元 k=1……..n I=1,2,….,n J=1,2,….,n,n+1 SUB INVERSE2(M as integer) for K=1 to M FOR I=1 TO M IF IK THEN FOR J=1,M+1 IF JK THEN A(I,J)=A(I,J)-A(I,K)*A(K,J)/A(K,K) ENDIF NEXT J End if Next i for I=1 to M IF (IK) then A(I,K)=-A(I,K)/A(K,K) endif Next i for J=1 to M+1 IF (jK) then A(K,J)=A(K,J)/A(K,K) ENDIF Next j A(K,K)=1/A(K,K) Next k END 列主元消去法求逆 Sub InvZuQu(M As Integer) Dim L1(5, 5) As Double Dim L2(5, 5) As Double Dim TT As Double For I = 1 To M For J = 1 To M If (I = J) Then L1(I, J) = 1# Else L1(I, J) = 0# End If Next J Next I For k = 1 To M For L = k To M If (Abs(A(k, k)) Abs(A(L, k))) Then For M = 1 To M t = A(k, M) A(k, M) = A(L, M) TT = L1(k, M) L1(k, M) = L1(L, M) L1(L, M) = TT A(L, M) = t Next M End If Next L For I = 1 To M If I k Or I k Then L2(I, k) = A(I, k) / A(k, k) For J = 1 To M A(I, J) = A(I, J) - L2(I, k) * A(k, J)