秩亏自由网平差
在前面介绍的经典平差中,都是以已知的起算数据为基础,将控制网固定在已知数据上。如水准网必须至少已知网中某一点的高程,平面网至少要已知一点的坐标、一条边的边长和一条边的方位角。当网中没有必要的起算数据时,我们称其为自由网,本节将介绍网中没有起算数据时的平差方法,即自由网平差。
在经典间接平差中,网中具备必要的起算数据,误差方程为
(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)式得
(8-2-8)
上式就是根据广义逆理论直接求解参数的唯一最小范数解的公式。由于广义逆计算较为复杂,下面将公式做进一步改化:
令
(8-2-9)
(8-2-10)
式中行满秩,即,于是有
(8-2-11)
而,所以为满秩方阵,按照降阶法求矩阵广义逆的方法,即:如果有矩阵
其中存在凯利逆,则有的g逆
(8-2-12)
根据上式可得
(8-2-13)
代入(8-2-8)式,得
(8-2-14)
或写成
(8-2-15)
未知参数的协因数阵为:
(8-2-16)
二、附加条件法(伪观测值法)
前面已提及,秩亏自由网平差就是在满足最小二乘和最小范数的条件下,求参数一组最佳估值的平差方法,实际上就是求相容方程组的最小范数解。附加条件法的基本思想:由于网中没有起算数据,平差时多选了d个未知参数,因此在u个参数之间必定满足d个附加条件式,即在原平差函数模型中需要加入d个未知参数间的限制条件方程,从而可以按附有条件的间接平差法求解。问题的关键是如何导出等价于的限制条件方程的具体形式。
为叙述方便,我们先给出该限制条件方程,然后再推导平差计算公式,最后证明,在给定的限制条件方程下所求得的解,就是相容方程组的最小范数解。
设等价于约束条件的限制条件方程为
(8-2-17)
式中且满足称为附加阵。故秩亏自由网平差的函数模型为
权阵为
按照附有条件的间接平差可得法方程
(8-2-18)
式中,且,唯一不同的是这里为秩亏阵。
为解决秩亏问题,将(8-2-18)中的第二式左乘矩阵后,再加到第一组中得:
(8-2-19)
式中,且
根据附有条件的间接平差原理,上式的解为
(8-2-20)
(8-2-21)
由于上述解是通过增加未知参数间满足的d个附加条件,按照附有条件的间接平差法而实现的,因此人们把此法称为附加条件法。但它又不同于经典的附有条件的间接平差法,其主要表现为:当阵满足时,必定有下式成立(证明从略)
(8-2-22)
将(8-2-22)式代入(8-2-21)式,可得参数的解为
(8-2-23)
现在只需证明,按(8-2-23)式求得的解就是法方程的最小范数解。为此只需证明是的最小范数g逆中的一个即可,即只需证明满足以下两式:
(8-2-24)
现证明如下:因为 ,所以有
右乘阵并展开,则有
而,所以有
(8-2-25)
由于,存在逆阵,则有
(8-2-26)
所以有
(8-2-27)
(8-2-28)
因此(8-2-24)第一式得到验证。
由(8-2-27)式得
考虑到(8-2-26)式,则上式为
(8-2-29)
(8-2-28)、(8-2-29)两式说明是的最小范数g逆中的一个,因此按(8-2-23)式求得的一定是相容方程组的最小范数解。
三、精度评定
单位权中误差估值的计算
(8-2-30)
式中可以直接计算,也可以按下式求得
(8-2-31)
未知参数的协因数阵为
(8-2-32)
实际计算时,通常要对进行标准化,设标准化后的阵用表示,即不仅要求满足,还要求满足,此时(8-2-26)式变成,转置后有,因此(8-2-32)式将变成如下形式
(8-2-33)
四、两点说明
①若将代入法方程,则法方程变为
上式相当于下列误差方程联合组成的法方程
上式的第一式为观测值的误差方程,第二式可以看作是为求最小范数解而人为增设的d个虚拟误差方程,因此附加条件法又叫伪观测值法。
②该方法的特点就是用求凯利逆替代了求广义逆,因此便于计算和计算机编程,但首要条件是必须知道附加阵,关于附加阵的确定问题,本教材不准备作详细讨论,下面直接给出常见控制网的附加阵及其标准化后的矩阵的具体形式:
水准网(设有u个点)
; (8-2-34)
测边网(设有m个点)
(8-2-35)
式中为第I点的近似坐标
(8-2-36)
式中是以中心坐标为原点的第I点的近似坐标,它们的计算如下:
,
测角网(设有m个点)
只需在(8-2-35)式中增加一行元素、在(8-2-36)式中增加一行元素即可得到相应的阵和阵。
例[8-3] 如图8-2水准网,点全为待定点,同精度独立高差观测值为,,,平差时选取三个待定点的高程平差值为未知参数,并取近似值
试分别用直接法和附加条件法求解参数的平差值及其协因数阵。
解:1.直接解法
误差方程为
法方程为
由法方程易知
, ,
所以有
未知参数的改正数为
未知参数的平差值为
未知参数的协因数阵为
2.附加条件法
解法一中已求得法方程为的具体形式为:
该水准网有3个待定点,所以附加阵为
则有
所以有
未知参数的的协因数阵为
结果与直接解法完全相同。