3 控制网平差

法。本节讨论独立三角网按角度进行条件平差
时，条件方程式列立、法方程式组成和解算的
详细步骤和方法。
条件平差时，关键是列出条件方程。独立 三角网的观测量主要是三角形的内角，这些角 在几何上应该满足一定的条件，这些条件就是 列立条件方程的基础。 根据几何条件的不同，独立三角网的条件 方程分为图形条件、圆周角条件、极条件、基 线条件四种类型。
差时就以这 t 个独立量为参数，模型中的所有
量都一定是这 t 个独立参数的函数，亦即每个
观测量都可表达成所选 t 个独立参数的函数。
• 选择几何模型中 t 个独立量为平差参数，将每
一个观测量表达成所选参数的函数，即列出 n
个这种函数关系式，以此为平差的函数模型，
称为间接平差法，又称参数平差法。
例如： △ABC中，观测量为其中的三个内角，选 定∠A和∠B为平差参数，设为X1和 X2，将 每一个观测量均表达为这两个平差参数的 函数，构成数学模型: C
p1 . P . .
. p2 . .
. . ... .
. . . pn
显然 P 是一个对角阵，其逆存在，且：
1 p1 . . . . 1 p2 . . . . ... . . . . 1 pn
则间接平差的函数模型可用以下矩阵形式表达：
L+V=BX+d 或： V=BX – l 此式称为间接平差误差方程。 式中，L 为观测值向量（ n 1 阶）； V 为改正数向量（ n 1 阶） ； B 为系数矩阵（ n t 阶） ； X 为未知数向量（ t 1 阶） ；
l =L – d 为常数矩阵（ n 1 阶） 。
bi ai
ci
an
bn
bi
cn
(a)
(b)
• 对于大地四边形，
可以列出7个图形条件， 但是只有 3 个是相互独 立的，其余几个可以由 这 3 个方程推导出来：
b2
a3
b3
a4
a2 b1
b4 a1
a1 b1 a2 b2 180 0 a2 b2 a3 b3 180 0 a3 b3 a4 b4 180 0
va1+vb1+va2+ vb2 +w1=0; va2+vb2+va3+ vb3 +w2=0; va3+vb3+va4+ vb4 +w3=0; w1= a1+ b1 +a2 +b2- 180º
a2 b2
a3
b3
a4
b4 b1 a1
w2= a2+ b2 +a3 +b3- 180º
w3= a3+ b3 +a4 +b4- 180º
是确定条件方程满足VTPV=min 的唯
一解。
根据计算函数的条件极值的拉格朗日乘 数法则组成新函数： Ф = VTPV－ 2KT（AV+W)
其中： K =（k1, k2，…，kr )T 是拉格朗日 乘数，测量平差中称之为联系数向量。 显然，只要令Ф对V的一阶导数等于零就 可以求出 VTPV 的极值。
（3）解法方程得到联系数
（4）计算改正数
（5）计算平差值
L L V
（6）精度评定（计算单位权方差、观测值中误 差、平差值函数的中误差等）
第三节 独立三角网条件平差
根据三角网中起算数据的多少，三角网有 独立三角网（网中仅有必要的起算数据）和非
独立三角网（网中具有多余的起算数据）之分。
三角网平差有按角度平差和按方向平差两种方
一、图形条件
图形条件通常是指平差后三角形内角应满足 的几何条件，所也称为三角形闭合条件。
1.平差值表达的图形条件
• 对于有 n 个三角形组成的中点多边形 (a) 或三角锁 (b)，可以列出 n 个图形条件：
ai bi ci 180 0 （i =1，2，… ，n）
a1
b1
ai
c1 ci
hi hi vi
代入(a)式得： 其中：
v1 v2 v3 w1 0 v1 v4 v6 w2 0, v2 v4 v5 w3 0
w1 h1 h2 h3 w2 h1 h4 h6 w3 h2 h4 h5
(b)
令：
1 A 1 0 1 0 1 1 0 0 0 1 1 0 0 1 0 w1 1 , W w2 w 0 3
• 为了确定一个几何模型就必须进行观测。如果 观测个数 n 少于必要元素的个数，即 n＜t，显 然无法确定该模型，出现了数据不足的情况； 若观测了 t 个独立量，n =t，则可唯一地确定 该模型。在这种情况下，如果观测结果中含有 错误，将无法发现。为了能及时发现错误，并 提高测量成果的精度，就必须使 n＞t，即必须 进行多余观测。多余观测的个数在测量中又称 “自由度”。令 r=n–t 显然， r 就是多余观测数。
• 测量平差通常是基于线性函数模型的， 当函数模型为非线性形式时，是将其用 泰勒公式展开，并取其一次项化为线性 形式。
• 对于一个实际平差问题，可建立不同形 式的函数模型，相应地就有不同的平差 方法。测量中常见的控制网平差方法有 条件平差和间接平差两种。
1、条件平差法
以观测量之间必须满足一定的条件方程为 函数模型的平差方法，称为条件平差法 。 例如：为了确定B、 C、D三点的高程， 其必要观测数 t =3， 实际观测了6 段高 差, 故多余观测数 r = n–t =3，应列出 3个线性无关条件 方程.
E
....... SOA sin a1 sin a2 sin an SOA sin b1 sin b2 sin bn
a2 ② c2 ci
b2 ai C ⓘ
bi D
即：
sin a1 sin a2 sin an 1 sin b1 sin b2 sin bn
h1 A B
h2
C h4 h5 h3 h6
D
这个水准网可以列出7个条件方程,其中只有 3个是相互独立的,我们取：
h1 h3 h2 0 h1 h6 h4 0 h2 h5 h4 0
式中：
（a）
hi
表示观测量 hi 的平差值。
这就是用平差值表达的条件方程。
由于平差值应该等于观测值与其改正数之和， 即：
第二节
条件平差原理
条件方程 AV +W=0 中，
A 为 r n 阶矩阵，
V 为 n 1 列阵， 即有 r 个方程，n 个未知数,且 r <n， 这样的方程组有无穷多组解。然而，根 据最小二乘准则，观测量的最或然值应 T 该满足V PV=min。
在 AV +W=0的条件下确定 VTPV 的最 小值，这在数学中是求函数Ф=VTPV 的条件极值问题。条件平差，实际上
w限 2m n
式中， m 为角度观测中误差； n 为圆周角的个数。
三、极条件
以中点多边形为例，若从OA边出发， 依次解算三角形①、②、…，最后解算 出的OA边长应与出发边OA相等。即：
A sin a1 SOB SOA a1 sin b1 sin a2 sin a1 sin a2 ① ⓝ SOC SOB SOA B b1 c1 O sin b2 sin b1 sin b2
L1 X 1 L2 X 2 L 3 X 1 X 2 180
A
X1 L1
L3
X2 L2
B
令： L1 L1 v1 L L 2 L2 v2 L V v L3 L 3 3 0 1 0 X1 B 0 1 , X , d 0 X 2 1 1 180
令其等于零，注意到 (PV )T = V T P，从而有： V T P =K T A 转置后左乘 P –1 得： V =P –1 ATK (1) 该公式表达了改正数 V 与联系数 K 的关系。
二、法方程式
将（1）式代入条件方程 AV +W=0 中得： AP –1 AT K+W=0 （ 2） 这就是条件平差的法方程式。式中，P为观测值 的权矩阵，设第 i 个观测值的权为 pi , 则
V = ( v1
v2
v3
v4
v5
v6)T （c）
则条件方程可表达为以下矩阵形式：
AV +W=0
这就是条件平差函数模型的一般形式。
条件方程
AV +W=0 中，
A －为r n 阶矩阵，称为系数矩阵；
V －为n 1列阵，称为改正数向量； W－为r 1列阵，称为闭合差向量。
2、间接平差法
• 一个几何模型中，只会有 t 个独立量，如果平
P
1
三、法方程的解
令 N = AP –1 AT （ 3）
则法方程式的形式为
N K+W =0 其中N 称为法方程式系数矩阵，是一个满 秩二次型方阵，其逆存在。从而可解出联系 数向量： K = －N –1 W （ 4）
四、条件平差的一般过程
（1）列出条件方程 （2）组成法方程系数矩阵 AV +W=0 N = AP –1 AT K = － N –1 W V=P –1 ATK
各角度值之和不等于360°的现象，平差时
除了要满足三角形闭合条件外，还必须使中
心点处的角度满足下列条件：
a1
c
i
360 0
b1
c1 ci ai bi
用改正数表达的圆周角条件为：
v
ci
wo 0 ； wo ci 360
其中，wo 称为圆周角条件闭合差。 对wo应用误差传播定律，并以2倍中误差作 为限差，则圆周角闭合差的限差为：