第04章 间接平差

1、第04章 间接平差 zqz99 zqz99 第四章12 3 4 间接平差 1 间接平差原理2 误差方程式 3 间接平差的法方程 4 间接平差的精度评定 5 5 间接平差实例 zqz99 教 学 目 的 理解间接平差的基本思想 理解间接平差的基本原理 把握间接平差的方法和精度评定 能用间接平差进行水准网平差 能用间接平差进行平面掌握网平差 zqz99 两点之间最短的距离不肯定是直线在人与人的相处以及做事情的过程中，我们有时很 难直截了当就把事情做好。我们有时需要等待，有时需 要合作，有时需要技巧。当我们遇到困难时，我们先要学会分析自己有没有这个力量去克服，假如临时还没有， 那我们不肯定要硬挺、硬

2、冲，我们可以选择绕过困难， 绕过障碍，这并不是躲避，更不是轻意放弃，而是换一 条路连续前行，或许这样，一切会变得更顺当。“通往 广场的路不止一条”告知我们的也是这个道理。 zqz99 1 1 间接平差原理1.间接平差的基本思想 在测量工作中，确定某些量或某个图形所需要的最少观测 个数 t ,称为必要观测数。 例如，为了确定三角形的外形，我们观测了其中的两个内 角，则三角形的外形可由这两个内角唯一确定，不需要进行平 差，该问题必要观测数 t =2。 但为了能准时发觉粗差并提高测量成果的精度，通常对三 个内角都进行观测。由于观测误差的存在，致使观测值之间产 生不符值，即三角形内角和不等于180,三

3、角形的外形不能唯 一确定，需要采纳测量平差的方法来进行处理，消退冲突，获 取最优结果，最终唯一确定三角形的外形。 zqz99 那么，当观测结果存在多余观测的状况时，如何 应用测量平差方法解决因多余观测而产生有冲突呢?在上一章中，我们介绍了条件平差方法 它是以 n 个观测值的平差值作为未知数，通过它们之间存在的r 个条件方程来消退观测值之间的不符值，同时运用 求条件极值的原理来解出唯一的那组改正数，从而求 得各观测量的最或是值 本章将介绍测量平差的另一基本方法间接平差 法。 zqz99 下面就来争论间接平差法怎样解决问题 如图4-1所示的单三角形，三角形的外形可由任意两个 内角唯一确定，我们不妨

4、选取两个内角(如1和2)作为未 知参数，这两个角的大小一经确定，这个三角形的外形便唯一确定了，第三个角 3=180-1- 2,也就是说， 其它各量都可以表示成这 t 个独立未知参数的函数。 zqz99 设同精度观测三个内角分别为l1,l2,l3,其改正数为 l v ,l l v ,l l v ,选择 v1,v2,v3,平差值 l 1 1 1 2 2 2 3 3 3 t =2个量作为未知参数(未知数),这里就以1 2的平差值 和x 可惟一确定该三角形的外形。 、x 为未知数，由 x x 1 2 1 2 而每一个观测值的平差值都可表示为两个未知数的函数，即 l1 v1 x 1 l v x2 2 (

5、4-1) 2 x 180 l3 v3 x 1 2 式(4-1)为观测量的平差值与未知参数间的函数关系，称平差 值方程。将l1,l2,l3移至等式右端，可得误差方程 l v1 x 1 1 l v x2 2 2 (4-2) x 180 l v3 x 1 2 3 zqz99 式(4-2)中未知数 、x x 1 2 的数值一经确定，与其相应 的一组改正数便随之确定。对任意的一组未知量，都可 求得一组对应的改正数，且不论未知数取何值，只要是 由式(4-2)求得的改正数，都满意下式(l1+v1)+(l2+v2)+(l3+v3)=180 这说明由式(4-2)求得的改正数，都可达到消退不符 值的目的。但由式(

6、4-2)所求得的改正数随未知数的取值不同而不同，有无穷多组解。而我们盼望最终得到是这 无穷多组解中是最优的那一组解。那么怎样来猎取这样 最优的那一组解呢?测量平差中，就是依据最小二乘法 原理来猎取最优解的。 zqz99 依据最小二乘法原理的要求，在这无穷多组解中应选取其中能使v tpv=min(最小)的那一组改正数，它是最优的。 考虑到观测值等精度，观测值的权为1,权阵为单位阵， vtpv=min可表达为2 2 2 v tv v1 v2 v3 min 将式(4-2)代入，得v tv=(x1-l1)2 + (x2-l2)2 + (-x1-x2 + 180-l3)2 =min 由于未知数x1、x2

7、之间不存在关系，是自由变量，要使 l 上式为最小，可按数学中自由极值的求解方法，使 v1 x 1 1 t v v v2 x 2 l2 0 (4 2) x 1 180 v3 x 1 t x l 2 3 v v x 2 0 zqz99 得法方程 2x1 + x2 -l1 + l3-180=0 x1 + 2x2-l2 + l3-180=0 上式称为法方程，解算可得未知数x1、x2的最优估量值 x1=l1-(l1 + l2 + l3-180)/3x2=l2-(l1 + l2 + l3-180)/3 将x1、x2代入式(4-2)求得改正数 v1= -(l1 + l2 + l3-180)/3 v2= -(

8、l1 + l2 + l3-180)/3 l v1 x 1 1 l v x2 2 2 v3= -(l1 + l2 + l3-180)/3 x 180 l v3 x 1 2 3 (4 2) 由改正数便可求出观测值的平差值。 zqz99 综上所述，间接平差法思想：针对详细的平差问题，选定 t 个未知量，建立观测值的平差值与未知数间 的函数关系(列平差值方程),进而转化为误差方程(将 观测值移到等号的右边),并依据最小二乘法原理，按 求自由极值的方法解算出未知量的最优估值。 间接平差法是以最小二乘为平差原则，以平差值 方程、误差方程作为函数模型的平差方法。 下面介绍间接平差法的一般公式。 zqz99

9、2.间接平差的基本原理设某平差问题中，有n个独立观测值l1,l2,ln,其相应的权为p1,p2,pn,权阵为 0 p1 0 0 p 0 2 p 0 0 p n ,x , 必要观测数为 t,选定 t 个未知参数，分别以 x 1 2未知数的函数，即列出 n 个平差值方程 ,x t 表示 t 个未知参数的最优估值，将每一个观测值的平差值表达为 b x t x d l1 v1 a1 x 1 1 2 1 t 1 b x t x d l2 v2 a2 x 2 2 2 2 t 2 b x t x d ln vn an x 1 n 2 n t n (4-3) zqz99 将式中观测值移至等式右端，并令 li

10、= li di (i=1,2,n) 得到误差方程的一般形式为 b x v1 a1 x 1 1 2 b x v a x l t1 x t 1 t 2 x t l2 2 2 2 2 2 b x t x l v n an x 1 n 2 n t n (4-4) (4-5) 令 v1 v v 2 n 1 vn v v1 n 1 v2 vn t zqz99 l1 l l 2 n 1 ln x 1 x 2 x t 1 x n n 1 l l1 l2 ln t x x 1 t 1 x 2 x n t a1 a b 2 n t an b1 b2 bn t1 t2 tn zqz99 则可得误差方程的矩阵形式 l

11、 v bxn 1 n t t 1 n 1 (4-6) 满意式(4-6)的解有无穷多组，为求解惟一一组最优估 必需满意 vtpv=min 值，依据最小二乘法原理，上式的 x的要求，由于 t 个参数独立，可按数学中求自由极值的方 法，使 v t pv v t 2v p 2v t pb 0 x x转置后得 btpv = 0 (4-7) 满意上式 t 个方程的改正数便能使 vtpv 取得微小值。 zqz99 l v bx bt pv 0 (4 6) (4 7) 以上所得的式(4-6)和式(4-7)中，待求量的个数是 n 个改正数v 加 t 个未知参数，方程的个数为n+ t,有惟一 的解。满意这两个方程

12、的解，可消退观测值之间的不符值，同时又是符合最小二乘法原理的最优估值，故称式 (4-6)和式(4-7)为间接平差的基础方程。 解算这组基础方程，通常是将式(4-6)代入式(4-7), 以便先消去v,得 bt pl 0 bt pbx (4-8) 令nbb=btpb,w=btpl,并代入式 (4-8),得： w 0 (4-9) n xbb zqz99 由式(4-8)或式(4-9)可解算出未知数，称该方程组为 间接平差的法方程。解之得t 1 n 1 w x bb 或 t 1 -1 t x (bt pb) b pl 代入误差方程可得改正数v 将x l v bxn 1 n t t 1 n 1 观测值的平

13、差值为 l v l 当权 阵p为对角阵，即观测值之间相互独立时，将误差方程的系数阵b、常数阵 l 及权阵p代入法方程(4-8)中， 绽开后，可得法方程的纯量形式为 zqz99 pab x paa x 1 2 pbb x pab x1 2 pbt x pat x 1 2 pal 0 pat x t pbl 0 pbt x t (4-11) ptt x t ptl 0 ,再将未 实际计算时，先由式(4-9)解算未知数 x 知数代入误差方程式(4-6),计算观测值的改正数v, ,这 由改正数v与观测值求和即得观测值的平差值 l 些平差值之间已经消退了冲突。至此，解决了间接平 差中求最或是值的问题。 zqz99 3.间接平差求平差值的步骤及示例 依据上述原理，实际间接平差法的计算步骤如下: (1) 依据平差问题的性质，确定必要观测个数 t , 选定 t 个独立量作为未知参数； (2) 将每一个观测量的平差值分别表达成所选未知参数的函数，即平差值方程，并列出误差方程； (3) 由误差方程的系数 b 与自由项 l 组成法