第五章 条件平差
第5章附有条件的条件平差
按求条件极值的方法组成新的函数: 按求条件极值的方法组成新的函数: T ˆ ˆ Φ = V T PV − 2 K T ( AV + Bx − W ) − 2 K S (Cx − W x ) 分别对
V
和
ˆ x
求一阶偏导数并令一阶偏导数为零,得 求一阶偏导数并令一阶偏导数为零,
n×n n×1
∂Φ = 2V T P − 2 K T A = 0 ∂V
− − N bb = B T N aa1 B We = B T N aa1W u×u
u ×1
ˆ N bb x − C T K s − We = 0
− ˆ x = N bb1 (C T K s + We )
(d )
6
第五章 附有条件的条件平差
§5-2 精度评定
一、单位权方差估值的计算公式
V T PV V T PV ˆ σ = = r c−u + s
u×c c×1 u ×s s×1
基础方 程
B T K+ C T K S= 0
u ×c c×1
u× s
s×u u ×1
ˆ C x − Wx = 0
s×1
s×1
由(3 )得: 改正数方 程 V = P −1 AT K = QAT K
n×1
法方程
法方程的矩阵形式: 法方程的矩阵形式:
代入(1): 代入(
c×n n×n n×c
2010-11-15
− ˆ V T PV = W T N aa1W − WeT x + W xT K s
7
第五章 附有条件的条件平差
§5-2 精度评定
二、各种向量的协因数阵
ˆ ˆ 基本向量: L,W,X,K,K s,V,L 基本向量:
第五章条件平差
令
a1 b A= 1 ⋅⋅⋅ r 1 a2 b2 ⋅⋅⋅ r2 ⋅⋅⋅ ⋅⋅⋅ ⋅⋅⋅ ⋅⋅⋅ an bn ⋅⋅⋅ rn
v1 Wa Wb v W = ,V = 2 M M W v r n
第五章
条件平差
第一节
条件平差原理 2
第二节 第三节
条件方程 16 精度评定 39
第四节
条件平差公式汇编和水准网平差示例
第一节
条件平差的函数模: 条件平差的函数模:
r ,n n,1
条件平差原理
或
A∆ −W = 0
~ A L+ A0 = O
r ,1 r ,1
~ ˆ 当 L 的估值为 L ,∆的估值为V时,则有 的估值为 时
v L an 1 wa v L bn 2 + w b = 0 M L cn wc v n
组成新函数: 组成新函数:
K与方程个数相同 与方程个数相同
2 2 Φ = [ p 1 v 12 + p 2 v 2 + L + p n v n ]− 2k a ( [av ] + w a )
一、条件平差原理
设有r个 设有 个平差值线性条件方程
ˆ ˆ ˆ a1L1 + a2 L2 +L+ anLn + a0 = 0 ˆ + b L +L+ b L + b = 0 ˆ ˆ b1L1 2 2 n n 0 L L L L L L L ˆ ˆ ˆ r1L1 + r2 L2 +L+ rnLn + r0 = 0
第五章 条件平差
cot
a3
va3
sin a1 sin b1
sin sin
a2 b2
sin sin
a3 b3
cot
b1
vb1
sin a1 sin b1
sin sin
a2 b2
sin sin
a3 b3
cot
b2
vb2
sin a1 sin b1
sin sin
a2 b2
sin a3 sin b3
cot
b3
vb3
0
cot a1va1 cot a2va2 cot a3va3 cot b1vb1 cot b2vb2 cot b3vb3
§5-2 条件方程
4.水准网中条件方程的列立方法
• 列条件方程的原则:足数、独立、最简 • (1)先列附合条件,再列闭合条件 • (2)附合条件按测段少的路线列立,附合条件的个
数等于已知点的个数减1 • (3)闭合条件按小环建立(保证最简),一个水准
网中有多少小环,就列多少个闭合条件。
§5-2 条件方程
• (2)极条件
A a4 b4
a1 B b1
AB AC AD 1 0
AC AD AB
D b3a3
sin aˆ2 sin(aˆ3 bˆ3) sin aˆ1 sin(aˆ1 bˆ1) sin bˆ2 sin bˆ3
1
0
b2a2 C
cot a2va2 cot(a3 b3 ) (va3 vb3 ) cot a1va1
• 原则:将复杂图形分解成典型图形
三角形
大地四边形
中心多边形
§5-2 条件方程
7.条件方程的类型
• 图形条件(内角和条件):三角形内角和等于180 • 圆周条件(水平条件):圆周角等于360 • 极条件(边长条件):由不同推算路线得到的同一
条件平差原理
§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 ˆˆˆˆ2111⨯n L nn P ⨯1⨯n V 1ˆ⨯n L 1ˆ⨯n L则(9-1)及(9-2)上两式的矩阵表达式为0ˆ0=+A LA (9-4) 0=+W AV (9-5)上改正数条件方程式中V 的解不是唯一的解,根据最小二乘原理,在V 的无穷多组解中,取PV V T = 最小的一组解是唯一的,V 的这一组解,可用拉格朗日乘数法解出。
测量平差第五章
1
1
v1
AP1AT K W 0
1v2
9
0
4Ka 12 0
vv43
Ka 3
Lˆ1 Lˆ2
Lˆ3 Lˆ4
L1
L2
L3 L4
v1 v2 v3 v4
A)dSb
(Sc
Sb
cos
A)dSc
]
由图5-10知: SbSc sin A Sbhb (2倍三角形面积) Saha
Sb Sc cos A Sa cos C,Sc Sb cos A Sa cos B
代入上式,得: dA
1 ha
(dSa
cos CdSb
cos
BdS c )
1
1
1
Lˆ1 Lˆ2 Lˆ3 Lˆ4
360
0
W AL A0 [1 1
573216
1
1]
730308
360
12
1265128
1023320
③计算改正数和平差值 1
一、测角网条件方程
圆周条件只有一个!
图形条件不只列 三个,但独立的 只有三个!
§5.2 条件方程
各边均与D有关,即以D为极,故称为极条件。 在中点多边形测角网中,ai、bi、ci通常按顺时针编号,且ci以 中点为顶点,称为圆周角或间隔角;ai称为求距角;bi称为传 距角。这样,在列极条件时规律性很强!
条件平差
得法方程: AQATK-W=0 T 1 T N AQA AP A 令 aa r .r r .nn.nn.r 则有: NaaK-W=0
法方程系数阵Naa是一个r阶的满秩方阵,且可逆
N11k1 N12k 2 N1r k r W1 0 N 21k1 N 22k 2 N 2 r k r W2 0 N r1k1 N r 2 k 2 N rn k r Wr 0
目标函数:f x min n1 x a h x min F a , x f 1k 约束条件: h x 0 k 1 n1 F a, x
0 a F a, x 0 x
L2
L4 L1 L3 L2
A
B
C
§6-2 条件方程
条件方程的个数等于多余观测数r。条件方程之间 不能线性相关,在一个平差问题中,条件方程的个 数是固定不变的.
一、r的确定: r=n-t 二、条件方程的列立: 原则:足数(r个),线性无关,形式简单,易 于列立
控制网常见几何模型
水准网 三角网(测角网) 三边网(测边网) GPS基线向量网 单一附合导线
由此可得联系数K的解:
r ,1
K ( AQA ) W
T
T 1
V QA K
条件平差的 最小二乘解:
n,1
ˆ L V L
三、条件平差计算步骤:
1.根据平差问题的具体情况,列出条件方程,条 件方程的个数等于多余观测数r。 2.组成法方程式,法方程的个数等于多余观测数r 3.解法方程,求出联系数K值。 4.将K代入改正数方程式,求出V值,并求出观测 值的平差值=L+V。 5.检验平差计算的正确性(可用平差值重新列出 平差值条件方程式,看其是否满足方程)。
第五章条件平差
二、法方程及改正数方程
将V T PV min的原则作用于条件方程 。
组成新函数:
V T PV-2k T AV W
式中
r 1
k k a , kb , k r 条件方程联系数
T
对新函数求导: T T 2V P 2A k ---改正数方程
dSCD ˆ f T dL SCD ˆ SCD T 2 T ˆ f D f f QL ˆL ˆ ˆL ˆ f 0 L S CD
得测边相对中误差为: 3、大地四边形测角网
2
ˆS
CD
SCD
=
ˆ 0 f T QL ˆL ˆ f
设
F ( f1 , f 2 , f m )
T T
G ( g1 , g 2 , g m ) 有
均为m维向量函数,且 f i、g i 均为x的函数, d F G dG F T dG T dF F G dx dx dx dx
注意:当N为满秩方阵时,才有 N 1唯一存在,法方程才有唯
测方向网
测角网
测角网
三角网
测边网
测边长
测边+测方向
边角网
(导线网) 测边+测角
三、三角网的布设--从高级到低级逐级布设 四、三角网平差的方法 1。严密平差 ----遵守VTPV=min原则 ; 2。近似平差
5.3 测角网条件平差
独立网(经典自由网)---只有必要起算数据d。
非独立网(附合网)---已知条件超过必要起算数据。
3 图形条件: n=12 t=2×2+4=8 r =4 1 极条件:
v2 v1 v6 v5 v11 v10 W1 0
测量程序设计_条件平差和间接平差
程序代码如下:
disp(‘-------水准网间接平差示例-------------’) disp(‘已知高程’) Ha = 5.015 % 已知点高程,单位m Hb = 6.016 % 已知点高程,单位m
A h2 D h1
C h6 E h7 B h4
h5
h3
disp(‘观测高差,单位m’)
L = [1.359; 2.009; 0.363; 1.012; 0.657; -0.357] disp(‘系数矩阵B’)
则: PV AT K
V P A K QA K
T
1 T
4、法方程: 将条件方程 AV+W=0代入到改正数方程V=QATK 中,则得到:
AQAT K W 0
r1 r1 r1
记作: 由于
N aa K W 0
rr
R( Naa ) R( AQAT ) R( A) r
Naa为满秩方阵, K Naa1W ( AQAT )1 ( AL A0 )
if H(1,1)+H(2,1)-H(3,1)+HA-HB==0 && H(2,1)H(4,1)==0 disp(‘检核正确') else disp(‘检核错误') end disp(‘平差后的高程值') HC = HA + H(1,1) HD = HA + H(1,1) + H(4,1)
二、间接平差的基本原理
其中l=L-d.
ˆ 设误差Δ和参数X的估计值分别为V 和 X
则有
ˆ V AX l
X0 为了便于计算,通常给参数估计一个充分接近的近似值
ˆ ˆ X X0 x
则误差方程表示为
测量平差 第五章 条件平差
北京建筑工程学院 测绘工程系
求解法方程,求的联系数K
N aa K + W = 0
− K = − ( A P − 1 A T ) − 1 W = − N aa1W
回代求解观测值改正数
V = P −1 A T K = QA T K
观测值平差值
ˆ L = L +V
误差理论与测量平差基础
北京建筑工程学院 测绘工程系
四、条件平差的计算步骤
1. 首先确定条件方程的个数 r=n-t
ˆ 列出平差值条件方程式, AL + Ao = 0
列出改正数条件方程 AV + W = 0 定观测值的权阵P 2. 组成法方程式 4. 求改正数V值 5. 求出平差值 6. 检核
误差理论与测量平差基础
( AP −1 AT ) r×r K r×1 + W = 0
误差理论与测量平差基础
北京建筑工程学院 测绘工程系
水准网平差例2 A
h1
t = 7 −1− 3 = 3
C
h3
r = n−t = 7−3= 4
E
h2 h6 h5
符合条件方程
G
h4
h7
F D
B
ˆ ˆ ˆ h1 + h2 − h3 − ( H C − H A ) = 0 ⎫ ⎪ ˆ + h − h −( H − H ) = 0⎪ ˆ ˆ h1 6 7 B A ˆ + h + h − ( H − H ) = 0 ⎬ 闭合条件方 h7 ˆ 5 ˆ 4 ⎪ D B ⎪ ˆ ˆ ˆ h2 − h5 − h6 = 0 ⎭
矩阵计算基础知识 1. 向量 矩阵 2. 矩阵转置 3. 矩阵相乘 4. 矩阵微分 5. 矩阵求逆 6. 特殊矩阵 7. Matlab 矩阵计算
误差理论与测量平差基础第五章条件平差ppt课件.pptx
5-2 条件方程的列立
故有:
dA
1 ha
(dSa
cos CdS b
cos BdSc
)
将微分换成改正数,并将弧度换
成角度,得:
vA
ha (vSa
cos CvSb
cos BvSc
)
上式称为角度改正数方程。它具有明显的规律:
任意角度的改正数,等于其对边的改正数分别减去两邻 边的改正数乘以其邻角的余弦,然后再除以该角至其对边的
3、几种非线性条件方程的线性形式
极条件: 在图5-4中,极条件为 线性化得:
sin aˆ1 sin aˆ2 sin aˆ3 sin bˆ1 sin bˆ2 sinbˆ3
1
sin(a1 va1 )sin(a2 sin(b1 vb1 )sin(b2
va2 )sin(a3 va3 ) vb2 )sin(b3 vb3 )
dV
dV
dV
VTP VTP
2V T P
5-1 条件平差原理
2.2 求偏导
2.3 法方程 改正数方程
d 2V T P 2K T A 0 dV
AP1 AT K W 0
V P1 AT K
举例
水准网如右图:观测值及其权阵如下:
L 0.023 1.114 1.142 0.078 0.099 1.216 T m
m1
yA yˆi yB 0 i 1
5-2 条件方程的列立
➢GIS数字化数据采集中,折角均为90度的N边形的条件 方程
1、观测值
观测值为N个顶点的坐标,其个数为n=2 N。
2、必要观测个数
t=N+1
h
3、多余观测个数
r=n-t=2N-N-1=N-1 4、条件方程的类型
6第五章条件平差资料
l2
A
B
C
22
解:(1)条件方程:
✓此例n=4,t=2,故r=2;可列两个条件方程:
lˆ 1
lˆ2 lˆ2
lˆ3
lˆ4
0
0
lˆi li vi
v 1
v2
v4
3
0
v2 v3 2 0
✓写成矩阵形式:
10
3
3 6
ka
kb
3 2
0
(2)定权:✓100米量距为单位权:Pi=100/Si
F
39
§2 条件方程
sin sin
L3 L2
sin sin
L5 L4
sin sin
L1 L6
cos L1
sin L1
v1
sin sin
L1 L4
sin sin
L3 L6
sin sin
L5 L2
cos L2
sin L2
v2
sin L1 sin sin L2 sin
L3 sin L5 L4 sin L6
cos L6
sin L6
v6
( sin L1 sin L2
sin sin
L3 sin L5 L4 sin L6
1)
1
sin L1 sin L2
sin sin
L3 sin L5 L4 sin L6
(cot
L1v1
cot
L2v2
cot
L3v3
cot
L4v4
cot
L5v5
cot
L6v6 )
( sin L1 sin sin L2 sin
② 求L(x,y)对x和y的一阶偏导数,令它们等于零,并 与附加条件联立,即 L'x(x,y)=ƒ'x(x,y)+λφ'x(x,y)=0, L'y(x,y)=ƒ'y(x,y)+λφ'y(x,y)=0, φ(x,y)=0
5条件平差4
1、数学模型:
AV W 0
2 2 D 0 Q 0 P 1
关于矩阵导数
2、平差准则:
V T PV min
3、按求函数极值的拉格朗日乘数法,构造 新的函数
V T PV 2K T ( AV W ) min
4、将对V求一阶导数,并令其等于零
4、将对V求一阶导数,并令其等于零,得:
例:下图测角网,n=35
t=?
r=?
类型?
∆
∆
非自由网条件
1.边长条件
ˆ ˆ ˆ ˆ S AB sin L1 sin L11 sin L4 sin L5 SCD 0 ˆ sin L sin L sin L ˆ ˆ ˆ sin L2 10 8 6
ˆ ˆ ˆ ˆ ˆ ˆ 2.方位角条件 L1 L2 L3 L4 L5 L6 AB DC n 180o 0
授课内容: 5-1 条件平差原理 5-2 条件方程 5-3 精度评定 5-4 水准网平差示例
§5-1 条件平差原理
一、算例
ˆ ˆ ˆ L1 L2 L3 1800 0
v1 v2 v3 w 0
[v ] v v v min
2 2 1 2 2 2 3 2 2 G v12 v2 v3 2k (v1 v2 v3 w)
§5-2 条件方程
一、水准网条件方程
自由网条件
非自由网条件
例题
二、平面测角网条件方程
基本图形
三角网是由以下三种基本图形构成的:
1、单三角形 2、大地四边形 3、中点多边形。
自由网条件
ˆ ˆ ˆ 1.图形条件(三角形闭合条件) Lai Lbi Lci 180 o 0 vai vbi vci w 0
误差理论和测量平差习题5(含答案)
第五章条件平差习题第五章思考题参考答案5.1(a)n=6,t=3,r=3(b)n=6,t=3,r=3(c)n=14,t=5,r=95.2(a)n=13,t=6,r=7共有7个条件方程,其中有5个图形条件,2个极条件。
(b)n=14,t=8,r=6共有6个条件方程,其中有3个图形条件,3个极条件。
(c)n=16,t=8,r=8共有8个条件方程,其中有6个图形条件,2个极条件。
(d)n=12,t=6,r=6共有6个条件方程,其中有4个图形条件,1个圆周条件,1个极条件。
5.3n=23,t=6,r=17共有17个条件方程,其中有9个图形条件,1个圆周条件,1个固定角条件,1个固定边条件,5个极条件。
5.4 (1)n=22,t=9,r=13:7个图形条件,1个圆周条件,2个极条件,2个边长条件,一个基线条件。
(2)12837941314121520111718195610166101119910111213510ˆˆˆ1800ˆˆˆ1800ˆˆˆ1800ˆˆˆ1800ˆˆˆˆ1800ˆˆˆˆ1800ˆˆˆˆ1800ˆˆˆˆˆ1800ˆˆˆsin sin sin L L L L L L L L L L L L L L L L L L L L L L L L L L L L L L L L ++-=++-=++-=++-=+++-=+++-=+++-=++++-=171961116203614184715192211151217121318124ˆsin 1()ˆˆˆˆsin sin sin sin ˆˆˆˆsin sin sin sin 1()ˆˆˆˆsin sin sin sin ˆˆ()ˆˆˆˆsin sin sin sin ˆˆ(ˆˆˆˆsin sin sin sin FG FG L L L L L L L L L L L L L S S S S L L L L S S L L L L ===→=以大地四边形中心为极以中点四边形D 点为极的边长条件1213611891719ˆˆ)ˆˆˆˆsin sin sin sin ˆˆˆˆsin sin sin sin FG AB S S L L L L S S L L L L →=的边长条件(基线条件)5.5 n=8,t=4,r=4;有多种条件方程的列法,其中之一为:1001000100110000120001001104000011014V ⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥-=⎢⎥⎢⎥--⎢⎥⎢⎥---⎣⎦⎣⎦(注意常数项单位为mm ) 5.6 (1)P=3/2,(2)P=15.7 (1)P B =1.6,P C =2.1,P D =2.1,P E =1.6(2)P hCD =1.85.8 []ˆ 2.4998 1.9998 1.3518 1.8515h=2P σ=0.32(mm)5.9 1234561110009100110900101016V V V V V V ⎡⎤⎢⎥⎢⎥-⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥+=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦⎢⎥⎢⎥⎢⎥⎣⎦ []045452TV mm =---[]ˆ 1.576 2.219 3.7950.867 2.443 1.352T h m =--- 5.10 (1)1ˆ10.3556h m = 2ˆ15.0028h m = 3ˆ20.3556h m = 4ˆ14.5008h m =5ˆ 4.6472h m = 6ˆ 5.8548h m = 7ˆ10.5020h m =(2)±2.2mm。
误差理论与测量平差基础习题集
第五章条件平差§5-1条件平差原理条件平差中求解的未知量是什么?能否由条件方程直接求得5. 1. 02 设某一平差问题的观测个数为n.必要观测数为t,若按条件平差法进行平差,其条件方程、法方程及改正数方程的个数各为多少?5. 试用符号写出按条件平差法平差时,单一附合水准路线中(如图5-1所示)各观测值平差值的表达式。
图5-15. 1. 04 在图5-2中,已知A ,B的高程为Ha = m , Hb=11. 123m,观测高差和线路长度为:图5-2S1=2km,S2=Ikm,S3=,h1=,h2= m,h3= m,求改正数条件方程和各段离差的平差值。
在图5-3的水准网中,A为已知点B、C、D为待定点,已知点高程HA=,观测了5条路线的高差:h1=,h2=0. 821 m,h3=,h4=,h5= m。
各观测路线长度相等,试求:(1)改正数条件方程;(2)各段高差改正数及平差值。
有水准网如图5-4所示,其中A、B、C三点高程未知,现在其间进行了水准测量,测得高差及水准路线长度为h1=1 .335 m,S1=2 km;h2= m,S2=2 km;h3= m,S3=3km。
试按条件平差法求各高差的平差值。
如图 5-5 所示,L1=63°19′40″,=30″;L2=58°25′20″,=20″;L3=301°45′42″,=10″.(1)列出改正数条件方程;(2)试用条件平差法求∠C的平差值(注: ∠C是指内角)。
5-2条件方程5. 对某一平差问题,其条件方程的个数和形式是否惟一?列立条件方程时要注意哪些问题?如何使得一组条件方程彼此线性无关?. 10 指出图5-6中各水准网条件方程的个数(水准网中P i表示待定高程点,h i表示观测高差)。
(a) (b)图5-65. 2. 11指出图5-7中各测角网按条件平差时条件方程的总数及各类条件的个数(图中P i 为待定坐标点)。
第5章 附有限制条件的条件平差
测量平差太原理工大学测绘科学与技术系第五章附有限制条件的条件平差附有限制条件的条件平差§5-1 基础方程和它的解§5-2 精度评定§5-3 各种平差方法的共性和特性§5-4 平差结果的统计性质§5-1 基础方程和它的解条件平差、附有参数的条件平差、间接平差、附有条件的间接平差等四种经典平差方法,除条件平差不增选参数外,其它三种方法都要增选数量不等的参数参与平差,其未知参数的个数分别是u<t,u=t,u>t,且要求参数间彼此独立,在u>t 的情况下,也要求必须包含t个独立参数,从函数模型上看,四种平差方法总共包含如下四类的方程:基础方程和它的解前三类方程中都含有观测量或同时含有观测量和未知参数,而最后一种方程则只含有未知参数而无观测量,为了便于区别起见,特将前三类方程统称为一般条件方程,而最后一类条件方程称为限制条件方程。
~0)~(0=+=A L A L F ,线性形式为:dX B L X F L +==~~)~(~,线性形式为:0~~0)~,~(0=++=A X B L A X L F ,线性形式为:0~0)~(0=+=ΦC X C X ,线性形式为:基础方程和它的解在第二章中介绍过附有条件的条件平差的模型建立方法,该方法也要增选u 个参数,方程的总数为r+u 个。
如果在u 个参数中有s 个是不独立的,或者说在这u 个参数中存在着s 个函数关系式,则建立平差模型时应列出s 个限制条件方程,除此之外再列出c=r+u-s 个一般条件方程,因此方程总数也可以认为是c+s 个,形成如下的函数模型若为线性形式,则为0)~,~(1=⨯X L F c 0)~(1=Φ⨯X S 0~~1011=++⨯⨯⨯⨯⨯c u u c n n c A X B L A 0~101=+⨯⨯⨯s u u s C X C基础方程和它的解无论线性模型还是非线性模型,按照第二章介绍的线性化方法和结论,并考虑到则可写出其线性化后的函数模型为以和的估值和代入上式,则∆+=L L ~x X X ~~0+=0~111=-+∆⨯⨯⨯⨯⨯c u u c n n c W x B A 0~11=-⨯⨯⨯s x u u s W x C ∆x ~V x ˆ0ˆ111=-+⨯⨯⨯⨯⨯c u u c n n c W x B V A 0ˆ11=-⨯⨯⨯s x u u s W xC基础方程和它的解式中以上式作为函数模型而进行的平差,称为附有限制条件的条件平差,有的文献也称其为概括平差函数模型。
第五章 条件平差
v1 v V 2 n ,1 vn
W AL A0
则相应方程的矩阵表达式分别为
AV W 0
第五章:条件平差
3. 基础方程
按求函数极值的拉格朗日乘数法,设乘数 联系数向量。组成函数 将 φ 对V 求一阶导数,并令其为零,得
r , n n ,1
A V W 0 ——改正数条件方程 r ,1
W AL A0 —改正数条件方程常数项(闭合差)计算式
第五章:条件平差
例题 :右图中L1、L2、L3为观测角度, 试列出该图形的条件方程和改正数条件方 程。 解:t=2, r=n-t=3-2=1 条件方程:
ˆ L ˆ L ˆ 180 0 L 1 2 3
试列出条件方程 解:t=2p-q-4=4,r=n-t=9-4=5 条件方程为:
ˆ L ˆ L ˆ 1800 0 L 1 2 3 ˆ L ˆ L ˆ 1800 0 L 4 5 6 ˆ L ˆ L ˆ 1800 0 L
7 8 9
ˆ L ˆ L ˆ 3600 0 L 3 6 9 ˆ sin L ˆ sin L ˆ sinL 1 4 7 1 ˆ ˆ ˆ sinL sin L sin L
第五章:条件平差
4.基础方程的解
将改正数方程代入改正数条件方程,得
AQAT K W 0
令 则有
N aa AQAT AP1 AT
N aa K W 0 ——联系数法方程
秩 RN aa RAQAT RA r ,即 N aa 是个r阶的满秩方阵,由此 可解出
试按条件平差法求C、D点高程的平差值。 解:此例 n=4,t=2,r=n-t=2,可列出两个条件方程。 (1)列条件方程:
误差平差:条件平差
0 ′ ′ L v 421217′ 1 1 0 ˆ ′ ′ L = L+V = L +v2 = 78 0906′ 2 3 0 ′ ′ L v3 59 3837′
例5-2.P75
ˆ ˆ ˆ •h +h −h +HA −HB =0 1 2 3 ˆ ˆ h −h =0
条件方程的个数: 条件方程的个数:
等于多余观测数=观测总数-必要观测数。 等于多余观测数=观测总数-必要观测数。
条件方程一般是依两种思路建立的, 条件方程一般是依两种思路建立的,即:
1)沿水准路线的闭合环确定几何关系; 沿水准路线的闭合环确定几何关系; 沿附合路线来确定几何关系。 2)沿附合路线来确定几何关系。
ˆ ˆ ˆ L +L +L −1800 =0 5 6 13 ˆ ˆ ˆ L +L +L −1800 =0
7 8 14
ˆ ˆ ˆ L +L +L −1800 =0 9 10 15
ˆ ˆ ˆ ˆ ˆ L +L +L +L +L −3600 =0 11 12 13 14 15
9
ˆ ˆ ˆ ˆ ˆ sinL sinL sinL sinL sinL 1 3 5 7 9 =1 ˆ ˆ ˆ ˆ ˆ sinL sinL sinL sinL sinL 2 4 6 8 10
(5-2-1)
例如:
a2 b1 b2 c1 a1 c2
0 图形条件式为: ˆ1 ˆ ˆ 图形条件式为: a +b +c −180 =0 1 1
(5-2-2) (5-2-3)
ˆ ˆ2 ˆ a2 +b +c2 −1800 =0
5-条件平差.
0 0 s3 0
0 2 0 0 0 0 s4 0
0 1 0 0 0 2 0 0 0 1.5 0 0
§5-1 条件平差原理
二、条件平差步骤及示例
5、组成法方程,求联系数K
2 0 0 1 1 1 0 0 1 0 1 T AP A 0 1 0 1 0 0 2 0 0 0 0 1 0 0 1 1 5 0 1 0 1 1.5 平差步骤及示例
4、确定观测值的权,令C=1,则由定权公式
C 1 pi Si Si
有
p11 0 P 1 0 0
0
1 p2
0 0
1 p3
0 0
0
0 s 0 1 0 0 s2 0 0 0 1 0 0 p4
解:此例n=4,t = 2,r = n-t= 2,可列出两个条件方程。 1、列条件方程:
ˆ h ˆ h ˆ H H 0 h 1 2 3 A B ˆ ˆ h h 0 2 4
h1
C
h4
D
h3
B
h2
A
§5-1 条件平差原理
二、条件平差步骤及示例
2、计算改正数条件方程闭合差
§5-1 条件平差原理
一、基础方程及其解
改正数条件方程
r ,n n ,1
AV W 0
r ,1
改正数条件方程(纯量方程)
a1v1 a2 v 2 ... an vn wa 0 b v b v ... b v w 0 1 1 2 2 n n b ... ... ... r1 v1 r2 v 2 ... rn vn wr 0
第五章-附有参数的条件平差2009
§5.1 平差原理
在附有参数的条件平差的函数模型中,有 c = r + u < n 个条件方程,即
c× n n×1
ˆ +W = 0 AV + B x
c×u u×1 c×1
c×1
(5.1.1)
ˆ 的个数总共为 n + u 个。由于 c = r + u < n + u ,因此利用 c 个条件方程不可能解出 而未知数 V 、 x ˆ 的唯一解,但可以按照最小二乘原理求 V 和 x ˆ 的最或然值,从而求得观测值向量的最或然值 V和x ˆ = L + V 和参数的最或然值 X ˆ = X0 + x ˆ ,它们统称为平差值。 L
−1 -1 T −1 V = P −1A T N aa (BNbb B N aa − I )W
ˆ = L +V L
且观测值向量 L 的协因素阵是 Q LL = Q ,那么根据协因素阵传播定律就有
Q WW = AQA T = N aa
其余随机向量都可表示为 W 的函数,它们的协因素阵的具体表述见表 5.1.1。 表 5.1.1:附有参数条件平差的协因素阵
下面将要证明,按最小二乘准则、利用附有参数的条件平差法求得的结果,满足评定一个统计 量具有的最优性质,即满足无偏性、一致性和有效性。
ˆ 具有无偏性 ˆ 、L 一、估计量 X
ˆ 满足无偏性,下面将从另一个角度来证明。由于 ˆ 、L 在上节第四问题中,实际上已经证明了 X
用真值和真误差表示的附有参数的条件平差的函数模型为
2 0
(5.1.8)
它与平差时是否选取参数 X 无关。关于上式的严格证明,稍后进行。 (三)协因素阵 由于附有参数的条件平差中各基本向量的表达式为
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
h3
B
§5-2 条件方程
4.水准网中条件方程的列立方法
• 列条件方程的原则:足数、独立、最简 • (1)先列附合条件,再列闭合条件 • (2)附合条件按测段少的路线列立,附合条件的个
数等于已知点的个数减1 • (3)闭合条件按小环建立(保证最简),一个水准
网中有多少小环,就列多少个闭合条件。
§5-2 条件方程
空间数据误差处理
Surveying Adjustment
第五章 条件平差
第五章 条件平差
❖§5-1 条件平差原理 ❖§5-2 条件方程 ❖§5-3 精度评定 ❖§5-4 条件平差公式汇编和水准网平差实例 ❖作 业
预备知识
❖拉格朗日乘数法
求函数Z =f(x,y)在满足附加条件(x, y) 0
个数等于多余观测数r。
ALˆ A0 0
AV W 0
组成法方程式,法方程的个数等于多余观测
数r。Naa K W 0 (式中Naa AQAT)
§5-1 条件平差原理
解算法方程,求出联系数K值。 K Naa1W 将K值代入改正数方程式,求出V值,并求出平
差值
V QAT K P1AT K
P为对角阵时
• 改正数方程: V P1AT K QAT K
vi
1 pi
(aika bi
n
• 法方程: i1
kb L
ai ai pi
ka
rikr ) Qii
n i 1
aibi pi
kb
L
(aika
n airi p i1 i
bi kb
kr
L
wa
0
ri kr
)
n i 1
aibi pi
的情况下的极值问题,首先构成辅助函数,
F f (x, y) (x, y)
其中λ为某一常数。即拉格朗日乘数。再来求上 式的极值
预备知识
❖函数向量关于向量的求导规则
(1)dC 0 dX m,n
(2)Z F G m,1 m,1 m,1
(3)F A X m,1 m,n n,1
(4)F X T A X , d ( X T AX ) 2 X T A
Lˆ L V
为了检查平差计算的正确性,常用平差值 Lˆ重
新列出平差值条件方程式,看其是否满足方程。
F (Lˆ) 0
§5-1 条件平差原理
❖ 例1.设对图中的三个内角做同精度观测,得观测值 L1=
42°12′20″ ,L2= 78°09′09″ ,
L3= 59°38′40″ 。
试按条件平差求三个内角的平差值。
平差值方程
A
r,n
Lˆ
n,1
A0
r ,1
0
r ,1
a1Lˆ1 b1Lˆ1
a2 Lˆ2 b2 Lˆ2
anLˆn a0 0 bnLˆn b0 0
r1Lˆ1 r2Lˆ2 rnLˆn r0 0
§5-1 条件平差原理
条件方程
AV W 0
a1v1 a2v2 anvn wa 0
L3
L1
L2
§5-1 条件平差原理
h4
C
D
h3
h2
❖ 例2.A,B为已知水准点,其高程为HA=12.013m, HB=10.013m,可视为无误差。为了确定C及D点的高程,
共观测了四个高差,高差观测值及相应水准路线的距离 为:h1 1.004m, S1 2km ; h2 1.516m, S2 1km
1,1 1,n n,n n,1
dX
1,n n,n
§5-1 条件平差原理
➢ 基础方程及其解 ➢ 条件平差的的计算步骤
§5-1 条件平差原理
❖函数模型
A
r,n
Lˆ
n,1
A0
r ,1
0
r ,1
或
AV W 0
❖随机模型
D
n,n
2 0
Q
n,n
2 0
P1
n,n
❖平差准则
V T PV min
§5-1 条件平差原理
的平差值 l$ 4 ,1 l4
l1
A
B
l3
l2 C
§5-2 条件方程
➢ 水准网条件方程 ➢ 测角网条件方程 ➢ 测边网条件方程
§5-2 条件方程
❖一、水准网
1.水准网的分类及水准网的基准
• 分为有已知点和无已知点两类 • 要确定各点的高程,需要1个高程基准
2.水准网中必要观测数t的确定
• 有已知点: t = 网中待定点数 • 无已知点: t = 网中待定点数 - 1
❖二、测角网
1.测角网的组成
• 由三角形、大地四边形和中点多边形等三种基本图 形互相邻接或互相重叠而成。
2.测角网的观测值
• 测角网的观测值很简单,全部是角度观测值
ka
n i 1
bibi pi
kb
L
n i 1
bi ri pi
kr
wb
0
LLLL
n i 1
ai ri pi
ka
n biri p i1 i
kb L
n
i 1
ri ri pi
kr
wr
0
§5-1 条件平差原理
❖ 二、条件平差的计算步骤
根据平差问题的具体情况,列出平差值条件
方程并转化为改正数条件方程,条件方程的
V T PV min
AV W 0
r ,n n,1 r ,1
问题:根据上述方程求解V值 方法:拉格朗日乘数法
§5-1 条件平差原理
基础方程
V Βιβλιοθήκη P1AT KQAT K
AV W 0
法方程:Naa K W 0
解向量:
(1)K
N
W 1
aa
(2)V P1AT K QAT K
(3)L$ L V
§5-1 条件平差原理
h3 2.512m, S3 2km ; h4 1.520m, S4 1.5km
试求C和D点高程的平差值。
§5-1 条件平差原理
❖ 例3.如图,A、B、C三点在一直线上,测出了AB、BC及 AC的距离,得4个独立观测值:l1=200.010m, l2=300.050m,l3=300.070m,l4=500.090m。若令100m量距 的权为单位权,试按条件平差法确定A、C之间各段距离
b1v1
b2v2 bnvn
wb
0
r1v1 r2v2 rnvn wr 0
方程闭合差
wa a1L1 a2L2 anLn a0
W
AL
A0
wbb1L1
b2L2 bn Ln
b0
wr r1L1 r2L2 rn Ln r0
§5-1 条件平差原理
❖一、基础方程及其解
§5-2 条件方程
3.水准网中条件方程的分类
• 分为附合条件和闭合条件两类
• 已知点个数大于1:存在闭合和附合两类条件
• 已知点个数小于等于1:只有闭合条件
A
例:水准网如图,A、B、C三点高程 h1
未知,观测值h1—h3,列出条件方程 C
n=3,t=2,r=1
h$1 h$2 h$3 0
h2