矿业大学环境与测绘学院测绘工程《测量平差》第三章条件平差-2521演示教学.ppt
《矿业大学环境与测绘学院测绘工程《测量平差》第三章条件平差-2521演示教学.ppt》由会员分享,可在线阅读,更多相关《矿业大学环境与测绘学院测绘工程《测量平差》第三章条件平差-2521演示教学.ppt(204页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、矿业大学环境与测绘学院测绘工程测量平差第三章条件平差-2521一、条件平差原理一、条件平差原理 设在某个测量作业中,有设在某个测量作业中,有n个观测值个观测值 ,均含有相,均含有相互独立的偶然误差,相应的权阵为互独立的偶然误差,相应的权阵为 ,改正数为,改正数为 ,平差值为,平差值为 ,表示为,表示为其中其中:为对角阵;为对角阵;即即(3-1-6)在这在这n个观测值中,有个观测值中,有t个必要观测数,多余观测个必要观测数,多余观测数为数为r。可以列出。可以列出r个平差值线性条件方程个平差值线性条件方程(3-1-7)式中,式中,(i=1,2,n)为各平差值条件方)为各平差值条件方程式中的系数,程
2、式中的系数,为各平差值条件方程式中为各平差值条件方程式中的常数项。的常数项。将(将(3-1-6)式代入)式代入(3-1-7)式,得相应的改正数式,得相应的改正数条件方程式条件方程式(3-1-8)式中式中 称为改正数条件方程的闭合差称为改正数条件方程的闭合差(或不符值),即(或不符值),即(3-1-9)若取若取 (3-1-7)、(、(3-1-8)和()和(3-1-9)式可分别表达)式可分别表达成矩阵形式如下成矩阵形式如下(3-1-10)(3-1-11)(3-1-12)按求函数极值的拉格朗日乘数法,引入乘系数按求函数极值的拉格朗日乘数法,引入乘系数(又称为联系数向量)(又称为联系数向量),构成函数
3、:构成函数:(3-1-13)为引入最小二乘法为引入最小二乘法,将将对对V求一阶导数求一阶导数,并令其为零并令其为零得得上式两端转置,得上式两端转置,得由于由于P是主对角线阵,则是主对角线阵,则 ,得,得 将上式两边左乘权逆阵将上式两边左乘权逆阵 ,得,得(3-1-14)此式称为改正数方程,其纯量形式为此式称为改正数方程,其纯量形式为(3-1-15)将(将(3-1-14)式代入()式代入(3-1-11)式,得)式,得 此式称为联系数法方程(简称法方程),其纯量此式称为联系数法方程(简称法方程),其纯量形式为形式为(3-1-16)(3-1-17)取法方程的系数阵取法方程的系数阵 ,由上式易知,由上
4、式易知N阵阵关于主对角线对称,得法方程表达式关于主对角线对称,得法方程表达式(3-1-18)上式两边左乘法方程系数阵上式两边左乘法方程系数阵N的逆阵的逆阵 ,得联系,得联系数数K的唯一解:的唯一解:(3-1-19)即,即,N是一个是一个r阶的满秩方阵,且可逆。将(阶的满秩方阵,且可逆。将(3-1-18)式移项,得式移项,得法方程数阵法方程数阵N的秩的秩 将(将(3-1-19)式代入()式代入(3-1-14)或()或(3-1-15)式,可计算出式,可计算出V,再将,再将V代入(代入(3-1-6),即可计算出即可计算出所求的观测值的最或然值所求的观测值的最或然值 。通过观测值的平差值通过观测值的平
5、差值 ,可以进一步计算一,可以进一步计算一些未知量(如待定点的高程、纵横坐标以及边些未知量(如待定点的高程、纵横坐标以及边的长度、某一方向的方位角等)的最或然值。的长度、某一方向的方位角等)的最或然值。由上述推导可看出,由上述推导可看出,K、V及及 都是由(都是由(3-1-11)和()和(3-1-14)式解算出的,因此我们把()式解算出的,因此我们把(3-1-11)和)和(3-1-14)式合称为条件平差的基础方程。式合称为条件平差的基础方程。二、精度评定二、精度评定 在第一个问题中已经阐述了计算未知量最或然值的在第一个问题中已经阐述了计算未知量最或然值的原理和公式,下面来论述测量平差的第二个任
6、务原理和公式,下面来论述测量平差的第二个任务,即即评定测量成果的精度。精度评定包括单位权方差评定测量成果的精度。精度评定包括单位权方差 和单位权中误差和单位权中误差 的计算、平差值函数的计算、平差值函数()的的协因数协因数 及其中误差及其中误差 的计算等。的计算等。在第二章中已经介绍过,当已知单位权方差在第二章中已经介绍过,当已知单位权方差 时,时,如果知道某量的权为如果知道某量的权为p,则该量的方差为,则该量的方差为 在实际工作中,由于观测值的个数在实际工作中,由于观测值的个数n是有限值,因此,是有限值,因此,只能求出只能求出 的估值的估值 和和 的估值的估值 。则有。则有 根据协因数的定义
7、,有了单位权方差根据协因数的定义,有了单位权方差 和某平和某平差值函数的验后协因数阵差值函数的验后协因数阵 ,也可按下式计算该,也可按下式计算该平差值向量的协方差阵。平差值向量的协方差阵。(3-1-22)估值形式为估值形式为(3-1-20)(3-1-21)例如,已知观测值的平差值例如,已知观测值的平差值 的协因数阵的协因数阵 ,则则 的协方差阵为的协方差阵为 下面,我们分别讨论单位权中误差下面,我们分别讨论单位权中误差 和平差值和平差值函数协因数阵函数协因数阵 的计算方法。的计算方法。1.计算单位权方差和中误差的估值计算单位权方差和中误差的估值 根据第二章中对中误差的定义,单位权中误差的根据第
8、二章中对中误差的定义,单位权中误差的计算公式为计算公式为 在一般情况下,观测值的真误差在一般情况下,观测值的真误差 是不知道的,是不知道的,也就不可能利用上式计算单位权中误差。但在条件也就不可能利用上式计算单位权中误差。但在条件平差中,可以通过观测值的改正数平差中,可以通过观测值的改正数V来计算单位权来计算单位权方差和中误差:方差和中误差:(3-1-23)(2)由()由(3-1-14)和()和(3-1-11)式导出)式导出 式中式中r为多余观测值个数,为多余观测值个数,r=n t。在(在(3-1-24)中,须先算出)中,须先算出VTPV的值,才能计算的值,才能计算单位权中误差。单位权中误差。V
9、TPV可用下列几种方法计算:可用下列几种方法计算:(1)直接利用定义式()直接利用定义式(3-1-23)计算。)计算。纯量形式为纯量形式为(3-1-25)(3-1-24)即即(3-1-26)其纯量形式为其纯量形式为(3-1-27)2、协因数阵、协因数阵条件平差的基本向量条件平差的基本向量L、W、K、V、都可以表都可以表达成随机向量达成随机向量L的函数的函数将向量将向量L、K、V、组成列向量,并以组成列向量,并以Z表示之表示之(3-1-28)式中等号右端第二项是与观测值无关的常数项阵,按式中等号右端第二项是与观测值无关的常数项阵,按协因数传播律,得协因数传播律,得Z的协因数阵为的协因数阵为整理后
10、得整理后得(3-1-29)由上式可见,平差值由上式可见,平差值 与闭合差与闭合差W、联系数、联系数K、改、改正数正数V是不相关的统计量,又由于它们都是服从正态是不相关的统计量,又由于它们都是服从正态分布的向量,所以分布的向量,所以 与与W、K、V也是相互独立的向也是相互独立的向量。量。3平差值函数的协因数平差值函数的协因数 在条件平差中,平差计算后,首先得到的是各个在条件平差中,平差计算后,首先得到的是各个观测量的平差值。例如,水准网中的高差观测值的观测量的平差值。例如,水准网中的高差观测值的平差值,测角网中的观测角度的平差值,导线网中平差值,测角网中的观测角度的平差值,导线网中的角度观测值和
11、各导线边长观测值的平差值等。而的角度观测值和各导线边长观测值的平差值等。而我们进行测量的目的,往往是要得到待定水准点的我们进行测量的目的,往往是要得到待定水准点的高程值、未知点的坐标值、三角网的边长值及方位高程值、未知点的坐标值、三角网的边长值及方位角值等,并且评定其精度。这些值都是关于观测值角值等,并且评定其精度。这些值都是关于观测值平差值的函数。平差值的函数。设有平差值函数设有平差值函数(3-1-30)对上式全微分得对上式全微分得(3-1-31)取全微分式的系数阵为取全微分式的系数阵为由协因数传播律得由协因数传播律得(3-1-32)(3-1-33)根据(根据(3-1-29)式,知)式,知代
12、入(代入(3-1-33)式得)式得 将(将(3-1-34)式代入()式代入(3-1-22)式,可求得该平)式,可求得该平差值函数的方差差值函数的方差(3-1-35)即即(3-1-34)此式即为平差值函数式(此式即为平差值函数式(3-1-30)的协因数表达式。)的协因数表达式。三、条件平差的计算步骤三、条件平差的计算步骤 综合以上所述,按条件平差的计算步骤可归结为综合以上所述,按条件平差的计算步骤可归结为以下几步:以下几步:(1)根据实际问题,确定出总观测值的个数)根据实际问题,确定出总观测值的个数n、必、必要观测值的个数要观测值的个数t及多余观测个数及多余观测个数r=n-t,进一步列,进一步列
13、出最或是值条件方程出最或是值条件方程(3-1-10)或改正数条件方程(或改正数条件方程(3-1-11););(2)根据()根据(3-1-16)式,组成法方程式;)式,组成法方程式;(3)依据()依据(3-1-19)式计算出联系数)式计算出联系数K;(4)由()由(3-1-14)式计算出观测值改正数)式计算出观测值改正数V;并依;并依据(据(3-1-6)式计算出观测值的平差值;)式计算出观测值的平差值;(5)根据()根据(3-1-23)和()和(3-1-24)计算单位权方)计算单位权方差差 和单位权中误差和单位权中误差 ;(6)列出平差值函数关系式()列出平差值函数关系式(3-1-30),并对其
14、),并对其全微分,求出其线性函数的系数阵全微分,求出其线性函数的系数阵f,利用(,利用(3-1-34)式计算出平差值函数的协因数)式计算出平差值函数的协因数 ,代入(,代入(3-1-22)计算出平差值函数的协方差)计算出平差值函数的协方差 。为了检查平差计算的正确性,可以将平差值为了检查平差计算的正确性,可以将平差值 代代入平差值条件方程式(入平差值条件方程式(3-1-10),看是否满足方程),看是否满足方程关系。关系。例例3-1 如图如图3-1所示,所示,A和和P点为等级三角点,点为等级三角点,PA方方向的方位角已知,在测站向的方位角已知,在测站P上等精度测得的各方向的上等精度测得的各方向的
15、夹角观测值如下:夹角观测值如下:T PA=482436L 1=573216L 2=730308L 3=1265128L 4=1043320试用条件平差法,计算各观测值的平差值、试用条件平差法,计算各观测值的平差值、PC方方向的方位角向的方位角TPC,及,及TPC的精度的精度 。解:本题中解:本题中n=4,t=3,则条件方程个数为,则条件方程个数为 r=n t=1。因为是等精度观测,取因为是等精度观测,取观测值权阵观测值权阵由由 ,列出平差值条件方程的纯量形式,列出平差值条件方程的纯量形式其矩阵形式为其矩阵形式为由,由,计算闭合差计算闭合差由由 ,写出改正数条件方程式,写出改正数条件方程式其纯量
16、形式为其纯量形式为根据根据 ,写出法方程,写出法方程4 ka -12=0纯量形式为纯量形式为4 ka +12=0由由 ,计算联系数,计算联系数 k a =-0.25 12=-3其纯量形式为其纯量形式为k a=-3由由 ,计算各改正数,计算各改正数由由 =+,计算观测值平差值,计算观测值平差值由由(3-1-24)式,计算单位权中误差式,计算单位权中误差其矩阵式为其矩阵式为其中系数阵为其中系数阵为 PC边的方位角边的方位角=482436+573213+730305=1785954计算计算PC边的协因数边的协因数则则PC边方位角的中误差为边方位角的中误差为 3-2 高程网条件平差高程网条件平差v一、
17、高程网条件方程的个数及条件方程式一、高程网条件方程的个数及条件方程式 高程网包括水准网和三角高程网。对高程网进行高程网包括水准网和三角高程网。对高程网进行条件平差时,一般以已知高程点的高程值作为起算条件平差时,一般以已知高程点的高程值作为起算数据,以各测段的观测高差值作为独立观测值,写数据,以各测段的观测高差值作为独立观测值,写出其满足的条件关系式,按照条件平差的原理解算出其满足的条件关系式,按照条件平差的原理解算各高差值的改正数和平差值,然后再计算出各待求各高差值的改正数和平差值,然后再计算出各待求点的高程平差值,并进行精度评定。点的高程平差值,并进行精度评定。进行条件平差时,首先要确定条件
18、方程的个数。进行条件平差时,首先要确定条件方程的个数。从上节内容可知道,在一般情况下,条件方程式的个从上节内容可知道,在一般情况下,条件方程式的个数与多余观测的个数数与多余观测的个数r相符。而要确定多余观测个数就相符。而要确定多余观测个数就必须先确定必要观测个数必须先确定必要观测个数t。高程测量(包括三角高程测量和水准测量)的主要高程测量(包括三角高程测量和水准测量)的主要目的是确定未知点的高程值。如图目的是确定未知点的高程值。如图3-2所示高程网中,所示高程网中,有有2个已知高程点个已知高程点A、B,3个未知高程点个未知高程点C、D、E和和8个高差观测值。从图中可以看出,要确定个高差观测值。
19、从图中可以看出,要确定3个未知点个未知点的高程值,至少需要知道其中的的高程值,至少需要知道其中的3个高差观测值(如个高差观测值(如h1、h2、h3,或,或h6、h7、h8,或,或h2、h4、h5 等等多种选择),即必要观测个数多种选择),即必要观测个数t=3 则多余观测个数则多余观测个数r=n t=8-3=5,可,可以写出这以写出这5个条件方个条件方程式程式相对应的改正数条件方程式形式相对应的改正数条件方程式形式 其中其中 这些条件方程式(或改正数条件方程式),大体这些条件方程式(或改正数条件方程式),大体上分为两类:其一是闭合路线情况,如条件方程式中上分为两类:其一是闭合路线情况,如条件方程
20、式中前四个条件方程式,可称为闭合条件方程式;其二是前四个条件方程式,可称为闭合条件方程式;其二是附合路线情况,如条件方程式中第五个,反应的是从附合路线情况,如条件方程式中第五个,反应的是从A点出发后测得的点出发后测得的B点的高程值是否与点的高程值是否与B点的已知高程值点的已知高程值相等的问题,可称为附合条件方程式。相等的问题,可称为附合条件方程式。再如图再如图3-3所示高程网中,有所示高程网中,有4个已知高程点、个已知高程点、4个个未知高程点和未知高程点和8个高差观测值,即个高差观测值,即n=8。则必要观测。则必要观测个数为个数为t=4,多余观测个数为,多余观测个数为r=n t=4。可以写出。
21、可以写出4个最或是值条件方程式或个最或是值条件方程式或4个改正数条件方程式,其个改正数条件方程式,其中有中有1个闭合条件方程式和个闭合条件方程式和3个附合条件方程式。个附合条件方程式。对应的改正数条件方程式为对应的改正数条件方程式为其中其中 如果高程网中没有必要的起算数据(即没有已知高如果高程网中没有必要的起算数据(即没有已知高程点,如图程点,如图3-4所示),就不存在某点高程值的已知所示),就不存在某点高程值的已知值与计算值的矛盾,也就不存在附合条件方程式,只值与计算值的矛盾,也就不存在附合条件方程式,只能写出有关的闭合条件方程式。为确定条件方程式个能写出有关的闭合条件方程式。为确定条件方程
22、式个数或者多余观测个数,可将高程网中一个点的高程值数或者多余观测个数,可将高程网中一个点的高程值假定为已知值,再按上文中所述方法操作。由此不难假定为已知值,再按上文中所述方法操作。由此不难判断图判断图3-4所示的高程网中有所示的高程网中有4个必要观测值,个必要观测值,4个多个多余观测值,可写出余观测值,可写出4个闭合条件方程式个闭合条件方程式相对应的改正数条件方程式形式相对应的改正数条件方程式形式 其中其中 二、高程网平差举例二、高程网平差举例 图图3-5为一水准网,为一水准网,A、B为两个高程已知点,为两个高程已知点,C、D、E、F分别为待定点。已知高程值和高差观测值如分别为待定点。已知高程
23、值和高差观测值如表表3-1中所示,计算各待定点的高程平差值。中所示,计算各待定点的高程平差值。解:水准网中总观测个数解:水准网中总观测个数n=8,必要观测,必要观测数数t=4,多余观测,多余观测r=n t=4。平差值条件方程式平差值条件方程式 为为改正数条件方程式改正数条件方程式 为为由条件方程得:由条件方程得:令令C=1,观测值的权倒数:,观测值的权倒数:3-3 导线网条件平差计算导线网条件平差计算v一单一附合导线条件平差一单一附合导线条件平差 导线网,包括单一附合导线、单一闭合导线和结导线网,包括单一附合导线、单一闭合导线和结点导线网,是目前较为常用的控制测量布设方式之点导线网,是目前较为
24、常用的控制测量布设方式之一,其观测值有长度观测值和角度观测值。在本节一,其观测值有长度观测值和角度观测值。在本节中我们主要讨论单一导线的平差计算,先讨论单一中我们主要讨论单一导线的平差计算,先讨论单一附合导线问题。附合导线问题。如图如图3-6所示,在这个导线中有四个已知点、所示,在这个导线中有四个已知点、n-1个未知点、个未知点、n+1个水平角观测值和个水平角观测值和n条边长观测值,条边长观测值,总观测值数为总观测值数为2n+1。从图中可以分析,要确定一个。从图中可以分析,要确定一个未知点的坐标,必须测一条导线边和一个水平角,即未知点的坐标,必须测一条导线边和一个水平角,即需要两个观测值;需要
25、两个观测值;要确定全部要确定全部n-1个未知点,则需观测个未知点,则需观测n-1个导线边和个导线边和n-1个水平角,即必要观测值数个水平角,即必要观测值数t=2n-2;则多余观测;则多余观测个数个数r=(2n+1)t=3。也就是说,在单一附合导线。也就是说,在单一附合导线中,只有三个条件方程。下面讨论其条件方程式及改中,只有三个条件方程。下面讨论其条件方程式及改正数条件方程式的写法。正数条件方程式的写法。设设AB边方位角已知值为边方位角已知值为TAB=T0,CD边方位角已边方位角已知值为知值为TCD、计算值为、计算值为Tn+1,B点坐标的已知值为点坐标的已知值为(,)或者或者(x1,y1),C
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 测量平差 矿业大学 环境 测绘 学院 工程 测量 第三 条件 2521 演示 教学
限制150内