commsub.bas

三角形闭合差计算程序

Attribute VB_Name = "CommSub"
Public Sub Transe(PotN As String, PCode As Integer, PNum As Integer)
    Dim I As Integer
    Dim arecord1 As Recordset
    Set arecord1 = g_d_Base.OpenRecordset("PointNameTable", dbOpenTable)
        With arecord1
              If .RecordCount > 0 Then
                 .MoveFirst
                 For I = 1 To .RecordCount
                    g_PN(I) = .Fields(1)
                    If I < .RecordCount Then
                        .MoveNext
                     End If
                 Next I
              End If
         End With
    arecord1.Close
    For I = 1 To PNum
          If PotN = g_PN(I) Then
                PCode = I
                GoTo L
          End If
     Next I
L:
End Sub

Public Sub GW(Ang As Double, kg As Integer)
        Dim t1 As Double
        Dim t2 As Double
        Dim T3 As Double
        Dim T4 As Double
        Dim T5 As Double
        Dim t As Double
        If Abs(Ang) < 1E-24 Then
            t = 1
        Else
             t = Ang / Abs(Ang)
        End If
        t1 = Abs(Ang)
        t2 = Int(t1 + 0.0000001)
        T3 = (t1 - t2) * 100#
        T4 = Int(T3 + 0.0000001)
        If kg = 1 Then
               T5 = (t2 + T4 / 60# + (T3 - T4) / 36#) * g_Pi / 180#
        Else
               T5 = (t2 + T4 / 100# + (T3 - T4) / 100#) * g_Pi / 200#
        End If
        Ang = T5 * t
End Sub

Public Sub WG(Ang As Double, kg As Integer)
        Dim t1 As Double
        Dim t2 As Double
        Dim T3 As Double
        Dim T4 As Double
        Dim T5 As Double
        Dim t As Double
        If Abs(Ang) < 1E-23 Then
            t = 1
        Else
             t = Ang / Abs(Ang)
        End If
        If kg = 1 Then
            t1 = Abs(Ang) * 180# / g_Pi
            t2 = Int(t1 + 0.0000001)
            T3 = (t1 - t2) * 60#
            T4 = Int(T3 + 0.0000001)
            T5 = t2 + T4 / 100# + (T3 - T4) * 0.006
        Else
            t1 = Abs(Ang) * 200# / g_Pi
            t2 = Int(t1)
            T3 = (t1 - t2) * 100
            T4 = Int(T3)
            T5 = t2 + T4 / 100# + (T3 - T4) * 0.01
        End If
        Ang = T5 * t
 End Sub

Public Sub ALFA(x As Double, y As Double, A As Double)
      Dim R As Double
      Dim t As Double
        If y < 0.0000000000001 Then
             t = 1
        Else
              t = y / Abs(y)
        End If
        If Abs(x) < 0.00000000001 Then
              R = t * g_Pi / 2#
        Else
              R = Atn(y / x)
        End If
        If x >= 0# Then
            If y >= 0# Then
                A = R
            Else
                A = R + 2# * g_Pi
            End If
        Else
            A = R + g_Pi
        End If
End Sub
        
Public Function ACos(x As Double) As Double
     Dim t As Double
     Dim t1 As Double
     Dim t2 As Long
     Dim T3 As Double
     Dim T4 As Double
     Dim T5 As Double
     t1 = 1#
     t2 = 3#
     T3 = -t1 / (t1 + 1#)
     T4 = x ^ t2 / t2
     t = T3 * T4
     T5 = g_Pi / 2# - x + t
     Do While Abs(t) > 0.000000001
        t1 = t1 + 2#
        t2 = t2 + 2
        T3 = T3 * t1 / (t1 + 1#)
        T4 = x ^ t2 / t2
        t = T3 * T4
        T5 = T5 + t
     Loop
     ACos = T5
End Function