附合导线按条件平差算例
按条件平差解算单导线和导线网

n
n
1
wy = ∑ ∆yi,i+1
i =1
n
i =2
ρ
5.6.1 单导线的平差解算和精度评定
2、 无定向导线和单定向导线的平差解算 、
β
T1
2
β
s
2
n
2
s
1
n
1
Tn
sn
n+1
引入方位角条件未知数T1,Tn,采用附有未 知数的条件平差模型计算。其条件方程如下: 方位角条件:
∑ vβ + δ T − δ T + wT = 0
1
J
n
J+1
每个导线段可列出3个条件:
坐标方位角的条件式:∑ v β + δ TJ − δ TJ +1 + wT
i =1
i
nj
J , J +1
=0
其中: wT
J , J +1
= TJ + ∑ v β − TJ0+1 − n J × 180
0 i =1
i
nj
5.6.2 导线网的平差解算和精度评定
纵坐标条件式: ∑ cos Ti v s − ∑
i
ρ
1
ρ
n
wy = yA + ∑ ∆yi0,i+1 − yc
i =1
n
5.6.1 单导线的平差解算和精度评定
如果只在导线的一端引入方位角未知数 T1,则只 能列出如下两个坐标条件:
0 0 ∑ cosTi vs − ∑ ( yn+1 − yi )vβ − ( yn+1 − y1 )δ T + wx = 0 0 0 0 i =1
附合导线平差计算过程说明

附合导线平差计算过程说明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改最后,推算出终止边的坐标,与原有设计值相等,以作检核。
附和导线、闭合导线的简易平差与支导线的计算

否 无法计算,
所
计算的类 型
跳出提示选 择的提醒框
示
是
点击计算
是否
否
保存
是
保存结果TXT
退出
•
程 序 设 计 流 程
图
基础理论
• 计算方位角闭合差fβ,判断Sβ<fβ容
• β=α’CD -αCD
• fβ容=60或40
n:角的个数
• 分配角闭合差(平均分配),计算改正后的角
度
• V=—β/n βi改=βi+V
政
关键技术和难点
• 度分秒与弧度间的互相转换
•
public static class hddfm
•
{
•
const double pi = 3.1415926;
•
//角度转化为弧度
•
public static double angle_r(double alfa)
•
{
度
•
double alfa1, alfa2;
弧 度
•
return (alfa1 / 180d * pi);
的
•
}
•
//弧度转化为角度
相
•
public static double redian_a(double alfa)
互
•
{
•
double alfa1, alfa2;
•
alfa = alfa * 180d / pi;
转 换
•
alfa1 = Math.Floor(alfa) + Math.Floor((alfa - Math.Floor(alfa)) * 60d) / 100d;
1-11单导线条件平差计算--附合导线坐标条件方程

平差值条件方程: ALˆ A0 0
改正数条件方程: AV W 0
STEP1
表达成角度和边
长平差值的形式
STEP2
表达成角度和边
长改正数的形式
附合导线坐标条件方程
1、纵坐标附合条件方程式
xˆN 1 xC 0
STEP1
xˆN 1 xB [xˆi ]1N
xˆi Sˆi cosTˆi
3
N
SN-1
SN
D
TCD
βN+1
C (N+1)
纵坐标符合条件方程式就是:xˆN 1 xC 0 横坐标符合条件方程式就是:yˆ N 1 yC 0
B点坐标的已知值 (xB , yB ) 或 (x1, y1)
C点坐标已知值
(xC , yC ) 或 (xN 1, yN 1)
C点坐标观测值
A
D
β2
β4
βN
2
4
β1 S1 S2 β3 S3
N
SN-1
SN βN+1
3
B(1)
附合导线示例图
C(N+1)
附合导线坐标条件方程
补充:具体计算时的单位选取
如果x、y 以米为单位,w、vS 以厘米为单位、vβ 以秒为单位。
[cos Ti
vSi
]1N
1
[( yN 1
yi )vi
]1N
STEP2
表达成角度和边 长改正数的形式
xˆN 1 xB [xˆi ]1N
xˆi Sˆi cosTˆi
非线性从这个式子着手!!
Tˆi TAB ˆ j
i
180
1-15条件平差精度评定--条件平差单位权中误差的计算

B
……
……
方位角T ……
条件平差精度评定
1.角度和边长的平差值(观测值平差值)的精度 2.未知点的坐标平差值(观测值平差值函数)的精度
已有的量
观测值的权阵
PLL
观测值的协因数阵
QLL
P 1 LL
待求的量
1.观测值平差值的协方差阵
DLˆLˆ
建立 联系
2.观测值平差值函数的协方差阵
测量平差的第二个任务
评定测量成果的精度
? 具体来说评定什么量的精度
导线条件平差实例
图示为一四等附合导线,《规范》测角 中误差 2.5 ,测边所用测距仪的标称精 度公式 S 5mm 5 ppmDkm 。 已知数据和观测值见表1、表2。 试按条件平差法求:1.角度和边长的平差值;
0 0 0 0 0 1 0 0 0 0.833
0
0
0
0
0
0
1
0
0
0.612
0 0 0 0 0 0 0 1 0 0.049
0 0 0 0 0 0 0 0 1 3.246
3.2460
V T PV (W )T K 3.9 -1.6 1.71.0518 20.2761
3.4906
ˆ0
V T PV
r
20.276132 2.6 3
后验单位权中误差估值
单位权中误差
? 0 2.5" 先验单位权中误差
理论
感谢聆听,批评指导
公式
思考
平差
算例
间接法 矩阵式
V T PV V T P(P1AT K ) V T AT K ( AV )T K (W )T K
附合导线平差计算(函数版)

22468
数据精度
由 两 点 坐 标 , 求 得 方 位 角 , 返 还 十 进 制 角 度
∑△y = 1464.520 fy=∑△y-(Yc-Yb)= 0.016 f=√(fx*fx+fy*fy) = 0.082 K=f/∑D≈1/22468<1/15000 符合精度要求
3.43 #### 22′ 57.8″ ####### -121.569 148.092 3.43 43° 27′ 32.2″ ####### 115.012 3.43 61° 25′ 16.7″ ####### 116.022 3.43 #### 31′ 57.1″ ####### -99.987 3.43 64° 25′ 23.5″ ####### 3.43 #### 20′ 41.0″ ####### -65.986 3.43 29° 19′ 45.4″ ####### 210.005 3.43 #### 01′ 57.8″ ####### -68.994 3.43 44° 57′ 49.3″ ####### 148.645 3.43 83° 28′ 20.7″ 148.456 200.002 117.990 225.004 78.012 162.994 140.009 212.987 108.986
验 算
1801° 30′ 21″ fβ = Σ fβ 测-Σ fβ 理= -34″ fβ 允= ±10√n = ±32″ fβ >fβ 允 , 不符合精度要求
∑△x = 311.160 fx=∑△x-(Xc-Xb)=0.080 f=√(Βιβλιοθήκη x*fx+fy*fy) =
K=f/∑D≈1/22468<1/15000
线 计 算 表
附合导线平差内业计算

附合导线平差内业计算一、主要公式(左角):如图:有:α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、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
§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。