【转】方差分量估计

§8-4 方差分量估计

2学时

我们知道,平差前观测值向量的方差阵一般是未知的,因此平差时随机模型都是使用观测值向量的权阵。而权的确定往往都是采用经验定权,也称为随机模型的验前估计,对于同类观测值可按第一章介绍的常用定权方法定权;对于不同类的观测值,就很难合理地确定各类观测值的权。为了合理地确定不同类观测值的权,可以根据验前估计权进行预平差,用平差后得到的观测值改正数来估计观测值的方差,根据方差的估计值重新进行定权,以改善第一次平差时权的初始值,再依据重新确定的观测值的权再次进行平差,如此重复,直到不同类观测值的权趋于合理,这种平差方法称为验后方差分量估计。此概念最早由赫尔默特(F.R.Helmert)在1924年提出,所以又称为赫尔默特方差分量估计。

 

一、赫尔默特方差分量估计公式

 

为推导公式简便起见,设观测值由两类不同的观测量组成,不同类观测值之间认为互不相关,按间接平差时的数学模型为

                                   (函数模型)              8-4-1

                   (随机模型)               8-4-2

其误差方程为

                                 权阵                       8-4-3

                            权阵                      8-4-4

作整体平差时,法方程为

                                                           8-4-5

式中

                     

                 

一般情况下,由于第一次给定的权是不恰当的,或者说它们对应的单位权方差是不相等的,设为,则有

                                                       8-4-6

但只有才认为定权合理。方差分量估计的目的就是根据事先初定的权进行预平差,然后利用平差后两类观测值的来求估计量,再根据(8-4-6)式求出,由这个方差估值再重新定权,再平差,直到为止。为此需要建立与估计量之间的关系式。

由数理统计知识可知,若有服从任一分布的q维随机变量,已知其数学期望为,方差阵为,则向量的任一二次型的数学期望可以表达为:

                       8-4-7

式中为任意q阶的对称可逆阵。

现用向量代替上式中的向量,则其中的应换为应换为阵可以换成权阵,于是有

                                 8-4-8

 前面已经证明,于是有:

                                                   8-4-9

                    

                      

对上式应用协因数传播律得

             

代入上式,整理后得

       

将上式代入(8-4-9)式,得

      

顾及矩阵迹的性质,上式可写为

同理可得

    

去掉上面两式的期望符号,相应的单位权方差也改用估值符号表示,整理顺序后得

                                                                   8-4-10

                                                                   8-4-11

其矩阵形式可写为

                                                             8-4-12

                                                               8-4-13

式中

 

            

            

    8-4-12)、(8-4-13)两式即为赫尔默特方差分量估计的严密公式。由此式可以求得两类观测值的单位权方差估值,从而可以根据(8-4-6)式求得观测值方差的估值,以此方差估值再次定权,再次平差,直至满足要求为止。

现将以上推导扩展至m组观测值。误差方程为

                       

                       

                   

                   

则得参数的估值为

                             

按照上述类似的推导,则有

去掉期望符号,相应的单位权方差也改为用估值符号,则有

                                                              8-4-14

式中

                     

                     

二、计算步骤

 

1.将观测值分类,并进行验前权估计,即确定各类观测值的权的初值

2.进行第一次平差,求得

3.按(8-4-14)式求各类观测值单位权方差估值

4.按(8-4-6)式计算各类观测值方差的估值;

5.依据定权公式再次定权,再次平差,如此反复,直到各类单位权方差的估值相等或接近相等为止。

原文地址:https://www.cnblogs.com/xiongyunqi/p/3735933.html