以平面三角形 为例,由于测量误差的存在,对三个角度进行多次测量后,得到的各角度观测值内角和并不等于180°,此时需要进行平差以求得各角度的最佳估计值
在测量学 中,测量平差 ,也称为平差 或最小二乘平差 (英语:Least-squares adjustment ),是指依据某类最优化 准则对带有观测误差 的测量数据进行调整,以求得测量对象的最佳估计值 的理论和方法。测量平差的问题来源于测量过程中的多余观测。[ 1] [ 2] 受到测量误差 的影响,通过多余观测得到的测量值必然无法精确满足测量对象之间应存在的数学关系,这些测量值之间的不一致性被称为不符值或闭合差。[ 3] [ 4] 包含多余观测的观测值在数学上组成了一个无精确解的超定系统 ,但根据所选取的最优化准则,可以从该系统中求得一个符合该准则的近似解。[ 5] [ 6]
测量平差的基本任务即是处理观测值之间的不符值,并依据选取的最优化准则,求得观测量的最佳估计值,并对其精度 进行评定。[ 7] 由于经典的测量平差方法通常选取最小二乘准则 作为最优化准则,这类经典的平差方法也被称为最小二乘平差 。[ 8] 依据经典测量平差求得的估计值是测量对象的最优线性无偏估计 。[ 9] [ 10]
平差的函数模型是描述观测量与未知量之间的数学关系的模型。这些数学关系即可是几何关系,也可是物理关系。例如,大地测量 中的测量控制网 和摄影测量 中的共线方程 描述的是几何关系,而在重力测量 、卫星定轨或是形变监测中使用的模型描述的则是物理关系。函数模型中的观测量在实际的观测过程中确定,而未知量可根据模型的需要进行选取,测量平差的目的即是对这些未知量做出最优估计。[ 1] [ 2]
函数模型的观测量指经测量得到的测量对象的测量值,可分为必要观测量和多余观测量。
必要观测量是唯一确定函数模型所需要的观测量,仅与函数模型有关,必要观测量之间相互独立 。[ 11] 以平面三角形为例,唯一确定其形状 需要两个角度观测量 ,必要观测数为2;而唯一确定其形状和大小 则需要两个角度及一条边长 、一个角度及两条边长,或是三条边长的观测量,必要观测数为3。
平差模型中的总观测数
n
{\displaystyle n}
与必要观测数
t
{\displaystyle t}
之间必须满足
t
≥
n
{\displaystyle t\geq n}
,才能保证该模型满秩 ,平差模型有解。[ 12] [ 13]
多余观测量则是在必要观测以外进行额外观测得到的观测量,其数目
r
{\displaystyle r}
亦称为自由度或冗余度(英语:Redundancy )[ 8] ,由总观测数
n
{\displaystyle n}
和必要观测数
t
{\displaystyle t}
确定:
r
=
n
−
t
{\displaystyle r=n-t}
由多余观测可列出与其数目相对应的条件方程,从而检核出条件方程中的不符值或闭合差,进而对测量对象的理论值和测量精度进行估计。[ 14]
函数模型的未知量通常特指模型的待求参数。其与观测量的区别在于,这些未知量既可以是被直接观测的对象(例如上述的角度和边长),亦可以是对观测量产生影响的因素(例如几何图形中点的坐标)。将观测量表达成未知量的函数的方程即为观测方程。函数模型中未知量(参数)的数目
u
{\displaystyle u}
与必要观测量的数目
t
{\displaystyle t}
之间的数量关系与函数模型的类型一一对应:[ 1] [ 15]
当
u
=
0
{\displaystyle u=0}
时,函数模型中只存在表达观测量间约束关系的条件方程,此时的平差方法为条件平差;
当
0
<
u
<
t
{\displaystyle 0<u<t}
时,除了表达观测量之间约束关系的条件方程外,还需在函数模型中增设
u
{\displaystyle u}
条对参数作出约束的条件方程,此时的平差方法为附有参数的条件平差;
当
u
=
t
{\displaystyle u=t}
时,函数模型中的所有观测量都可通过
u
{\displaystyle u}
个相互独立的参数所组成的观测方程进行表达,此时的平差方法为间接平差;
当
u
>
t
{\displaystyle u>t}
时,除了
u
{\displaystyle u}
个相互独立的参数,函数模型中还存在
s
=
u
−
t
{\displaystyle s=u-t}
个可表达为其他参数的函数的多余参数,因此需增设
s
{\displaystyle s}
个条件方程对其进行约束,此时的平差方法为附有限制条件的间接平差。
测量平差中的条件方程是指描述各观测量之间,以及各观测量与参数之间应当满足的数学关系的方程式 ,其数学形式可统一为[ 1] [ 15]
F
(
L
^
,
X
^
)
=
0
{\displaystyle \operatorname {F} ({\hat {L}},{\hat {X}})=0}
式中,
L
^
{\displaystyle {\hat {L}}}
表示经测量得到的观测量的理论值,
X
^
{\displaystyle {\hat {X}}}
表示选取的参数的理论值。根据函数模型中选取参数的不同,条件方程有以下特殊形式:
当平差模型中没有另选参数时,条件方程的形式化为
F
(
L
^
)
=
0
{\displaystyle \operatorname {F} ({\hat {L}})=0}
,该条件方程是间接平差的函数模型;
当条件方程仅表示对参数的约束条件,方程当中没有观测量时,条件方程的形式化为
Φ
(
X
^
)
=
0
{\displaystyle \operatorname {\Phi } ({\hat {X}})=0}
,被特称为约束条件方程或限制条件方程,而其他形式的条件方程亦被统称为一般条件方程。
由两个已知点(
A
{\displaystyle A}
、
B
{\displaystyle B}
)和两个未知点(
C
{\displaystyle C}
、
D
{\displaystyle D}
)组成的测角三角网,共观测了
a
i
{\displaystyle a_{i}}
、
b
i
{\displaystyle b_{i}}
和
c
i
{\displaystyle c_{i}}
等九个角度
在大地测量 中,条件方程选取的依据通常是观测量之间的几何关系 。[ 16]
以右图的平面测角三角网 为例,其被测量的
a
i
{\displaystyle a_{i}}
、
b
i
{\displaystyle b_{i}}
和
c
i
{\displaystyle c_{i}}
等九个角度需要满足内角和条件、圆周条件和边长条件三类条件方程:[ 1] [ 17]
内角和条件,指同一三角形内观测的三个内角之和的理论值之和应等于180度,即
a
^
i
+
b
^
i
+
c
^
i
−
180
∘
=
0
{\displaystyle {\hat {a}}_{i}+{\hat {b}}_{i}+{\hat {c}}_{i}-180^{\circ }=0}
;
圆周条件,指在中间点
D
{\displaystyle D}
处观测的三个角度
c
1
{\displaystyle c_{1}}
、
c
2
{\displaystyle c_{2}}
和
c
3
{\displaystyle c_{3}}
的理论值之和应等于360度,即
c
^
1
+
c
^
2
+
c
^
3
−
360
∘
=
0
{\displaystyle {\hat {c}}_{1}+{\hat {c}}_{2}+{\hat {c}}_{3}-360^{\circ }=0}
;
边长条件,指经由不同方式[ 注 1] 求得的同一条边的长度理论值应当相等,对于右图中的测角三角网为
sin
a
^
1
sin
b
^
1
sin
a
^
2
sin
b
^
2
sin
a
^
3
sin
b
^
3
−
1
=
0
{\displaystyle {\frac {\sin {{\hat {a}}_{1}}}{\sin {{\hat {b}}_{1}}}}{\frac {\sin {{\hat {a}}_{2}}}{\sin {{\hat {b}}_{2}}}}{\frac {\sin {{\hat {a}}_{3}}}{\sin {{\hat {b}}_{3}}}}-1=0}
上式可由正弦定理 导出,将未知边长
C
D
¯
{\displaystyle {\overline {CD}}}
通过两种方式表示由已知边长
A
B
¯
{\displaystyle {\overline {AB}}}
的求得:
C
D
¯
=
A
B
¯
sin
a
^
1
sin
c
^
1
sin
a
^
2
sin
b
^
2
=
A
B
¯
sin
b
^
1
sin
c
^
1
sin
b
^
3
sin
a
^
3
{\displaystyle {\overline {CD}}={\overline {AB}}{\frac {\sin {{\hat {a}}_{1}}}{\sin {{\hat {c}}_{1}}}}{\frac {\sin {{\hat {a}}_{2}}}{\sin {{\hat {b}}_{2}}}}={\overline {AB}}{\frac {\sin {{\hat {b}}_{1}}}{\sin {{\hat {c}}_{1}}}}{\frac {\sin {{\hat {b}}_{3}}}{\sin {{\hat {a}}_{3}}}}}
消去作为系数的已知边长
A
B
¯
{\displaystyle {\overline {AB}}}
,即可得到仅由角度测量值表示的边长条件方程。
测量平差中的观测方程是指将观测量表达为参数的函数的方程式,其数学形式可统一成:[ 1] [ 15]
L
^
=
F
(
X
^
)
{\displaystyle {\hat {L}}=\operatorname {F} ({\hat {X}})}
在观测方程中,观测量和参数分别列于等号的两侧,且每一个观测量都对应一条观测方程。[ 17] [ 18] 相较于条件方程,观测方程无需考虑观测量之间应满足的数学关系,更为简单且便于电算。[ 19]
平差的随机模型描述的是平差模型中各类随机变量 自身以及随机量之间的统计相关性质。在经典平差方法中,观测量被假设为仅包含偶然误差的随机变量,服从于正态分布 。因此,常以方差 和中误差 描述这类随机变量自身的精度,并以协方差 描述随机变量之间的相关性,而方差-协方差矩阵即是这类随机变量的随机模型:[ 12] [ 13]
D
=
σ
0
2
Q
=
σ
0
2
P
−
1
{\displaystyle D=\sigma _{0}^{2}Q=\sigma _{0}^{2}P^{-1}}
式中,
D
{\displaystyle D}
是随机变量的方差-协方差阵,通常直接简称为方差阵或协方差阵;
Q
{\displaystyle Q}
和
P
{\displaystyle P}
分布是随机变量的协因数阵和权阵,且互为对方的逆 ;
σ
0
2
{\displaystyle \sigma _{0}^{2}}
是该随机变量的先验单位权方差。
假设独立观测值
L
i
{\displaystyle L_{i}}
的方差是
σ
i
2
{\displaystyle \sigma _{i}^{2}}
,则其有方差阵
D
=
[
σ
1
2
0
⋯
0
0
σ
2
2
⋯
0
⋮
⋮
⋱
⋮
0
0
⋯
σ
n
2
]
{\displaystyle D={\begin{bmatrix}\sigma _{1}^{2}&0&\cdots &0\\0&\sigma _{2}^{2}&\cdots &0\\\vdots &\vdots &\ddots &\vdots \\0&0&\cdots &\sigma _{n}^{2}\end{bmatrix}}}
选定先验单位权中误差为
σ
0
2
{\displaystyle \sigma _{0}^{2}}
,按下式确定观测值的权
p
i
{\displaystyle p_{i}}
:
p
i
=
σ
0
2
σ
i
2
{\displaystyle p_{i}={\frac {\sigma _{0}^{2}}{\sigma _{i}^{2}}}}
因此,各观测值之间权的比例仅取决于观测值的方差(即观测值的精度指标 ),即有
p
1
:
p
2
:
⋯
:
p
n
=
σ
1
−
2
:
σ
2
−
2
:
⋯
:
σ
n
−
2
{\displaystyle p_{1}:p_{2}:\cdots :p_{n}=\sigma _{1}^{-2}:\sigma _{2}^{-2}:\cdots :\sigma _{n}^{-2}}
,按如下方式组成权阵
P
=
[
p
1
0
⋯
0
0
p
2
⋯
0
⋮
⋮
⋱
⋮
0
0
⋯
p
n
]
=
σ
0
2
D
−
1
{\displaystyle P={\begin{bmatrix}p_{1}&0&\cdots &0\\0&p_{2}&\cdots &0\\\vdots &\vdots &\ddots &\vdots \\0&0&\cdots &p_{n}\end{bmatrix}}=\sigma _{0}^{2}D^{-1}}
又观测值的协因数阵与权阵互逆,则有
Q
=
P
−
1
{\displaystyle Q=P^{-1}}
.
若相关观测值 有协方差阵
D
=
[
σ
1
2
σ
12
⋯
σ
1
n
σ
21
σ
2
2
⋯
σ
2
n
⋮
⋮
⋱
⋮
σ
n
1
σ
n
2
⋯
σ
n
2
]
{\displaystyle D={\begin{bmatrix}\sigma _{1}^{2}&\sigma _{12}&\cdots &\sigma _{1n}\\\sigma _{21}&\sigma _{2}^{2}&\cdots &\sigma _{2n}\\\vdots &\vdots &\ddots &\vdots \\\sigma _{n1}&\sigma _{n2}&\cdots &\sigma _{n}^{2}\end{bmatrix}}}
式中
σ
i
2
{\displaystyle \sigma _{i}^{2}}
为各观测量自身的方差,
σ
i
j
=
σ
j
i
(
i
≠
j
)
{\displaystyle \sigma _{ij}=\sigma _{ji}\,(i\neq j)}
为两相关观测量之间的协方差。
选定先验单位权中误差为
σ
0
2
{\displaystyle \sigma _{0}^{2}}
,相关观测值的协因数阵为
Q
=
1
σ
0
2
D
{\displaystyle Q={\frac {1}{\sigma _{0}^{2}}}D}
,相关观测量的权阵则被定义为
P
=
Q
−
1
{\displaystyle P=Q^{-1}}
.
以最小二乘为最优化准则的平差方法,可根据函数模型的不同分为条件平差、间接平差和混合平差等类型。这些模型的共同点在于,其模型中的误差均被假设为仅包含偶然误差,且通过最小二乘准则得到的观测量和参数的估计值具有相同的统计性质。因此,这些平差模型亦被统称为经典平差模型,其在形式上可以统一为概括平差模型。[ 20] [ 21] [ 22]
各类平差模型的计算步骤类似,首先建立函数模型并对其进行线性化,其次组成法方程并求解,再次由解得的联系数或参数求出观测值的改正数和观测值的估计值,最后进行精度评定。
条件平差方法是以仅给出观测量之间约束条件的条件方程作为函数模型的平差方法,其函数模型的一般形式为
F
(
L
^
)
=
0
{\displaystyle \operatorname {F} ({\hat {L}})=0}
,所有的观测值和非零常数都位于方程组的同一侧。方程组中,条件方程的个数与多余观测量
r
{\displaystyle r}
相等,观测量的个数则为
n
{\displaystyle n}
。
取
A
=
d
F
d
L
^
|
L
{\displaystyle A=\left.{\operatorname {d} \!\operatorname {F} \over \operatorname {d} \!{\hat {L}}}\right\vert _{L}}
作为系数矩阵,
A
0
=
F
(
L
^
)
−
A
L
^
{\displaystyle A_{0}=\operatorname {F} ({\hat {L}})-A{\hat {L}}}
作为常数向量,该函数模型可线性化为
A
L
^
+
A
0
=
0
{\displaystyle A{\hat {L}}+A_{0}=0}
;
又由
L
^
=
L
+
V
{\displaystyle {\hat {L}}=L+V}
,因此取闭合差为
W
=
A
L
+
A
0
{\displaystyle W=AL+A_{0}}
,将上式表述为观测值改正数
V
{\displaystyle V}
与观测值闭合差
W
{\displaystyle W}
之间的关系式,即
A
V
+
W
=
0
{\displaystyle AV+W=0}
以最小二乘准则
V
T
P
V
=
min
{\displaystyle V^{T}PV=\min }
作为该函数模型的附加条件进行解算,最终观测值改正数和观测量估计值的解分别为
V
=
P
−
1
A
T
K
=
Q
A
T
K
{\displaystyle V=P^{-1}A^{T}K=QA^{T}K}
L
^
=
L
+
V
{\displaystyle {\hat {L}}=L+V}
由协因数传播律,观测值改正数和观测量估计值的协因数阵分别为
Q
V
V
=
Q
A
T
N
a
a
−
1
A
Q
{\displaystyle Q_{VV}=QA^{T}N_{aa}^{-1}AQ}
Q
L
^
L
^
=
Q
−
Q
V
V
{\displaystyle Q_{{\hat {L}}{\hat {L}}}=Q-Q_{VV}}
解算过程中出现的各矩阵和向量的大小、计算公式和含义如下表所示:
条件平差中各符号的大小含义和计算公式
符号
行数
列数
计算公式
含义
L
{\displaystyle L}
n
{\displaystyle n}
1
{\displaystyle 1}
经实际测量获得
观测量的观测值
V
{\displaystyle V}
n
{\displaystyle n}
1
{\displaystyle 1}
V
=
Q
A
T
K
{\displaystyle V=QA^{T}K}
观测值的改正数
L
^
{\displaystyle {\hat {L}}}
n
{\displaystyle n}
1
{\displaystyle 1}
L
^
=
L
+
V
{\displaystyle {\hat {L}}=L+V}
观测量的估计值,具有最优无偏估计的性质
A
{\displaystyle A}
r
{\displaystyle r}
n
{\displaystyle n}
A
=
d
F
d
L
^
|
L
{\displaystyle A=\left.{\operatorname {d} \!\operatorname {F} \over \operatorname {d} \!{\hat {L}}}\right\vert _{L}}
条件方程的系数矩阵,表达观测量之间的限制关系
A
0
{\displaystyle A_{0}}
n
{\displaystyle n}
1
{\displaystyle 1}
A
0
=
F
(
L
^
)
−
A
L
^
{\displaystyle A_{0}=\operatorname {F} ({\hat {L}})-A{\hat {L}}}
条件方程的常数向量,表达与观测量的取值无关的部分
W
{\displaystyle W}
n
{\displaystyle n}
1
{\displaystyle 1}
W
=
A
L
+
A
0
{\displaystyle W=AL+A_{0}}
条件方程的常数向量,表达观测量之间的不符值或闭合差
P
{\displaystyle P}
n
{\displaystyle n}
n
{\displaystyle n}
由随机模型确定
对称矩阵,观测值的权阵
Q
{\displaystyle Q}
n
{\displaystyle n}
n
{\displaystyle n}
Q
=
P
−
1
{\displaystyle Q=P^{-1}}
对称矩阵,观测值的协因数阵
N
a
a
{\displaystyle N_{aa}}
r
{\displaystyle r}
r
{\displaystyle r}
N
a
a
=
A
Q
A
T
{\displaystyle N_{aa}=AQA^{T}}
法方程的系数矩阵,仅当其可逆时该平差模型有唯一解
K
{\displaystyle K}
r
{\displaystyle r}
1
{\displaystyle 1}
K
=
−
N
a
a
−
1
W
{\displaystyle K=-N_{aa}^{-1}W}
拉格朗日函数中的联系数向量
间接平差方法是以观测方程为函数模型的平差方法,亦称参数平差,其函数的一般形式为
L
^
=
F
(
X
^
)
{\displaystyle {\hat {L}}=\operatorname {F} ({\hat {X}})}
,方程组中一侧仅有观测值,另一侧则是所选参数和非零常数,观测方程的个数与观测量的个数
n
{\displaystyle n}
相等,参数的个数则与必要观测数
t
{\displaystyle t}
相等。
类似于条件平差,取
B
=
d
F
d
X
^
|
X
0
{\displaystyle B=\left.{\operatorname {d} \!\operatorname {F} \over \operatorname {d} \!{\hat {X}}}\right\vert _{X^{0}}}
作为系数矩阵,
d
=
L
^
−
B
X
^
{\displaystyle d={\hat {L}}-B{\hat {X}}}
作为常数向量,该函数模型可线性化为
L
^
=
B
X
^
+
d
{\displaystyle {\hat {L}}=B{\hat {X}}+d}
;
设参数有近似值
X
0
{\displaystyle X^{0}}
,即有
X
^
=
X
0
+
x
^
{\displaystyle {\hat {X}}=X^{0}+{\hat {x}}}
;又由
L
^
=
L
+
V
{\displaystyle {\hat {L}}=L+V}
,因此取约化观测量为
l
=
L
−
F
(
X
0
)
=
L
−
(
B
X
0
+
d
)
=
L
−
L
0
{\displaystyle l=L-\operatorname {F} (X^{0})=L-(BX^{0}+d)=L-L^{0}}
,线性化后的函数模型亦可表示为
V
=
B
x
^
−
l
{\displaystyle V=B{\hat {x}}-l}
该式也称误差方程,式中由参数改正数
x
^
{\displaystyle {\hat {x}}}
表达的观测值改正数(即观测值的后验误差)需满足最小二乘准则
V
T
P
V
=
min
{\displaystyle V^{T}PV=\min }
,以此求得自由极值下的参数改正数和参数估计值分别为[ 23] [ 24]
x
^
=
N
B
B
−
1
W
{\displaystyle {\hat {x}}=N_{BB}^{-1}W}
X
^
=
X
0
+
x
^
{\displaystyle {\hat {X}}=X^{0}+{\hat {x}}}
在误差方程中,由参数改正数
x
^
{\displaystyle {\hat {x}}}
表达的观测量、最终观测值改正数和观测量估计值的解分别为
V
=
B
x
^
−
l
{\displaystyle V=B{\hat {x}}-l}
L
^
=
L
+
V
{\displaystyle {\hat {L}}=L+V}
由协因数传播律,参数估计值、观测值改正数和观测量估计值的协因数阵分别为
Q
X
^
X
^
=
N
B
B
−
1
{\displaystyle Q_{{\hat {X}}{\hat {X}}}=N_{BB}^{-1}}
Q
V
V
=
Q
−
B
N
B
B
−
1
B
T
{\displaystyle Q_{VV}=Q-BN_{BB}^{-1}B^{T}}
Q
L
^
L
^
=
B
N
B
B
−
1
B
T
{\displaystyle Q_{{\hat {L}}{\hat {L}}}=BN_{BB}^{-1}B^{T}}
解算过程中出现的各矩阵和向量的大小、计算公式和含义如下表所示:
间接平差中各符号的大小含义和计算公式
符号
行数
列数
计算公式
含义
L
{\displaystyle L}
n
{\displaystyle n}
1
{\displaystyle 1}
经实际测量获得
观测量的观测值
V
{\displaystyle V}
n
{\displaystyle n}
1
{\displaystyle 1}
V
=
B
x
^
−
l
{\displaystyle V=B{\hat {x}}-l}
观测值的改正数
L
^
{\displaystyle {\hat {L}}}
n
{\displaystyle n}
1
{\displaystyle 1}
L
^
=
L
+
V
{\displaystyle {\hat {L}}=L+V}
观测量的估计值,具有最优无偏估计的性质
X
0
{\displaystyle X^{0}}
t
{\displaystyle t}
1
{\displaystyle 1}
根据实际情况估计
参数的近似值
x
^
{\displaystyle {\hat {x}}}
t
{\displaystyle t}
1
{\displaystyle 1}
x
^
=
N
B
B
−
1
W
{\displaystyle {\hat {x}}=N_{BB}^{-1}W}
参数的改正数
X
^
{\displaystyle {\hat {X}}}
t
{\displaystyle t}
1
{\displaystyle 1}
X
^
=
X
0
+
x
^
{\displaystyle {\hat {X}}=X^{0}+{\hat {x}}}
参数的估计值,具有最优无偏估计的性质
B
{\displaystyle B}
r
{\displaystyle r}
n
{\displaystyle n}
B
=
d
F
d
X
^
|
X
0
{\displaystyle B=\left.{\operatorname {d} \!\operatorname {F} \over \operatorname {d} \!{\hat {X}}}\right\vert _{X^{0}}}
观测方程的系数矩阵,表达观测量和参数的函数关系
d
{\displaystyle d}
n
{\displaystyle n}
1
{\displaystyle 1}
d
=
L
^
−
B
X
^
{\displaystyle d={\hat {L}}-B{\hat {X}}}
观测方程的常数向量,表达与参数的取值无关的部分
L
0
{\displaystyle L^{0}}
n
{\displaystyle n}
1
{\displaystyle 1}
L
0
=
B
X
0
+
d
{\displaystyle L^{0}=BX^{0}+d}
观测值的近似值,由参数的近似值确定
l
{\displaystyle l}
n
{\displaystyle n}
1
{\displaystyle 1}
l
=
L
−
L
0
{\displaystyle l=L-L^{0}}
误差方程的常数向量,表达观测量与近似观测量的差异
P
{\displaystyle P}
n
{\displaystyle n}
n
{\displaystyle n}
由随机模型确定
对称矩阵,观测值的权阵
Q
{\displaystyle Q}
n
{\displaystyle n}
n
{\displaystyle n}
Q
=
P
−
1
{\displaystyle Q=P^{-1}}
对称矩阵,观测值的协因数阵
N
B
B
{\displaystyle N_{BB}}
t
{\displaystyle t}
t
{\displaystyle t}
N
B
B
=
B
T
P
B
{\displaystyle N_{BB}=B^{T}PB}
法方程的系数矩阵,仅当其可逆时该平差模型有唯一解
W
{\displaystyle W}
t
{\displaystyle t}
1
{\displaystyle 1}
W
=
B
T
P
l
{\displaystyle W=B^{T}Pl}
法方程的常数向量
概括平差方法的提出旨在同一条件平差和间接平差等经典平差方法的函数模型,以对其共性和特性进行研究。从概括平差方法出发可以证明各经典平差模型之间的等价性,从而将通过概括平差方法证明的数学性质推广到所有经典平差方法中。附有限制条件的间接平差模型和附有限制条件的条件平差模型均已被证明可作为概括模型的函数模型。[ 20] [ 21] [ 22]
在此以附有限制条件的条件平差为例,其具体形式为
{
F
(
L
^
,
X
^
)
=
0
Φ
(
X
^
)
=
0
{\displaystyle {\begin{cases}\operatorname {F} ({\hat {L}},{\hat {X}})=0\\\operatorname {\Phi } ({\hat {X}})=0\end{cases}}}
模型中共有
n
{\displaystyle n}
个观测量和
u
{\displaystyle u}
个参数,其中独立观测量的个数为
r
{\displaystyle r}
,独立参数的个数为
u
−
s
{\displaystyle u-s}
;因此,一般条件方程的个数为
c
=
r
+
u
−
s
{\displaystyle c=r+u-s}
,限制条件方程的个数为
s
{\displaystyle s}
,且满足
c
+
s
=
r
+
u
{\displaystyle c+s=r+u}
的关系。
取
A
=
∂
F
∂
L
^
|
L
{\displaystyle A=\left.{\operatorname {\partial } \!\operatorname {F} \over \operatorname {\partial } \!{\hat {L}}}\right\vert _{L}}
、
B
=
∂
F
∂
X
^
|
X
0
{\displaystyle B=\left.{\operatorname {\partial } \!\operatorname {F} \over \operatorname {\partial } \!{\hat {X}}}\right\vert _{X^{0}}}
和
C
=
∂
Φ
∂
X
^
|
X
0
{\displaystyle C=\left.{\operatorname {\partial } \!\Phi \over \operatorname {\partial } \!{\hat {X}}}\right\vert _{X^{0}}}
作为系数矩阵,
W
=
F
(
L
,
X
0
)
{\displaystyle W=\operatorname {F} \left(L,X^{0}\right)}
和
W
x
=
Φ
(
X
0
)
{\displaystyle W_{x}=\Phi \left(X^{0}\right)}
作为常数项,将该函数模型线性化为
{
A
V
+
B
x
^
+
W
=
0
C
x
^
+
W
x
=
0
{\displaystyle {\begin{cases}AV+B{\hat {x}}+W=0\\C{\hat {x}}+W_{x}=0\end{cases}}}
类似于条件平差模型,以最小二乘准则
V
T
P
V
=
min
{\displaystyle V^{T}PV=\min }
作为该函数模型的附加条件进行解算,最终参数改正数和观测值改正数的解分别为[ 1] [ 20]
x
^
=
−
(
N
b
b
−
1
−
N
b
b
−
1
C
T
N
c
c
−
1
C
N
b
b
−
1
)
W
e
−
N
b
b
−
1
C
T
N
c
c
−
1
W
x
{\displaystyle {\hat {x}}=-(N_{bb}^{-1}-N_{bb}^{-1}C^{T}N_{cc}^{-1}CN_{bb}^{-1})W_{e}-N_{bb}^{-1}C^{T}N_{cc}^{-1}W_{x}}
V
=
−
Q
A
T
N
a
a
−
1
(
W
+
B
x
^
)
{\displaystyle V=-QA^{T}N_{aa}^{-1}(W+B{\hat {x}})}
参数估计值和观测量估计值分别为
X
^
=
X
0
+
x
^
{\displaystyle {\hat {X}}=X^{0}+{\hat {x}}}
L
^
=
L
+
V
{\displaystyle {\hat {L}}=L+V}
参数估计值、观测值改正数和观测量估计值的协方差阵分别为
Q
X
^
X
^
=
N
b
b
−
1
−
N
b
b
−
1
C
T
N
c
c
−
1
C
N
b
b
−
1
{\displaystyle Q_{{\hat {X}}{\hat {X}}}=N_{bb}^{-1}-N_{bb}^{-1}C^{T}N_{cc}^{-1}CN_{bb}^{-1}}
Q
V
V
=
Q
A
T
(
N
a
a
−
1
−
N
a
a
−
1
B
Q
X
^
X
^
B
T
N
a
a
−
1
)
A
Q
{\displaystyle Q_{VV}=QA^{T}\left(N_{aa}^{-1}-N_{aa}^{-1}BQ_{{\hat {X}}{\hat {X}}}B^{T}N_{aa}^{-1}\right)AQ}
Q
L
^
L
^
=
Q
−
Q
V
V
{\displaystyle Q_{{\hat {L}}{\hat {L}}}=Q-Q_{VV}}
解算过程中出现的各矩阵和向量的大小、计算公式和含义如下表所示:
概括平差中各符号的大小含义和计算公式
符号
行数
列数
计算公式
含义
L
{\displaystyle L}
n
{\displaystyle n}
1
{\displaystyle 1}
经实际测量获得
观测量的观测值
V
{\displaystyle V}
n
{\displaystyle n}
1
{\displaystyle 1}
V
=
−
Q
A
T
N
a
a
−
1
(
W
+
B
x
^
)
{\displaystyle V=-QA^{T}N_{aa}^{-1}(W+B{\hat {x}})}
观测值的改正数
L
^
{\displaystyle {\hat {L}}}
n
{\displaystyle n}
1
{\displaystyle 1}
L
^
=
L
+
V
{\displaystyle {\hat {L}}=L+V}
观测量的估计值,具有最优无偏估计的性质
X
0
{\displaystyle X^{0}}
t
{\displaystyle t}
1
{\displaystyle 1}
根据实际情况估计
参数的近似值
x
^
{\displaystyle {\hat {x}}}
t
{\displaystyle t}
1
{\displaystyle 1}
x
^
=
−
(
N
b
b
−
1
−
N
b
b
−
1
C
T
N
c
c
−
1
C
N
b
b
−
1
)
W
e
−
N
b
b
−
1
C
T
N
c
c
−
1
W
x
{\displaystyle {\hat {x}}=-(N_{bb}^{-1}-N_{bb}^{-1}C^{T}N_{cc}^{-1}CN_{bb}^{-1})W_{e}-N_{bb}^{-1}C^{T}N_{cc}^{-1}W_{x}}
参数的改正数
X
^
{\displaystyle {\hat {X}}}
t
{\displaystyle t}
1
{\displaystyle 1}
X
^
=
X
0
+
x
^
{\displaystyle {\hat {X}}=X^{0}+{\hat {x}}}
参数的估计值,具有最优无偏估计的性质
A
{\displaystyle A}
c
{\displaystyle c}
n
{\displaystyle n}
A
=
∂
F
∂
L
^
|
L
{\displaystyle A=\left.{\operatorname {\partial } \!\operatorname {F} \over \operatorname {\partial } \!{\hat {L}}}\right\vert _{L}}
一般条件方程的系数矩阵,表达一般条件方程和观测量之间的函数关系
B
{\displaystyle B}
c
{\displaystyle c}
u
{\displaystyle u}
B
=
∂
F
∂
X
^
|
X
0
{\displaystyle B=\left.{\operatorname {\partial } \!\operatorname {F} \over \operatorname {\partial } \!{\hat {X}}}\right\vert _{X^{0}}}
一般条件方程的系数矩阵,表达一般条件方程和参数之间的函数关系
C
{\displaystyle C}
s
{\displaystyle s}
u
{\displaystyle u}
C
=
∂
Φ
∂
X
^
|
X
0
{\displaystyle C=\left.{\operatorname {\partial } \!\Phi \over \operatorname {\partial } \!{\hat {X}}}\right\vert _{X^{0}}}
限制条件方程的系数矩阵,表达参数之间的限制关系
W
{\displaystyle W}
c
{\displaystyle c}
1
{\displaystyle 1}
W
=
F
(
L
,
X
0
)
{\displaystyle W=\operatorname {F} \left(L,X^{0}\right)}
一般条件方程的常数向量,表达一般条件方程中与观测量和参数的取值无关的部分
W
x
{\displaystyle W_{x}}
c
{\displaystyle c}
1
{\displaystyle 1}
W
x
=
Φ
(
X
0
)
{\displaystyle W_{x}=\Phi \left(X^{0}\right)}
限制条件方程的常数向量,表达限制条件方程中与参数的取值无关的部分
P
{\displaystyle P}
n
{\displaystyle n}
n
{\displaystyle n}
由随机模型确定
对称矩阵,观测值的权阵
Q
{\displaystyle Q}
n
{\displaystyle n}
n
{\displaystyle n}
Q
=
P
−
1
{\displaystyle Q=P^{-1}}
对称矩阵,观测值的协因阵
N
a
a
{\displaystyle N_{aa}}
c
{\displaystyle c}
c
{\displaystyle c}
N
a
a
=
A
Q
A
T
{\displaystyle N_{aa}=AQA^{T}}
法方程的系数矩阵,仅当其可逆时该平差模型有唯一解
N
b
b
{\displaystyle N_{bb}}
u
{\displaystyle u}
u
{\displaystyle u}
N
b
b
=
B
T
N
a
a
−
1
B
{\displaystyle N_{bb}=B^{T}N_{aa}^{-1}B}
对称矩阵
N
c
c
{\displaystyle N_{cc}}
s
{\displaystyle s}
s
{\displaystyle s}
N
c
c
=
C
N
b
b
−
1
C
T
{\displaystyle N_{cc}=CN_{bb}^{-1}C^{T}}
对称矩阵
W
e
{\displaystyle W_{e}}
u
{\displaystyle u}
1
{\displaystyle 1}
W
e
=
B
T
N
a
a
−
1
W
{\displaystyle W_{e}=B^{T}N_{aa}^{-1}W}
常数向量
根据概括平差方法和各类经典平差方法的函数模型,可以得出如下转换关系:[ 1] [ 20]
概括平差方法与各类经典平差方法间的转换关系
概括平差的函数模型
概括平差的法方程
转换关系
转换后的平差方法
转换后的平差方法的函数模型
转换后的平差方法的法方程
{
A
V
+
B
x
^
+
W
=
0
C
x
^
+
W
x
=
0
{\displaystyle {\begin{cases}AV+B{\hat {x}}+W=0\\C{\hat {x}}+W_{x}=0\end{cases}}}
{
N
a
a
K
+
B
x
^
+
W
=
0
B
T
K
+
C
T
K
s
=
0
C
x
^
+
W
x
=
0
{\displaystyle {\begin{cases}N_{aa}K&+&B{\hat {x}}+W&=0\\B^{T}K&+&C^{T}K_{s}&=0\\&&C{\hat {x}}+W_{x}&=0\end{cases}}}
A
=
−
I
{\displaystyle A=-I}
附有限制条件的间接平差方法
{
V
=
B
x
^
−
l
C
x
^
+
W
x
=
0
{\displaystyle {\begin{cases}V=B{\hat {x}}-l\\C{\hat {x}}+W_{x}=0\end{cases}}}
{
N
B
B
x
^
+
C
T
K
s
−
W
=
0
C
x
^
+
W
x
=
0
{\displaystyle {\begin{cases}N_{BB}{\hat {x}}+C^{T}K_{s}&-\,W&=0\\C{\hat {x}}&+\,W_{x}&=0\end{cases}}}
C
=
0
{\displaystyle C=0}
附有未知参数的条件平差方法
A
V
+
B
x
^
+
W
=
0
{\displaystyle AV+B{\hat {x}}+W=0}
{
N
a
a
K
+
B
x
^
+
W
=
0
B
T
K
=
0
{\displaystyle {\begin{cases}N_{aa}K+B{\hat {x}}+W&=0\\B^{T}K&=0\end{cases}}}
C
=
0
{\displaystyle C=0}
且
B
=
0
{\displaystyle B=0}
条件平差方法
A
V
+
W
=
0
{\displaystyle AV+W=0}
N
a
a
K
+
W
=
0
{\displaystyle N_{aa}K+W=0}
C
=
0
{\displaystyle C=0}
且
A
=
−
I
{\displaystyle A=-I}
间接平差方法
V
=
B
x
^
−
l
{\displaystyle V=B{\hat {x}}-l}
N
B
B
x
^
−
W
=
0
{\displaystyle N_{BB}{\hat {x}}-W=0}
值得注意的是,当系数矩阵
A
=
−
I
{\displaystyle A=-I}
时,该方法转化为附有限制条件的间接平差方法,且该方法中亦包含有一般观测方程和限制条件方差。因此附有限制条件的间接平差方法是与附有限制条件的条件平差方法等价的概括平差方法。[ 21] [ 22]
^ 这里的不同方式是指,通过不同方式列出的数学表达式中至少包含一个不同的独立变量
^ 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 武汉大学测绘学院测量平差学科组.误差理论与测量平差基础(第三版).武汉:武汉大学出版社,2014.ISBN 978-7-307-12922-1 .
^ 2.0 2.1 宁津生,陈俊勇,李德仁,等.测绘学概论(第三版).武汉:武汉大学出版社,2016.ISBN 978-7-307-18610-1 .
^ Fan, Huaan. Theory of Errors and Least Squares Adjustment . Tekniska högskolan. 1997. ISBN 978-91-7170-200-5 (英语) .
^ Ghilani, Charles D. Adjustment Computations: Spatial Data Analysis . John Wiley & Sons. 2017-10-23. ISBN 978-1-119-38598-1 (英语) .
^ Xue, S. Q.; Yang, Y. X. Gauss–Jacobi combinatorial adjustment and its modification . Survey Review. 2013-12-19, 46 (337): 298–304. ISSN 0039-6265 . doi:10.1179/1752270613y.0000000086 (英语) .
^ Awange, Joseph; Paláncz, Béla. Geospatial Algebraic Computations: Theory and Applications . Springer. 2016-01-29. ISBN 978-3-319-25465-4 (英语) .
^ Sneeuw, Nico, Friedhelm Krumm, and Matthias Roth. "Adjustment Theory." University Stuttgart: Stuttgart, Germany (2014). (英语)
^ 8.0 8.1 Okwuashi, Onuwa and Inemesit E Asuquo. “Basics of Least Squares Adjustment Computation in Surveying.” (页面存档备份 ,存于互联网档案馆 ) (2014). (英语)
^ Teunissen, Peter. Adjustment theory: an introduction . VSSD Press. 2000 [2020-05-21 ] . ISBN 978-90-407-1974-5 . (原始内容存档 于2019-04-14) (英语) .
^ Kennie, T. J. M. Engineering Surveying Technology . CRC Press. 2014-04-21. ISBN 978-1-4822-6937-6 (英语) .
^ 姚宜斌,邱卫宁.测量平差问题中必要观测数的确定 .测绘通报,2007(03):14-15+18.
^ 12.0 12.1 夏春林主编.误差理论与测量平差.北京:清华大学出版社,2015.ISBN 978-7-302-38074-0 .
^ 13.0 13.1 王穗辉.误差理论与测量平差.上海:同济大学出版社,2010.ISBN 978-7-560-84209-7 .
^ 王侬,过静珺主编.现代普通测量学 (页面存档备份 ,存于互联网档案馆 ).北京:清华大学出版社,2001.ISBN 978-7-302-04717-9 .
^ 15.0 15.1 15.2 孙海燕. 测量平差函数模型的若干讨论 . 武汉大学学报信息科学版. 2003, 28 (S1): 62–64. ISSN 1671-8860 .
^ Ghilani, Charles D. Adjustment Computations: Spatial Data Analysis . John Wiley & Sons. 2017-10-23 [2020-06-05 ] . ISBN 978-1-119-38598-1 . (原始内容存档 于2020-06-08) (英语) .
^ 17.0 17.1 吴究. 測量平差 . 国立中央大学出版中心. 2012-09-16 [2020-06-05 ] . ISBN 978-986-87944-7-4 . (原始内容存档 于2020-06-08).
^ Kahmen, Heribert; Faig, Wolfgang. Surveying . Walter de Gruyter. 1988 [2020-06-02 ] . ISBN 978-3-11-008303-3 . (原始内容存档 于2020-06-11) (英语) .
^ 潘正风等.数字地形测量学.武汉:武汉大学出版社,2015.ISBN 978-7-307-15677-7
^ 20.0 20.1 20.2 20.3 于宗俦; 王正林. 各种平差方法的统一与特性 . 武汉测绘科技大学学报. 1986, 26 (04): 30-43 [2020-06-11 ] . doi:10.13203/j.whugis1986.04.004 . (原始内容存档 于2020-06-11).
^ 21.0 21.1 21.2 王新洲. 论经典测量平差模型的内在联系 . 测绘通报. 2004, (2) [2020-06-09 ] . ISSN 0494-0911 . (原始内容存档 于2020-06-09).
^ 22.0 22.1 22.2 赵超英; 张勤. 再论经典测量平差模型间的内在联系 . 测绘通报. 2006, (03): 26-27 [2020-06-11 ] . ISSN 0494-0911 . (原始内容存档 于2020-06-11).
^ 観測方程式が . nippo1.co.jp. [2020-06-05 ] . (原始内容存档 于2016-11-09) (日语) .
^ Brinker, Russell C.; Minnick, Roy. The Surveying Handbook 2. Springer US. 1995 [2020-06-09 ] . ISBN 978-1-4613-5858-9 . (原始内容存档 于2020-06-09) (英语) .
崔希璋.广义测量平差(第二版).武汉:武汉大学出版社,2009.ISBN 978-7-307-07268-8 .
王新洲.高等测量平差.北京:测绘出版社,2006.ISBN 978-7-503-01396-6 .
武汉大学测绘学院测量平差学科组.误差理论与测量平差基础(第三版).武汉:武汉大学出版社,2014.ISBN 978-7-307-12922-1 .
吴究.测量平差 (页面存档备份 ,存于互联网档案馆 ).国立中央大学出版中心,2012.ISBN 978-986-87944-7-4 .
Fan, Huaan. Theory of Errors and Least Squares Adjustment . Royal Institute of Technology (KTH), Division of Geodesy and Geoinformatics, Stockholm, Sweden, 2010, ISBN 91-7170-200-8 .
Ghilani, Charles D. Adjustment Computations: Spatial Data Analysis. John Wiley & Sons, 2017. ISBN 978-1-119-38598-1 .
Teunissen, Peter. Adjustment theory: an introduction . VSSD Press, 2000, ISBN 978-9-040-71974-5 .
Thomas Wallace Wright.The Adjustment of Observations by the Method of Least Squares: With Applications to Geodetic Work. (页面存档备份 ,存于互联网档案馆 ) D. Van Nostrand Company, 1906.