1、Sub madd(mtxA() As Double, mtxB() As Double, mtxC() As Double) 矩阵相加Dim i As Integer, j As IntegerDim m As Integer, n As Integerm = UBound(mtxA, 1) - LBound(mtxA, 1) + 1n = UBound(mtxA, 2) - LBound(mtxA, 2) + 1For i = 1 To nFor j = 1 To mmtxC(i, j) = mtxA(i, j) + mtxB(i, j)Next jNext iEnd SubSub mcut
2、(mtxA() As Double, mtxB() As Double, mtxC() As Double) 矩阵相减Dim i As Integer, j As IntegerDim m As Integer, n As Integerm = UBound(mtxA, 1) - LBound(mtxA, 1) + 1n = UBound(mtxA, 2) - LBound(mtxA, 2) + 1For i = 1 To nFor j = 1 To mmtxC(i, j) = mtxA(i, j) - mtxB(i, j)Next jNext iEnd SubSub Mtrans(mtxA(
3、) As Double, mtxAT() As Double) 矩阵转置Dim i As Integer, j As IntegerDim m As Integer, n As Integerm = UBound(mtxA, 1) - LBound(mtxA, 1) + 1n = UBound(mtxA, 2) - LBound(mtxA, 2) + 1For i = 1 To nFor j = 1 To mmtxAT(i, j) = mtxA(j, i)Next jNext iEnd SubSub mmul(mtxA() As Double, mtxB() As Double, mtxC()
4、 As Double) 矩阵相乘Dim m As IntegerDim n As IntegerDim l As IntegerDim i As Integer, j As Integer, K As Integerm = UBound(mtxA, 1) - LBound(mtxA, 1) + 1n = UBound(mtxA, 2) - LBound(mtxA, 2) + 1l = UBound(mtxB, 2) - LBound(mtxB, 2) + 1For i = 1 To mFor j = 1 To lmtxC(i, j) = 0#For K = 1 To nmtxC(i, j) =
5、 mtxC(i, j) + mtxA(i, K) * mtxB(K, j)Next KNext jNext iEnd SubFunction MRinv(mtxA() As Double) As Boolean 矩阵求逆Dim n As Integern = UBound(mtxA, 1) - LBound(mtxA, 1) + 1ReDim nIs(0 To n) As Integer, nJs(0 To n) As IntegerDim i As Integer, j As Integer, K As IntegerDim D As Double, p As DoubleFor K = 1
6、 To nD = 0#For i = K To nFor j = K To np = Abs(mtxA(i, j)If (p D) ThenD = pnIs(K) = inJs(K) = jEnd IfNext jNext iIf (D + 1# = 1#) ThenMRinv = FalseExit FunctionEnd IfIf (nIs(K) K) ThenFor i = 1 To np = mtxA(i, K)mtxA(i, K) = mtxA(i, nJs(K)mtxA(i, nJs(K) = pNext iEnd IfmtxA(K, K) = 1# / mtxA(K, K)For j = 1 To nIf (j K) ThenFor j = 1 To nIf (j K) Then mtxA(i, K) = -mtxA(i, K) * mtxA(K, K)Next iNext KFor K = n To 1 Step -1If (nJs(K) K) ThenFor i = 1 To np = mtxA(i, K)mtxA(i, K) = mtxA(i, nIs(K)mtxA(i, nIs(K) = pNext iEnd IfNext KMRinv = TrueEnd Function
