《vb矩阵运算函数.doc》由会员分享,可在线阅读,更多相关《vb矩阵运算函数.doc(3页珍藏版)》请在金锄头文库上搜索。
1、Sub madd(mtxA() As Double, mtxB() As Double, mtxC() As Double) 矩阵相加 Dim i As Integer, j As Integer Dim m As Integer, n As Integer m = UBound(mtxA, 1) - LBound(mtxA, 1) + 1 n = UBound(mtxA, 2) - LBound(mtxA, 2) + 1 For i = 1 To n For j = 1 To m mtxC(i, j) = mtxA(i, j) + mtxB(i, j) Next j Next iEnd Su
2、bSub mcut(mtxA() As Double, mtxB() As Double, mtxC() As Double) 矩阵相减 Dim i As Integer, j As Integer Dim m As Integer, n As Integer m = UBound(mtxA, 1) - LBound(mtxA, 1) + 1 n = UBound(mtxA, 2) - LBound(mtxA, 2) + 1 For i = 1 To n For j = 1 To m mtxC(i, j) = mtxA(i, j) - mtxB(i, j) Next j Next iEnd S
3、ubSub Mtrans(mtxA() As Double, mtxAT() As Double) 矩阵转置 Dim i As Integer, j As Integer Dim m As Integer, n As Integer m = UBound(mtxA, 1) - LBound(mtxA, 1) + 1 n = UBound(mtxA, 2) - LBound(mtxA, 2) + 1 For i = 1 To n For j = 1 To m mtxAT(i, j) = mtxA(j, i) Next j Next iEnd SubSub mmul(mtxA() As Doubl
4、e, mtxB() As Double, mtxC() 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 m For j = 1 To l mtxC(i, j) =
5、 0# For K = 1 To n mtxC(i, j) = mtxC(i, j) + mtxA(i, K) * mtxB(K, j) Next K Next 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
6、D As Double, p As DoubleFor K = 1 To nD = 0#For i = K To n For j = K To n p = Abs(mtxA(i, j) If (p D) Then D = p nIs(K) = i nJs(K) = j End IfNext jNext iIf (D + 1# = 1#) Then MRinv = False Exit FunctionEnd IfIf (nIs(K) K) Then For j = 1 To n p = mtxA(K, j) mtxA(K, j) = mtxA(nIs(K), j) mtxA(nIs(K), j
7、) = p Next jEnd IfIf (nJs(K) K) Then For i = 1 To n p = mtxA(i, K) mtxA(i, K) = mtxA(i, nJs(K) mtxA(i, nJs(K) = p Next iEnd IfmtxA(K, K) = 1# / mtxA(K, K)For j = 1 To n If (j K) Then mtxA(K, j) = mtxA(K, j) * mtxA(K, K)Next jFor i = 1 To n If (i K) Then For j = 1 To n If (j K) Then mtxA(i, j) = mtxA
8、(i, j) - mtxA(i, K) * mtxA(K, j) Next j End IfNext iFor i = 1 To n If (i K) Then mtxA(i, K) = -mtxA(i, K) * mtxA(K, K)Next iNext KFor K = n To 1 Step -1 If (nJs(K) K) Then For j = 1 To n p = mtxA(K, j) mtxA(K, j) = mtxA(nJs(K), j) mtxA(nJs(K), j) = p Next jEnd IfIf (nIs(K) K) ThenFor i = 1 To n p = mtxA(i, K) mtxA(i, K) = mtxA(i, nIs(K) mtxA(i, nIs(K) = pNext iEnd IfNext KMRinv = TrueEnd Function