条件平差详解

条件平差

平差方法

基本的线代基础
$$\begin{array}{rl}& (AB{)}^T=B^T{A}^T\\ & [(AB{)}^{-1}{]}^T=[(AB{)}^T{]}^{-1}\\ & \frac{\partial A{B}^T}{\partial X}=A^T\frac{\partial B}{\partial X}+B^T\frac{\partial A}{\partial X}\end{array}$$
条件平差的函数模型：
$$\begin{array}{l}AV-W=0\\ W=-(AL+A_0)\end{array}$$
随机模型：
$$D={\sigma }_0^{2}Q={\sigma }_0^2{P}^{-1}$$
估计准则：
$$V^{T}PV=\min$$
根据拉格朗日乘数法：
$$\Phi =V^{T}PV-2K^T(AV-W)$$
对其求导：
$$\frac{\partial \Phi }{\partial V}=2V^{T}P-2K^{T}A=0$$
即：
$$P^{T}V=A^{T}K\to V=P^{-1}{A}^{T}K=Q{A}^{T}K$$
其中权阵$P$的转置等于其自身

将式代入： $AV-W=0$

则有：
$$AQ{A}^{T}K-W=0\to K=(AQ{A}^T{)}^{-1}W=N^{-1}W$$
再将$K$的解代回：$V=Q{A}^{T}K$

则有：
$$V=Q{A}^T{N}^{-1}W=Q{A}^T(AQ{A}^T{)}^{-1}W$$

精度评定

按照单位权中误差的定义，有：
$${\hat{\sigma }}_0=\sqrt{\frac{V^{T}PV}{n-t}}=\sqrt{\frac{V^{T}PV}{r}}$$
上述参数都可根据平差过程中求得，在此不再赘述。

按照协因数阵与方差阵之间的关系，有：
$$D_{\hat{L}\hat{L}}={\hat{\sigma }}_0^2{Q}_{\hat{L}\hat{L}}$$
令
$$Z=\left[\begin{array}{c}L\\ W\\ K\\ V\\ \hat{L}\end{array}\right]=\left[\begin{array}{c}I\\ -A\\ -N^{-1}A\\ -P^{-1}{A}^T{N}^{-1}A\\ I-P^{-1}{A}^T{N}^{-1}A\end{array}\right]L+\left[\begin{array}{c}0\\ -A_0\\ -N^{-1}{A}_0\\ -P^{-1}{A}^T{N}^{-1}{A}_0\\ -P^{-1}{A}^T{N}^{-1}{A}_0\end{array}\right]$$
此处推导 $\hat{L}$的计算过程：
$$\hat{L}=L+V=L+Q{A}^T{N}^{-1}W=L+Q{A}^T{N}^{-1}[-(AL+A_0)]=(I-Q{A}^T{N}^{-1}A)L-Q{A}^T{N}^{-1}{A}_0$$
则：
$$Q_{\hat{L}\hat{L}}=(I-P^{-1}{A}^T{N}^{-1}A)Q_{LL}(I-P^{-1}{A}^T{N}^{-1}A{)}^T$$
接下来是变换的几个关键内容：

由于：$N=AQ{A}^T$为可逆矩阵，则 $(N^{-1}{)}^T=(N^T{)}^{-1}=[(AQ{A}^T{)}^T{]}^{-1}=(AQ{A}^T{)}^{-1}=N^{-1}$

将上述关系式进行化简：
$$\begin{array}{rl}{Q}_{\hat{L}\hat{L}}& =(Q_{LL}-Q_{LL}{A}^T{N}^{-1}A{Q}_{LL})(I-Q_{LL}{A}^T{N}^{-1}A{)}^T\\ & =(Q_{LL}-Q_{LL}{A}^T{N}^{-1}A{Q}_{LL})(I-A^T{N}^{-1}A{Q}_{LL})\\ & =Q_{LL}-Q_{LL}{A}^T{N}^{-1}A{Q}_{LL}-Q_{LL}{A}^T{N}^{-1}A{Q}_{LL}+Q_{LL}{A}^T{N}^{-1}(A{Q}_{LL}{A}^T)N^{-1}A{Q}_{LL}\\ & =Q_{LL}-Q_{LL}{A}^T{N}^{-1}A{Q}_{LL}-Q_{LL}{A}^T{N}^{-1}A{Q}_{LL}+Q_{LL}{A}^T{N}^{-1}N{N}^{-1}A{Q}_{LL}\\ & =Q_{LL}-Q_{LL}{A}^T{N}^{-1}A{Q}_{LL}-Q_{LL}{A}^T{N}^{-1}A{Q}_{LL}+Q_{LL}{A}^T{N}^{-1}A{Q}_{LL}\\ & =Q_{LL}-Q_{LL}{A}^T{N}^{-1}A{Q}_{LL}\end{array}$$
其余结果可参考下表：
$$\begin{array}{rl}{Q}_{ZZ}& =\left[\begin{array}{ccccc}{Q}_{LL}& Q_{LW}& Q_{LK}& Q_{LV}& Q_{\hat{L}\hat{L}}\\ Q_{WL}& Q_{WW}& Q_{WK}& Q_{WV}& Q_{W\dot{L}}\\ Q_{KL}& Q_{KW}& Q_{KK}& Q_{KV}& Q_{K\dot{L}}\\ Q_{VL}& Q_{VW}& Q_{VK}& Q_{VV}& Q_{V\hat{L}}\\ Q_{\hat{L}L}& Q_{\hat{L}W}& Q_{\hat{L}K}& Q_{\hat{L}V}& Q_{\hat{L}\hat{L}}\end{array}\right]\\ & =\left[\begin{array}{ccccc}Q& -Q{A}^T& -Q{A}^T{N}^{-1}& -Q{A}^T{N}^{-1}AQ& Q-Q{A}^T{N}^{-1}AQ\\ -AQ& N& I& AQ& 0\\ -N^{-1}AQ& I& N^{-1}& N^{-1}AQ& 0\\ -Q{A}^T{N}^{-1}AQ& Q{A}^T& Q{A}^T{N}^{-1}& Q{A}^T{N}^{-1}AQ& 0\\ Q-Q{A}^T{N}^{-1}AQ& 0& 0& 0& Q-Q{A}^T{N}^{-1}AQ\end{array}\right]\end{array}$$

条件平差应用

应用	必要起算数据个数d	必要观测数据个数t	测角网中条件方程的列立
水准网	1
测角网	4	$2p-q-4$	图形条件：三角形内角和=$180^{\circ }$ 极条件：边长条件 圆周角条件：圆周角为 $360^{\circ }$
侧边网	3	$2p-q-3$	图形条件：三角形内角和=$180^{\circ }$ 极条件：边长条件 圆周角条件：圆周角为 $360^{\circ }$
边角网	3	$2p-q-3$	方位角条件：从一个已知方位角出发，推算到另一个已知方位角，方位角的推算值要与已知方位角相等。方位角条件的个数等于已知方位角的个数-1。 纵横坐标条件 ： 从一个已知点的纵横坐标出发，推算到不同已知点组中另一个已知点时，纵横坐标的推算值要与已知值相等。纵横坐标的条件个数等于多余已知点的个数乘以2-2。 图形条件：三角形内角和=$180^{\circ }$ 边长条件：极条件，余弦条件 圆周角条件：圆周角为 $360^{\circ }$
单一导线网		$2p^{\prime }$	方位角条件：从一个已知方位角出发，推算到另一个已知方位角，方位角的推算值要与已知方位角相等。方位角条件的个数等于已知方位角的个数-1。 纵横坐标条件 ： 从一个已知点的纵横坐标出发，推算到不同已知点组中另一个已知点时，纵横坐标的推算值要与已知值相等。纵横坐标的条件个数等于多余已知点的个数乘以2-2。
GPS基线向量网		$3p^{\prime }$	方位角条件 直角条件 距离条件 面积条件

其中t为必要观测数，p为所有三角点的个数，q为多余起算数据的个数，d为基准个数

一弧度等于多少秒：
$${\rho }^{\prime \prime }=\frac{(180\ast 3600{)}^{\prime \prime }}{\pi }$$