控制测量学水准网按条件平差算例(新)
在图 表9-1
试求： （１）1P 、2P 及3P 点高程之最或然值； （２）1P 、2P 点间平差后高差的中误差。
解:(1)列条件方程式,不符值以“mm ”为单位。
已知3,7==t n ，故437=-=r ，其条件方程式为
⎪⎪⎭
⎪
⎪
⎬⎫
=--+=-+--=-+--=++-01030707742643765521v v v v v v v v v v v v
（2）列函数式：
555v h x F +==
故 15=f 0764321======f f f f f f
（3）组成法方程式。
1）令每公里观测高差的权为1，按1/i i s p =,将条件方程系数及其与权倒数之乘积填于表9-2中。
2)由表9-2数字计算法方程系数,并组成法方程式:
⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----------5221251021411013⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡d c b a k k k k +⎥⎥⎥⎥⎦
⎤⎢⎢⎢⎢⎣⎡---1377=0 表9-2 条件方程系数表
（4）法方程式的解算。
1）解算法方程式在表9-3中进行。
2）[]pvv 计算之检核。
[][]wk pvv -= []467.35=-wk
由表9-3中解得[]47.35-=pvv ，两者完全一致，证明表中解算无误。
（5）计算观测值改正数及平差值见表9-4。
（6）计算321,,P P P 点高程最或然值。
359.3611=+=x H H A P m 012.3722=+=x H H A P m
表9-4 改正数与平差值计算表
（7）精度评定。
1）单位权（每公里观测高差）中误差
2）21,P P 点间平差后高差中误差 mm 0.34
47
.35±=±=μ。
控制测量学第九章控制网平差
r nt
3
4
L1
p1 0 0
5 6
记：观测值
L n1
L2
相应权阵
p
nn
0
p2
0
7 8
Ln
0
0
pn
9
10
4 /7
v1
平差值改正数
V
n1
v2
vn
平差值
Lˆ
n1
Lˆ1 Lˆ2
Lˆn
L1
L2
Ln
v1 v2
vn
4
控制网平差
1 2 3 4 5 6 7 8 9
10
13 /7 4
2）平差值函数的权倒数
设有平差值函数为 f Lˆ1, Lˆ2,, Lˆn
它的权函数式为：
d
(
Lˆ1
)dLˆ1
(
Lˆ2
)dLˆ2
(
Lˆn
)dLˆn
f1dLˆ1 f2dLˆ2 fndLˆn
令
f T f1, f2, fn
则 d f T dLˆ
9
2） 改正数条件方程
10
以 Lˆi L（i v=i 1、i 2、… ）代入n (1)得纯量形式为：
5 /7 4
a1v1 a2v2 anvn wa 0
b1v1 b2v2 bnvn wb 0
r1v1 r2v2 rnvn wr 0
（2）
1 2 3 4 5 6 7 8 9
控制网平差
第九章 控制网平差
1
2
3 4 5
[本章提要] 9.1 条件平差数学模型和公式
6
9.2 水准网按条件平差算例
7 8
9.3 附合导线按条件平差算例
水准网按条件平差算例
在图 表9-1试求： （１）1P 、2P 及3P 点高程之最或然值； （２）1P 、2P 点间平差后高差的中误差。
解:(1)列条件方程式,不符值以“mm ”为单位。
已知3,7==t n ，故437=-=r ，其条件方程式为⎪⎪⎭⎪⎪⎬⎫=--+=-+--=-+--=++-01030707742643765521v v v v v v v v v v v v（2）列函数式：555v h x F +==故 15=f 0764321======f f f f f f（3）组成法方程式。
1）令每公里观测高差的权为1，按1/i i s p =,将条件方程系数及其与权倒数之乘积填于表9-2中。
2)由表9-2数字计算法方程系数,并组成法方程式:⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----------5221251021411013⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡d c b a k k k k +⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---1377=0 表9-2 条件方程系数表（4）法方程式的解算。
1）解算法方程式在表9-3中进行。
2）[]pvv 计算之检核。
[][]wk pvv -= []467.35=-wk由表9-3中解得[]47.35-=pvv ，两者完全一致，证明表中解算无误。
（5）计算观测值改正数及平差值见表9-4。
（6）计算321,,P P P 点高程最或然值。
359.3611=+=x H H A P m 012.3722=+=x H H A P m表9-4 改正数与平差值计算表（7）精度评定。
1）单位权（每公里观测高差）中误差2）21,P P 点间平差后高差中误差mm 0.3447.35±=±=μmmP m FF 2.252.00.31±=±=±=μ。
水准网平差程序设计.ppt
2.2 数据存储方案设计
这里的数据存储是指数据在内存中的存储。
测量平差程序中用到的数据较多，一般先将磁盘
中的数据读到内存中，按照不同的类别有序地放
在变量或数组中，然后再进行平差计算，最后将 计算结果写到磁盘文件中。
三、水准网平差界面设计
要整洁、简练、实用、方便有关数据的输入 或输出，必要的话要考虑和用户的交互，尽量不 要画蛇添足。
第三章
水准网平差程序设计
本章难点： 1、近似高程计算 2、最短路线的计算
3、误差方程及法方程的构建
4、直接计算出法方程的系数矩阵BTPB和常 数矩阵BTPL
一、水准网间接平差算法概述
该课程中所采用的平差模型为间接平差，即所 选的独立参数的个数等于必要观测数，这样可以 将每个观测值表示成这t个参数的函数，组成观测 方程。 L B X d 间接平差的函数模型： n ,1 n ,t t ,1 n ,1 0 平差时一般对所选参数取近似值 X X x 代入上式 并令 l L ( BX 0 d ) L L0
1）基本信息部分：该部分仅占一行，其内容为
已知点数N1，未知点数N2，高差观测值个数NS。
2）已知点高程信息：该部分占N1行，每行格式为
已知点编号，该点的高程（单位：米）。 注意：在给控制点编号时，先给待定点编号，然后 给已知点编号。编号从1开始顺序编号。 3）高差观测值信息：该部分占NS行，每行格式为 测段编号，测段起点编号，测段终点编号，测段观 测高差（单位：米），测段长度（单位：千米）。 4）点名信息：该部分占N（N=N1+N2）行，每行 格式为 水准点编号，该点名称 其中点名长度不超过8个字符，即4个汉字。但最好 用字符串命名。
式中，正负号取决于高差起始点到终点的方向与 推算路线是否一致。
水准网按条件平差算例
§ 9.3 水准网按条件平差算例在图(9-5)所示水准网中，A,B两点高程及各观测高差和路线长度列于表(9-1)中。
•.一h s/ I \ ：'1 h5 「6A丈' \ 』4\ \ I\ \ x J、\rh2、丿P sJL ___ ■--P2 h7图9-5试求：(1) R、P2及P s点高程之最或然值；(2) P i、F2点间平差后高差的中误差。
解：(1)列条件方程式，不符值以“ mm”为单位。
已知n =7,t =3，故r =7 -3 =4，其条件方程式为w -V2 +V5 +7 =0-V5 -V6 ■ V7 _ 7 = 0 _V3 1V4 ■ V6 - 3 = 0 IV2 ■ V4 -V7 -1 = 0 I(2)列函数式：F = x5 = h5 V5故f5 -1 f^f^f^f^f^f^0(3 )组成法方程式。
1) 令每公里观测高差的权为1,按1/ P i =s，将条件方程系数及其与权倒数之乘积填于表9-2中。
2) 由表9-2数字计算法方程系数，并组成法方程式：表9-2条件方程系数表观测号 abcdsfF s1 111 2-113-1-1-1 4-1151-111 6-110 0 71-1z1-1-11" jg*1观测号%/P%dP %fP %1 1 1112 1-113 2-2-2-24 2-225 11-1116 1-11722-2z1-311（4）法方程式的解算。
1 ）解算法方程式在表 9-3中进行。
2） Ipvv 计算之检核。
pvv ] = -Wk I -Wk I-35.467由表9-3中解得Pvvl--35.47，两者完全一致，证明表中解算无误。
（5） 计算观测值改正数及平差值见表 9-4。
（6） 计算R,P 2,P 3点高程最或然值。
H P l= H A X ! =36.359 mH p 2 =H A x 2 =37.012 m-1 0 -14 -1 -2 -15 -2 -2 -2 5_7 -7 —3=0k a k b k cR 36 =H B+X4 =35.360 m5 4A(7) 精度评定。
高程控制网平差
i
i
i
h h V 改厕厕短的改正数， 代入上式，得：
i
i
i
V1 V2 V3 V4 W 0
W H A h1 h2 h3 h4 H B
1.附合水准路线的条件数和条件方程式组成
观测值5个，待定水准点2 个，所以条件有3个，可 以列出3个条件方程：
h1
H B h1 h2 H A 0
V 1 V 3 V 2 W a 0 V 2 V 4 V 6 W b 0 V 4 V 5 V 3 W c 0
(二)观测值权的确定：
1.各水准路线都进行了往返观测，每公里水准路线的观测中误差为 ，
则m：i
R mi2
1 4n
n i
2 i
i
式中，为测往返测高程不符值，以mm为单位；R为测段长度，以km为单位；n
H A h2 h3 h5 H D 0
H B h1 h3 h4 H C 0
一般以1个已知点为起点，其它已知点为终点，所构成的附合 水准路线为已知点数减1，这样可以列出的条件方程式为已知 水准点个数减1.
2.闭合水准路线的条件数和条件方程式的组成
从一个水准点出发，经过若干水准测段，又回到该 水准点，这样的水准路线称为闭合水准路线。
V 1 V 7 V 8 W b 0
V 2 V 8 V 7 W c 0
V 3 V 5 V 8 W d 0
V 4 V 6 V 5 W e 0
2.闭合水准路线的条件数和条件方程式的组 成
图(c)是四边形状水准网，网中有4个待定点，没有已知点， 在平差计算时，只能确定个待定水准点之间的相互关系，如 果确定一个水准点的高程，就可以确定其他点的高程。因此， 该网的必要观测是3个，观测值总数是6个，又3个多余观测， 可以列出3个条件方程。为了让所列立的条件方程式互相独 立，没个条件方程都要求有一个其他方程没有用到的观测值， 即：
水准网的条件平差
目录目录 (1)观测误差 (2)摘要： (2)关键词： (2)引言 (3)1水准测量 (4)1.1水准测量的原理 .............................................................................. •1.2水准网...................................................................................... •2条件平差 .. (6)2.1 衡量精度的指标 (6)2.2条件平差的原理 .............................................................................. •3水准网的平差 .. (14)3.1必要观测与多余观测 (14)3.2条件方程 (14)3.3条件平差法方程式 (14)3.4条件平差的精度评定 (15)3.5水准网的条件平差 (18)致谢 (20)参考文献 (21)观测误差由观测者、外界环境引起的偶然误差学生：xxx 指导教师：xxx摘要:对一系列带有偶然误差的观测值，采用合理的的方法消除它们间的不符值，得出未知量的最可靠值；以及评定测量成果的精度。
关键词:偶然误差；观测值；精度引言测量工作中，要确定地面点的空间位置，就必须进行高程测量，确定地面点的高程。
几何水准测量是高程测量中最基本、最精密的一种方法。
通过测量仪器，工具等任何手段获得的以数字形式表示的空间信息，即观测量。
然而，测量是一个有变化的过程，受仪器、观测值、外界环境因素的影响，观测的结果与客观上存在的一个能反映其真正大小的数值，即真值（理论值），有一定的差异。
可以说在测量中产生误差是不可避免的。
所以，观测值不能准确得到，在测量上称这种差异为观测误差。
根据其对观测结果影响的性质，可将误差分为系统误差和偶然误差两种。
水准网的条件平差
目录目录 (1)观测误差 (2)摘要: (2)关键词: (2)引言 (3)1水准测量 (4)1.1水准测量的原理 (4)1.2水准网 (5)2条件平差 (6)2.1衡量精度的指标 (6)2.2条件平差的原理 (8)3水准网的平差 (14)3.1必要观测与多余观测 (14)3.2条件方程 (14)3.3条件平差法方程式 (14)3.4条件平差的精度评定 (15)3.5水准网的条件平差 (18)致 (21)参考文献 (21)观测误差—由观测者、外界环境引起的偶然误差学生： xxx 指导教师：xxx摘要:对一系列带有偶然误差的观测值，采用合理的的方法消除它们间的不符值，得出未知量的最可靠值；以及评定测量成果的精度。
关键词:偶然误差；观测值；精度引言测量工作中，要确定地面点的空间位置，就必须进行高程测量，确定地面点的高程。
几何水准测量是高程测量中最基本、最精密的一种方法。
通过测量仪器，工具等任何手段获得的以数字形式表示的空间信息，即观测量。
然而，测量是一个有变化的过程，受仪器、观测值、外界环境因素的影响，观测的结果与客观上存在的一个能反映其真正大小的数值，即真值（理论值），有一定的差异。
可以说在测量中产生误差是不可避免的。
所以，观测值不能准确得到，在测量上称这种差异为观测误差。
根据其对观测结果影响的性质，可将误差分为系统误差和偶然误差两种。
前者可以通过在观测过程中采取一定的措施和在观测结果中加入改正数，消除或减弱它的影响，使其达到忽略不计的程度。
但是，观测结果中，不可避免地包含了后者，它是不可消除的，但可以选择较好的观测条件或采用适当的数据处理方法减弱它。
现在我们要讨论的就是采用适当的数据处理方法来减弱其对水准测量中的影响。
1 水准测量1.1水准测量的原理1.1.1 水准测量的基本原理水准测量是利用水准仪提供的水平视线在水准尺上读数，直接测定店面上两点的高差，然后根据已知点高程及测得的高差来推算待定点的高程。
(完整word版)控制测量学水准网按间接平差算例
+58.00
+17.00
-1.8889
-5.00
-1.22
+0.1414
+3.00
+6.45
-0.8206
+58.00
35.45
+1
-0.1111
-1
-0.78
+0.0911
0
0
0
0
-0.18
+6
-0.6667
+4
+5.34
-0.6239
+5
+7.86
-1.0000
0
+6
-0.6667
1.114
1.142
0.078
0.099
1.216
5
5
5
2
2
2
图9-10
（2）列误差方程，得
式中
mm mm
mm mm
mm mm
（3）组成法方程，将误差方程系数、常数项以及权填入表9-10。
表9-10
编号
（mm）
-1.0918
0.9074
0.1794
1
2
3
4
5
6
-1
0
-1
+1
0
0
+1
-1
0
0
+1
0
0
+1
+4
+5.33
-0.6238
+5
+7.86
-1.0000
0
-0.18
-1
+0.1111
[整理]10控制网平差计算
§9.1 条件平差原理在条件观测平差中，以n 个观测值的平差值1ˆ⨯n L作为未知数，列出v 个未知数的条件式，在min =PV V T 情况下，用条件极值的方法求出一组v 值，进而求出平差值。
9.1.1基础方程和它的解设某平差问题，有n 个带有相互独立的正态随机误差的观测值 ，其相应的权阵为 ， 它是对角阵，改正数为 ，平差值为 。
当有r 个多余观测时，则平差值 应满足r 个平差值条件方程为：⎪⎪⎭⎪⎪⎬⎫=++++=++++=++++0ˆˆˆ0ˆˆˆ0ˆˆˆ221122112211οοοr L r L r L r b L b L b L b a L a L a L a n n n n n n （9-1） 式中i a 、i b 、…i r （i =1、2、…n ）——为条件方程的系数；0a 、0b 、…0r ——为条件方程的常项数以ii i v L L +=ˆ（i =1、2、…n ）代入(9-1)得条件方程（9-2）式中a w 、b w 、……r w 为条件方程的闭合差，或称为条件方程的不符值，即（9-3） 令⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⨯n n n n r r r r b b b a a a A212121⎪⎪⎭⎪⎪⎬⎫++⋅⋅⋅++=++⋅⋅⋅++=++++=022110221102211r L r L r L r w b L b L b L b w a L a L a L a w n n n n n b n n a ⎪⎪⎭⎪⎪⎬⎫=++⋅⋅⋅++=++⋅⋅⋅++=++⋅⋅⋅++000221122112211r n n b n n a n n w v r v r v r w v b v b v b w v a v a v a ⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⨯n n L L L L 211⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⨯n n L L L L ˆˆˆˆ211⎪⎪⎪⎪⎫ ⎛=⨯b a r w w W 1⎪⎪⎪⎪⎫ ⎛=⨯n v v V 211⎪⎪⎪⎫⎛=⨯οοb a A o r 11⨯n L nn P ⨯1⨯n V 1ˆ⨯n L 1ˆ⨯n L⎪⎪⎪⎫⎛=⨯n n p p P 000021则(9-1)及（9-2）上两式的矩阵表达式为0ˆ0=+A LA （9-4） 0=+W AV （9-5）上改正数条件方程式中V 的解不是唯一的解，根据最小二乘原理，在V 的无穷多组解中，取PV V T = 最小的一组解是唯一的，V 的这一组解，可用拉格朗日乘数法解出。
条件平差算例
一、水准网条件平差示例 范例：有一水准网（如图8-3所示），已知点A ，B 的高程为: HA=50.000m ， HB=40.000 m ，观测高差及路线长度见表8-1。
试用条件平差求:(1) 各观测高差的平差值；(2) 平差后P 1到P 2点间高差的中误差。
图8-3【解】1）、求条件方程个数；由图易知：n=7，t=3，条件式r=4。
故应列4个平差值条件方程，三个闭合环，一个附和路线2）、列平差值条件方程； 所列4个平差值条件方程为：⎪⎪⎭⎪⎪⎬⎫=-+-=--=-+=+-0ˆˆ0ˆˆˆ0ˆˆˆ0ˆˆˆ31643765521BA H H h h h h h h h h h h h 3)、转换成改正数条件方程；以ii i V L L +=ˆ代入上式可得： ⎪⎪⎭⎪⎪⎬⎫=-+-+-=--+--=-++-+=+-++-00003131643643765765521521B A H H h h v v h h h v v v h h h v v v h h h v v v 化简可得：⎪⎪⎭⎪⎪⎬⎫=--=+--=+-+=++-0403070731643765521mm mm mm mm v v v v v v v v v v v 可知条件方程系数阵为：⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-----000101010110011100000010011⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=2101001000210000210000010000001称对P ⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=2010010002000020000010000001称对Q ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡-----⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-----=320125100141101300100110001101001100000110010002010102200211000000100114)、组成法方程； 先求权阵P ；以1km 观测高差为单位权观测高差，则： 11=P ，12=P ，213=P ，214=P ，15=P ，16=P ，217=P ，而各观测高差两两相互独立，所以权阵为：，则协因数阵为：则，法方程的系数阵Naa 为：⎥⎥⎦⎤⎢⎢⎣⎡-----⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-----==-=00010101011001110000001001120100100020000200000100000010001010101100111000000100111TT AQA T A AP aa N 称对所以，法方程为：043773212510014110134321=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----k k k k 5）、解算法方程，求出联系数K⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--------=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎦⎤⎢⎢⎣⎡34831.213483.002247.177528.2437758427.025843.012360.023596.025843.032584.011236.012360.012360.011236.031461.014607.023596.012360.014608.046067.04377320125100141101314321k k k k 6）、求V 及高差平差值Lˆ 所以4210.212.118.3213.0214.418.214.0ˆ22222220⨯+⨯-+⨯-+⨯-+⨯-+⨯+⨯-==）（）（）（）（）（r PV V T σ⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡⎥⎥⎦⎤⎢⎢⎣⎡-----=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡-----=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-------⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡==045.2157.1798.3270.0427.4775.2427.034831.213483.002247.177528.2002001100011020022000001100134831.213483.002247.177528.200001010101100111000000100112010010002000020000010000001m m T K T QA V 称对mmmm v v v v v v v h h h h h h h h h h h h h h L ⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡-----+⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=5020.108548.56472.45007.143556.200028.153556.100.22.18.33.04.48.24.0500.10856.5651.4501.14360.20000.15356.10ˆˆˆˆˆˆˆˆ7654321765432176543217）、精度评定1）、单位权方差估值计算mm 98.24605.35±==2）、建立所求精度的平差值函数的算式，并按误差传播律求平差值函数的精度 依题意列平差值函数为： 5ˆh =ϕ 则:[]Tf 0010000=[][][][]51687.048313.01)16853.3146.0(1001111236.001124.016853.03146.0100110011111ˆˆ=-=+-=⨯---=-=-=--TTT T T aaaa N AQf N QA f Qf fQ ϕϕ所以：mm Q 14.251687.098.2ˆˆ0ˆ±=⨯==ϕϕϕσσ【答】：各观测高差的平差值为：}{m m m m m m m5020.108548.56472.45007.143556.200028.153556.10平差后P1到P2点间高差的中误差为：±2.14mm987654321ACPB 图8-11二、测角网条件平差 范例：有一测角网（如图8-11所示），A 、B 、C 三点为已知三角点，P 为待定点。