一些重要的通用程序.docVIP

求AB的坐标方位角，输入是两点坐标，输出的是弧度值 Public Function DirectAB(Xa#, Ya#, Xb#, Yb#) As Double Dim detX#, detY#, tana# detX = Xb - Xa detY = Yb - Ya If Abs(detX) 0.000001 Then If detY 0 Then DirectAB = PI / 2 Else DirectAB = PI * 3 / 2 End If Else tana = detY / detX DirectAB = Atn(tana) If detX 0 Then DirectAB = PI + DirectAB ElseIf detX 0 And detY 0 Then DirectAB = PI * 2 + DirectAB End If End If End Function 弧度化为度.分秒的形式：输入弧度值，输出度.分秒（各占两位） Public Function HuToDo(ByVal Hu As Double) As Single Dim du%, fen%, miao% Hu = Hu * 180 / PI du = Fix(Hu) Hu = (Hu - du) * 60 fen = Fix(Hu) Hu = (Hu - fen) * 60 miao = Fix(Hu + 0.5) If miao = 60 Then fen = fen + 1 miao = 0 End If HuToDo = du + fen / 100 + miao / 10000 End Function 将度.分秒形式化为弧度：输入为度.分秒形式，输出为弧度 Public Function DoToHu(ByVal DoFenMiao As Double) As Single Dim du%, fen%, miao%, angle# du = Fix(DoFenMiao) DoFenMiao = (DoFenMiao - du) * 100 fen = Fix(DoFenMiao) miao = (DoFenMiao - fen) * 100 angle = du + fen / 60 + miao / 3600 DoToHu = angle * PI / 180 End Function 通用的间接平差解算过程：输入系数矩阵A、权矩阵P、常数向量L和解向量X，求出X，并通过参数传出去 Public Sub InAdjust(B, P, L, X) …… MatrixTrans B, Bt ‘系数矩阵转置 Matrix_Multy BtP, Bt, P ‘求解出BTP Matrix_Multy Nbb, BtP,B 法方程系数矩阵 Matrix_Multy W, BtP, L 法方程常数向量 MajorInColGuass Nbb, W, X ‘求解平差值 End Sub 通用的条件平差解算过程：输入系数矩阵A、权矩阵P、常数向量L和解向量X，求出X，并通过参数传出去 Public Sub CondiAdjust(A, P, W, V) …… For i = 1 To p1 求Q矩阵 Q(i, i) = 1 / P(i, i) Next i MatrixTrans A, At 转置系数阵A Matrix_Multy QAt, Q, At 求解出P-1AT Matrix_Multy Naa, A, QAt 法方程系数矩阵 MajorInColGuass Naa, W, K 解法方程 Matrix_Multy V, QAt, K 求改正数 End Sub 条件平差和间接平差中的代码说明：MatrixTrans（