【正文】 rDim n1 As IntegerDim m2 As IntegerDim n2 As IntegerDim mystring() As StringDim l() As Stringmystring() = Split(, vbCrLf)m1 = UBound(mystring, 1) LBound(mystring, 1) + 1 39。第一行數(shù)據(jù)n1 = UBound(l, 1) LBound(l, 1) + 1 39。重新定義矩陣AFor i = 1 To m1 l() = Split(mystring(i 1), ) n1 = UBound(l, 1) LBound(l, 1) + 1 39。獲取矩陣A Next jNext imystring() = Split(, vbCrLf)m2 = UBound(mystring, 1) LBound(mystring, 1) + 1 39。第一行數(shù)據(jù)n2 = UBound(l, 1) LBound(l, 1) + 1 39。重新定義矩陣BFor i = 1 To m2 l() = Split(mystring(i 1), ) n2 = UBound(l, 1) LBound(l, 1) + 1 39。獲取矩陣B Next jNext iReDim C(m1, n1) 39。調(diào)用矩陣相減的函數(shù) = For i = 1 To m1 For j = 1 To n1 = + Str(C(i, j)) + Next j = + vbCrLfNext iEnd Sub39。獲取行數(shù)l() = Split(mystring(0), ) 39。獲取列數(shù)ReDim A(m1, n1) 39。獲取列數(shù) For j = 1 To n1 A(i, j) = l(j 1) 39。獲取行數(shù)l() = Split(mystring(0), ) 39。獲取列數(shù)ReDim B(m2, n2) 39。獲取列數(shù) For j = 1 To n2 B(i, j) = l(j 1) 39。重新定義矩陣CCall mmul(A, B, C) 39。矩陣轉(zhuǎn)置提取數(shù)據(jù)Private Sub Cmd_zhuanzhi_Click()Dim n1 As IntegerDim m1 As IntegerDim n2 As IntegerDim m2 As IntegerDim AT() As DoubleDim mystring() As StringDim l() As Stringmystring() = Split(, vbCrLf)m1 = UBound(mystring, 1) LBound(mystring, 1) + 1 39。第一行數(shù)據(jù)n1 = UBound(l, 1) LBound(l, 1) + 1 39。重新定義矩陣AFor i = 1 To m1 l() = Split(mystring(i 1), ) n1 = UBound(l, 1) LBound(l, 1) + 1 39。獲取矩陣A Next jNext iReDim AT(n1, m1) 39。調(diào)用矩陣轉(zhuǎn)置的函數(shù) = For i = 1 To n1 For j = 1 To m1 = + Str(AT(i, j)) + Next j = + vbCrLf Next iEnd Sub39。獲取行數(shù)l() = Split(mystring(0), ) 39。獲取列數(shù)ReDim A(m1, n1) 39。獲取列數(shù) For j = 1 To n1 A(i, j) = l(j 1) Next jNext iCall MRinv(A) 39。Cmd_qc_Click()事件完成數(shù)據(jù)的清除Private Sub Cmd_qc_Click() = = = End SubSub madd(mtxA() As Double, mtxB() As Double, mtxC() As Double) 39。矩陣相減 Dim i As Integer, j As Integer Dim m As Integer, n As Integer m = UBound(mtxA, 1) LBound(mtxA, 1) + 1 n = UBound(mtxA, 2) LBound(mtxA, 2) + 1 For i = 1 To m For j = 1 To n mtxC(i, j) = mtxA(i, j) mtxB(i, j) Next j Next iEnd SubSub mmul(mtxA() As Double, mtxB() As Double, mtxC() As Double) 39。矩陣轉(zhuǎn)置 Dim i As Integer, j As Integer Dim m As Integer, n As Integer m = UBound(mtxA, 1) LBound(mtxA, 1) + 1 n = UBound(mtxA, 2) LBound(mtxA, 2) + 1 For i = 1 To n For j = 1 To m mtxAT(i, j) = mtxA(j, i) Next j Next iEnd SubFunction MRinv(mtxA() As Double) As Boolean 39。利用角度與弧度的相互轉(zhuǎn)換函數(shù)，根據(jù)坐標正算公式XB=XA+DAB*Cos(αAB)，YB=YA+DAB* Sin(αAB)計算出B點坐標。Form_Load()Private Sub Form_Load()opt_fwj(0) = TrueEnd SubPrivate Sub Cmd_js_Click()Dim dab, αab, xa, ya, xb, yb, dw%xa = Val(txt_xa)ya = Val(txt_ya)dab = Val(txt_dab)αab = Val(txt_αab)dw = 0 39。輸入角度單位為度dw = 1ElseIf opt_fwj(2) Then 39。調(diào)用角度與弧度相互轉(zhuǎn)換函數(shù)xb = xa + dab * Cos(αab) 39。利用公式計算B點的縱坐標xb = Round(xb, 3) 39。對yb進行四舍五入，小數(shù)點后保留3位txt_xb = Trim(Str(xb)) 39。去掉字符串兩邊的空格End SubPrivate Sub Cmd_tc_Click()Unload frm_zbzsEnd Sub39。定義常量piDim d%, f%, m, fh%fh = Sgn(jd) 39。取輸入角度的絕對值39。d = Int(jd) 39。提取分m = ((jd d) * 100 f) * 100 39。得到角度總和，單位為度Case 1 39。輸入單位為弧度jdzh = jd * 180 / piEnd Select39。輸出單位轉(zhuǎn)換為弧度jdzh = jdzh * pi / 180 * fhCase 1 39。輸出單位為分jdzh = jdzh * 60 * fhCase 3 39。d = Int(jdzh) 39。提取分m = ((jdzh d) * 60 f) * 60 39。利用函數(shù)jddw（）將弧度轉(zhuǎn)換為不同的角度單位及表現(xiàn)形式，根據(jù)兩點之間的距離公式求出水平距離，再運用坐標反算公式計算出坐標方位角。自定義函數(shù)jddw（）的功能是將弧度轉(zhuǎn)換為不同的角度單位及表現(xiàn)形式Function jddw(hd, n%) As VariantConst pi = Dim jd, d%, f%, m, fh%If n = 3 Then 39。jddw = Format((d + f / 100 + m / 10000) * fh, )ElseIf n = 1 Then 39。分’秒”jddw = Format(d, 0) amp。 amp。 ′ amp。 ″ElseIf n = 2 Then 39。Form_Load（）事件完成組合框（Cmb_fwj）項目的加載Private Sub Form_Load() 度176。Cmd_js_Click()事件完成坐標反算Private Sub Cmd_js_Click()Const pi = Dim xa, ya, xb, yb, dab, aabDim dx, dyxa = Val()ya = Val()xb = Val()yb = Val()dx = xb xady = yb yadab = Sqr(dx ^ 2 + dy ^ 2) 39。求方位角 = Format(dab, ) = jddw(aab, )End Sub39。Cmd_tc_Click()事件退出界面Private Sub Cmd_tc_Click()Unload Frm_zbfsEnd Sub（4）運行結(jié)果實驗四 前方交會（1） 原理前方交會分角度前方交會和距離前方交會兩種。其中，角度交會按式（41）進行計算；距離交會先按式（42）計算出α和β角度值，再按式（41）計算出P點