重慶分公司,新征程啟航
為企業提供網站建設、域名注冊、服務器等服務
為企業提供網站建設、域名注冊、服務器等服務
矩陣求逆的VB程序
成都創新互聯公司是一家專業提供合肥企業網站建設,專注與成都做網站、成都網站建設、H5網站設計、小程序制作等業務。10年已為合肥眾多企業、政府機構等服務。創新互聯專業網站設計公司優惠進行中。
Private Function MRinv(N As Integer, mtxA() As Double) As Boolean
'****************************************************************************************
' 功能: 實現矩陣求逆的全選主元高斯-約當法
' 參數: n - Integer型變量,矩陣的階數
' mtxA - Double型二維數組,體積為n x n。存放原矩陣A;返回時存放其逆矩陣A-1。
' 返回值:Boolean型,失敗為False,成功為True
'****************************************************************************************
ReDim nIs(N) As Integer, nJs(N) As Integer
Dim i As Integer, j As Integer, k As Integer
Dim D As Double, p As Double
' 全選主元,消元
For k = 1 To N
D = 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 If
Next j
Next i
' 求解失敗
If (D + 1# = 1#) Then
MRinv = False
Exit Function
End If
If (nIs(k) k) Then
For j = 1 To N
p = mtxA(k, j)
mtxA(k, j) = mtxA(nIs(k), j)
mtxA(nIs(k), j) = p
Next j
End If
If (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 i
End If
mtxA(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 j
For i = 1 To N
If (i k) Then
For j = 1 To N
If (j k) Then mtxA(i, j) = mtxA(i, j) - mtxA(i, k) * mtxA(k, j)
Next j
End If
Next i
For i = 1 To N
If (i k) Then mtxA(i, k) = -mtxA(i, k) * mtxA(k, k)
Next i
Next k
' 調整恢復行列次序
For 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 j
End If
If (nIs(k) k) Then
For i = 1 To N
p = mtxA(i, k)
mtxA(i, k) = mtxA(i, nIs(k))
mtxA(i, nIs(k)) = p
Next i
End If
Next k
' 求解成功
MRinv = True
End Function
來源:
1采用二維數組存放矩陣各元素值
2根據數矩陣運算規則進行相應計算
求矩陣的逆常用的有如下三種做法。經濟數學團隊幫你解答,請及時采納。謝謝!
一、公式法:A的逆陣=(1/|A|)A*,其中A*是A的伴隨陣。
二、初等變換法:對分塊矩陣(A,E)做行初等變換,前半部分A化成單位陣E時,后半部分E就化成了A的逆陣。
三、猜測法:如果能通過已知條件得出AB=E或BA=E,則B就是A的逆矩陣。