【转】秩亏自由网平差

秩亏自由网平差

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

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

                                         (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个待定点,所以附加阵为

则有

所以有

未知参数的的协因数阵为

结果与直接解法完全相同。

时间: 2024-11-03 03:28:47

【转】秩亏自由网平差的相关文章

【转】GPS网平差

进行GPS网平差的目的主要有三个: (1)消除由观测量和已知条件中存在的误差所引起的GPS网在几何上的不一致.包括闭合环闭合差不为0:复测基线较差不为0:通过由基线向量所形成的导线,将坐标由一个已知点传算到另一个已知点的符合差不为0等.通过网平差,可以消除这些不一致. (2)改善GPS网的质量,评定GPS网的精度.通过网平差,可得出一系列可用于评估GPS网的精度指标,如观测值改正数.观测值验后方差等等.结合这些精度指标,还可以设法确定出可能存在粗差或质量不佳的观测值,并对它们进行相应的处理,从而

空间谱专题11:子阵平滑与秩亏缺

作者:桂. 时间:2017-09-29  21:20:18 链接:http://www.cnblogs.com/xingshansi/p/7612984.html 未完待续... 前言 主要分析在解相干算法中,子阵平滑的有效性. 一.前向平滑 以均匀线阵(ULA)为例,第l个接收阵元的信号为: 其中,M为阵元数,N为信号个数. 以前向平滑为例: 令β = ,对于第k个子阵有 其中D为: 假设信号与噪声不相关,且噪声为白噪声,计算相关矩阵: 前向平滑修正的协方差矩阵: 证明:当子阵阵元数m>=N,

空间谱专题10:MUSIC算法

作者:桂. 时间:2017-09-20  07:43:12 链接:http://www.cnblogs.com/xingshansi/p/7553746.html 未完待续… 前言 MUSIC(Multiple Signal Classification)算法通常用来进行到达角(DOA,Direction of arrival)估计. 一.MUSIC原理简介 根据前文的分析,模型依然建立在窄带信号的基础上: X为接收阵元,F为入射信号,a为对应的导向矢量,W为噪声.可直接记作矩阵形式 通常借助相

【转】GPS基线的精化处理

影响基线解算结果的因素主要有: (1)基线解算时所设定的起点坐标不准确. 起点坐标不准确,会导致基线出现尺度和方向上的偏差,造成的影响目前还没有较容易的方法来加以判别,因此,在实际工作中,只有尽量提高起点坐标的准确度,以避免这种情况的发生. (2)少数卫星的观测时间太短,导致这些卫星的整周未知数无法准确确定. 当卫星的观测时间太短时,会导致与该颗卫星有关的整周未知数无法准确确定,而对于基线解算来讲,对于参与计算的卫星,如果与其相关的整周未知数没有准确确定,就将影响整个基线解算的结果. 对于卫星观

【转】GPS基线解算模式

GPS基线向量是利用2台或2台以上GPS接 收机所采集的同步观测数据形成的差分观测值,通过参数估计得方法所计算出的两两接收机间的三维坐标差.与常规地面测量中所测定的基线边长不同,基线向量是 既具有长度特性又具有方向特性的矢量,而基线边长则是仅具有长度特性的标量.基线向量主要采用空间直角坐标的坐标差的形式. 在一个基线解算结果中,可能包含很多项内容,但其中最主要的只有两项,即基线向量估值及其验后方差-协方差阵. 对于一组具有一个共同端点的同步观测基线来说,由于在进行基线解算时用到了一部分相同的观测

空间谱专题06:宽带信号处理思路

作者:桂. 时间:2017-09-09  20:04:22 链接:http://www.cnblogs.com/xingshansi/p/7413139.html 前言 目前分析的问题,仍然限定在布阵的环节,暂不涉及后处理及硬件实现. 一.宽带处理的一般方式 前面分析的阵列信号模型,都是建立在窄带信号的基础上,对于宽带阵列信号,有两个思路:1)信道化为窄带信号,按窄带的思路进行处理:2)利用聚焦变换,聚焦到特定频率,进行处理. 聚焦变换的思路对相干信号同样有效,这样便可以省去子阵平滑的步骤,需要

自适应滤波:矩阵求逆

作者:桂. 时间:2017-04-02  10:36:09 链接:http://www.cnblogs.com/xingshansi/p/6658655.html 声明:欢迎被转载,不过记得注明出处哦~  [读书笔记09] 前言 西蒙.赫金的<自适应滤波器原理>第四版第八章:最小二乘法.因为最小二乘涉及到矩阵求逆,因为通常对于秩缺矩阵其逆是不可求的,这就需要借助广义逆矩阵.而广义逆矩阵可以借助奇异值分解(SVD,Singularly Valuable Decomposition)进行求解. 有

空间谱专题03:相干信号的处理方式

作者:桂. 时间:2017-09-09  16:44:12 链接:http://www.cnblogs.com/xingshansi/p/7413133.html 主要理论空间平滑算法,首先分析相干信号特性. 〇.相干信号特性 以两个信号为例 求相关矩阵 如果两个信号的相关系数ρ满足: 1)ρ=0,则认为两信号不相关: 2)0<ρ<1,则认为两信号相关: 3)ρ = 1,则两信号相干. 当两信号相干时,ρ=1,对于相关矩阵: 秩为1,这就造成了秩亏,对于子空间等空间谱估计算法便不再适用. 也可

论文学习-sparse methods for direction of arrival estimation1.

翻译自Sparse Methods for Direction-of-Arrival Estimation(Zai Yang∗†, Jian Li‡, Petre Stoica§, and Lihua Xie†) direction of arrival(DOA) 1.引言 DOA(direction of arrival)estimation指接收一些电磁波的方向信息的过程,这些电磁波来自许多形成阵列传感器的接收雷达的输出. 传统的波束形成器(beamformer)仅仅使用了对空间采样数据的傅