【转】秩亏自由网平差

  

秩亏自由网平差

    在前面介绍的经典平差中，都是以已知的起算数据为基础，将控制网固定在已知数据上。如水准网必须至少已知网中某一点的高程，平面网至少要已知一点的坐标、一条边的边长和一条边的方位角。当网中没有必要的起算数据时，我们称其为自由网，本节将介绍网中没有起算数据时的平差方法，即自由网平差。

在经典间接平差中，网中具备必要的起算数据，误差方程为

                                                          （8-2-1）

式中系数阵为列满秩矩阵，其秩为 。在最小二乘准则下得到的法方程为

                                                             （8-2-2）

由于其系数阵的秩为，所以为满秩矩阵，即为非奇异阵，具有凯利逆，因此具有唯一解，即

                                                               （8-2-3）

当网中无起算数据时，网中所有点均为待定点，设未知参数的个数为u，误差方程为

                              （8-2-4）

式中

                               

为必要的起算数据个数。尽管增加了个参数，但的秩仍为必要观测个数，即

其中为不满秩矩阵，称为秩亏阵，其秩亏数为。

组成法方程

                                                              （8-2-5）

式中 ，且，所以也为秩亏阵，秩亏数为：   

                                  （8-2-6）

由上式知，不同类型控制网的秩亏数就是经典平差时必要的起算数据的个数。即有：

在控制网秩亏的情况下，法方程有解但不唯一。也就是说仅满足最小二乘准则，仍无法求得的唯一解，这就是秩亏网平差与经典平差的根本区别。为求得唯一解，还必须增加新的约束条件，来达到求唯一解的目的。秩亏自由网平差就是在满足最小二乘和最小范数的条件下，求参数一组最佳估值的平差方法。

下面将推导自由网平差常用两种解法的有关计算公式。

 一、直接解法

根据广义逆理论，相容方程组虽然具有无穷多组解，但它有唯一的最小范数解，即：

                                                           （8-2-7）

式中，称为矩阵的最小范数g逆。称为矩阵的g逆。代入（8-2-7）式得