可靠性度量
可靠性概念是从研究模型误差的角度而提出来的,在大量的工程实践中,人们往往倾向于注重平差系统的精度,而忽略了系统的可靠性。实际上,可靠性与精度同等重要,且两者并不等价。GPS精密单点定位中,高精度的解算结果并不代表其可靠性也好。根据Baarda (1968)提出的测量系统“可靠性”理论,可靠性指标是衡量一个平差系统是否具有良好的抗御模型误差品质的参数,它可分为内部可靠性和外部可靠性。
内部可靠性
内部可靠性是指在一定显著水平
α
α
α和检验功效
β
β
β的条件下,设计的检验统计量能够探测到的最小模型误差(Minimal Detectable Bias,MDB),它表征了系统发现模型误差的能力。根据Baarda的数据探测理论,一组观测值
L
i
L_i
Li的内部可靠性可由其边界值
ε
i
ε_i
εi来表征:
ε
i
=
k
i
σ
0
ε_i=k_i σ_0
εi=kiσ0
表示第
i
i
i个模型误差的边界值为
k
i
k_i
ki倍的单位权中误差,而
k
i
k_i
ki取决于非中心化参数
λ
λ
λ和平差系统的结构:
k
i
=
λ
/
(
c
i
T
P
i
Q
V
i
V
i
P
i
c
i
)
1
/
2
k_i={λ/(c_i^T P_i Q_{V_i V_i } P_i c_i)}^{1/2}
ki=λ/(ciTPiQViViPici)1/2
当P为对角阵,即观测不相关时,
R
=
P
Q
V
V
R=PQ_{VV}
R=PQVV为平差系统的多余观测分量,
r
i
r_i
ri为矩阵
R
R
R中对角线第i行i列元素。顾及
σ
i
=
σ
0
√
(
Q
i
i
)
=
σ
0
/
√
(
P
i
)
σ_i=σ_0 √(Q_ii )=σ_0/√(P_i )
σi=σ0√(Qii)=σ0/√(Pi),可进一步表示为:
ε
i
=
(
√
λ
σ
0
)
/
(
P
i
√
(
Q
V
i
V
i
)
)
=
(
√
λ
σ
i
)
/
√
(
P
i
Q
V
i
V
i
)
=
√
(
λ
/
r
i
)
σ
i
ε_i=(√λ σ_0)/(P_i √(Q_{V_i V_i } ))=(√λ σ_i)/√(P_i Q_{V_i V_i } )=√(λ/r_i ) σ_i
εi=(√λσ0)/(Pi√(QViVi))=(√λσi)/√(PiQViVi)=√(λ/ri)σi
上式表明,第i个观测量的可发现最小模型误差
ε
i
ε_i
εi与该观测中误差
σ
i
σ_i
σi成正比,与
√
(
r
i
)
√(r_i )
√(ri)成反比。当观测精度越高或多余观测分量
R
i
R_i
Ri越大,系统可发现模型误差的能力越强,可靠性越好。因此,要提高系统的可靠性就要从提高观测精度和增加多余观测分量这两方面着手。
外部可靠性
外部可靠性是指未被检测到的且与MDB数值相同的模型误差对参数(或部分参数)估计的影响(Mininal Detectable Effect,MDE)。
ε
i
ε_i
εi对参数x的影响可表示为:
∇
x
i
=
(
A
T
Q
y
−
1
A
)
−
1
A
T
Q
y
−
1
c
y
i
ε
i
∇x_i=(A^T Q_y^{-1} A)^{-1} A^T Q_y^{-1} c_{yi} ε_i
∇xi=(ATQy−1A)−1ATQy−1cyiεi
式中,KaTeX parse error: Expected 'EOF', got '̂' at position 4: ∇x ̲̂_i是一个有量纲的指标,该向量包含较多的数值。对于精密单点定位而言,它包含了三维坐标、接收机钟差、对流层天顶湿延迟、模糊度等参数,难以对可靠性进行解释。因此,在此基础上构造一个无量纲指标,即偏差-噪声率(Bias to Noise Ratio,BNR)来表征外部可靠性:
G
(
x
i
)
=
√
(
(
∇
x
i
)
T
Q
x
−
1
(
∇
x
i
)
)
G_(x_i)=√((∇x_i )^T Q_x^{-1} (∇x_i))
G(xi)=√((∇xi)TQx−1(∇xi))
该指标描述了所有参数的总体可靠性,其数值越小,表明模型误差对参数估值的影响越小,即外部可靠性越强。
可以证明:
G
(
x
i
)
=
√
(
(
∇
x
i
)
T
Q
x
−
1
(
∇
x
i
)
)
=
√
(
λ
(
(
1
−
r
i
)
/
r
i
)
)
G_(x_i)=√((∇x_i )^T Q_x^{-1} (∇x_i))=√(λ_ ((1-r_i)/r_i ))
G(xi)=√((∇xi)TQx−1(∇xi))=√(λ((1−ri)/ri))
实际应用中,我们可能只关心未被探测的模型误差对参数的某个线性函数的影响,例如,对于精密单点定位,普通用户更加关注未被模型化的误差对位置参数(点位误差)的影响,因此,可将
∇
x
i
∇x_i
∇xi向量分为位置参数和其他参数(如接收机钟差、对流层天顶湿延迟、模糊度等)两部分,即
σ
F
=
σ
0
√
(
Q
F
F
)
=
σ
0
√
(
f
T
Q
x
x
f
)
σ_F=σ_0 √(Q_{FF} )=σ_0 √(f^T Q_{xx} f)
σF=σ0√(QFF)=σ0√(fTQxxf)