控制测量中三角锁推算元素的精度估算

《控制测量中三角锁推算元素的精度估算》由会员分享，可在线阅读，更多相关《控制测量中三角锁推算元素的精度估算（12页珍藏版）》请在金锄头文库上搜索。

1、12.3 三角锁推算元素的精度估算在 1.1.2 小节中已经提到，控制测量工作的第一阶段就是控制网的设计阶段。论述控制网的精度是否能满足需要是技术设计报告的主要内容之一。虽然对于评定控制网的优劣、费用的高低也是一项重要的指标，但是，通常首先考虑的是精度，只有在精度指标满足要求的情况下，才考虑选择费用较低廉的布设方案。本节着重介绍估算三角锁边长精度的方法。近 20 年来，随着电子计算机的广泛应用，以近代平差理论为基础的控制网优化设计理论获得了迅速地发展。例如，仅在表达控制网质量的指标方面，无论在广度和深度上，均非过去所能比。2.3.1 精度估算的目的和方法精度估算的目的是推求控制网中边长、方位角

2、或点位坐标等的中误差，它们都是观测量平差值的函数，统称为推算元素。估算的方法有两种。1.公式估算法此法是针对某一类网形导出计算某种推算元素（例如最弱边长中误差）的普遍公式。由于这种推算过程通常相当复杂，需经过许多简化才能得出有价值的实用公式，所以得出的结果都是近似的。而对另外一些推算元素，则难以得出有实用意义的公式。公式估算法的好处是，不仅能用于定量地估算精度值，而且能定性地表达出各主要因素对最后精度的影响，从而为网的设计提供有用的参考。推导估算公式的方法以最小二乘法中条件分组平差的精度计算公式为依据，现列出公式如下。设控制网满足下列两组条件方程式02121rnbanwvrvba （）021v

3、vn（）推算元素 是观测元素平差值的函数，其一般形式为F ),(21nvllvF式中， 为观测值， 为其权， 为其相应的改正数。实际上 的数值很小，可将上iliPi iv式按台劳级数展开，并舍去二次以上各项，得到其线性式（2-nvffvF21021）式中 ),(210nllF， ， 1lflf nlf根据两组平差的步骤，首先按第一组条件式进行平差，求得第一次改正后的观测值，然后改化第二组条件方程式。设改化后的第二组条件方程式为 021AnwvAv BB则 的权倒数为F（2- 111 2222 PfAfPbfafPfF 2）如果平差不是按克吕格分组平差法进行的，即全部条件都是第一组，没有第二组条

4、件，则在计算权倒数时应将上式的后两项去掉。的中误差为F（2-FFPm13）式中， 为观测值单位权中误差。2.程序估算法此法根据控制网略图，利用已有程序在计算机上进行计算。在计算过程中，使程序仅针对所需的推算元素计算精度并输出供使用。通常这些程序所用的平差方法都是间接平差法。设待求推算元素的中误差、权（或权系数）分别为 、 ，后者与网形和边角观测值权的比例有关（对边角网iM)(iQP而言） ，不具有随机性。至于单位权中误差 ，对验后网平差来说，是由观测值改正数求出的单位权标准差的估值，具有随机性。但对于设计的控制网来说，用于网的精度估算，可取有关规范规定的观测中误差或经验值。这时需要计算的主要是

5、 或iP1，所用程序最好具有精度估算功能。否则，应加适当修改，以使其自动跳过用观iQ测值改正数计算 的程序段，而直接由用户将指定值赋给 。如此计算出的 即为所 iM需结果。在这种情况下，运行程序开始时应输入由网图量取的方向和边长作为观测值，3各观测值的精度也应按设计值给出。输入方式按程序规定进行。2.3.2 三角锁推算边长的精度估算1.单三角形中推算边长的中误差图 2-8 中，设 为三角形的起算边， 为推算边， 、0ssA、 C为角度观测值，于是由 推算 的函数式为B0sBAin单三角形中有下列图形条件018C按角度平差时，条件方程式的系数为， ，AaBa对角度 、 、 C的偏导数（各角以弧度

6、为单位）如下sABAsAssf cotininco00，BtCf设角 、 、 为等精度观测，中误差为 ，代入（2-2）式（去掉后两项）得m)cot(t22BAsfttaf3于是 )cotcot(t32)cot(t1)cot(t1 22222 BAAsBAsBAsP 将上式的结果代入（2-3）式，并注意上式在求导数时角度是以弧度为单位的，因而相应的测角中误差也应化成以弧度为单位，即为 ，于是可得m)cotcot(t32122BAAsPmss 写成相对中误差的形式为（2-)cotcot(t322Bs4）过去经常使用边长对数的中误差，为此可利用微分式图 2-84sdlg式中， =0.434 29 为

7、常用对数的模，将上式换成中误差的形式有（2-6lg10sms5）式中的 是以对数第 6 位为单位的。于是（2-5）式又可改写为smlg（2-)cotcot(t321026lg BAms 6）将上式右端的 乘以根号内的 和 可得610AcotBt（2-)(322lg BAsm7）式中 （2-BBAcot10t68）、 的含义可以这样理解，因为AB（ 以秒为单位）dAdcotsinlg当 =1时d tsil左端为正弦对数每秒的增量，在对数表上即为相应每 1的正弦对数表列值之差，简称为正弦对数每秒表差。若以对数第 6 位为单位，则上式可写为 AAd610cotsinlg由此可见， 等于角 的正弦对数

8、每秒表差（以对数第 6 位为单位） 。A若令（2-RBAA29）则（2-7）式可写为5（2-Rms32lg10）6表 2-5 （以对数第六位为单位）BAAR27如果已知的不是测角中误差 ，而是方向中误差 （有关方向和方向观测的概念见m r第三章） ，则利用关系 2代入（2-10）式可得 Rrms34lg或（2-Ps1lg11）由（2-9）和（2-8）式可知 与三角形的内角有关，亦即与三角形的形状有关。通常R将 称为三角形的图形权倒数，也就是以方向的权为单位权，三角形推算边（一般sPlg1是指精度最差的边，即最弱边）边长对数的权倒数称为三角形的图形权倒数。关于图形权倒数的这个定义不仅适用于三角形

9、，也适用于下面讲述的大地四边形等其他图形。为了便于计算图形权倒数，已将 列成数表，以角度 为引数查取（见表 2-5） 。RBA,2.三角形的最有利形状以上导出了三角形的图形权倒数公式，并说明了它同三角形的形状有关。由此，我们自然会提出什么样的三角形图形权倒数最小，亦即推算出的边长精度最高的问题。为了便于研究，选取（2-6）式进行分析。令 。欲BAQcotcott22使 最小，亦即 最小，则应使 最sPlg1smlgQ小。表面看来这是个多元极值问题，但应注意， 三个角为三角形的内角，CBA,此外由图 2-8，从已知边 推求任一边0s或 应使它们精度相等，则应使 。s CA于是考虑这两个条件，可写

10、出 AB2180)(180因而使 最小变成了一元极值问题。首先求出Q图 2-98AABcot21ct)2180cot(t 将上式代入 表达式内，得到QA22222 tan41cot3)ct1(ot)ct1(cot 为了求 的极小值，将上式对 取一阶导数，并令其为零，则 0sectan21csot232AdQ经整理得方程 3t4A因此，CA652 8274180B这个结果说明，以 为底边，角度 的等腰三角形，对推算边长的精度最0s652为有利。然而上述结果只是从推算边长精度最高这一要求得出的。如果用这种等腰三角形布设三角锁，则三角形的边长将越来越短（见图 2-9） ，因而将无法扩展下去。这说明实

11、际布网时不能只从精度考虑，而必须顾及各方面的条件。若按正三角形布网，则不仅点位密度均匀而且正三角形的 值（=4.4）与上述最有利图形（=4.0）也比较接近。R因此从两个方面的要求综合考虑，可以认为正三角形是布网的理想图形。3.三角形锁中推算边长的中误差图 2-10 代表一段三角形单锁，其中 为起算边， 为传距边。在每个0snss,321三角形中与传距边相对的角为传距角，用 和 表示。三角形中另一个角用 表示，iAiBiC称为间隔角，与之相对的边称为间隔边。设三角形单锁是按角度观测和按角度平差的，也就是所有 等,21BA角都是等精度独立观测值并按此参加平差。现在导出计算 的边长对数中误差的公式。

12、ns由图 2-10 可以看出 是由 依次经过第 1，第 2，第 个三角形推算而得的，ns0由于在平差时只是将第 个三角形的角度i闭合差平均分配在三个内角 、 、iAiB上，因此平差后只有这三个内角是相iC图 2-109关的，而不同三角形之间各角是互不相关的。于是每个三角形对推算边长，所产生的误差可以认为是互相独立的。因而根据协因数传播律可知，由起始边 通过各三角形0s推算最末边 的权倒数将是各三角形图形权倒数之和，即ns（2-12）nsRPn1lg344.大地四边形和中点多边形推算边长的中误差在两相邻三角形内加测一条对角线所构成的图形，称为大地四边形，如图 2-11、2-12 所示。这种图形在

13、工程控制网中应用颇广，例如桥梁三角网，通常就采用一个或几个大地四边形构成。图 2-13 所示的图形为中点多边形。大地四边形和中 点多边形都是构成三角网的主要图形。图 2-11、2-12 和 2-13 中的 是已知0s边， 是推算边。图 2-11 和 2-12 两种图形中既含有若干图形条件（前者有 3 个）又s含有一个极条件因此不易推出边长中误差的普遍公式。图 2-12对于大地四边形，此处只给出两种典型情况的图形权倒数公式。一种是图 2-11 （a）所示的矩形大地四边形和图 2-12（a）所示的菱形大地四边形（由两个等边三角形加测对角线所构成的图形） 。按方向平差时它们的图形权倒数如下：矩形大地

14、四边形 （2-RP75.01矩13）菱形大地四边形 （2-2.菱14）式中 （见图 2-11（b）和图 2-12（b） 、 （c） ） 。21R图 2-1110在图 2-12（a）中，如果不加测长对角线-，而按图 2-12（b）计算三角形单锁的图形权倒数，则得 （见（2-12）式） 。与（2-14）式比较，可见加测RP3.1长对角线后， 前面的系数仅由 1.33 降低为 1.25，这说明图形强度增强很少。但R长对角线给观测带来困难，如在平地还须增加觇标高度。由此可见，在两个近似等边的三角形内一般不宜加测长对角线。虽然对于任意角度的大地四边形计算图形权倒数的普遍公式不易求得，但是在实际作业中所选

15、出的大地四边形通常总是介于矩形与菱形大地四边形之间，因此可近似地取（2-13）式和（2-14）式中系数的平均值，作为计算任意角度大地四边形图形权倒数的系数，即（2-RP四115）按上式计算大地四边形权倒数时有两个不同的推算路线（见图 2-11（b）和图 2-12 （c） ） ，应取其中较小的 。 较小的那条推算路线又称最佳推算路线。R对于中点多边形，现给出三种图形的最弱边边长对数的权倒数如下：中点五边形 RP75.01中点六边形 2.中点七边形 05.1图 2-13可见采用中点五边形或中点六边形较为有利。实际作业时所选定的中点多边形一般不符合等边情况，因此计算权倒数时常采用近似公式（2-RP中116）用上式计算中点多边形图形权倒数，同样存在两条推算路线（见图