附合导线严密平差算法复习总结.docVIP

编制：测试科/黄QN 附合导线严密平差算法总结 图1 如图的单一附合导线，有4个已知点A、B、C、D，2个未知点TP1、TP2。设观测边数为n, 则未知点数为n-1, 观测角数为n+1。以上图为例，n = 3。观测边为：S1 = B-TP1, S2 = TP1-TP2, S3 = TP2-C 思路： 由于A、B坐标已知，则可以算出起始方位角，再根据B点坐标和每个观测角（夹角，左角）推算出TP1、TP2、C点的近似坐标值。如果是用全站仪进行测量，则用盘左盘右重复观测求平均的方式，直接测出TP1、TP2、C点的近似坐标值以及CD的方位角。再根据c点的已知坐标与近似坐标求坐标闭合差，由CD的已知方位角和近似方位角求角度闭合差，两个闭合差联立求得边长和角度的改正数，最后求得未知点的坐标平差值。 条件平差过程： 建立条件方程，求得条件系数 求法方程系数 求权阵 计算出联系数K 解算出观测值改正数V 由观测值和改正数计算平差值 详细步骤如下： 建立条件方程 在单一附合导线中，只需要三个条件方程即：方位角附合条件、纵坐标附合条件和横坐标附合条件方程。（省略了条件方程的推导过程，详细过程请查看参考资料：《测量平差.pps》） （1）方位角附合条件 [Vβi]n+11 + Wβ = 0 式中， Wβ = - （T0 – TCD+ [βi]n+11 - （n+1）*180°）（角度闭合差） βi —— 角度观测值（夹角，左角） Vβi —— 各观测角的改正数。 如果是用全站仪观测，则 Wβ = - （TCD – TCD） 式中，TCD ——CD的方位角观测值， TCD ——CD的已知方位角 （2）纵坐标X附合条件方程 [Cos TI * VSi]1n - （1/ ρ”） * [(Yn+1 - Yi) * Vβi]1n - Wx = 0; 式中，TI —— 各方位角观测值（近似值） VSi —— 边长改正数 Yn+1 —— C点即终止点的横坐标Y的观测值（近似值） Yi —— 待定点的横坐标Y的观测值 Wx = - （Xn+1 - XC） XC —— C点即终止点的纵坐标X的已知值 ρ” = 2062.65 （3）横坐标Y附合条件方程 [Sin TI * VSi]1n + （1/ ρ”） * [(Xn+1 - Xi) * Vβi]1n – WY = 0; 式中，TI —— 各方位角观测值（近似值） VSi —— 边长改正数 Xn+1 —— C点即终止点的纵坐标X的观测值（近似值） Xi —— 待定点的纵坐标X的观测值 WX = - （Yn+1 - YC） YC —— C点即终止点的横坐标Y的已知值 ρ” = 2062.65 求条件方程的系数矩阵 联立3个方程得改正数条件方程组： [Vβi]n+11 + Wβ = 0 [Cos TI * VSi]1n - 1/ ρ” * [(Yn+1 - Yi) * Vβi]1n - Wx = 0; [Sin TI * VSi]1n + 1/ ρ” * [(Xn+1 - Xi) * Vβi]1n – WY = 0; 其系数矩阵arrA为：（即改正数V的系数，此处以图1为例， n = 3） VSi  (边长改正数) Vβi (角度改正数) 0 0 0 1 1 1 1 CosT1 CosT2 CosT3 （-1/ ρ” ）* (Y4 – Y1) （-1/ ρ”） * (Y4 – Y2) （-1/ ρ”） * (Y4 – Y2) 0 SinT1 SinT2 SinT3 （-1/ ρ”）* (X4 – X1) （-1/ ρ”） * (X4 – X1) （-1/ ρ”） * (X4 – X1) 0 联系数法方程（简称法方程） AP-1ATK – W = 0 A——系数矩阵arrA K —— 乘系数 P —— 权阵 W —— 闭合差矩阵 由上得法方程的系数阵N： N = AP-1AT （权的推导见参考资料：《全站仪观测导线测量平差方法的研究.pdf》） 角度权： Pβi = 1；（因为角度的标称精度是固定的，各观测角权值相等） 边长权： PSi = (μ0 *μ0 ) / （MD * MD）（误差比例系数固定，边长的误差与距离有关，因此不一致） 式中，μ0 ——先验测角中误差，以秒为单位 MD —— 距离观测中误差 若 导线边长为Si （米），e1 为仪器的边长标称固定误差（mm），e2为仪器的边长比例误差系数（无单位），则 MD = ± （e1 + e2 * Si * 0.001