测量平差程序设计7.pptVIP

测量平差程序设计 Private Sub equation(B() As Double, p As Double, l As Double) 组法方程通用过程 n = 2 * (zds - yds) For i = 1 To n 按未知数循环 If Abs(B(i)) 0 Then 第i个未知数的系数不为0 UX(i) = UX(i) + p * B(i) * l 组法方程常数项 测量平差程序设计 For j = 1 To i If Abs(B(j)) 0 Then h = (i - 1) * i / 2 + j 计算xi、xj的互乘系数在一维按列上三角存储法方程系数阵中的序号 NX(h) = NX(h) + B(j) * B(i) * p 组法方程系数阵 End If Next j End If Next i End Sub 测量平差程序设计 Case 3 高斯约化法解法方程 n = 2 * (zds - yds) ReDim nb(n), nc(n) 数组清零 For i = 1 To n nb(i) = -UX(i) 将法方程常数项赋予一工作数组,保留原值用于求pvv Next i 测量平差程序设计 For k = 1 To n - 1 从法方程的第一行到倒数第二行, k实际上控制约化次数 m1 = k * (k + 1) / 2 ‘m1是k行的自乘元素在法方程系数阵中序号（一维按列上三角阵） For j = k + 1 To n 在j循环内完成k次约化 m2 = (j - 1) * j / 2 + k ‘m2在j循环内依次是k行自乘系数右边各元素序号 nb(j) = nb(j) - nb(k) * NX(m2) / NX(m1)”约化常数项 测量平差程序设计 For i = j To n i循环内完成法方程系数阵j行元素的k次约化 d = (i - 1) * i / 2 n1 = d + k 与约化元素同列的 k行元素 n2 = d + j 约化元素在法方程系数阵中编号 NX(n2) = NX(n2) - NX(n1) * NX(m2) / NX(m1) Next I Next j Next k 测量平差程序设计 For k = n To 1 Step -1 在k循环内回代求解未知数 u = nb(k) If k n Then For i = k + 1 To n nub = (i - 1) * i / 2 + k u = u - NX(nub) * nb(i) Next i nb(k) = u End If nub = k * (k + 1) / 2 nb(k) = nb(k) / NX(nub) Next k 测量平差程序设计 pvv2 = ll For k = 1 To n 按[pvv]=[pll]+W*X求[[pvv] pvv2 = pvv2 + UX(k) * nb(k) Next k 测量平差程序设计 For i = 1 To zds 求坐标平差值 cha = charact(i, k) 自定义函数，查点i是否已知点，如不是,用k返回i前面有几个已知点 If cha = n Then 序号i的点不是已知点 h = i - k d = 2 * (h - 1) + 1 计算i点x坐标未知数在未知数点集中的序号 x(i) = x(i) + nb(d) / 100: y(i) = y(i) + nb(d + 1) / 100 End If Next i 测量平差程序设计 ma = Sqr(pvv1 / (nl(cds) + ns(cds) - n - cds)) 求单位权中误差 Call inversion(NX(), nb()) 将约化后的法方程系数阵送入通用过程inversion()求逆 ReDim nb(n) 工作数组清零 测量平差程序设计 For i = 1 To zds 按点号循环,求点位中误差 cha = charact(i, k) 自定义函数，查点i是否已知点，如不是,用k返回i前面有几个已知点 If cha = n Then 不是已知点 n1 = 2 * (i - k - 1) + 1 计算i点x坐标未知数在未知数点集中的序号 n2 = n1 * (n1 + 1) / 2 i点x坐标未知数协因数在协因数阵中序号 n3 = (n1 + 1) * (n1 + 2) / 2 i点y坐标未知数协因数在协因数阵中序号