Excel VB如何应用在附合导线内业计算中

Excel VB如何应用在附合导线内业计算中摘要：Excel是大家经常使用的办公软件，相信对大家来说并不陌生。

如果你是从事测量工作，你是否感觉到Excel有很不方便的时候？比如，当你要计算一个方位角的三角函数值，而角度的单位是60进制的，此时，你一定感到很无奈，因为，Excel本身无法直接计算60进制的角度的三角函数；还有，在附合导线内业计算中要计算角度的三角函数值以及坐标增量、以传统的°′″形式显示观测角值和方位角值，要做到这些较专业的工作，Excel电子表的一些常用功能是无法办到的。

当你遇到这种情况，你只有拐弯抹角了，通过配合使用计算器一步一步地计算，很麻烦！同时也失去了Excel电子表的强大的、高效的计算能力。

但你只要对Excel进行一些挖掘，就可以发现Excel的功能我们还没有充分的利用。

Excel本身提供了强大的二次开发功能，只要我们仔细的研究，就可以解决工作中的一些难题。

下面，是本人在测量工作中积累的一些工作经验和方法，希望能和大家一起交流。

关键词：Excel VB 导线计算
初识VB语言，首先，你必须懂得一些简单的VB编程常识。

如果不懂就只有通过其他的途径去学习了。

但用不着深入的研究，只要静下心来，几个小时就可以了。

一、60进制的角度的三角函数计算
打开Excel，按Alt+F11或从菜单栏中的（工具->宏->Visual Basic编辑器）进入VB编程状态，学过VB的人一看就知道那就是熟悉的VB 界面。

下面看看如何定义一个函数，然后利用它来解决60进制的角度
的三角函数计算问题。

在菜单上依次点击[插入]->[模块]，然后输入如下代码：
Public Const pi = 3.14159265359
Public Function DEG(n As Double)
Dim A As Double, B As Double, C As Double, D As Double, E As Double, F As Double, G As Double, KA As Double
D = Abs(n) + 0.000000000000001
F = Sgn(n)
A = Int(D)
B = Int((D - A) * 100)
C =
D - A - B / 100
DEG = F * (A + B / 60 + C / 0.36) * pi / 180
End Function
这样，就定义了一个名字叫DEG的函数，它的作用就是转换60进制的角度为Excel认识的弧度。

编辑完后按Alt+Q即返回Excel。

当在Excel中某一单元格输入=sin(deg(A1))(A1既可以是单元格的角度值,也可以是输入的角度值),回车后结果便出来了(或者在编辑栏中点击fx弹出函数选择框，并选择用户定义函数中的DEG())。

你可以用计算器检验一下计算结果是否正确。

如果出现#NAME？那就要设置一下安全设置。

依次点[工具]->[宏]->[安全性]，在安全级选项卡上选择“中”或者“低”，然后关闭后重新打开就可以了，以后只要是60进制的角度，就用它转换，非常方便快捷了。

二、以度分秒的形式显示弧度值
在Excel电子表中，进入VB宏编辑状态，在菜单上依次点击[插入]->[模块]，然后输入如下代码：
Public Function RAD(Nu As Double) As Double
Dim a As Double, b As Double, c As Double, d As Double, E As Double, F As Double, G As Double, p As Double
d = Abs(Nu)
F = Sgn(Nu)
p = 180# / pi
G = p * 60#
a = Int(d * p)
b = Int((d - a / p) * G)
W = b
c = (
d - a / p - b / G) * 20.62648062
RAD = (c + a + b / 100) * F
End Function
这样，便定义了一个名字叫RAD的函数，它的作用就是转换Excel 所认识的弧度值为我们所熟悉的60进制的角度值。

编辑完后按Alt+Q 即返回Excel。

在Excel中某一单元格输入=RAD(deg(A1))(A1既可以是单元格的角度值,也可以是输入的角度值),回车，这样就完成了以度分秒（°′″）的形式计算的角度值。

你可以用计算器检验一下计算结果是否正确。

三、在Excel电子表中计算两点间的方位角
以同样的方法进入VB宏编辑状态，输入以下代码：
Public Function FWJ(a As Double, b As Double, c As Double, d As Double) E = c - a: F = d - b
If (E > 0 And F > 0) Then FWJ = RAD((Atn(F / E)))
If (E < 0 And F > 0) Then FWJ = RAD((pi + Atn(F / E)))
If (E < 0 And F < 0) Then FWJ = RAD((pi + Atn(F / E)))
If (E > 0 And F < 0) Then FWJ = RAD((2 * pi + Atn(F / E)))
End Function
这样，也完成了一个名叫FWJ的函数的编辑，它的作用就是计算两点间的方位角，并以度分秒的形式显示在Excel电子表的单元格中。

编辑完后按Alt+Q即返回Excel。

在Excel中某一单元格输入=FWJ(a,b,c,d)或点击编辑栏的fx弹出函数参数输入框，在完成参数的输入后回车。

这样就计算出了两点间的方位角。

（其中：第一点坐标为（a,b）,第二点坐标为（c,d）。

输入顺序依次a->b->c->d）
四、下面是一段导线内业简易平差计算例子（见附表1—1）
其中，点号T102、T103、T145、T112为已知的控制点。

计算时只需输入点号、转角（度分秒形式，357.3708即：357°37′08″）、距离、测站数n以及已知点的坐标，剩余的计算工作由Excel通过拖、拉或复制自动完成。

在打印时第B列可隐藏。

计算表中数值的计算：
第B列的值=DEG (第C列的值)，第B列最后一行为观测角的总和（弧度）：B30=SUM(B6：B28)。

②第C列转角（°′″）：测站数C31=12，C32=RAD(B30)；
C33=C32-C31*180，其它单元格为输入的观测角的平均值。

③第D列中转角改正数D6~D28=(－$D$33/$C$31)，$D$33、$C$31为绝对引用，拖、拉或复制不会引起单元格数值的改变。

D33=RAD(DEG(C33)-DEG(E33))。

④第E列中方位角E5=FWJ(K4,L4,K6,L6)、E29=FWJ(K28,L28,K30,L30) 即方位角T102->T103和T145->T112；
E7=RAD(DEG(E5)+DEG(C6)+DEG(D6)-DEG(180)),E9=(RAD(DEG(E7) +DEG(C8)+DEG(D8)-DEG(180))-360),以此类推并通过拖拉来实现Excel 的自动计算。

⑤第F列中F7~F27为距离观测值的平均值，导线总长
F33=SUM(F7:F27)。

⑥第G列中X坐标增量G7=F7*COS( DEG(E7))，即G列的值= F列的值*COS( DEG(E列的值))。

坐标增量之和G33=SUM(G7:G27)。

⑦第H列中坐标增量改正数H7=-F7/$F$33*$H$33，再拖拉操作计算H9到H27的值。

⑧第i列中Y坐标增量和第j列中改正数的计算类似于第G、H列。

⑨第K、L列为点的最终坐标值，K8=K6+G7+H7、L8=L6+I7+J7，再拖拉操作计算得出K9~K26、L9~L26的坐标。

K33=K28-K6、L33=L28-L6为理论坐标增量。

X坐标增量误差fx：H33=G33-K33、Y坐标增量误差fy：J33=G33-K33。

⑩相对闭合差K的计算：F34=SQRT(H33*H33+J33*J33)
->H34=INT((F33/F34)) ->K=1/ H34；角度闭合差fβ=D33。

以上所述是测量工作过程中经常遇到的计算问题，涉及角度、弧度、方位角等较为专业的领域，通过Excel电子表并结合前面提到的DEG
函数、RAD函数、FWJ函数就方便、高效地解决了。

Excel的功能是非常强大的，如果有兴趣，你还可以在AutoCAD中直接与Excel通讯，比如一条三维多段线的所有结点的三维坐标直接导入到Excel，比在AutoCAD中用列表的方法要方便的多，同样，在Excel中利用曲线要素表数据来计算公路路线的逐桩坐标或逐桩高程也极为有
效的方法。

限于篇幅，在此不作详细叙述了。

如果读者有兴趣，可以深入的学习和探讨。

附表1—1。