1、空间数据误差处理空间数据误差处理Surveying Adjustment第1页第1页第四章第四章 平差数学模型与平差数学模型与最小二乘原理最小二乘原理第2页第2页回回顾顾测量平差测量平差求平差值求平差值精度评估精度评估何谓平差值何谓平差值如何去求?如何去求?精度指标精度指标精度指标求法精度指标求法方差、协方差阵方差、协方差阵权、协因数阵权、协因数阵定义定义广义传播律广义传播律单位权方差单位权方差相应权倒数相应权倒数第3页第3页第四章第四章平差数学模型与最小二乘原理平差数学模型与最小二乘原理v4-1 测量平差概述测量平差概述v4-2函数模型函数模型v4-3 函数模型线性化函数模型线性化v4-4测
2、量平差数学模型测量平差数学模型v4-5参数预计与最小二乘原理参数预计与最小二乘原理v小小结结第4页第4页4-1测量平差概述测量平差概述 几种主要概念几何模型、几何量、函数模型、必要观测、多出观测、条件方程 测量平差第5页第5页4-1测量平差概述测量平差概述v一、主要概念一、主要概念w1.几何模型几何模型在测量工作中,为了拟定待定点高程,需要建立水在测量工作中,为了拟定待定点高程,需要建立水准网,为了拟定待定点平面坐标,需要建立平面控准网,为了拟定待定点平面坐标,需要建立平面控制网(包括测角网、测边网、边角网),我们常把制网(包括测角网、测边网、边角网),我们常把这些网称为这些网称为几何模型几何
3、模型。第6页第6页4-1测量平差概述测量平差概述v2.几何量几何量v每种几何模型都包含有不同几何元素,如水准每种几何模型都包含有不同几何元素,如水准网中包含点高程、点网中包含点高程、点间间高差,平面网中包含角高差,平面网中包含角度、度、边长边长、边边坐坐标标方位角以及点二方位角以及点二维维或三或三维维坐坐标标等元素。等元素。这这些元素都被称些元素都被称为为几何量。几何量。第7页第7页4-1测量平差概述测量平差概述w3.函数模型函数模型要拟定一个几何模型,并不需要知道其中所有元素要拟定一个几何模型,并不需要知道其中所有元素大小,只需知道其中一部分就能够了,其它元素能大小,只需知道其中一部分就能够
4、了,其它元素能够通过它们之间函数描述而拟定出来,这种描述所够通过它们之间函数描述而拟定出来,这种描述所求量与已知量之间关系式称为求量与已知量之间关系式称为函数模型函数模型。第8页第8页4-1测量平差概述测量平差概述w4.必要观测必要观测(1 1)必要元素必要元素:能够唯一地拟定一个几何模型所必能够唯一地拟定一个几何模型所必要元素,称为要元素,称为必要元素必要元素。必要观测元素个数用。必要观测元素个数用t 表示,表示,称为称为必要观测个数必要观测个数。(2)必要元素选取与性质)必要元素选取与性质 能能唯一唯一拟定该模型拟定该模型至少至少需要需要元素间元素间不存在不存在任何拟定关系任何拟定关系第9
5、页第9页4-1测量平差概述测量平差概述(3 3)独立量)独立量一个几何模型必要观测元素之间是不存在任何拟定一个几何模型必要观测元素之间是不存在任何拟定函数关系,即其中任何一个必要观测元素不也许表函数关系,即其中任何一个必要观测元素不也许表示为其余必要观测元素函数。这些彼此不存在函数示为其余必要观测元素函数。这些彼此不存在函数关系量称为关系量称为函数独立量函数独立量。第10页第10页4-1测量平差概述测量平差概述w5.多出观测多出观测假设对模型中几何量总共观测假设对模型中几何量总共观测n 个,个,n t,能够拟定模型,还能够发觉粗差能够拟定模型,还能够发觉粗差多出观测多出观测就是观测值个数就是观
6、测值个数n 与必要观测个数与必要观测个数t 之之差,用差,用r 表示。表示。r=n t 在统计学中也叫在统计学中也叫自由度自由度。第11页第11页4-1测量平差概述测量平差概述v例例1 1v例例2 2n=6,t=3,r=3n=7,t=4,r=3ADCBh1h4h3h5h6h2第12页第12页4-1测量平差概述测量平差概述w6.条件方程条件方程每增长一个多出观测,在它们中间就必定增长且只每增长一个多出观测,在它们中间就必定增长且只增长一个拟定函数关系式,有多少个多出观测,就增长一个拟定函数关系式,有多少个多出观测,就会增长多少个这样关系式。这种函数关系式,在测会增长多少个这样关系式。这种函数关系
7、式,在测量平差中称为量平差中称为条件方程条件方程。第13页第13页4-1测量平差概述测量平差概述v二、测量平差二、测量平差误差存在误差存在三角形图形三角形图形条件闭合差条件闭合差第14页第14页4-1测量平差概述测量平差概述求更正数求更正数V消除矛盾消除矛盾产生矛盾产生矛盾多出观测多出观测平差值平差值一个测量平差问题,首先要由观测值和未知量间一个测量平差问题,首先要由观测值和未知量间构成函数模型,然后采用一定平差原则对未知量进构成函数模型,然后采用一定平差原则对未知量进行预计,最后分析结果精度。行预计,最后分析结果精度。评估精度评估精度函数函数模型模型VTPV=min第15页第15页4-2函数
8、模型函数模型 条件平差函数模型 附有参数条件平差函数模型 间接平差函数模型 附有限制条件间接平差函数模型第16页第16页4-2函数模型函数模型测量数据函数模型:测量数据函数模型:几何模型、物理模型几何模型、物理模型或或几何、物理模型几何、物理模型。(测量控制网如水准网、三角网、测量控制网如水准网、三角网、GPS网等都属于几何模型网等都属于几何模型)建立不同函数模型,就有不同平差方法,测量中惯用有:1.条件平差函数模型2.附有参数条件平差函数模型3.间接平差函数模型4.附有限制条件间接平差函数模型第17页第17页4-2函数模型函数模型v一、条件平差函数模型一、条件平差函数模型w1.条件平差法条件
9、平差法:以条件方程为函数模型平差办以条件方程为函数模型平差办法,称为条件平差法。法,称为条件平差法。ADCBh1h4h3h5h6h2第18页第18页4-2函数模型函数模型第19页第19页4-2函数模型函数模型w2.特点特点(1)一个平差问题,条件形式不唯一!)一个平差问题,条件形式不唯一!选取形式最选取形式最简为宜简为宜。(2)各条件式之间是)各条件式之间是独立独立。找出观测值真值之间应。找出观测值真值之间应当满足当满足r个不相关关系式个不相关关系式。第20页第20页w3.模型模型4-2函数模型函数模型ADCBh1h4h3h5h6h2第21页第21页4-2函数模型函数模型第22页第22页4-2
10、函数模型函数模型第23页第23页 假如有假如有n个观测值个观测值 ,必要观测个数为,必要观测个数为 t,则应,则应列出列出 r=n-t 个条件方程个条件方程:线性形式:线性形式:令令条件平差自由度为多出观测数条件平差自由度为多出观测数r,即条件方,即条件方程个数。程个数。4-2函数模型函数模型第24页第24页4-2函数模型函数模型v例例1 1:已知点:已知点A、B高程高程;观测值:;观测值:h1h5t=2,n=5,r=3第25页第25页4-2函数模型函数模型v二、附有参数条件平差函数模型二、附有参数条件平差函数模型w1.附有参数条件平差法附有参数条件平差法:设在平差问题中,观设在平差问题中,观
11、测值个数为测值个数为n,t为必要观测数,则可列出为必要观测数,则可列出r=n-t个条件方程,既有增设了个条件方程,既有增设了u个独立量作为参数,个独立量作为参数,而而0ut,每增设一个参数应增长一个条件方程,每增设一个参数应增长一个条件方程。以含有参数条件方程作为平差函数模型,称。以含有参数条件方程作为平差函数模型,称为附有参数条件平差法。为附有参数条件平差法。第26页第26页4-2函数模型函数模型w2.模型模型(1)仍然按条件平差列仍然按条件平差列r 个方程个方程(2)然后增长然后增长一个一个参数,就会增长参数,就会增长一个一个条件方程:条件方程:(3)联立)联立写成矩阵形式写成矩阵形式第2
12、7页第27页4-2函数模型函数模型(1)按条件平差列按条件平差列函数模型函数模型(2)选)选C点高程为参数,则点高程为参数,则增长一个增长一个条件方程:条件方程:(3)写成矩阵形式写成矩阵形式t=2,n=5,r=3第28页第28页4-2函数模型函数模型普通而言,其普通形式为:普通而言,其普通形式为:假如条件方程是线性,其形式为:假如条件方程是线性,其形式为:令令 此平差问题,由于选择了此平差问题,由于选择了u u个独立参数,方个独立参数,方程总数由程总数由r r 个增长到个增长到 c=r+uc=r+u 个,故平差自由度个,故平差自由度为为 r=c-ur=c-u。第29页第29页4-2函数模型函
13、数模型v3.特点特点v(1)“特殊特殊”条件平差条件平差v(2)不同于)不同于间间接平差,参数接平差,参数u选择选择:0ut个参数,其中包括个参数,其中包括t个独立参数个独立参数,则多选,则多选s=u-t个参个参数必定是数必定是t个独立参数函数,即在个独立参数函数,即在u个参数之间存在着个参数之间存在着s个个函数关系式。方程总数函数关系式。方程总数c=n+s=r+t+s=r+u个,建立模型时个,建立模型时,除了列立,除了列立n个个观测方程(间接平差)观测方程(间接平差)外,还要增长参外,还要增长参数之间满足数之间满足s个条件方程,以此作为平差函数模型平差个条件方程,以此作为平差函数模型平差办法
14、称为附有限制条件间接平差法。办法称为附有限制条件间接平差法。第41页第41页4-2函数模型函数模型w2.模型模型(1)t=2(2)选选u=2+1 个参数:个参数:(3)产生约束条件产生约束条件(4 4)函数模型:)函数模型:(5 5)写成矩阵形式写成矩阵形式间接平差间接平差约束条件约束条件第42页第42页4-2函数模型函数模型间接平差间接平差约束条件:约束条件:第43页第43页4-2函数模型函数模型函数模型普通形式为:函数模型普通形式为:线性形式函数模型为:线性形式函数模型为:令:令:第44页第44页4-2函数模型函数模型w3.3.特点特点(1 1)“特殊特殊”间接平差,即仍要选参数,参数个数
15、间接平差,即仍要选参数,参数个数 ut(2 2)多选参数个数)多选参数个数 s=u-t,参数之间不独立,产生,参数之间不独立,产生 s 个函数个函数式式(3 3)函数模型由两部分构成函数模型由两部分构成:间接平差构成方程间接平差构成方程 参数之间条件方程参数之间条件方程(4 4)函数模型总数)函数模型总数 c=n+s=r+t+s=r+u(5 5)自由度)自由度 r=n-t=n-(u-s)第45页第45页4-2函数模型函数模型v条件条件w特殊条件平差特殊条件平差w需要选参数,且参数间独立需要选参数,且参数间独立w参数个数参数个数u:0utw函数模型个数:函数模型个数:c=r+u=r+t+s=n+
16、sw函数模型类型函数模型类型 1.1.按间接平差观测方程按间接平差观测方程 2.2.参数之间条件方程(限制条件参数之间条件方程(限制条件式)式)附有参数条件平差与附有限制条件间接平差比较附有参数条件平差与附有限制条件间接平差比较第46页第46页4-2函数模型函数模型四种平差办法与参数关系四种平差办法与参数关系参数个数参数个数平差办法函数模型普通式u=0条件平差条件平差0ut附有限制条件间接平差第47页第47页4-2函数模型函数模型v例例5.5.如图所表示,如图所表示,A为已知点,为已知点,B、C、D、E为待定点,为待定点,观测了观测了9条路线高差条路线高差h1h9,列出下面四种情况函数模型,列
17、出下面四种情况函数模型,试指出所列情况属于基本平差办法哪一类函数模型,并指试指出所列情况属于基本平差办法哪一类函数模型,并指出方程个数。出方程个数。ACB9876521E43D第48页第48页4-2函数模型函数模型(1 1)选取)选取B、C、D三三点高程平差值为参数点高程平差值为参数(2 2)选取)选取h1h5高差平差值为参数高差平差值为参数(3 3)选取)选取h5h8平差值为参数平差值为参数(4 4)选取)选取C、E两点间高差为参数两点间高差为参数ACB9876521E43u=3t,附有限制条件间接平差模型,附有限制条件间接平差模型,c=r+u=5+5=10 u=4=t,间接平差函数模型,间
18、接平差函数模型,c=n=9 u=10,有,有,则称则称是是 一致预计量。一致预计量。第68页第68页4-5参数预计与最小二乘原理参数预计与最小二乘原理v引例引例1.为了拟定刀具磨损速度,做下列试验:通过一定期为了拟定刀具磨损速度,做下列试验:通过一定期间(如每隔一小时),测量一次刀具厚度,得到一组试验间(如每隔一小时),测量一次刀具厚度,得到一组试验数据下列:数据下列:编号编号0 01 12 23 34 45 56 67 7时间时间ti0 01 12 23 34 45 56 67 7刀具厚度刀具厚度yi27.027.026.826.826.526.526.326.326.126.125.725
19、.725.325.324.824.8第69页第69页4-5参数预计与最小二乘原理参数预计与最小二乘原理f(t)=at+b问题:拟定问题:拟定a,b选取选取a,b,使得使得 f(t)在在t0,t1,t2t7函数值与试验数据函数值与试验数据y0,y1,y2y7 相差很小,即使偏差相差很小,即使偏差 yi-f(ti)最小。最小。最小二乘法最小二乘法第70页第70页4-5参数预计与最小二乘原理参数预计与最小二乘原理v引例引例2.匀速运动质点在时刻匀速运动质点在时刻 位置位置y表示为:表示为:第71页第71页4-5参数预计与最小二乘原理参数预计与最小二乘原理写成矩阵:写成矩阵:间接平差函数模型间接平差函
20、数模型第72页第72页4-5参数预计与最小二乘原理参数预计与最小二乘原理最小二乘法最小二乘法第73页第73页4-5参数预计与最小二乘原理参数预计与最小二乘原理v二、最小二乘原理二、最小二乘原理w1.定义定义在满足在满足条件下解出参数估值,这种求预计量办法就称为最条件下解出参数估值,这种求预计量办法就称为最小二乘法小二乘法或或第74页第74页4-5参数预计与最小二乘原理参数预计与最小二乘原理w2.权阵权阵独立观测值:独立观测值:相关观测值:相关观测值:第75页第75页4-5参数预计与最小二乘原理参数预计与最小二乘原理v例例1 1:设对某物理量:设对某物理量 进行了进行了n次同精度独立观测,得观测次同精度独立观测,得观测值值 ,试按最小二乘准则求该量预计值。,试按最小二乘准则求该量预计值。第76页第76页4-5参数预计与最小二乘原理参数预计与最小二乘原理v例例2.2.一个三角形,等精度观测了三个内角一个三角形,等精度观测了三个内角Li,求三内角平,求三内角平差值。差值。解:(解:(1 1)函数模型)函数模型 (2 2)随机模型)随机模型 (3 3)平差准则)平差准则第77页第77页平差办法函数模型函数模型随机模型随机模型备注备注条件平差条件平差附有参数条件平差间接平差间接平差附有限制条件间接平差小小结结第78页第78页