zoukankan      html  css  js  c++  java
  • 【转】方差分量估计

    §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.依据定权公式再次定权,再次平差,如此反复,直到各类单位权方差的估值相等或接近相等为止。

  • 相关阅读:
    解决windows上安装TortoiseSVN后不能使用命令行问题
    Python里Pure paths、PurePosixPath、PureWindowsPath的区别
    PHP数组运算符
    global,local,static的区别
    echo和print的区别
    PHP中foreach循环传值问题
    Matlab入门学习(文件读写)
    Matlab入门学习(程序设计)
    IDEA中使用Maven下载依赖时报错:unable to find valid certification path to requested target
    全国县市区编码表
  • 原文地址:https://www.cnblogs.com/xiongyunqi/p/3735933.html
Copyright © 2011-2022 走看看