lecpxgaussmainmodule.bas

编程中经常用到的一些算法,算法已模块化,来自于<科学工程与数值算法大全>

Attribute VB_Name = "MainModule"
Sub Main()
    Dim dblAR(4, 4) As Double, dblAI(4, 4) As Double
    Dim dblBR(4) As Double, dblBI(4) As Double
    Dim sCpx1 As String, sCpx2 As String, sCpx3 As String, sCpx4 As String
    
    '系数矩阵
    dblAR(1, 1) = 1
    dblAR(1, 2) = 3
    dblAR(1, 3) = 2
    dblAR(1, 4) = 13
    dblAR(2, 1) = 7
    dblAR(2, 2) = 2
    dblAR(2, 3) = 1
    dblAR(2, 4) = -2
    dblAR(3, 1) = 9
    dblAR(3, 2) = 15
    dblAR(3, 3) = 3
    dblAR(3, 4) = -2
    dblAR(4, 1) = -2
    dblAR(4, 2) = -2
    dblAR(4, 3) = 11
    dblAR(4, 4) = 5
    
    dblAI(1, 1) = 3
    dblAI(1, 2) = -2
    dblAI(1, 3) = 1
    dblAI(1, 4) = 6
    dblAI(2, 1) = -2
    dblAI(2, 2) = 7
    dblAI(2, 3) = 5
    dblAI(2, 4) = 8
    dblAI(3, 1) = 9
    dblAI(3, 2) = -3
    dblAI(3, 3) = 15
    dblAI(3, 4) = 1
    dblAI(4, 1) = -2
    dblAI(4, 2) = -2
    dblAI(4, 3) = 7
    dblAI(4, 4) = 6
    '常数向量
    dblBR(1) = 2
    dblBR(2) = 7
    dblBR(3) = 3
    dblBR(4) = 9
    
    dblBI(1) = 1
    dblBI(2) = 2
    dblBI(3) = -2
    dblBI(4) = 3
    
    '求解
    If LECpxGauss(4, dblAR, dblAI, dblBR, dblBI) = True Then
        sCpx1 = dblBR(1) & " + " & dblBI(1) & "j"
        sCpx2 = dblBR(2) & " + " & dblBI(2) & "j"
        sCpx3 = dblBR(3) & " + " & dblBI(3) & "j"
        sCpx4 = dblBR(4) & " + " & dblBI(4) & "j"
        MsgBox "求解成功！" & Chr$(13) & Chr$(13) & "x1 = " & sCpx1 & Chr$(13) & "x2 = " & sCpx2 & Chr$(13) & "x3 = " & sCpx3 & Chr$(13) & "x4 = " & sCpx4
    Else
        MsgBox "求解失败！"
    End If
    
End Sub