九章节控制网平差.ppt

《九章节控制网平差.ppt》由会员分享，可在线阅读，更多相关《九章节控制网平差.ppt（74页珍藏版）》请在三一办公上搜索。

1、第九章 控制网平差,本章提要9.1 条件平差数学模型和公式9.2 水准网按条件平差算例 9.3 附合导线按条件平差算例 9.4 参数平差数学模型和公式9.5 高程网参数平差及算例 9.6 三角网网参数平差及算例 习题,本章提要,本章讲述条件平差与参数平差的原理及基本数学模型，两种方法计算结果是完全相同的。还介绍了高程网条件平差，三角网条件平差，附合导线条件平差。高程网参数平差，三角网参数平差，并给出了算例。,1条件平差与参数平差原理 2条件差的步骤及相应数学模型；3能分别采用条件平差与参数平差解决高程控制网，平面控制网平差。,知识点及学习要求,难点在本章学习过程中，伴随有大量的公式推导与应用。

2、特别是控制网条件方程与误差方程列立，法方程解算为本章的突破点。,9.1 条件平差数学模型和公式,设某一平差问题中有个 误差独立的观测值，个函数独立的未知数（必要观测数），多余观测数为,记：观测值,相应权阵,平差值改正数,平差值,1、条件平差的数学模型和公式 1）平差值方程（）,（1）,式中、（=1、2、）为条件方程的系数；、为条件方程的常项数,2）改正数条件方程,以（=1、2、）代入(1)得纯量形式为：,（2）,式中、为条件方程的闭合差，或称为条件方程的不符值，即,令,矩阵形式为：,（4）,（3）,3）改正数方程 上改正数条件方程式中 的解不是唯一的解，根据最小二乘原理，在 的无穷多组解中，取

3、=最小的一组解是唯一的，的这一组解，可用拉格朗日乘数法解出。为此，设，称为联系数向量，它的唯数与条件方程个数相等，按拉格朗日乘数法解条件极值问题时，要组成新的函数：,将对 求一阶导数，并令其为零得：,（5）,上式称为改正数方程，其纯量形式为：,（i=1,2,n）（6）,4）法方程,将 代入 得,矩阵形式为：,（7）,上式称为联系数法方程，简称法方程。式中法方程系数距阵，为,因 故，是 阶的对称方阵。,法方程的纯量形式为,（8）,从法方程解出联系数K后，将 值代入改正数方程，求出改正数 值，再求平差值，这样就完成了按条件平差求平差值的工作。,2、条件平差法求平差值的步骤,根据平差问题的具体情况，

4、列出平差值条件方程式（1），并转化为改正出数的条件方程（2），条件方程的个数等于为多余观测的个数r；根据条件方程的系数，闭合差及观测值的权组成法方程式（8）；法方程的个数等于多余观测的；解法方程，求出联系数K,并代入法方程检验；将K代入改正数方程（6），求改正数 值；将V代入平差值方程 求平差值；将平差值 代入平差值方程，检验是否满足条件；精度评定。,2、精度评定,1）单位权中误差,从中误差计算公式可知，为了计算，关键是计算。下面将讨论 的计算方法。,由 直接计算,由联系数 及常数项 计算,直接在高斯杜力特表格中解算,2）平差值函数的权倒数,设有平差值函数为 它的权函数式为：,令则,这就是高斯

5、约化表中 的计算公式，其规律与 计算规律完全相同。,9.2 水准网按条件平差算例,在如图1所示水准网中,两点高程及各观测高差和路线长度列于（表1）中。,试求：（）、及 点高程之最或然值；（）、点间平差后高差的中误差。,解:列条件方程式,不符值以“mm”为单位。,已知，故，其条件方程式为,列函数式。,故,组成法方程式。,1）令每公里观测高差的权为1，按1/,将条件方程系数及其与权倒数之乘积填于（表2）中。,2)由下表2数字计算法方程系数,并组成法方程式:,+,=0,表2 条件方程系数表,续表,4）法方程式的解算。1）解算法方程式在（表3）中进行。2）计算之检核。,由表3中解得，两者完全一致，证明

6、表中解算无误。5）计算观测值改正数及平差值(见表4)6）计算 点高程最或然值。,m,m,m,表3 高斯-杜力特表格,表4 改正数与平差值计算表,7）精度评定。单位权（每公里观测高差）中误差,mm,点间平差后高差中误差,mm,9.3 附合导线按条件平差算例,1.附合导线的条件平差方程式,如图1所示，符合在已知，之间的单一符合导线有 条 与 是已知方位角。,设观测角为、，测角中误差为，观测边长为、，测边中误差为（1、2、）。,此导线共有 个观测值，有 个未知数，故 则。因此，应列出三个条件方程，其中一个是坐标方位角条件，另两个是纵、横坐标条件。,1）坐标方位角条件 设观测角 的改正数为（1、2、1

7、），观测边 的改正数为（1、2、）。由图1知,式中 方位角条件的不符值，按,若导线的A 点B 与点重合，则形成一闭合导线，由此坐标方位角条件就成了多边形的图形闭合条件。,2）纵、横坐标条件 设以、表示（图1）中各导线边的纵坐标增量之平差值；、表示（图1）中各导线边的横坐标增量之平差值；由图可写出以坐标增量平差值表示的纵、横坐标条件。,（1）,令,则,（2）,将上式代入式（2）得纵坐标条件式，且同理已可得横坐标的条件式即,（3）,上式就是单一符合导线的纵、横坐标条件方程、为条件式的不符值，按,式中、是由观测值计算的各导线点的近似坐标。,计算时一般 以秒为单位，、以cm为单位；若、以m为单位，则，

8、从而使全式单位统一。若单一导线的 与 点重合形成闭合导线，则纵、横坐标条件成为多边形各边的坐标增量闭合条件，以增量平差值表示为,2.符合导线的精度评定1）单位权中误差：单一符合导线计算单位权中误差公式与边角网相同，按,2）平差值的权函数式：为了平定平差值函数的精度，必须要列出权函数式。一般有下列三种函数式。边长平差值权函数式 由导线边 故其权函数式为,坐标平差值的权函数式 由（23页图1）得点坐标平差值的权函数式为,坐标方位角平差值权函数式 由（23页图1）得单一符合导线的任一边的坐标方位角的计算式为,3.附合导线按条件平差算例,在下图2所示附合导线中A,B 为已知点，其坐标为,方位角，应用红

9、外测距仪观测导线的转折角 和边长 列入下表1。试按条件平差法，求各观测值及平差后 边的边长相对中误差。,表1 近似坐标计算,解：（1）确定观测值的权。测角中误差,边长中误差按仪器给定公式为,（cm）,式中,以cm为单位。,由上式算得,cm,cm,cm,cm,以角度观测的权为单位权，即,表2 条件方程及权函数式系数表,续表,则边长的权为,边长权倒数为,（2）计算条件方程式不符值。由表1得,（3）计算条件方程式系数及权函数式系数，列于上表2中。,（4）组成法方程式并解算。根据上表2中系数组成法方程系数，然后填于下表3中相应行内。法方程式的解算在下表3中进行。,cm,cm,表3 法方程式解算表,（5

10、）计算改正数和平差值。由法方程解算表解得的联系数k 和观测边加相应改正数，即得角度和边长平差值。计算见下表4。,表4 观测值之平差值计算,（6）计算边的精度。1）单位权中误差，按,计算 边的中误差,2）,边长相对中误差,9-4 参数平差数学模型和公式,1.参数平差数学模型,（1）平差值方程,设平差问题中，有 个不等精度的独立观测，相应权为（1，2,），并设需 个必要观测，用 表示选定的未知数，按题列出 个平差值方程为,(2)误差方程,令 则1式为,其中（1，2，,）,则 2式的矩阵形式为,若设,（3）法方程式中有 个待定的改正数和 个未知数，共 个待定量，而方程只有 个，所以有无穷多组解。为了

11、寻求一组唯一的解，根据最小二乘原理按 的准则求，按数学上求函数自由极值的理论，即,转置后得,代3式代入4式得法方程,令,5式可表示为,其纯量形式为,将上式算得的 代入式2求出改正数向量V，进而求出观测平差值。,2、按间接平差法求平差值的计算步骤,根据平差问题的性质，选择 个独立量作为未知数；将每一个观测量的平差值分别表达成所选参数的函数。若函数为非线性要将其线性化，列出误差方程；由误差方程系数 和自由项 组成法方程，法方程个数等于未知数个数；解算法方程，求出未知数改正数，计算未知数；由误差方程计算V，求出观测量平差。,3、精度评定,（1）单位权中误差和 的计算同条件平差一样，间接平差单位权中误

12、差公式为,的计算方法为,或,由,直接计算,由未知数改正数 及法方程常数项 及计算,在高斯杜力特表中解算,（2）未知数函数的权倒数和中误差设某平差问题的未知数的函数为,它的权函数式为,则上式的矩阵形式为,根据权逆阵的传播律，得未知数的权倒数,则 的中误差为：,（3）未知数的权逆阵和中误差,由法方程,得,由权逆阵的传播律得,即法方程系数距阵的逆阵就是知数向量的权逆阵，令,因为法方程系数阵是一个对称方阵，故它的逆阵也为对称方阵未知数的权逆阵，通常又称为权系数阵，其对角线上的元素 为未知数的权倒数，非对角线上的元素 称为未知数 关于 的相关权倒数，而所有的元素又称为权系数。权系数的计算除了用矩阵求逆的

13、方法以外，还可以用高斯约化法求权系数的方法。则任一未知数 的中误差为：,9.5 高程网参数平差及算例,1、未知数个数的确定在间接平差中，未知数个数就等于必要观测个数。2、未知数的选取 在水准网中，即可以选取待定点高程作为未知数，也可选取高差作为未知数，但一般是选高程为未知数，例如在下图1中，就是选取待定点高程作为未知数的。3、算例在下图1所示水准网，已知水准点A的高程为 m，为求B，C，D三点的高程，进行了水准测量，及观测高差及水准路线长度载于下表1中，试按间接平差法求定B、C、D三点高程的平差值。,解：按题意知必要观测数 t=3，选取B、C、D三点高程、为参数 根据图示水准路线写出5个平差值

14、方程,则误差方程为,将观测高差和已知点高程代入上式，即可计算误差方程的常数项，此时，这些常数项具有56个数字，这对后续计算是不利的。为了便于计算，选取参数的近似值,这样，后续计算求定的只是未知数近似值改正数、，它们的关系为,将上二式代入误差方程，得,取10km的观测高差为单位权观测，即按,定权，得各观测值的权分别为,，,，,，,，,组成法方程为,解法方程，得,mm,mm，,mm，,代人误差方程得,mm，,mm，,mm，,mm，,最后得平差值,m，,m，,m，,m，,m，,m，,mm,m,m,9.6 三角网参数平差及算例,1.未知数的选定 平面控制网参数平差总是选择未知点x,y的坐标为平差参数。

15、,2、测角网坐标平差误差方程列立,这里讨论测角网中选择待定点坐标为未知数时，误差方程列立及线性化问题。如下图1为某一测角网的任一角 为三个待定点，它们的近似坐标为 改正数为，则平差值分别为,由图1可得 的平差值方程为,令,误差方程为,式中：,或写为：,讨论：（1）若测站点 为已知点时，则=0 有，得,若照准点 为已知点，则有=0，得,（3）若某边的两个端点均为已知点，则,=,=,=,=0,，,（4）同一边的正反坐标方位角的改正数相等，它们与坐标改正数的关系也一样。,即,因为,顾及,综上所述，对于角度观测的角网，采用间接平差，选择待定点的坐标为未知数时，列误差方程的步骤为：计算各待定点的近似坐标

16、；由待定点的近似坐标和已知点坐标计算各边的近似坐标方位角和近似边长；计算各边的坐标方位角改正数，并计算其系数列出误差方程。,据此，实际计算时，只要对每条待定边计算一个方向的坐标方位角改正数方程即可。,3、测边网坐标平差误差方程列立,如下图2为某一测边网中的任意一条边，为两个待定点，它们的近似坐标为，改正数为 则，的坐标平差值为,由右图2可写出 的平差值方程为,按台劳公式展开，得,讨论：（1）若 为已知点，则,若 为已知点，则,（2）若，均为知点，则该边为固定边，不需要列误差方程。,（3）某边的误差方程，按 方向列立与按 方向列立结果完全相同。,若按 方向，则,按 方向，则,顾及,4、三边网参数

17、平差算例,同精度测得如右图3中的三个边长,其结果为,m,m,m,已知点A，B，C 的起算数据列于下表1。试列出误差方程并求平差值。,表1,解:（1）本题,选择待定点的坐标和为未知数,其近似值和由已知点、和观测边、交会计算而得。图4中，设为三角形底边上的高，为在上的投影。得,按此,计算待定点的近似坐标,m,m,m,m,图4,（2）根据近似坐标和已知点坐标计算误差方程的系数和常数项，计算数值列在表2中。,表2,由表2的最后三列数值,按式写出误差方程,（3）求平差值,根据误差方程组成法方程,解算法方程,得m,m。待定点坐标平差值为,将 代入误差方程计算观测值的改正数并从而求得平差值,（m），,习 题

18、,1、在图1的水准网中，观测高差及路线长度见下表：,已知A、B点高程为：,=50.000m,=40.000m,试用条件平差法求,(1)、各观测高差的平差值；,(2)、平差后到点间高差的中误差。,2、有闭合导线如下图2，观测4条边长和5个左转折角，已知测角中误差，边长中误差 mm,计算（以m为单位），起算数据为：,观测值如下：,试按条件平差:,(1)、列条件方程和法方程；,(2)、求改正数和平差值；,(3)、求导线点2、3、4的坐标平差值。,3、在下图3的水准网中，为已知点，m，观测各点间的高差为：,m、,m、,m、,m、,设，试按间接平差法求：,（）待定点、的高程平差值；,（）平差后至点间高差平差值及中误差。,4、如图4的测角网中，为已知三角点，点为待定点，起算数据如下表：,起 算 数 据,试用坐标平差法求待定点的坐标值及其点位中误差。,5、在图5中测边网中，点为已知点，点为待定点。已知点坐标为，同精度测得边长观测值为：,试按间接平差求点坐标平差值及其协因数阵。,设待定点的近似坐标为：,