测量平差程序设计6.ppt

测量平差程序设计 Public Function angle(l As Double) As Double ‘弧度转角度函数过程 l = l * 180 / pi b1 = Int(l) b3 = (l - b1) * 0.6 b2 = Int(b3 * 100) l = b1 + b2 / 100 + (b3 - b2 / 100) * 0.6 angle = l End Function 测量平差程序设计 Private Function side(n As Integer, m As Integer) 反算已知点间边长函数过程 For i = 1 To yds If ym(i) = dm(n) Then k1 = i If ym(i) = dm(m) Then k2 = i Next i dx = xo(k1) - xo(k2): dy = yo(k1) - yo(k2) side = Sqr(dx ^ 2 + dy ^ 2) End Function 测量平差程序设计 Private Function azimuth(x1 As Double, y1 As Double, x2 As Double, y2 As Double) 反算坐标方位角 dx = x2 - x1: dy = y2 - y1 azimuth = Atn(dy / dx) If dx 0 Then azimuth = azimuth + pi Else If dy 0 Then azimuth = azimuth + 2 * pi End If End Function 测量平差程序设计 Private Function seqn(str As String) As Integer 由点名查计算序号函数 For i = 1 To zds If str = dm(i) Then seqn = i 将查到的序号赋给函数名,返回调用处 Next i End Function 测量平差程序设计 Private Function sid(n As Integer, m As Integer) As Double 提取边长函数，n是测站号，m是方向号 ss = 0 If ns(n) - ns(n - 1) - 1 = 0 Then 该测站有观测边，确定m方向的边长 测量平差程序设计 For k = ns(n - 1) + 1 To ns(n) If sb(k) = lb(m) Then ss = s(k) 依次将边的照准点名与m方向照准点名对比，相等则赋予ss。 Next k End If If ns(n) - ns(n - 1) - 1 0 Or ss 1 Then 如测站n未找到m方向的边，则转到lb(m)站去找 测量平差程序设计 h = seqn(lb(m)) 由点名查计算编号 If ns(h) - ns(h - 1) - 1 = 0 Then For k = ns(h - 1) + 1 To ns(h) If sb(k) = dm(n) Then ss = s(k) 依次将k边的照准点名与n测站点名对比，相等则赋予ss。 Next k End If End If sid = ss 将查到的边长赋给函数名,返回调用处 End Function 测量平差程序设计 Private Function charact(n1 As Integer, n2 As Integer) As String 根据序号判断是否已知点,小于该序号的点中有几个已知点 n2 = 0: charact = n For k = 1 To yds 查那些点是已知点 If n1 = m(k) Then m(k)存的是已知点序号，此判断检查序号为i的点是否为已知点 charact = y 是已知点 测量平差程序设计 Else If n1 m(k) Then n2 = k n1不是已知点,用n2记下该序号前有几个已知点,m()数组已经排过序 End If Next k End Function 测量平差程序设计 Private Sub order(m() As Integer, n As Integer) 排序通用过程 For i = 1 To n - 1 For j = i + 1 To n If m(j) m(i) Then A = m(j): m(j) = m(i): m(i) = A End If Next j Next i End Sub‘由于作为测站出现的顺序，与作为已