武汉大学测量平差[第2部分-2]

（1）
式中V为观测值L的改正数，xˆ为参数近似值 X 0的改正数。其系数
矩阵的秩分别为 rk ( A) = c, rk (B) = u 。其随机模型为：
D LL
=
σ
2 0
Q
LL
=
σ
2 0
P
−1
（1）式中的未知数为n个观测值的改正数V 和u个参数近似值的
改正数样 xˆ ，即未知数的个数为m = n + u，而方程的个数为
Bxˆ
+
W
)
（9）
（8）式和（9）式就是附有参数的条件平差的最终解。
第五讲 平差方法——附有参数的条件平差（续）
2、附有参数的条件平差的计算步骤
由以上推导，可总结出附有参数的条件平差的计算步骤如下：
（1）、根据具体的平差问题，选取u个独立的参数，并列出附有参数的条件
方程（1）式。 A V + B xˆ + W = 0 c×n n×1 c×u u×1 c×1 c×1
L,X 0
=
F (L,
X
0)
则上式变为：A Δ+ B ~x +W = 0
c.n n.1 c.u u.1 c.1 c.1
第五讲 平差方法——附有参数的条件平差（续）
二、附有参数的条件平差原理
一般地，附有参数的条件平差的函数模型为：
A V + B xˆ + W = 0
c×n n×1 c×u u×1 c×1 c×1
N aa −1左乘（4）式的第一式，可得：
K
=
−
N
−1 aa
(
Bxˆ
+W)
(5)
再以
BT
N
−1 aa
左乘（4）式的第一式，并减去第二式，得：
B
T
N
−1 aa
Bxˆ
+
B
T
N
W −1
aa
=
0
第五讲 平差方法——附有参数的条件平差（续）
令 则有
N bb
=
BT
N
−1 aa
B
Nbb xˆ +
B
T
N
W −1
aa
第五讲 平差方法——附有参数的条件平差
一、问题的提出
由条件平差知，对于n个观测值，t个必要观测（n>t）的条件 平差问题，可以列出r=n-t个独立的条件方程，且列出r个独立的条 件方程后就可以进行后继的条件平差计算。然而，在实际工作 中，有些平差问题的r个独立的条件方程很难列出。例如，在下图 所示的测角网中，A、B为已知点，AC为已知边。观测了网中的9 个角度，即n=9。要确定C、D、E三点的坐标，其必要观测数为 t=5，故条件方程的个数为r=n-t=9-5=4，即必须列出4个独立的条 件方程。由图知，三个图形条件很容易列出，但第四个条件却不 容易列出。
（2）、组成法方程（4）式。
Naa K + Bxˆ + W = 0
BT K = 0
（3）、按（8）式和（9）式计算参数近似值的改正数 xˆ 和观测值L的改正数
V。
xˆ
=
−
N
−1 bb
BT
N
W −1
aa
V
=
−
P
−1
AT
N
−1 aa
(
Bxˆ
+
W
)
（4）、计算观测值和参数的平差值。Lˆ = L +V , Xˆ = X 0 + xˆ （5）、用平差值重新列平差值条件方程，检核整个计算的正确性。F(L~, X~) = 0
由此可得观测值和参数的平差值为：
( ) LˆT = 60D00′05.4′′ 59D59′56.3′′ 59D59′58.3′′ 60D00′05.0′′ 59D59′56.0′′ 59D59′59.0′′
Xˆ = 30D00′05.4′′
检核略。
第五讲 平差方法——附有参数的条件平差
三、精度评定 1、单位权方差的估值
F (L~, X~) = 0
第五讲 平差方法——附有参数的条件平差（续）
1）如果条件方程为线性形式，则有
A L~+
c.n n.1
B
c.u
X~
u.1
+ A0
c.1
=
0
c.1
A0 : 常数向量 将L~ = L + Δ, X~ = X 0 + ~x代入上式，并令
W = AL + BX 0 + A0 则上式变为：A Δ+ B ~x +W = 0
x A = 1000 .00m, y A = 0.00m, xB = 1000 .00m, y B = 1732 .00m , S BD = 1000 .00m
各角的同精度独立观测值见表1。现选 ∠BAD 的最或 是值为参数，试按附有参数的条件平差求观测值的平 差值和参数的平差值。
第五讲 平差方法——附有参数的条件平差（续）
第五讲 平差方法——附有参数的条件平差（续）
为了解决这个问题，可以选择某个（或某几个）非观测量作 为参数。例如图中选择 X 作为参数。设选择了u个参数，则原来 的r个条件方程就变为c = r+u个了。如图中，由于选择了X 作为参 数，则条件方程的个数就变为c = r+u = 4+1=5个，即除了三个图 形条件外，还可以列出1个极条件和1个固定边条件。如下图，若 以A点为极，则极条件(以A点为极）为：
Lˆ2
或
S AC sin(Lˆ6 + Lˆ8 ) sin Lˆ2 S AB sin Xˆ sin Lˆ3
=1
第五讲 平差方法——附有参数的条件平差（续）
设在平差问题中，观测值个数为n,必要观测数为t,则 可列出r=n-t个条件方程。现又增设了u个独立量作为参 数，而0<u<t，每增设一个参数应增加一个条件方程。 以含有参数的条件方程作为平差的函数模型的平差方 法，称为附有参数的条件平差法。
c = r + u。由于m – c = n – r = t > 0，所以（1）式是一组具有无
穷多组解的相容方程组。必须根据最小二乘原理，求出能使
V T PV = min 的一组解。为此，下面就来求解这组解。
第五讲 平差方法——附有参数的条件平差（续）
1、基础方程及其解 为了求得能使 V T PV = min 的一组解，按求函数之条件 极值的方法，组成新函数： Φ = V T PV − 2K T ( AV + Bxˆ + W ) 式中K是对应（6-1）式的联系数向量。
表1
角号
观测值
角号
观测值
1
60D00′03′′
4
59D59′57′′
2
60D00′02′′
5
59D59′56′′
3
60D00′04′′
6
59D59′59′′
第五讲 平差方法——附有参数的条件平差（续）
解：本例中n = 6，t = 3，r = 3，u = 1，故c = r + u = 4
由图知，可列2个图形条件，1个极条件和1个固定边条件。 这4个条件如下：
⎟ ⎟ ⎟⎟⎠
xˆ
+
⎜ ⎜ ⎜⎜⎝
−8 ⎟ −56.1.09561⎟⎟⎟⎠
=
0
由于为同精度独立观测，故 P = I 。于是由（4）式得法方程为：
⎜⎛ 3.000 0 1.732 0.577 ⎟⎞ ⎜⎛ 0 ⎟⎞ ⎜⎛ 9 ⎟⎞
⎜0 ⎜⎜⎜⎝10..753727
3.000 0
0.577
0 6.334 − 0.999
⎜ ⎜
QVL
QXˆW QKW QVW
QXˆXˆ QKXˆ QVXˆ
QXˆK QKK QVK
QXˆV QKV QVV
QXˆLˆ
⎟ ⎟
QKLˆ QVLˆ
⎟ ⎟ ⎟
⎜⎝QLˆL QLˆW QLˆXˆ QLˆK QLˆV QLˆLˆ ⎟⎠
⎜⎛
QLL
QLL AT
−QLLAT Na−a1BQXˆXˆ
− QLL AT QKK
−000..56.9769769 ⎟⎟⎟⎟⎠ K
+
⎜ ⎜ ⎜⎜⎝
−13.7.043624 ⎟⎟⎟⎟⎠ xˆ
+
⎜ ⎜ ⎜⎜⎝
−8 ⎟ −56.1.09561⎟⎟⎟⎠
=
0
(0 0 − 3.464 1.732)K = 0
第五讲 平差方法——附有参数的条件平差（续）
解得： xˆ = 5.4009′′, V T = (2.4′′ − 5.7′′ − 5.7′′ 8.0′′ 0.0 0.0)
v1 + v2 + v3 + wa = 0
v4 + v5 + v6 + wb = 0
sin Lˆ4 sin(Lˆ1 − Xˆ )sin(Lˆ3 + Lˆ5 sin Lˆ5 sin(Lˆ2 + Lˆ4 )sin Xˆ
)
=
1(以B点为极）
SAB sin SBD sin(Lˆ3
Xˆ +
Lˆ5 )
=
1
第五讲 平差方法——附有参数的条件平差（续）
第五讲 平差方法——附有参数的条件平差
N −1 aa
N −1 aa
AQLL AQLL
⎟ ⎟ ⎟
N −1 aa
(I
−
BQXˆXˆ
BT
N −1 aa
)
QKK AQLL
QLL AT QKK
QLL AT QKK AQLL
0
⎟
0