作者silent1017 (......)
看板Visual_Basic
标题[VB6 ]如何呼叫使用模组
时间Fri Aug 24 15:51:46 2007
如题
我都是被这些变数给搞混了！！
不知有没有人使用过这个模组
知道的高手说一下
感激不尽！！
我把我的模组PO上来
Public Function InverseLU(ByRef A() As Double, _
ByRef Pivots() As Long, _
ByVal N As Long) As Boolean
Dim Result As Boolean
Dim WORK() As Double
Dim I As Long
Dim IWS As Long
Dim J As Long
Dim JB As Long
Dim JJ As Long
Dim JP As Long
Dim JP1 As Long
Dim V As Double
Dim i_ As Long
Result = True
'
' Quick return if possible
'
If N = 0# Then
InverseLU = Result
Exit Function
End If
ReDim WORK(1# To N)
'
' Form inv(U)
'
If Not InvTriangular(A, N, True, False) Then
Result = False
InverseLU = Result
Exit Function
End If
'
' Solve the equation inv(A)*L = inv(U) for inv(A).
'
For J = N To 1# Step -1
'
' Copy current column of L to WORK and replace with zeros.
'
For I = J + 1# To N Step 1
WORK(I) = A(I, J)
A(I, J) = 0#
Next I
'
' Compute current column of inv(A).
'
If J < N Then
JP1 = J + 1#
For I = 1# To N Step 1
V = 0#
For i_ = JP1 To N Step 1
V = V + A(I, i_) * WORK(i_)
Next i_
A(I, J) = A(I, J) - V
Next I
End If
Next J
'
' Apply column interchanges.
'
For J = N - 1# To 1# Step -1
JP = Pivots(J)
If JP <> J Then
For i_ = 1# To N Step 1
WORK(i_) = A(i_, J)
Next i_
For i_ = 1# To N Step 1
A(i_, J) = A(i_, JP)
Next i_
For i_ = 1# To N Step 1
A(i_, JP) = WORK(i_)
Next i_
End If
Next J
InverseLU = Result
End Function
''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
'Inversion of a general matrix.
'
'Input parameters:
' A - matrix. Array whose indexes range within [1..N, 1..N].
' N - size of matrix A.
'
'Output parameters:
' A - inverse of matrix A.
' Array whose indexes range within [1..N, 1..N].
'
'Result:
' True, if the matrix is not singular.
' False, if the matrix is singular.
'
' -- ALGLIB --
' Copyright 2005 by Bochkanov Sergey
'
''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
Public Function Inverse(ByRef A() As Double, ByVal N As Long) As Boolean
Dim Result As Boolean
Dim Pivots() As Long
Call LUDecomposition(A, N, N, Pivots)
Result = InverseLU(A, Pivots, N)
Inverse = Result
End Function
我有一个3*3的矩阵，要将它做反(逆)矩阵（inv）
但是我不知怎麽去呼叫它
我的矩阵如下
1 5 6
7 4 2
-3 6 7
做完反（逆）矩阵如下
0.2462 0.0154 -0.2154
-0.8462 0.3846 0.6154
0.8308 -0.3231 -0.4769
这是用matlab做出来的，但是现在要用vb写...就一直卡在这里
请各位高手帮帮忙
感激不尽！！
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 59.113.139.87
※ 编辑: silent1017 来自: 59.113.139.87 (08/24 16:21)
※ 编辑: silent1017 来自: 61.223.244.68 (08/25 20:47)
※ 编辑: silent1017 来自: 61.223.232.178 (08/26 15:13)