以平面三角形 為例,由於測量誤差的存在,對三個角度進行多次測量後,得到的各角度觀測值內角和並不等於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.