Skip to content

Commit ae4bbba

Browse files
committed
测边网间接平差
1 parent 11e9b89 commit ae4bbba

4 files changed

Lines changed: 352 additions & 0 deletions

File tree

Lines changed: 186 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,186 @@
1+
# 10 间接平差的应用
2+
3+
## 测边网间接平差
4+
5+
以边长为观测值、以未知点纵横坐标为参数,可列立方程。
6+
7+
### 单测边误差方程
8+
9+
#### 单未知点
10+
11+
![](./images/sigle-unknown-segment.svg)
12+
13+
- 参数:坐标 $X_k,Y_k$
14+
- 近似值:坐标 $X_k^0,Y_k^0$,边长 $L_{jk}^0$
15+
16+
有观测方程
17+
18+
$$
19+
\hat L_{jk}=\sqrt{(\hat X_k-X_j)^2+(\hat Y_k-Y_j)^2}
20+
$$
21+
22+
记 $\hat X_k=X_k^0+\hat x_k$,$\hat Y_k=Y_k^0+\hat y_k$,有线性化
23+
24+
$$
25+
\begin{aligned}
26+
\hat L_{jk}
27+
&=L_{jk}^0+\frac{\partial \hat L_{jk}}{\partial\hat X_k}\hat x_k
28+
+\frac{\partial \hat L_{jk}}{\partial\hat Y_k}\hat y_k \\
29+
&=L_{jk}^0+\frac{X_k^0-X_j}{L_{jk}^0}\hat x_k
30+
+\frac{Y_k^0-Y_j}{L_{jk}^0}\hat y_k \\
31+
&=L_{jk}^0+\hat x_k\cos\alpha_{jk}^0+\hat y_k\sin\alpha_{jk}^0
32+
\end{aligned}
33+
$$
34+
35+
#### 双未知点
36+
37+
![](./images/double-unknown-segment.svg)
38+
39+
- 参数:坐标 $X_j,Y_j,X_k,Y_k$
40+
- 近似值:坐标 $X_j^0,Y_j^0,X_k^0,Y_k^0$,边长 $L_{jk}^0$
41+
42+
有观测方程
43+
44+
$$
45+
\hat L_{jk}=\sqrt{(\hat X_k-\hat X_j)^2+(\hat Y_k-\hat Y_j)^2}
46+
$$
47+
48+
记 $\hat X_j=X_j^0+\hat x_j$,$\hat Y_j=Y_j^0+\hat y_j$,$\hat X_k=X_k^0+\hat x_k$,$\hat Y_k=Y_k^0+\hat y_k$,有线性化
49+
50+
$$
51+
\begin{aligned}
52+
\hat L_{jk}
53+
&=L_{jk}^0+\frac{\partial \hat L_{jk}}{\partial X_j}\hat x_j
54+
+\frac{\partial \hat L_{jk}}{\partial\hat Y_j}\hat y_j
55+
+\frac{\partial \hat L_{jk}}{\partial\hat X_k}\hat x_k
56+
+\frac{\partial \hat L_{jk}}{\partial\hat Y_k}\hat y_k \\
57+
&=L_{jk}^0+\frac{X_k^0-X_j^0}{L_{jk}^0}(\hat x_k-\hat x_j)
58+
+\frac{Y_k^0-Y_j^0}{L_{jk}^0}(\hat y_k-\hat y_j) \\
59+
&=L_{jk}^0+(\hat x_k-\hat x_j)\cos\alpha_{jk}^0+(\hat x_k-\hat x_j)\sin\alpha_{jk}^0
60+
\end{aligned}
61+
$$
62+
63+
### 测边网求解
64+
65+
![](./images/edge-net.svg)
66+
67+
$$
68+
\begin{aligned}
69+
{\color{blue}\hat L_1}
70+
&=\sqrt{(X_A-{\color{red}\hat X_2})^2
71+
+(Y_A-{\color{red}\hat Y_2})^2} \\
72+
73+
{\color{blue}\hat L_2}
74+
&=\sqrt{(X_B-{\color{red}\hat X_2})^2
75+
+(Y_B-{\color{red}\hat Y_2})^2} \\
76+
77+
{\color{blue}\hat L_3}
78+
&=\sqrt{({\color{red}\hat X_1}-{\color{red}\hat X_2})^2
79+
+({\color{red}\hat Y_1}-{\color{red}\hat Y_2})^2} \\
80+
81+
{\color{blue}\hat L_4}
82+
&=\sqrt{(X_A-{\color{red}\hat X_1})^2
83+
+(Y_A-{\color{red}\hat Y_1})^2} \\
84+
85+
{\color{blue}\hat L_5}
86+
&=\sqrt{(X_B-{\color{red}\hat X_1})^2
87+
+(Y_B-{\color{red}\hat Y_1})^2} \\
88+
\end{aligned}
89+
$$
90+
91+
线性化得到 $\hat{\boldsymbol L}=\boldsymbol B\hat{\boldsymbol X}+\boldsymbol d$
92+
93+
$$
94+
\begin{bmatrix}
95+
\hat L_1\\
96+
\hat L_2\\
97+
\hat L_3\\
98+
\hat L_4\\
99+
\hat L_5
100+
\end{bmatrix}
101+
=
102+
\begin{bmatrix}
103+
0 & 0 & \cos\alpha^0_{A2} & \sin\alpha^0_{A2}\\
104+
0 & 0 & \cos\alpha^0_{B2} & \sin\alpha^0_{B2}\\
105+
\cos\alpha^0_{21} & \sin\alpha^0_{21} & -\cos\alpha^0_{21} & -\sin\alpha^0_{21}\\
106+
\cos\alpha^0_{A1} & \sin\alpha^0_{A1} & 0 & 0\\
107+
\cos\alpha^0_{B1} & \sin\alpha^0_{B1} & 0 & 0
108+
\end{bmatrix}
109+
\begin{bmatrix}
110+
\hat x_1\\
111+
\hat y_1\\
112+
\hat x_2\\
113+
\hat y_2
114+
\end{bmatrix}
115+
+
116+
\begin{bmatrix}
117+
L^0_{A2}\\
118+
L^0_{B2}\\
119+
L^0_{21}\\
120+
L^0_{A1}\\
121+
L^0_{B1}
122+
\end{bmatrix}
123+
$$
124+
125+
边长误差方程的特点:
126+
127+
- 两待定点的坐标参数的系数等值反号
128+
- 已知点坐标系数不出现在设计矩阵中
129+
- $j,k$ 均为已知点时,无误差方程
130+
- $j,k$ 边的误差方程,$j\to k$ 或 $k\to j$ 列立,其结果相同
131+
132+
### 精度评定
133+
134+
$$
135+
\hat{\sigma}_0^2
136+
=
137+
\frac{\boldsymbol V^{\rm T}\boldsymbol{PV}}{n-t}
138+
=
139+
\frac{\boldsymbol V^{\rm T}\boldsymbol{PV}}{r}
140+
\qquad
141+
\boldsymbol{Q}_{\hat{x}\hat{x}}=\boldsymbol{N}^{-1}
142+
$$
143+
144+
- 函数式
145+
146+
$$
147+
\hat{F}
148+
=
149+
\hat{L}_{jk}
150+
=
151+
\sqrt{
152+
(\hat{X}_k-\hat{X}_j)^2
153+
+
154+
(\hat{Y}_k-\hat{Y}_j)^2
155+
}
156+
$$
157+
158+
- 权函数式
159+
160+
$$
161+
\begin{gathered}
162+
\mathrm d\hat{F}
163+
=
164+
-\frac{\Delta X^0_{jk}}{L^0_{jk}}\mathrm d\hat{x}_j
165+
-\frac{\Delta Y^0_{jk}}{L^0_{jk}}\mathrm d\hat{y}_j
166+
+\frac{\Delta X^0_{jk}}{L^0_{jk}}\mathrm d\hat{x}_k
167+
+\frac{\Delta Y^0_{jk}}{L^0_{jk}}\mathrm d\hat{y}_k
168+
=
169+
\boldsymbol{f}^{T}\mathrm d\hat{\boldsymbol{x}} \\
170+
171+
Q_{\hat{F}\hat{F}}
172+
=
173+
\boldsymbol{f}^{T}\boldsymbol{Q}_{\hat{x}\hat{x}}\boldsymbol{f} \\
174+
175+
\hat{\sigma}_{\hat{F}}
176+
=
177+
\hat{\sigma}_0\sqrt{Q_{\hat{F}\hat{F}}}
178+
\end{gathered}
179+
$$
180+
181+
- 相对精度
182+
$$
183+
\frac{\hat{\sigma}_{\hat{F}}}{\hat{F}}
184+
=
185+
\frac{\hat{\sigma}_0\sqrt{Q_{\hat{F}\hat{F}}}}{\hat{L}_{jk}}
186+
$$

0 commit comments

Comments
 (0)