土木工程測量学教程(下)教亲 第3讲 教学目标:掌握三角网条件评差方法和程序。 重点难点:条件方程列立,闭合差检核 第3章控制网平差 3-1概述 在测量工作中,常要确定某些几何量的大小。由几何量组成的模型称为几何模型。 为了确定一个几何模型,并不需要知道该模型中所有元素的大小,而只需要知道其中部 分元素的大小就行了,其它元素可以通过它们来确定。能够唯一地确定一个几何模型所必要 的元素,简称必要元素:必要元素的个数用t来表示。必要元素不仅要考虑其个数,而且要 考虑它的类型。由此可知,当某个几何模型给定之后,能够唯一确定该模型的必要元素的个 数t及其类型,t只与几何模型有关,与实际观测无关 约定:一一真值,L一一平差值, L一一观测值 在一个几何模型中,除了t个独立量以外,若再增加一个量,则必然产生一个相应函数 关系式。例,必要量选为L1、L2、S1,若增加一个量3 则存在 L1+L2+L3=180°, 若再增加一个量S2,则有S2=52 由此可知,一个几何模型的独立量个数最多为t个,除此之外,增加一个量必然要产生一个 相应的函数关系式,这种函数关系式,在测量平差中称为条件方程。 在测量工程中,为了求得一个几何模型中各量的大小就必须进行观测。如果总共观测了 该模型中n个量的大小,若观测个数少于必要元素的个数,即n<t,显然它无法确定该模型, 即出现了数据不足的情况:若观测了个独立量,n=t,则可唯一地确定该模型。由于它们都 是独立量,故不存在任何条件方程,在这种情况下,如果观测结果中含有粗差甚至错误,都 将无法发现,在测量工作中是不允许这样做的。为了能及时发现粗差和错误,并提高测量成 果的精度,就必须使n>t,若令 (3-1) 式中n为观测值的个数,称为必要观测数,r称为多余观测数。多余观测数在测量中又称“自 由度”。一个几何模型如果有r个多余观测,就产生r个条件方程 由于观测下不可避免地存在偶然误差,当nt时,几何模型中应该满足r=n-1个条件方
土木工程测量学教程(下)教案 3-1 第3讲 教学目标:掌握三角网条件评差方法和程序。 重点难点:条件方程列立,闭合差检核 第 3 章 控制网平差 3—1 概述 在测量工作中,常要确定某些几何量的大小。由几何量组成的模型称为几何模型。 为了确定一个几何模型,并不需要知道该模型中所有元素的大小,而只需要知道其中部 分元素的大小就行了,其它元素可以通过它们来确定。能够唯一地确定一个几何模型所必要 的元素,简称必要元素;必要元素的个数用 t 来表示。必要元素不仅要考虑其个数,而且要 考虑它的类型。由此可知,当某个几何模型给定之后,能够唯一确定该模型的必要元素的个 数 t 及其类型,t 只与几何模型有关,与实际观测无关。 约定: L ~ ――真值, L ――平差值, L――观测值 在一个几何模型中,除了 t 个独立量以外,若再增加一个量,则必然产生一个相应函数 关系式。例,必要量选为 1 ~ L 、 2 ~ L 、 1 ~ S ,若增加一个量 3 ~ L , 则存在 180 ~ ~ ~ L1 + L2 + L3 = , 若再增加一个量 2 ~ S ,则有 1 2 2 1 ~ sin ~ ~ ~ sin L L S = S 由此可知,一个几何模型的独立量个数最多为 t 个,除此之外,增加一个量必然要产生一个 相应的函数关系式,这种函数关系式,在测量平差中称为条件方程。 在测量工程中,为了求得一个几何模型中各量的大小就必须进行观测。如果总共观测了 该模型中 n 个量的大小,若观测个数少于必要元素的个数,即 n<t,显然它无法确定该模型, 即出现了数据不足的情况;若观测了 t 个独立量,n=t,则可唯一地确定该模型。由于它们都 是独立量,故不存在任何条件方程,在这种情况下,如果观测结果中含有粗差甚至错误,都 将无法发现,在测量工作中是不允许这样做的。为了能及时发现粗差和错误,并提高测量成 果的精度,就必须使 n>t,若令 r=n–t (3-1) 式中 n 为观测值的个数,t 称为必要观测数,r 称为多余观测数。多余观测数在测量中又称“自 由度”。一个几何模型如果有 r 个多余观测,就产生 r 个条件方程。 由于观测下不可避免地存在偶然误差,当 n>t 时,几何模型中应该满足 r=n-t 个条件方
土木工程測量学教程(下)教亲 程,实际存在闭合差而并不满足,如何调整观测值,即对观测值合理地加上改正数,使其达 到消除闭合差的目的,这是测量平差的主要任务 个测量平差问题,首先要由观测值和待求量间组成数学模型,然后采用一定的平差原 则对待求量进行估计,这种估计要求是最优的,最后计算和分析成果的精度 3-2测量平差的数学模型 在测量工程中,涉及的是通过观测量来确定某些几何量或物理量大小等有关的数量问 题,因而考虑的模型总是数学模型。平差的数学模型与一般数学只考虑函数模型不同,它还 要考虑随杋模型,因为观测量是一种随机变量。所以平差的数学模型同时包含函数模型和随 机模型两种,在研究任何平差方法时必须予以考虑。 函数模型是由描述观测量和待求量间的数学函数关系的模型,是确定客观实际的本质或 特征的模型。随机模型是描述观测量及其相互间统计相关性质的模型。对于一个实际平差问 题,可建立不同形式的函数模型,与此相应,就产生了不同的平差方法 、条件平差法 以条件方程为函数模型的平差方法,称为条件平差法 般而言,如果有n个观测值L,t个必要观测,则应列出r=n-1个条件方程 +1.=g (3-8) A为常数向量,将L=L+y代入,并令 W=AL+a (3-9) 则(3-8)式为 A+W=0 (3-8)或(3-10)式为条件平差的函数模型 条件平差的自由度极为多余观测数r,即条件方程的个数 、间接平差法 在一个几何模型中,最多只能选出t个独立量,如果在进行平差就选定t个独立量为参 数,那么通过这1个独立参数就能唯一地确定几何模型了。换言之,模型中的所有量都一定 是这t个独立参数的函数,亦即每个观测量都可表达成所选t个独立参数的函数。 选择几何模型中1个独立量为平差参数,将每一个观测量表达成所选参数的函数,即列 出n个这种函数关系式,以此为平差的函数模型,称为间接平差法,又称参数平差法
土木工程测量学教程(下)教案 3-2 程,实际存在闭合差而并不满足,如何调整观测值,即对观测值合理地加上改正数,使其达 到消除闭合差的目的,这是测量平差的主要任务。 一个测量平差问题,首先要由观测值和待求量间组成数学模型,然后采用一定的平差原 则对待求量进行估计,这种估计要求是最优的,最后计算和分析成果的精度。 3—2 测量平差的数学模型 在测量工程中,涉及的是通过观测量来确定某些几何量或物理量大小等有关的数量问 题,因而考虑的模型总是数学模型。平差的数学模型与一般数学只考虑函数模型不同,它还 要考虑随机模型,因为观测量是一种随机变量。所以平差的数学模型同时包含函数模型和随 机模型两种,在研究任何平差方法时必须予以考虑。 函数模型是由描述观测量和待求量间的数学函数关系的模型,是确定客观实际的本质或 特征的模型。随机模型是描述观测量及其相互间统计相关性质的模型。对于一个实际平差问 题,可建立不同形式的函数模型,与此相应,就产生了不同的平差方法。 一、条件平差法 以条件方程为函数模型的平差方法,称为条件平差法 。 一般而言,如果有 n 个观测值 n,1 L ,t 个必要观测,则应列出 r=n–t 个条件方程, r n A , ,1 , ,1 ~ r r n n L+ A = O (3-8) A 为常数向量, 将 L = L + v 代入,并令 W = AL + A (3-9) 则(3-8)式为 AV +W = 0 (3-10) (3-8)或(3-10)式为条件平差的函数模型。 条件平差的自由度极为多余观测数 r,即条件方程的个数。 二、间接平差法 在一个几何模型中,最多只能选出 t 个独立量,如果在进行平差就选定 t 个独立量为参 数,那么通过这 t 个独立参数就能唯一地确定几何模型了。换言之,模型中的所有量都一定 是这 t 个独立参数的函数,亦即每个观测量都可表达成所选 t 个独立参数的函数。 选择几何模型中 t 个独立量为平差参数,将每一个观测量表达成所选参数的函数,即列 出 n 个这种函数关系式,以此为平差的函数模型,称为间接平差法,又称参数平差法
土木工程測量学教程(下)教亲 尽管间接平差法是选了t个独立参数,但多余观测数不随平差不同而异,其自由度仍是 =n-t。 3-3条件平差原理 某一平差问题,如必要观测数为t,则只能选出t个独立的未知数,每增加一个未知数, 就会产生一个应满足的条件方程。在条件平差法中选n个观测量的平差值L作为未知数,由 于多选了r=n-1个未知数,因而就产生了r个条件方程。条件平差法就是在满足r个条件方 程的需求下,求函数P=min的V值,在数学中是求函数的条件极值问题。 基础方程和它的解 设有r个平差值线性条件方程 1L1 L bL1+b2L2+…+bnLn+b=0 (3-23) +rl+ 式中a、b、…、n(i=1,2,…;n)为条件方程系数,a0、b、…、m为条件方程常数项。将(3-22) 式代入上式,得条件方程为 a,v,+a2v2+.+a,v-W=o bV1+b2V2+…+b 0 (3-24) rv,+rv2+.+r,,-Wr=0 式中w、wb、…、W称为条件方程的闭合差,或称不符值,其值为 (a,,L +a,L W,=-(b, L+b,L2+.+b, Ln+bo) (3-25 (1L1+r2L2 Ln +ro) Li 6. b bn W V2
土木工程测量学教程(下)教案 3-3 尽管间接平差法是选了 t 个独立参数,但多余观测数不随平差不同而异,其自由度仍是 r=n-t。 3—3 条件平差原理 某一平差问题,如必要观测数为 t,则只能选出 t 个独立的未知数,每增加一个未知数, 就会产生一个应满足的条件方程。在条件平差法中选 n 个观测量的平差值 L ˆ 作为未知数,由 于多选了 r=n–t 个未知数,因而就产生了 r 个条件方程。条件平差法就是在满足 r 个条件方 程的需求下,求函数 V TPV=min 的 V 值,在数学中是求函数的条件极值问题。 一. 基础方程和它的解 设有 r 个平差值线性条件方程: + + + + = + + + + = + + + + = ˆ ˆ ˆ 0 ˆ ˆ ˆ 0 ˆ ˆ ˆ 0 1 1 2 2 0 1 1 2 2 0 1 1 2 2 0 r L r L r L r b L b L b L b a L a L a L a n n n n n n (3-23) 式中 ai、bi、…、ri(i=1,2,…,n)为条件方程系数,a0、b0、…、r0 为条件方程常数项。将(3-22) 式代入上式,得条件方程为 + + + − = + + + − = + + + − = 0 0 0 1 1 2 2 1 1 2 2 1 1 2 2 n n r n n b n n a r v r v r v w b v b v b v w a v a v a v w (3-24) 式中 wa 、 wb 、…、 wr 称为条件方程的闭合差,或称不符值,其值为 = − + + + + = − + + + + = − + + + + ( ) ( ) ( ) 1 1 2 2 0 1 1 2 2 0 1 1 2 2 0 w r L r L r L r w b L b L b L b w a L a L a L a r n n b n n a n n (3-25) 设 = n n n r n r r r b b b a a a A 1 2 1 2 1 2 , , = r b a r w w w W,1 , = n n v v v V 2 1 ,1 , = 0 0 0 ,1 r b a A r , = n n L L L L 2 1 ,1
土木工程測量学教程(下)教亲 则(3-24)、(3-25)式的矩阵表达式为 A+W=0 (3-26) 按求函数极值的拉格朗日乘数法,设其乘数为K=[。kb j,称为联系数向 量。组成函数 Φ=IP-2K(A+W) 将Φ对V求一阶导数,并令其为零,得 dv= 2VTP-2KA=O 将上式两边左乘权逆阵P1,得改正数V的计算公式为 V=PA K 上式称为改正数方程 将n个改正数方程(3-28)和r个条件方程(3-26)联立求解,就可以得一组唯一的解 n个改正数和r个联系数。 将(3-26)和(3-28)式合称为条件平差的基础方程。显然,由基础方程解出的一组改 正数V不仅能消除闭合差,也必能满足VPI=mn的要求。 解算基础方程式,是先将(3-28)式代入(3-26)式,得 AOAK+M 令 则有 (3-31) 上式称为联系数法方程,简称法方程。法方程系数阵N的秩R(N)=R(AQA)=R(A)=F, 即N是一个r阶的满秩方阵,且可逆。由此可得联系数K的唯一解。 解出联系K后,将K值代入改正数方程式,求出改正数Ⅴ值,再求平差值L=L+V 这样就完成了按条件平差求平差值的工作 三、三角网的条件方程式 角网有独立三角网(网中仅有必要的起算数据)和非独立三角网(网中具有多余的起 算数据)之分。三角网平差有按角度平差和按方向平差两种方法。以下仅讨论独立三角网的 条件和条件方程式
土木工程测量学教程(下)教案 3-4 则(3-24)、(3-25)式的矩阵表达式为 AV +W = 0 (3-26) W = AL + A0 (3-27) 按求函数极值的拉格朗日乘数法,设其乘数为 T a b r r K = k k k ,1 ,称为联系数向 量。组成函数 V PV 2K (AV W) T T = − + 将 对 V 求一阶导数,并令其为零,得 = 2 − 2 = 0 V P K A dV d T T 即 PV A K T = 将上式两边左乘权逆阵 P -1,得改正数 V 的计算公式为 V P A K −1 T = (3-28) 上式称为改正数方程 将 n 个改正数方程(3-28)和 r 个条件方程(3-26)联立求解,就可以得一组唯一的解: n 个改正数和 r 个联系数。 将(3-26)和(3-28)式合称为条件平差的基础方程。显然,由基础方程解出的一组改 正数 V 不仅能消除闭合差,也必能满足 V PV = min T 的要求。 解算基础方程式,是先将(3-28)式代入(3-26)式,得 AQA K +W = 0 T 令 T n r T n n r r r n N A Q A AP A 1 , , , , − = = (3-30) 则有 0 , ,1 ,1 + = r r r r N K W (3-31) 上式称为联系数法方程,简称法方程。法方程系数阵 N 的秩 R N R AQA R A r T ( ) = ( ) = ( ) = , 即 N 是一个 r 阶的满秩方阵,且可逆。由此可得联系数 K 的唯一解。 解出联系 K 后,将 K 值代入改正数方程式,求出改正数 V 值,再求平差值 L ˆ = L +V , 这样就完成了按条件平差求平差值的工作。 三、三角网的条件方程式 三角网有独立三角网(网中仅有必要的起算数据)和非独立三角网(网中具有多余的起 算数据)之分。三角网平差有按角度平差和按方向平差两种方法。以下仅讨论独立三角网的 条件和条件方程式
土木工程測量学教程(下)教亲 1.图形条件 图3-4所示的三角网中,有五个独立的三角形图形条件,每个三角形都应列出其图形条 件方程式。对于第个三角形而言,设a、b、c1为各三角形内角的平差值,则其图形条 件方程为 a1+b2+c1-180=0。 3-33) 设各三角形内角的观测值为a、b、c,则有 ai=ai tai bi=bi+vbi (3-34) +vci 上式代入(3-33)式,整理后可得 Vaitbitvatwi=0 W1=a1+b1+c1-180° (3-35) 式中w1是三角形图性条件的闭合差,简称三角形闭合差。 大地四边形每个三角形都可以组成一个图形条件方程式。但是,大地四边形中独立的图 形条件方程式只有三个,第四个图形条件方程式可由三个独立的图形条件方程式导出 图3-4 图3-5 图形条件闭合差的限值: 设各角度观测中误差为m,对式(3-35)应用误差传播定律可得闭合差w的中误差为 取闭合差中误差的两倍作为闭合差的限值,即 m=2m3 同理可得,n条边的多边形图形条件闭合差的限值为 w=2m√n 3-37)
土木工程测量学教程(下)教案 3-5 1. 图形条件 图 3-4 所示的三角网中,有五个独立的三角形图形条件,每个三角形都应列出其图形条 件方程式。对于第 i 个三角形而言,设 i a 、 i b 、 i c 为各三角形内角的平差值,则其图形条 件方程为 + + −180 = 0 i i i a b c 。 (3–33) 设各三角形内角的观测值为 ai、bi、ci ,则有 = + = + = + i i ci i i bi i i ai c c v b b v a a v (3–34) 上式代入(3–33)式,整理后可得 = + + − + + + = 180 0 i i i i ai bi ci i w a b c v v v w (3-35) 式中 wi 是三角形图性条件的闭合差,简称三角形闭合差。 大地四边形每个三角形都可以组成一个图形条件方程式。但是,大地四边形中独立的图 形条件方程式只有三个,第四个图形条件方程式可由三个独立的图形条件方程式导出。 图形条件闭合差的限值: 设各角度观测中误差为 m,对式(3-35)应用误差传播定律可得闭合差 w 的中误差为 mw = m 3 取闭合差中误差的两倍作为闭合差的限值,即 w限 = 2m 3 (3-36) 同理可得,n 条边的多边形图形条件闭合差的限值为 w限 = 2m n (3-37) 图 3–4 图 3–5