文档详情

三角网条件平差计算

jin****ng
实名认证
店铺
DOCX
130.48KB
约16页
文档ID:191580043
三角网条件平差计算_第1页
1/16

§3-4 三角网条件平差计算2学时三角网测量的目的,是通过观测三角形的各角度或边长,计算三角网中各未知点的坐标、 边的长度及方位角等三角网按条件平差计算时,首要的问题是列出条件方程因此了解三 角网的构成,总结其条件方程的种类及各种条件方程的组成规律是十分重要的三角网的种类比较多,网的布设形式也比较复杂根据观测容的不同,有测角网、测边 网、边角同测网等;根据网中起始数据的多少,有自由三角网和非自由三角网自由三角网 是指仅具有必要起算数据的三角网,网中没有多余的已知数据如果测角三角网中,只有两 个已知点(或者已知一个已知点的坐标、一条已知边的长度和一个已知的方位角),根据数 学理论,以这两个已知点为起算数据,再结合必要的角度测量值,就能够解算出网中所有未 知点的坐标如果三角网中除了必要的起算数据外还有其它的已知数据,或者说已知数据有 冗余,就会增加对网形的约束,从而增强其可靠性,这种三角网称之为非自由三角网无论 多么复杂的三角网,都是由单三角形、四边形和中点多边形组合而成的在本节,我们先讨论三角网条件平差中条件方程个数的确定问题,然后主要讨论测角三 角网的条件方程的形式问题一、网中条件方程的个数三角网平差的目的,是要确定三角点在平面坐标系中的坐标最或然值。

如图3-9 所示, 根据前面学到的测量基础知识,我们知道,必须事先知道三角网中的四个数据,如两个三角 点的4个坐标值,或者一个三角点的2 个坐标值、一条边的长度和一个方位角,这4个已知 数据我们称之为三角网的必要起算数据有了必要起算数据,就可以确定三角网在平面坐标 系中的位置、网的大小及其方位,就可以计算三角网中未知点的坐标要对三角网进行平差计算,还必须先知道网中的总观测数n、判定必要观测数t, 从而 确定了多余观测数:r = n - t由条件平差原理知,多余观测数与条件方程数是相等的,有了多余观测数,也就确定出了条件方程的个数因此,问题的关键是判定必要观测数t1. 网中有2个或2个以上已知点的情况三角网中有2个或2 个以上已知三角点,就一定具备了4个必要起算数据无论是测 角网、测边网还是边角同测网,如果有2个已知点相邻,要确定一个未知点的坐标,需要观 测两个观测值(2个角,或者1条边和1个角,或者2条边)也就是说,确定1个未知点 要有2个必要观测值;那么如果网中有p个未知点,必要观测数应等于未知点个数的两倍t = 2 • p (3-4-1)(1) 测角网图3-9 所示,三角网中有2个已知点,待定点个数为p = 6。

如果三角网中观测量全部 是角度时总观测值个数:n = 23必要观测数:t = 2 - p =12则多余观测数,即条件平差条件方程个数:r = n - t = 11(2) 测边网在图 3-9 中,如果三角网中观测量全部是边的长度时:总观测值个数:n = 14必要观测数:t = 2 - p =12则多余观测数,即条件平差条件方程个数:r = n - t = 2G15E16H101417ISAP11201921KB(3) 边角同测网在图3-9中,如果三角网中的所有的角度值和所有的边长值都进行观测时:总观测值个数:n = 37必要观测数:t = 2 • p =12则多余观测数,即条件平差条件方程个数:r = n - t =252. 网中已知点少于2 个的情况有些情况下,三角网中已知点可能少于2个,只有1 个已知点、1个已知边和1个已知 方位角,或者没有已知点和已知方位角只有1个已知边但是,不管怎样说, 1条已知边是 必须已知的,或者需要进行观测的如果没有已知点,可以假定网中的1 个未知点;如果没 有已知方位角,可以取网中的1个方向的方位角为某一假定值这样也就间接地等价于网中 有2 个相邻点的坐标是已知的。

1) 测角网三角网中共有P个三角点、1个已知方位角(也可以没有)、1个已知点(也可以没有 已知点)和1个已知边长S(或者也是观测得到的),并观测了所有的角度如果已知点和 已知方位角都没有,就要进行必要的假设则在进行条件平差时,必要观测数为:t = 2 - ( p - 2) (3-4-2)如图3-10所示,三角网中观测了所有角度值(如果没有已知边时,也观测1条边长作 为起算数据)网中三角点个数:P = 6角度观测值个数:n = 12必要观测数:t = 2 - ( p 一 2) = 8(2) 测边网或边角同测网若三角网中,共有p个三角点和1个已知点(或者也是假定的),并对所有的边长,或 者角度和边长进行了观测,观测值总个数为n在进行条件平差时,由于要加上必须的起算 边长,则必要观测(边或者边和角)的个数为t = 2 • ( p — 2)+1 (3-4-3)如图3-10 所示,网中三角点个数:P = 6如果是测边网,则总观测值个数:n = 9必要观测数:t = 2 • ( p 一 2) +1=9多余观测数,即条件平差条件方程个数:r = n - t = 0如果是边角同测网,则总观测值个数:n = 21必要观测数:t = 2 • ( p - 2) +1=9多余观测数,即条件平差条件方程个数:r = n - t = 12以上我们仅对几种三角网,讨论了条件平差时必要观测数及多余观测数和条件平差方程 数的确定方法,还有很多情况没有涉及到。

在实际平差计算中,应针对不同情况进行具体分 析二、条件方程的形式三角网中的条件方程主要有以下几种形式:1. 图形条件方程图形条件,又叫三角形角和条件,或三角形闭合差条件在三角网中,一般对三角形的 每个角都进行了观测根据平面几何知识,三角形的三个角的平差值的和应为180°,如图 3-12中的三角形ABP,其角平差值的和应满足下述关系:L L + L3 -180 = 0 D此即为三角形角和条件方程由于三角形是组成三角网的最基本的几何图形,因此, 通常称三角形角和条件为图形条件因此图形条件也是三角网的最基本、最常见的条件方程 形式与(3-4-4)式相对应的改正数条件方程为v + v + v 一 w = 0 (3-4-5)123w = —(L + L + L -180) (3-4-6)1232. 水平条件方程 水平条件,又称圆周条件,这种条件方程一般见于中点多边形中如图3-12 所示,在中点P上设观测站时,周围的五个角度都要观测这五个观测值的平差值之和应等于360, 即L + L + L + L + L — 360 0 (3-4-7)3 6 9 12 15相应的改正数条件方程为v + v + v + v + v — w = 0 (3-4-8)3 6 9 12 15w = —(L + L + L + L + L — 360 °) (3-4-9)3 6 9 12 153. 极条件方程极条件是一种边长条件,一般见于中点多边形和四边形中。

先看中点多边形的情况如 图 3-12 所示,中心 P 点为顶点,有五条边,从其中任一条边开始依次推算其它各边的长度 最后又回到起始边,则起始边长度的平差值应与推算值的长度相等在图3-12所示的三角网中,我们应用正弦定理,以BP边为起算边,依次推算AP、EP、DP、CP,最后回到起算边BP、,得到下式S = S 哑• bp bp sin L2sin Lsin L5sin L• 7八sin L8sin L• 卅sin L11八sin L• *sin L14整理得/\ /\ /\sin L sin L sin Lsin Lsin L(3-4-10)14 71013—1 二 0sin L sin L sin L sin L sin L2 5 8 11 14(3-4-10)式即为平差值的极条件方程为得到其改正数条件方程形式,可用泰勒级数对上式左边展开并取至一次项:sin L sin L sin L sin L sin L1 4 7 10 13 — 1 =/\ /\ /\ /\ -1.sin L sin L sin L sin L sin L2 5 8 11 14sin L sin L sin L sin L sin L1 4 7 10 13 — 1sin L sin L sin L sin L sin L2 5 8 11 14sin L sin L sin L sin L sin L v1 4 7 10 13 COt L 4 —sin L sin L sin L sin L sin L 4 7 10 13 COt L 2-sin L sin L sin L sin L sin L 2 p" p"2 5 8 11 14sin L sin L sin L sin L sin L vsin L sin L sin L sin L sin L v sin L sin L sin L sin L sin L v+ 1 4 7 10 13 COt L 4 — 1 4 7 10 13 COt L 5-sin L sin L sin L sin L sin L 4 p" sin L sin L sin L sin L sin L 5 p"2 5 8 11 14 2 5 8 11 14sin Lsin Lsin Lsin Lsin L+710-13sin Lsin Lsin Lsin Lsin L2581114sin Lsin Lsin Lsin Lsin L+71013sin Lsin Lsin Lsin Lsin L2581114sin Lsin Lsin Lsin Lsin L+1471013COt LCOt Lsin L sin L sin L sin L sin L2 5 8 11 14v1010 p""vCOt L 7—7 p〃sin L sin L sin L sin L sin L v1 4 7 10 13 COt L 8-sin L sin L sin L sin L sin L8 11 148 p ""sin L sin L sin L7 10 13 COt Lsin L sin L sin L sin L sin L2 5 8 11 14sin L sin L1 4v1111 p"v1313 p"sin L sin L sin L sin L sin L1 4 7 10 13 COt Lsin L sin L sin L sin L sin L 142 5 8 11 14v14 = 0 P"化简,即得极条件的改正数条件方程:ctgL v — ctgL v + ctgL v — ctgL v + ctgL v — ctgL v1 1 2 2 4 4 5 5 7 7 8 8 (3-4-11)+ ctgL v — ctgL v + ctgL v — ctgL v — w = 010 10 11 11 13 13 14 14,/ sin L sin L sin L sin L sin L '(3-4-12)W = —p 1 2 5 8 H 14( sin L sin L sin L sin L sin L ,1 4 7 10 13在四边形中的极条件方程与中点多边形稍有不同。

如图 3-11所示,可以取 D 点为极点,以 BD 为起始边,依次推算 AD、CD 再回到 BD 边仿照中点多边形的极条件方程,由正弦 定理,得四边形的极条件平差值方程sin(L + L )sin L sin L7 8 2 4 — 1 = 0sin L sin(L + L )sin L1 3 4 7整理得八 八 八 八sin L sin L sin(L + L )T 护 A 以—1 = 0 (3-4-13)sin L sin(L + L )sin L1 3 4 7相应的改正数条件方程ctgLv + ctgL v — ctg (L + L )v + (ctgL — ctg (L + L ))v1 1 2 2 3 4 3 4 3 4 4 (3-4-14)+ (ctg (L + L ) — ctgL )v + ctg (L + L )v — w = 07 8 7 7 7 8 8(3-4-15),/ sin L sin(L + L )sin L 'W = 一 p 1 1 3 4 7\ sin L sin L sin( L + L ),2 4 7 8前面讨论的三种条件方程在三角网中比较常见如果三角网中的起始数据有了变化,起 算数据不相邻,或者已知数据有冗余,还会增加一些限制条件,产生其它类型的条件方程 如方位角条件方程、边长条件方程、坐标条件方程等。

这些类型的条件方程常见于非自由三 角网中如图3-13所示,为一个非自由三角网,有 4个已知点、2 个未知点和12 个角度观测值必要观测个数t = 2 x 2 = 4,多余观测数r = n - t = 12- 4 = 8, 即共有8个条件方程,其中 图形条件方程有4 个,没有极条件,也没有水平角条件,那么另4 个是什么类型的呢?由于 三角网中有4 个已知点,每个已知点有2 个坐标值,共计8个已知数据,超过了4 个必要起 算数据,从而产生4 个冗余的已知数据这4 个多余的已知数据必然会导致4 个矛盾,进而 产生4个条件方程方位角条件,严格地说是方位角附合条件,是指从一个已知方位角出发,推算至另一个 已知方位角后,所得推算值应与原已知值相等如从4个已知点可以反算岀AB和EF两边的边长值和方位角值,这些值也可看作是已知 值,作为起算数据用设AB边的方位角T,EF边的已知方位角为T 如果从 AB向 E F推算,推算路线如AB EF图中所示,设EF方位角的推算值的最或然值为亍,近似值为T则方位角附合条件方程EF EF为八 —T — T = 0 (3-4-16)EF EF其中八 八 八 八 八T = T — L + L + L — L 土 3 -180。

EF AB 3 6 9 12代入(3-4-16)后,整理得八 八 八 八 —L + L + L — L + T — T 土 3 -180 03 6 9 12 AB EF(3-4-17)其相应的改正数条件方程—v + v + v — v — w = 0 3 6 9 12 T其中w =—(—L + L + L — L + T — T 土 3 -180) (3-4-18)T 3 6 9 12 AB EF5. 边长条件方程边长条件,严格地说是边长附合条件,是指从一个已知边长出发,推算至另一个已知边长后,所得推算值应与原已知值相等图3-13三角网中,设AB边的已知长度为S , EF边的已知长度为S 如果沿图中AB EF所示的推算路线,从AB向E F推算,得EF边长推算值的最或然值为S ,近似值为Sef则 EF EF边长附合条件方程为/\ S - S = 0EF EF(3-4-19)其中八SEFsin L sin L sin L sin L=S 1 4 7 10ab sin L sin L sin L sin L2 5 8 11将上式代入(3-4-19)式,并将边长条件整理为 八 八 八 八S sin L sin L sin L sin LAB 1 4 7 10 — 1 = 0 __ 八 八 八 jl vyS sin L sin L sin L sin LEF 2 5 8 11(3-4-20)仿照极条件式,将上式左边用泰勒级数展开,取至一次项,整理后得其改正数条件方程:ctgL v — ctgL v + ctgL v — ctgL v + ctgL v — ctgL v1 1 2 2 4 4 5 5 7 7 8 8+ ctgL v — ctgL v — w 二 010 10 11 11 S(3-4-21)二—p,( S sin L sin L sin L sin L 'rr -11 — 2 5 8 11. S sin L sin L sin L sin L 丿AB 1 4 7 10(3-4-22)6. 坐标条件方程 坐标条件方程,是指从一个已知点出发,推算至另一个已知点后,所得推算值应与 该点的已知坐标值相等。

图3-13 三角网中,设B点的已知坐标为(x,y ),E点的已知坐标为(x,y )B B E E如果沿图中所示的路线,从BfCfE进行推算,得E点坐标推算值的最或然值为(x,y ),EE近似值为(xE yE)则坐标条件方程为(3-4-23)(3-4-24)八 A 八 A 八 I Ix = x + Ajt + Ax = xE B BC CE B+ SBC C0S TBC + S CE C0S TCE(3-4-25)X E - XE = 0yE - yE = 0而其中八SBCsin L=S 1AB sin L2八 _S = SCE ABsin L sin L sin L sin L1 4 7 12sin L sin L sin L sin L2 5 8 11(3-4-26)(3-4-27)、 八、=T — L + 180BC AB 3(3-4-28)、 八 八 八 八' =T — L + L + L + L 土 2 -180CE AB 3 6 9 10(3-4-29)将上述(3-4-26) ~ (3-4-29)式代入(3-4-25)式,然后用泰勒级数展开,取至一次项,整理后得:( x — x )(ctgL v — ctgL v ) + (x — x )(ctgL v — ctgL v )E B 1 1 2 2 E C 4 4 5 5+ ( x — x )(ctgL v — ctgL v ) + ( x — x )(ctgL v — ctgL v )E C 7 7 8 8 E C 12 12 11 11—(y —y )(—v )—(y —y )(+v )—(y —y )(+v )—(y — y )(+v )E B 3 E C 6 E C 9 E C 10—w = 0x(3-4-30)为不使闭合差项w过大,影响平差结果的精度,在计算坐标条件方程时,可以考虑X、Xy以公里(km)为单位,而w中的坐标差项以米(m )为单位。

即X—P(xe — xE)=—206.265(x — x )1000 E E(3-4-31)其中同理可写出横坐标改正数条件方程(y —y )(ctgL v —ctgL v )+(y —y )(ctgL v — ctgL v )E B 1 1 2 2 E C 4 4 5 5+(y —y )(ctgL v — ctgL v )+(y —y )(ctgL v —ctgL vE C 7 7 8 8 E C 12 12 11 11— (x — x )(—v ) — ( x — x )(+v ) — (x — x )(+v ) — (x — x )(+vE B 3 E C 6 E C 9 E C—w = 0y(3-4-32)10)10w = -206.265(y — y ) (3-4-33)y E E坐标附合条件方程,尤其是改正数条件方程,形式上虽然比较复杂,但也非常具有规律 性这一点,请同学们结合图3-13认真地分析,看能否总结出其概括形式以上八种条件方程及其改正数条件方程的类型和形式,基本上涵盖了测角型三角网条件 方程的基本形式需要说明的是,三角网布设形式极其多样,条件方程的形式也较为繁杂, 但关键是要掌握其基本形式,并能融会贯通灵活运用。

三、例题如图3-14是一个三角网,A、B、E、F是已知点,C、D是待定点,等精度观测了所有角值,已知数据和观测数据如表3-4所示试列出用条件平差法时的改正数条件方程表3-4已知坐标(m)B (183120.420,29502560.540)E (181740.210,29505455.940)已知方位角T = 32°20‘ 14.9〃ABTEF355°53‘ 42.6〃已知边长(m)S = 2501.118AB%= 2582.529角度观测值B = 46°211156.1〃B =462°21/B = 74°59'241.4〃42.4〃B = 58°38‘317.2〃B =567O39Z43.6〃B =549°58‘3&9〃B = 58°03‘ 46.6〃1B = 910431 54.0〃1B = 53°15' 16.1〃B = 470211 49.9〃22B = 68°40‘ 54.3〃B = 400541 08.1〃33解:这是一个非自由测角三角网观测值总数n = 12必要观测数t = 4多余观测数r = n - t = 8即,有8个条件方程,而网中共有:4个图形条件、4个坐标附合条件、1 个方位角条件和1 个边长附条件可选。

由于坐标条件较为复杂,为计算方便,选4个图形条件、1 个方 位角条件、1 个边长条件和2个坐标条件运算线路如图中所示因为角度的观测精度相同,取Q = P-1=E 首先,根据观测值,利用余切公式计算有关近似坐标:C(181440.319 , 29503390.921)D(183084.184 , 29504111.735)E(181740.109 , 29505456.041)图形条件方程v + v + v 一 w 二 01 2 3 1v + v + v 一 w 二 04 5 6 2v + v + v 一 w 二 07 8 9 3v + v + v 一 w 二 010 11 12 4方位角条件方程一 v + v 一 v + v 一 w 二 03 6 9 12 5边长条件方程S sin L sin L sin L sin LAB 1 4 7 舟一1 = 0 八 八 八 /\ JL vy2 5 8 11其改正数形式cot L v11- cot L v + cot L v - cot L v + cot L v - cot L v2 2 4 4 55 7 7 88+ cot L v - cot L v - w = 010 10 11 11 6纵横坐标条件方程(yE)(cot L v - cot L v ) + ( x - x )(cot L v - cot L v )B 1 1 2 2 E C 4 4 5 5+ (x - x )(cot L v - cot L v ) - ( y - y )(-v )E D 7 7 88 E B 3-(y - y )(+v ) - (y - y )(-v ) - w = 0E C 6 E D 9 7- y )(cot L v - cot L v ) + ( y - y )(cot L v - cot L v )1 1 2 2 E C 4 4 5 5+(y -y )(cot L v -cotLv)-(x -x )(-v ) E D 7 7 8 8 E B 3-(x - x )(+v ) - (x - x )(-v ) - w = 0E C 6其中闭合差项为:=-(L1+L2+L3-180。

) = 5.3=-(L4+ L5+L6-180) = - 4.9=-(L7+L9-180) = 2.5=-(L10+L11+ L12 - 180°)= 8.0=-(-L3+ L - L + T - T9 12 AB EF土 3 ・ 180 ) = 一 7.8=—p— 、S sin L sin L sin L sin L 〜——―n = - 23.2701 4 7 10 丿/ —〃1-亠「 S sin L sin L sin L sin L ABw =-206.2657x - x ) = 20.8337 E Ew =-206・265"(y - y ) = 一 20.8338 E E以上改正数条件方程写成矩阵形式AV — W = 0V =[ vvvvvvvvvv v v ]T12 3 4 5 6 78 9 10 1112W=[wwwwwwww]T=[ 5.3 -4.9 2.5 8.0 -7.8 -23.270 20.833 -20.833 ]T12345678S sin L sin L sin L sin L EF计算改正数条件方程系数,并以矩阵形式表示为:111111111111-11-110.953 - 0.2680.524 - 0.4110.623- 0.747- 0.030 - 0.921- 1.315 0.3702.895 0.157 - 0.123- 2.065- 0.8371.0041.3442.759 - 0.7761.380 1.082 - 0.8490.3000.837 - 1.0041.3443.000-1.0000.6851.9503.3633.0001.0000.113- 2.0310.5333.000-1.000- 0.1241.5111.1773.0001.000- 0.951-1.0001.000-1.0001.0004.000- 6.304- 2.4240.6850.113- 0.124- 0.9513.219- 2.4915.0251.950- 2.0311.511- 6.304- 2.49118.066- 0.1683.3630.5331.177- 2.4245.025- 0.16815.615N=AQAT=-3.7869 -10.7199 - 2.1844 0.8237]tK = N-iW = [3.4485 -1.5924 —.0.0843 0.5308V=P -iAtK=[-1.6 4.9 2.0 -6.7 2.4 -0.6 -4.3 4.9 1.9 0.9 10.4 -3.3]tVtPV—P- =±5.6〃r2 5 8 11 14。

下载提示
相关文档
正为您匹配相似的精品文档