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

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

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

3、为式中，为观测值，为其权，为其相应的改正数。实际上的数值很小，可将上式按台劳级数展开，并舍去二次以上各项，得到其线性式 （2-1）式中， 根据两组平差的步骤，首先按第一组条件式进行平差，求得第一次改正后的观测值，然后改化第二组条件方程式。设改化后的第二组条件方程式为则的权倒数为 （2-2）如果平差不是按克吕格分组平差法进行的，即全部条件都是第一组，没有第二组条件，则在计算权倒数时应将上式的后两项去掉。 的中误差为 （2-3）式中，为观测值单位权中误差。 2.程序估算法 此法根据控制网略图，利用已有程序在计算机上进行计算。在计算过程中，使程序仅针对所需的推算元素计算精度并输出供使用。 通常这些程

4、序所用的平差方法都是间接平差法。设待求推算元素的中误差、权（或权系数）分别为、，后者与网形和边角观测值权的比例有关（对边角网而言），不具有随机性。至于单位权中误差，对验后网平差来说，是由观测值改正数求出的单位权标准差的估值，具有随机性。但对于设计的控制网来说，用于网的精度估算，可取有关规范规定的观测中误差或经验值。这时需要计算的主要是或，所用程序最好具有精度估算功能。否则，应加适当修改，以使其自动跳过用观测值改正数计算的程序段，而直接由用户将指定值赋给。如此计算出的即为所需结果。在这种情况下，运行程序开始时应输入由网图量取的方向和边长作为观测值，各观测值的精度也应按设计值给出。输入方式按程序规

5、定进行。图2-82.3.2 三角锁推算边长的精度估算1.单三角形中推算边长的中误差 图2-8中，设为三角形的起算边，为推算边，、为角度观测值，于是由推算的函数式为 单三角形中有下列图形条件 按角度平差时，条件方程式的系数为，对角度、的偏导数（各角以弧度为单位）如下，设角、为等精度观测，中误差为，代入（2-2）式（去掉后两项）得于是将上式的结果代入（2-3）式，并注意上式在求导数时角度是以弧度为单位的，因而相应的测角中误差也应化成以弧度为单位，即为，于是可得写成相对中误差的形式为 （2-4）过去经常使用边长对数的中误差，为此可利用微分式式中，=0.434 29为常用对数的模，将上式换成中误差的形

6、式有 （2-5）式中的是以对数第6位为单位的。于是（2-5）式又可改写为 （2-6）将上式右端的乘以根号内的和可得 （2-7）式中 （2-8）、的含义可以这样理解，因为 （以秒为单位）当=1"时左端为正弦对数每秒的增量，在对数表上即为相应每1"的正弦对数表列值之差，简称为正弦对数每秒表差。若以对数第6位为单位，则上式可写为由此可见，等于角的正弦对数每秒表差（以对数第6位为单位）。 若令 （2-9）则（2-7）式可写为 （2-10）14表2-5 （以对数第六位为单位）如果已知的不是测角中误差，而是方向中误差（有关方向和方向观测的概念见第三章），则利用关系代入（2-10）式可得

7、或 （2-11）由（2-9）和（2-8）式可知与三角形的内角有关，亦即与三角形的形状有关。通常将称为三角形的图形权倒数，也就是以方向的权为单位权，三角形推算边（一般是指精度最差的边，即最弱边）边长对数的权倒数称为三角形的图形权倒数。关于图形权倒数的这个定义不仅适用于三角形，也适用于下面讲述的大地四边形等其他图形。 为了便于计算图形权倒数，已将列成数表，以角度为引数查取（见表2-5）。2.三角形的最有利形状 以上导出了三角形的图形权倒数公式，并说明了它同三角形的形状有关。由此，我们自然会提出什么样的三角形图形权倒数最小，亦即推算出的边长精度最高的问题。图2-9 为了便于研究，选取（2-6）式进行

8、分析。令。欲使最小，亦即最小，则应使最小。表面看来这是个多元极值问题，但应注意，三个角为三角形的内角，此外由图2-8，从已知边推求任一边或应使它们精度相等，则应使。于是考虑这两个条件，可写出因而使最小变成了一元极值问题。首先求出将上式代入表达式内，得到为了求的极小值，将上式对取一阶导数，并令其为零，则经整理得方程因此 ，这个结果说明，以为底边，角度的等腰三角形，对推算边长的精度最为有利。 然而上述结果只是从推算边长精度最高这一要求得出的。如果用这种等腰三角形布设三角锁，则三角形的边长将越来越短（见图2-9），因而将无法扩展下去。这说明实际布网时不能只从精度考虑，而必须顾及各方面的条件。若按正三

9、角形布网，则不仅点位密度均匀而且正三角形的值（=4.4）与上述最有利图形（=4.0）也比较接近。因此从两个方面的要求综合考虑，可以认为正三角形是布网的理想图形。3.三角形锁中推算边长的中误差 图2-10代表一段三角形单锁，其中为起算边，为传距边。在每个三角形中与传距边相对的角为传距角，用和表示。三角形中另一个角用表示，称为间隔角，与之相对的边称为间隔边。 设三角形单锁是按角度观测和按角度平差的，也就是所有等角都是等精度独立观测值并按此参加平差。现在导出计算的边长对数中误差的公式。图2-10 由图2-10可以看出是由依次经过第1，第2，第个三角形推算而得的，由于在平差时只是将第个三角形的角度闭合

10、差平均分配在三个内角、上，因此平差后只有这三个内角是相关的，而不同三角形之间各角是互不相关的。于是每个三角形对推算边长，所产生的误差可以认为是互相独立的。因而根据协因数传播律可知，由起始边通过各三角形推算最末边的权倒数将是各三角形图形权倒数之和，即 （2-12）4.大地四边形和中点多边形推算边长的中误差图2-11 在两相邻三角形内加测一条对角线所构成的图形，称为大地四边形，如图2-11、2-12所示。这种图形在工程控制网中应用颇广，例如桥梁三角网，通常就采用一个或几个大地四边形构成。图2-13所示的图形为中点多边形。大地四边形和中 点多边形都是构成三角网的主要图形。图2-11、2-12和2-1

11、3中的是已知边，是推算边。图2-11和2-12两种图形中既含有若干图形条件（前者有3个）又含有一个极条件因此不易推出边长中误差的普遍公式。图2-12 对于大地四边形，此处只给出两种典型情况的图形权倒数公式。一种是图2-11 （a）所示的矩形大地四边形和图2-12（a）所示的菱形大地四边形（由两个等边三角形加测对角线所构成的图形）。按方向平差时它们的图形权倒数如下： 矩形大地四边形 （2-13） 菱形大地四边形 （2-14）式中（见图2-11（b）和图2-12（b）、（c）。 在图2-12（a）中，如果不加测长对角线-，而按图2-12（b）计算三角形单锁的图形权倒数，则得（见（2-12）式）。与

12、（2-14）式比较，可见加测长对角线后，前面的系数仅由1.33降低为1.25，这说明图形强度增强很少。但长对角线给观测带来困难，如在平地还须增加觇标高度。由此可见，在两个近似等边的三角形内一般不宜加测长对角线。 虽然对于任意角度的大地四边形计算图形权倒数的普遍公式不易求得，但是在实际作业中所选出的大地四边形通常总是介于矩形与菱形大地四边形之间，因此可近似地取（2-13）式和（2-14）式中系数的平均值，作为计算任意角度大地四边形图形权倒数的系数，即 （2-15） 按上式计算大地四边形权倒数时有两个不同的推算路线（见图2-11（b）和图2-12 （c），应取其中较小的。较小的那条推算路线又称最佳

13、推算路线。对于中点多边形，现给出三种图形的最弱边边长对数的权倒数如下：中点五边形 中点六边形 中点七边形 图2-13可见采用中点五边形或中点六边形较为有利。实际作业时所选定的中点多边形一般不符合等边情况，因此计算权倒数时常采用近似公式 （2-16）用上式计算中点多边形图形权倒数，同样存在两条推算路线（见图2-13（b）和图2-13（c），应取其中较小的。5.混合锁段图形权倒数的计算 实际作业时，由于受地形条件限制等原因，所选定的三角锁段常常是由几种图形混合组成的三角锁（见图2-14）。估算这种锁段的图形权倒数时，先按下列各式计算出每种图形的图形权倒数：三角形 （2-17）大地四边形 （2-18

14、）中点多边形 （2-19）式中的根据传距角由表2-5查出。但应注意，对后两种图形应取最佳推算路线求。然后取锁段传算路线上各图形权倒数之和，即为推算边长的图形权倒数 （2-20）对图2-14所示的锁段，推算边长的图形权倒数为图2-14应强调的是（2-17）（2-19）各式的单位权中误差应为方向中误差，于是最弱边边长对数的中误差为 （2-21）又由（2-5）式可知 （2-22） 以上的分析均未顾及起算边边长误差的影响。然而三角网平差后推算边长的精度不仅受水平角观测误差的影响，而且还受起算边边长误差的影响。在独立三角网中，这两种误差影响是彼此无关的、互相独立的。设起算边长相对中误差为，其边长对数中误

15、差为，当顾及起算边长误差的影响时，按误差传播定律，（2-21）式应改写为 （2-23）6.两端有起算边的三角形单锁最弱传距边边长的中误差 建立控制网时，为了提高精度，常在三角锁的两端布设起算边。当锁两端有起算边和时，最弱传距边大体上在锁的中央，即（见图2-15）。图2-15 设该锁按角度观测和按角度平差，可设想把全锁分为互相独立、大体相等的两个分段，分界边为（见图2-15）。由两端起算边和分别推算边长，可以得到两个互相独立的数值和，然后取其带权平均值。 （2-24）式中分别为、的权，则的权为。 设、分别为、的中误差，单位权中误差为，则有， （2-25）此外，由于权与误差的平方成反比，故有 （2-26）已知的权，再由（2-25）式可写出的中误差的平方为将（2-26）代入上式，化简后可得 （2-27）上式虽然是针对两端有起算边的三角形单锁导出的，但是从另一方面看，它也是在已知两个分量的中误差的情况下求其加权平均值的中误差的一般公式。 作为一个算例，现应用上式估算国家一等锁锁段中最弱边的中误差。设图2-15表示一等锁的一个锁段，两端有起算边。若两个分段的三角形形状及个数大致相同