精确度是精度和准确度的合成,是指观测结果与其真值的接近程度, 包括观测结果与其数学期望接近程度和数学期望与其真值的偏差。 精确度的衡量指标为均方误差:
MSE(X)=E(X−X¯)

测量不确定度

协方差传播律及权

数学期望的传播

协方差传播律

观测值线性函数的方差

设有观测值Xn1,其数学期望为μXn1,协方差阵为DXXnn,即
X=X1μ2…μn,μX=μ1μ2…μn=E(X1)E(X2)…E(Xn)=E(X)
DXX=E[(X−μX)(X−μX)T]=σ21σ21…σn1σ12σ22σn2………σ1nσ2n…σ2n(3-2-2)
其中 σi 为Xi 的方差 , σij 为Xi与Xj的协方差 , 又设有 X 的线性函数为
Zt1=K1nXn1+k011
式中:
K1n=[k1,k2,…,kn]
公式(3-2-2)的纯量形式为
Z=k1X1+k2X2+⋯+knXn+k0
则:
Z 的数学期望:
E(Z)=KμX+k0
Z 的方差:
DZZ11=σ2Z=KDXXKT(3-2-4)
将(3-2-4)写成纯量形式:
DZZ11==σ2Z=k21σ21+k22σ22+⋯+k2nσ2n+2k1k2σ12+2k1k3σ13+⋯+2k1knσ1n+…+2kn−1knσn−1,n
即协方差传播率
### 多个观测值线性函数的协方差阵
若令
Zt1=Z1Z2…Zt,Ktn=k11k21…kt1k12k22kt2………k1nk2n…ktn,K0t1=k10k20…kt0
Zt1=KtnXn1+K0t1
纯量展开式为:
Z1=k11X1+k12X2+⋯+k1nXn+k10Z2=k22X2+k22X2+⋯+k2nXn+k20………Zt=kttXt+kt2X2+⋯+ktnXn+kt0
Z 的数学期望为
E(Z)=KμX+K0
Z 的协方差阵
DZZtt=KtnDXXnnKTnt
### 非线性函数的协方差阵
设有 Xn1的非线性函数:
Z=f(X)(也写为Z=f(X1,X2,…,Xn))
则Z的方程 DZZ为:
DZZ=KDXXKT
其中:
K=[k1k2…kn]=[(∂f∂X1)0(∂f∂X2)0…(∂f∂Xn)0]

t个非线性函数的协方差阵

若有t个非线性函数
Z1=Z2=…Zt=f1(X1,X2,…,Xn)f2(X1,X2,…,Xn)ft(X1,X2,…,Xn)⎫⎭⎬⎪⎪⎪⎪⎪⎪⎪⎪
则Zt1的协方差阵
DZZ=KDXXKT
其中:
Ktn=⎡⎣⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢(∂f1∂X1)0(∂f2∂X1)0⋮(∂ft∂X1)0(∂f1∂X2)0(∂f2∂X2)0(∂ft∂X2)0………(∂f1∂Xn)0(∂f2∂Xn)0⋮(∂ft∂Xn)0⎤⎦⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥

协方差传播律的应用

权与定权的常用方法

协因数和协因数传播律

协因数与协因数阵

观测值的权与它的方差成反比, 设有观测值L_{i}和L_{j},它们的方差分别为σi和σj,它们之间的协方差为σij,令
Qii=1pi=σ2iσ20
Qjj=1pj=σ2jσ20
Qij=σijσ20
称Qii 和Qjj分别称为 Li和Lj的协因数或权倒数,称Qij为Li关于 Lj的协因数或相关权倒数。
假定有观测值向量(或者观测值向量函数) Xn1和Yr1,他们的方差阵分别为 DXX和DYY,X关于Y的互互协方差阵为 DXY,令
QXXnn=1σ20DXX
QYYrr=1σ20DYY
QXYnr=1σ20DXY
则称QXX和QYY分别为X和Y的协因数阵,QXY为X关于Y的互协因数阵(相关权逆阵)。

协因数传播律

设有观测值 X, 已知它的协因数阵为QXX, 又设有 X 的函数 Y 和Z:
Y=FX+F0
Z=KX+K0
按协方差传播律导出协因数阵为:
QYY=FQXXFT
QZZ=KQXXKT
QYZ=FQXXKT
这就是观测值的协因数阵与其线性函数的协因数阵的关系式,通常称之为协因数传播律,或称之为权逆阵传播律。
由真误差计算中误差及其实际应用

系统误差的传播

前几节所讨论的问题,是以观测值只含有偶然误差为前提的.
本节讨论观测值中同时含有系统误差时,观测值综合误差的方差(简称综合方差) 估计及系统误差的传播规律。

观测值的系统误差与综合方差

设有观测值 L,观测量的真值为Ln1˜,则L的综合误差可以定义为:
Ω=L˜−L
如果综合误差 Ω中只包含有偶然误差Δ, 由偶然误差的特性可知其数学期望应为 E(Ω)= E(Δ) = 0, 如果 Ω中除包含偶然误差Δ外, 还包含有系统误差 ε,即
Ω=Δ+ϵ=L˜−L
此时,由于系统误差 ε不是随机变量, 所以 Ω的数学期望为
E(Ω)=E(Δ)+ϵ=ϵ≠0
又因为
ϵ=E(Ω)=E(L˜−L)=L˜−E(L)
- ε也就是观测值上的数学期望对于观测量真值的偏差值,
- L 包含的系统误差愈小(ε愈小),则L的数学期望对于真值的偏差值愈小
当观测值 L 中既存在偶然误差Δ,又存在系统误差时, 其观测值的综合方差 DLL是用均方误差表示的
DLL=MSE(L)=E(L−L˜)=E(Ω2)=σ2+ε2
观测值的综合方差等于它的方差 σ2与系统误差的平方 ε2之和

系统误差的传播

线性函数
设已知观测值Li ( i = 1,2, ⋯, n) 的系统误差为
ϵi=E(Ωi)=Li˜−E(Li)(i=1,2,…,n)
其中,Li˜和Ωi是Li所对应的观测量的真值和综合误差.
设有线性函数
Z=k1L1+k2L2+⋯+knLn+k0
Z的综合误差ΩZ 与各个Li的综合误差Ωi之间的关系式为
ΩZ=k1Ω1+k2Ω2+⋯+knΩn
线性函数的系统误差的传播公式:
εZ=E(ΩZ)=p=∑i=1nkiεi
非线性函数
函数 Z 是非线性形式:
Z=f(L1,L2,…,Ln)
可以用它们的微分关系代替它们的误差之间的关系,即有
ΩZ=∂Z∂L1Ω1+∂Z∂L2Ω2+⋯+∂Z∂LnΩn

系统误差与偶然误差联合传播

系统误差为常熟、常系差的情况
对于线性函数
Z=k1L1+k2L2+⋯+knLn
与综合误差之间的关系为
ΩZ=k1Ω1+k2Ω1+⋯+knΩn
Z 的综合方差为:
DZZ=∑i=1nk2iσ2i+∑i=1nk2iε2i
随机系统误差情况
设观测值 L的综合误差是
Ω=Δ+ϵ
ϵ为随机性系统误差,一般与Δ相互独立,若已知Δ的偶然方差为σ2,ϵ的系统方差为σ2ϵ , 则可按协方差传播律,得
σ2L=σ2Ω=σ2+σ2ϵ
有观测值线性函数
Z=k1L1+k2L2+⋯+knLn
与综合误差之间的关系为
ΩZ=k1Ω1+k2Ω1+⋯+knΩn
Z 的综合方差为:
DZZ=∑i=1nk2iσ2i+DZZ=∑i=1nk2iσ2ϵi
平差数学模型与最小二乘原理

测量平差概述

函数模型

条件平差–条件方程

简接平差–误差方程

附有参数的条件平差

附有限制条件的间接平差

函数模型的线性化

测量平差的数学模型

参数估计与最小二乘原理

函数模型 (一般)函数模型 (线性)数学模型
条件平差F(L˜)=0ArnL˜n1+A0r1=0r1或者AΔ+W=0
简接平差L˜n1=F(X˜)L˜n1=BntX˜t1+dn1或者l+Δ=BX
附有参数的条件平差Ft1(L˜X˜)=0AtnL˜n1+BtuX˜u1+A0t1=0或者AΔ+BX˜+W=0
附有限制条件的间接平差L˜u1=F(X˜u1)Φs1(X˜)=0⎫⎭⎬⎪⎪L˜n1=BnuX˜u1+dn1或l+Δ=BX˜CsuX˜u1+Wxs1=0s1⎫⎭⎬⎪⎪

条件平差

条件平差原理

条件方程

精度评定

条件平差公式汇编和水准网平差示例

附有参数的条件平差

附有参数的条件平差原理

精度评定

公式汇编和示例

间接平差

间接平差原理

误差方程

精度评定

间接平差公式汇编和水准网平差示例

间接平差特例— — —直接平差

三角网坐标平差

测边网坐标平差

导线网间接平差

GPS 网平差

附有限制条件的间接平差

附有限制条件的间接平差原理

精度评定

公式汇编和示例

概括平差函数模型

基本平差方法的概括函数模型

附有限制条件的条件平差原理

精度评定

各种平差方法的共性与特性

平差结果的统计性质

误差椭圆

概述

点位误差

误差曲线

误差椭圆

相对误差椭圆

点位落入误差椭圆内的概率

平差系统的统计假设检验

统计假设检验概述

统计假设检验的基本方法

误差分布的假设检验

平差模型正确性的统计检验

平差参数的统计检验和区间估计

粗差检验的数据探测法

近代平差概论

序贯平差

附加系统参数的平差

秩亏自由网平差

最小二乘配置原理