附合导线按条件平差算例
- 格式:doc
- 大小:499.00 KB
- 文档页数:7
附合导线平差计算过程说明1)道路观测左角∑β测左=308°2.'38"+70°35'41"+156°56'39"+185°39'2"+205°21'59"+174°36'43"+197°31'46"+157°36'36"+135°14'40"+167°38'50"=1759°14'34"ƒβ测=a始边- a终边=-15"。
ƒβ容=± 40√n =±126"。
ƒβ测<ƒβ容,测角精度符合要求。
2)改正角:β=β测- ƒβ测/N。
3)坐标方位角的推算:根据起始边的坐标方位角及改正角,依据公式a下一边’= a始边+180°+转角(观测转左角)依次计算各边的坐标方位角。
4)坐标增量的计算及闭合差的调整坐标增量计算根据已经推算出的导线各边的坐标方位角和相应边的边长,按下面公式计算各边的坐标增量。
△ X AB=D AB*COS a AB,△ Y AB=D AB*SIN a AB,按附合导线要求,各边的坐标增量代数和的理论植,等于终起两点的已知坐标之差,所以,纵、横坐标增量闭合差按公式计算,Fx=∑△x测-(X终-X起)FY=∑△Y测-(Y终-Y起)导线全长闭合差f=√(ƒx2+ƒy2)=0.102m,k=f/∑D=1/38370<1/2000.满足精度要求。
5)根据后一点的坐标及改正后的坐标增量,按公式推算前一点坐标。
X前=X后+△x改Y前=Y后+△Y改最后,推算出终止边的坐标,与原有设计值相等,以作检核。
附合导线平差内业计算一、主要公式(左角):如图:有:αB1=αA B+βB±180°α12=αB1+β1±180°α23=α12+β2±180°α3C=α23+β3±180°αCD=α3C+βC±180°计算终边坐标方位角的一般公式为:α终边′=α始边+Σβ测±n×180°角度闭合差:?β测=α终边′-α终边,?β允=±10n(一级导线技术要求),式中n为导线观测角个数。
如果?β测>?β允,则说明测角误差超限,应停止计算,重新检测角度。
如果?β测>?β允,说明测角精度符合要求,此时需要进行角度闭合差的调整。
调整是应注意:当用左角计算α终边′时,改正数的符号与?β测符号相反;当用右角计算α终边′时,改正数的符号与?β测符号相同。
可将闭合差按相反符号平均分配给各观测角,而得出改正角:V改=-?β测/n式中n内角的个数,计算的改正数,取位至秒。
根据公式β=β测+V改得出改正后的观测角,继而算出新的方位角。
为了检核,最后应重新推算结束边的坐标方位角,它应与已知数值相等。
否则,应重新推算。
坐标增量闭合差:ΔX AB=D AB×COS(αAB),ΔY AB=D AB×SIN(αAB)ΔX B1=D B1×COS(αB1),ΔY B1=D B1×SIN(αB1)ΔX12=D12×COS(α12),ΔY12=D12×SIN(α12)ΔX23=D23×COS(α23),ΔY23=D23×SIN(α23)ΔX3C=D3C×COS(α3C),ΔY3C=D3C×SIN(α3C)按附合导线的要求,各边坐标增量代数和的理论值ΣΔX i、ΣΔY i,应等于终、起两点的已知坐标值之差。
附和导线平差计算详细教程,留着慢慢学习!本篇继续讲解附和导线的内业平差。
平差顾名思义就是把总误差进行平均分配,让每个点的误差都控制在允许的范围内。
平差有两种方式,一是手动平差,二是软件平差。
本文讲解手动平差,这个过程能让新手测量员们掌握平差的原理,和相关的基础知识。
本文还是以上篇的实例来讲解,开始前先来看看上篇文章中我们外业观测的记录。
第1步,制作平差计算表并填入已知数据在Excel中按适当格式制作一个《附和导线平差计算表》(表格我已经做好了,需要的可以给我发私信),然后按要求输入起始边和附和边的起、终点坐标并计算方位角和边长。
再参照观测记录表在”测点“栏中依次填入各个测点,在”观测左角值”栏中填入每个测站测得的平均角值,在“距离”栏中填入各导线边的平均边长。
填入后的效果如下图:注:已知边的方位角和边长的计算方式很多,比如用5800计算器的Pol函数,道路之星的测站、CASS查询等。
第2步,计算角度闭合差计算角度闭合差,是为了检验外业角度观测的精度是否满足相应等级导线的技术要求。
如果实测的角度闭合差<>相关计算公式:实测角度闭合差=实测附和边方位角-理论附和边方位角实测附和边方位角=起始边方位角-N*180+实测左角值之和容许角度闭合差各等级导线有相应规定(各等级导线的技术要求在上篇文章中)。
注:N为测量站数,方位角取值范围是0度(含)到360度(不含),大于360度的减去360度,小于0度的加上360度。
本案例经计算:角度闭合差=7.1秒,容许闭合差=22秒,观测精度合格。
第3步,计算左角改正数经过角度闭合差的计算,确定外业成果合格后,就要计算左角改正数。
左角改正数=角度闭合差的相反数/测站个数改正后左角值=观测左角值+左角改正数本案例角度闭合差=7.1秒,那么左角改正数=-7.1/5=-1.42秒。
分配说明:为尽量平均分配误差,我们可将改正数保留1位小数,所以案例中每个测站分配-1.4秒,这时还有0.1秒未得到分配。
附合导线严密平差算法总结图1如图的单一附合导线,有4个已知点A、B、C、D,2个未知点TP1、TP2。
设观测边数为n, 则未知点数为n-1, 观测角数为n+1。
以上图为例,n = 3。
观测边为:S1 = B->TP1,S2 = TP1->TP2, S3 = TP2->C思路:由于A、B坐标已知,则可以算出起始方位角,再根据B点坐标和每个观测角(夹角,左角)推算出TP1、TP2、C点的近似坐标值。
如果是用全站仪进行测量,则用盘左盘右重复观测求平均的方式,直接测出TP1、TP2、C点的近似坐标值以及CD的方位角。
再根据c点的已知坐标与近似坐标求坐标闭合差,由CD的已知方位角和近似方位角求角度闭合差,两个闭合差联立求得边长和角度的改正数,最后求得未知点的坐标平差值。
条件平差过程:1.建立条件方程,求得条件系数2.求法方程系数3.求权阵4.计算出联系数K5.解算出观测值改正数V6.由观测值和改正数计算平差值详细步骤如下:1、建立条件方程在单一附合导线中,只需要三个条件方程即:方位角附合条件、纵坐标附合条件和横坐标附合条件方程。
(省略了条件方程的推导过程,详细过程请查看参考资料:《测量平差.pps》)(1)方位角附合条件[Vβi]n+11+ Wβ= 0式中,Wβ= - (T0– T CD+ [βi]n+11 - (n+1)*180°)(角度闭合差)βi ——角度观测值(夹角,左角)Vβi ——各观测角的改正数。
如果是用全站仪观测,则Wβ= - (T CD– T CD)式中,T CD ——CD的方位角观测值,T CD ——CD的已知方位角(2)纵坐标X附合条件方程[Cos TI * VSi]1n - (1/ ρ”)* [(Yn+1- Yi) * Vβi]1n - W x = 0;式中,TI——各方位角观测值(近似值)VSi——边长改正数Yn+1—— C点即终止点的横坐标Y的观测值(近似值)Yi——待定点的横坐标Y的观测值Wx = - (Xn+1- XC)XC—— C点即终止点的纵坐标X的已知值ρ” = 2062.65(3)横坐标Y附合条件方程[Sin TI * VSi]1n + (1/ ρ”)* [(Xn+1- Xi) * Vβi]1n– W Y = 0;式中,TI——各方位角观测值(近似值)VSi——边长改正数Xn+1—— C点即终止点的纵坐标X的观测值(近似值)Xi——待定点的纵坐标X的观测值WX = - (Yn+1- YC)YC—— C点即终止点的横坐标Y的已知值ρ” = 2062.652、求条件方程的系数矩阵联立3个方程得改正数条件方程组:[Vβi]n+11+ Wβ= 0[Cos TI * VSi]1n - 1/ ρ”* [(Yn+1- Yi) * Vβi]1n - W x = 0;[Sin TI * VSi]1n + 1/ ρ”* [(Xn+1- Xi) * Vβi]1n– W Y = 0;其系数矩阵arrA为:(即改正数V的系数,此处以图1为例, n = 3)3、联系数法方程(简称法方程)AP-1A T K – W = 0A——系数矩阵arrAK ——乘系数P ——权阵W ——闭合差矩阵由上得法方程的系数阵N:N = AP-1A T(权的推导见参考资料:《全站仪观测导线测量平差方法的研究.pdf》)角度权:P βi = 1;(因为角度的标称精度是固定的,各观测角权值相等) 边长权:P Si = (μ0 *μ0 ) / (M D * M D )(误差比例系数固定,边长的误差与距离有关,因此不一致)式中,μ0 ——先验测角中误差,以秒为单位 M D —— 距离观测中误差若 导线边长为S i (米),e1 为仪器的边长标称固定误差(mm ),e2为仪器的边长比例误差系数(无单位),则M D = ± (e1 + e2 * S i * 0.001)mm 需转化为厘米:M D = M D ** 0.1 (cm )由于此处是要P 的逆矩阵P -1,因此要求P 的各元素的倒数(P 是对角矩阵,对角矩阵的逆矩阵就是原矩阵元素的倒数) Psi = 1/ Psi ;由上可得,P -1 矩阵如下:4、求改正数由于N 已经在前面的步骤中求出,求N 的逆矩阵。
§9.4 附合导线按条件平差算例
9.4.1附合导线的条件平差方程式
如图9-6所示,符合在已知),(A A y x A ,),(C C y x C 之间的单一符合导线有n 条AB α与CD α是已知方位角。
设观测角为
β、β、… …、β,测角中误差为 ,观测边长为s 、s 、… …、s ,
故t 1为v 1=i
i BA CD 01
1
=+∑+=a i n i v ω (9-2)
式中a ω—方位角条件的不符值,按
180)1(ˆ1
1+-∑+-=+=n i
n i CD BA a βααω (9-3) 若导线的A 点与C 点重合,则形成一闭合导线,由此坐标方位角条件就成了多边形的图形闭合条件。
2、纵、横坐标条件 设以1ˆx
∆、2ˆx ∆、…、n x ˆ∆表示图中各导线边的纵坐标增量之平差值;1ˆy
∆、2y ∆、…、n y ˆ∆表示图中各导线边的横坐标增量之平差值;由图可写出以坐标增量平差值表示的纵、横坐标条件。
⎪⎪⎭
⎪⎪⎬⎫∑+∆∑+=∆∑+=∑+∆∑+=∆∑+=∆∆yi n
i n A i n A C xi n
i n A i n A C v y y y y y v x x x x x 1111
11ˆˆ (9-4) βσ
令 ⎪⎪⎭
⎪⎪
⎬⎫--∆∑=--∆∑=)()(11
A C i n
y A C i n
x y y y x x x ωω (9-5) 则 ⎪⎪⎭
⎪
⎪
⎬⎫
=+∑=+∑∆∆0011y yi n
x xi n
v v ωω (9-6) 以微分量代替改正数,则有
)()()(211n xi n
x d x d x d v ∆++∆+∆=∑∆
{}ρ
α1
23121
1
)()()(cos v y y y y y y v v n C
si i n
xi n -'++-+--∑=∑∆
将上式代入式9-6得纵坐标条件式,且同理已可得横坐标的条件式即
⎪⎪⎭
⎪
⎪⎬⎫
=+-'∑+∑=+-'∑-
∑====0)(1sin 0)(1
cos 1111y i i C
n
i si i n i x i i C
n
i si i n
i v x x v v y y v ωραωρα (9-7)
上式就是单一符合导线的纵、横坐标条件方程x ω、y ω为条件式的不符值,按
⎪⎪⎭
⎪⎪⎬⎫
-'=-∆∑+=-'=-∆∑+=C C
C i n A y C C
C i n
A x y y y y y x x x x x 11
ωω (9-8) 式中i x 、i y 是由观测值计算的各导线点的近似坐标。
计算时一般i v 以秒为单位,si v 、x ω、y ω以cm 为单位;若x 、y 以m 为单位,则65.2062100206265==''ρ,从而使全式单位统一。
若单一导线的A 与C 点重合形成闭合导
线,则纵、横坐标条件成为多边形各边的坐标增量闭合条件,以增量平差值表示为
(9-9)
9.4.2符合导线的精度评定
⎪⎪⎭⎪⎪⎬⎫=∆∑=∆∑0ˆ0ˆ11
i n i n
y x {}ρ
ρn n C n C v y y v
y y y y y y )()()()(23423-'---'++-+--
1、单位权中误差:单一符合导线计算单位权中误差公式与边角网相同,按
(9-10)
2、平差值的权函数式:为了平定平差值函数的精度,必须要列出权函数式。
一般有下列三种函数式。
(1)边长平差值权函数式由导线边si i i v s s
+=ˆ 故其权函数式为 si Fsi v v = (9-11)
(2) 坐标方位角平差值权函数式
由图9-6得单一符合导线的任一边的坐标方位角的计算式为 180ˆ1n i n
BA i -∑+=βαα
n
i i F v v 1
∑=α (9-12)
(3) 坐标平差值的权函数式
由图9-6得j 点坐标平差值的权函数式为
⎪⎪
⎭
⎪
⎪
⎬⎫-∑+∑=-∑-∑=-=-=-=-=i i
j j i si i j i Fyi
i i
j j i si i j i Fxi v y x v v v y y v v ραρα11111
11
1sin cos (9-13)
9.4.3附合导线按条件平差算例
在图9-7所示附合导线中,B A ,为已知点,其坐标为 947.6556=A x m 735.4101=A y m 155.8748=B x m 647.6667=A y m
方位角4.130349'''= AB α,应用红外测距仪观测导线的转折角β和边长s 列入表9-5。
试按条件平差法,求各观测值及平差后3s 边的边长相对中误差。
[][][]r
v v P v v P r pvv s s s +±=±=βββσ0ˆ
表9-5 近似坐标计算
解:(1)确定观测值的权。
测角中误差0.3''
±
=
m
β
边长中误差按仪器给定公式为
26222)105()5.0()(i i c s s ppm s m m i ⨯⨯+±=+±=-(cm ) 式中i s 以cm 为单位。
由上式算得 96.01±=s m cm 82.02±=s m cm 79.03±=s m cm 92.04±=s m cm
以角度观测的权为单位权,即0.3''±==βμm
表9-6 条件方程及权函数式系数表
则边长的权为
76.9221
1==
s
m
s m m p β 38.132=s p
42.143=s p 63.104=s p 边长权倒数为 101.011=s p 074.012
=s p
070.013=s p 091.01
4
=s p (2)计算条件方程式不符值。
由表9-4-1得
0.54.1303494.080349''-='''-'''=-'=
AB AB a w αα
9.4155.8746204.8748+=-=-'=B B x x x w cm 9.2647.6667676.6667+=-=-'=B B y y y w cm
(3)计算条件方程式系数及权函数式系数列于表9-6中。
(4)组成法方程式并解算。
根据表9-6中系数组成法方程系数,然后填于表9-7中相应行内。
法方程式的解算在表9-7中进行。
表9-7 法方程式解算表
(5)计算改正数和平差值。
由法方程解算表解得的联系数和观测边加相应改正数,即得角度和边长平差值。
计算见表9-8。
(6)计算3s 边的精度。
1)单位权中误差,按 59.23
080
.26][''±=±=±
=r pvv μ 2) 计算3s 边的中误差。
cm p m s s 75.0064.095.21
3
3±=±=±=μ 3s 边边长相对中误
164000
112294275.033==s m s。